Python数据分析之锁具装箱问题问题重述问题分析建模与求解

问题重述

某厂生产一种弹子锁,其槽数高度可以用1到6中取5个来表示。其限制条件是:至少在5个中有3个不同的数;相邻槽的高度相差不能为5。在实际试验中,发现若二锁对应5个槽的高度中有4个相同,另一个差1则可能互开,否则,不可能互开。如果60个锁具装一箱,求一批锁的多少及装箱数,并要求提出一种方案,使团体顾客减少或不再抱怨,并对于所提出的方案,求出其最大无互开的箱数,并衡量原来随机装箱时,顾客抱怨互开的程度。

问题分析

  1. 锁具个数 首先把锁具及装箱问题抽象成数学概念,以5个数字的一个符合条件的组合或组成一个列表数据结构来代表一个相应的锁具,例如[1,2,3,4,5]代表一个锁具。利用排除法的思想,通过Python语言,把问题分为所有可能的组合A6^5,存储到列表结构中;再通过集合的概念,剔除列表中相同槽高,只保留一个,如果个数小于3,则不符合要求,进行剔除;最后,剔除减去相邻差为5的列表。
  2. 装箱方案设计 归类法:把一批锁具分成两类或按其指数之和分为8~27共20类,其中组合,[1,1,1,2,3]为最小,其数字和为8;[6,6,6,5,4]为最大,其数字和为27;所以共有(27-8+1)=20个数组来表示,在装箱时可依: d9→d11→…→d27→d 8→d10→d12→… →d27 或 d8→d10→…→d26→d9→ …→d26 二者类似,其中di表示5槽数字和为i的所有组合。 事实上,对于一数组来说,其各位数之和相同,所以其本身显然不可能有互开现象。对于其它情况可用下法来鉴别:利用Python中的Numpy第三方库进行列表向量化运算,两个列表相减,再进行求和,绝对值为1则二者可能互开;否则不可能互开,记数字和为i的数组di,另一组数据为d_i^',则有以下结论: d_i^'∩d_((i+1))≠∅ d_i^'∩d_((i))=∅ d_i^'∩d_((i+j))=∅ 可知: d_26^'∩d_27≠∅ d_8^'∩d_27≠∅ 而对数字和全为偶或全为奇的数组间一定不互开,所以以上设计符合要求。
  3. 顾客抱怨定量度量 对于原来随机装箱时,求其顾客的不满意程度。定量表示时,我们认为可用下列量来表示: A:一定量中出现的互开锁具数; B:一定量中可能互开的概率; C:用计算机模拟仿真随机装箱过程,随机抽点进行计算。

建模与求解

1.锁具个数 利用排除法的思想,通过Python语言,逐步剔除不符合要求的锁具,可得有5880个锁具,每60个一箱,可装98箱,代码如下:

lists = [1,2,3,4,5,6]
list_alls = []
for list_1 in lists:
    for list_2 in lists:
        for list_3 in lists:
            for list_4 in lists:
                for list_5 in lists:
                    list_alls.append([list_1,list_2,list_3,list_4,list_5])
print(len(list_alls))
list_deletes_1 = []
for list_all in list_alls:
    counts = list(set(list_all))
    if len(counts) <= 2:
        pass
    else:
        list_deletes_1.append(list_all)
print(len(list_deletes_1))
list_deletes_2 = []
for list_all_1 in list_deletes_1:
    if abs(list_all_1[0] - list_all_1[1]) == 5:
        pass
    elif abs(list_all_1[1] - list_all_1[2]) == 5:
        pass
    elif abs(list_all_1[2] - list_all_1[3]) == 5:
        pass
    elif abs(list_all_1[3] - list_all_1[4]) == 5:
        pass
    else:
        list_deletes_2.append(list_all_1)
print(len(list_deletes_2))

2.装箱方案设计 首先证明d9→d11→…→d27→d 8→d10→d12→… →d27无互开现象。设其中某列表为[a,b,c,d,e](a+b+c+d+e=i为奇数),则(a±1)bcde,a(b±1)cde,ab(c±1)de,abc(d±1)e,abcd(e±1),当筛选其含相邻差为5或有位为0的不合条件项外,则可知 a+b+c+d+e±1=i(±1) 则其为偶数,与d9→d27中任一元素(组合)都不可能相同,对于d8→d10→…→d26类似。 由此可证明49箱为最大无互开箱数,也就是奇偶分箱方案。 3.顾客抱怨定量度量 利用Python计算所有可能互开对为22778,平均互开系数α=0.001317,对于一箱互开个数约2.33,二箱互开个数约9.42,代码如下:

import random
i = 0
xiangjian_list = []
hukai_counts = []
list_randoms = random.sample(list_deletes_2,60)
for list_randoms_1 in list_randoms:
    i = i+1
    for list_randoms_2 in list_randoms[i:]:
        list_1 = [list_randoms_1[0] - list_randoms_2[0],list_randoms_1[1] - list_randoms_2[1],list_randoms_1[2] - list_randoms_2[2],list_randoms_1[3] - list_randoms_2[3],list_randoms_1[4] - list_randoms_2[4]]
        xiangjian_list.append(list_1)
for list_2 in xiangjian_list:
    if (list_2.count(0) == 4):
        if (list_2.count(1) == 1 | list_2.count(1) == 1):
            hukai_counts.append('a')
print(len(hukai_counts))

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏大数据风控

R中如何用ifelse进行数据分组

数据分组,根据数据分析对象的特征,按照一定的数值指标,把数据分析对象划分为不同的区间部分来研究,以揭示内在的联系和规律性; 在R中,我们常用ifelse函数来进...

2728
来自专栏tkokof 的技术,小趣及杂念

Sweet Snippet系列 之 随机选择

  平日工作学习时总会遇到一些令人欣喜的代码段子(Snippet),虽然都很短小,但是其间所含的道理都颇有意味,遂而觉得应该不时的将她们记下,一来算作复习整理,...

1092
来自专栏郭耀华‘s Blog

剑指offer第七天

30.连续子数组的最大和 HZ偶尔会拿些专业问题来忽悠那些非计算机专业的同学。今天测试组开完会后,他又发话了:在古老的一维模式识别中,常常需要计算连续子向量的最...

2629
来自专栏ml

由判断三一点是否在三角形内部而引发的思考.....

判断一个点是否在三角形里面(包括边界上),这个问题对于许多初学者来说,可谓是一头雾水,如何判断呢? 假如有四个点A(x0,y0),B(x1,y1),C(x2,y...

2948
来自专栏章鱼的慢慢技术路

笔试常考题型之时间复杂度

1473
来自专栏Python小屋

Python使用爬山算法寻找序列“最大值”

爬山算法是人工智能算法的一种,特点在于局部择优,所以不一定能够得到全局最优解,尽管效率比较高。使用爬山算法寻找序列最大值的思路是:在能看得到的局部范围内寻找最大...

3686
来自专栏潇涧技术专栏

Problem: Delete Number Problem

这题可以使用贪心策略,每次从高位向低位数,删除高位比低位数字小的那位上的数字,直到删除了k位之后,得到的数字肯定是最大值。

812
来自专栏潇涧技术专栏

Python Algorithms - C3 Counting 101

原书主要介绍了一些基础数学,例如排列组合以及递归循环等,但是本节只重点介绍计算算法的运行时间的三种方法

834
来自专栏WD学习记录

牛客网 连续子数组的最大和

HZ偶尔会拿些专业问题来忽悠那些非计算机专业的同学。今天测试组开完会后,他又发话了:在古老的一维模式识别中,常常需要计算连续子向量的最大和,当向量全为正数的时候...

842
来自专栏华章科技

程序员必须知道的十大基础实用算法及其讲解

快速排序是由东尼·霍尔所发展的一种排序算法。在平均状况下,排序 n 个项目要Ο(nlogn) 次比较。在最坏状况下则需要Ο(n2) 次比较,但这种状况并不常见。...

892

扫码关注云+社区