【项目总结】2018年全国大学生数学建模大赛B题简要分析（附代码）

杨利霞

【项目总结】2018年全国大学生数学建模大赛B题简要分析（附代码）
. T* C5 u- l2 q5 z' U, m6 A- }3 K9 K
, M/ n8 f; ^1 i; j  l* Z今天早上跟学姐室友去复旦把论文答辩做掉了，虽然整个项目基本上是我承担了主要的思路与代码部分，但是今天答辩我跟室友竟然连一句有用的话都没说出来，全场都靠学姐开场流畅地引入，临场随机应变，从而得以与答辩教授欢聚喜散。主要原因是教授竟然准确地问到了我代码里一个细节却相当致命的问题（是一个随机初始化问题，我下面代码部分会详细提到），正好学姐室友都不是特别熟悉我的随机初始化方法，我又不能当场跟他们两个解释这个随机初始化的问题。我差点当场就要以“这样随机初始化能够减少代码量”这种蹩脚的理由跟教授争辩了。好在姜还是老的辣，辩论队队长出身的学姐一顿 Speech Art 操作成功忽悠掉了两位教授，最终两位答辩教授还是认可了我们的模拟仿真方法[捂脸]。事后细想以后我成功也好，失败也罢，恐怕也是成也言语，败也言语。也许我确实能够成为一个有能力的人，但是说话艺术确实是一门很大的学问。不过看我运气一直这么差，大概率还是凡人一个落入俗套吧[摊手]。
$ W- C  E8 P/ G+ s
7 x- v1 t$ j% h# h6 }4 \8 e言归正传，本文主要介绍我们小组解决2018年全国大学生B题的思路分析，不代表标准答案。当然我还是有自知之明，本身水平不是很高，再加上三天时间限制，自己做出来的模型以及算法肯定是比较差的。这里仅仅从我个人的思考角度出发写一些参考思路作为分享讨论，希望各位读者朋友轻喷。

问题分析

今年的B题确实与往年有很大的不同。往年的数学建模问题往往具有比较好的开放性，问题解决存在较大的建模空间。今年的B题的题干本身就几乎是一个明确的模型（8台CNC+1台RGV+CNC定址），加上第二道任务要求我们根据给定三组数据完成八小时内的RGV详细调度方案，并写入四张Excel表格，给人的感觉就是要求我们去完成一道填空题，然后附带写一篇论文[尴尬]。

为了方便各位读者对赛题的阅读，这里给出链接：https://download.csdn.net/download/cy19980216/10708725
/ L- y$ P0 Y4 y5 K( S# R1 }7 N( p
/ ]! k! c4 \& _* m7 k( Y问题一共有四种不同的情况：一道工序无故障，一道工序有故障，两道工序无故障，两道工序有故障。

* y8 r$ L1 W. q4 U一道工序无故障

: u( J' Z  q: [5 O: p第一种情况是最简单的，直观上直接不停地1234567812345678……按顺序上料差不多就是最优了。但严谨地来说，虽然题目中给的三组数据确实都是用这种最幼稚的策略能够达到最优，但是如果对于一般的情况而言，比如最极端的情况下，RGV移动时间无穷大，那RGV显然就只会不停地在121212121212……这样原地上下料了。
2 w6 y* l- v1 Q6 {: _- S
3 y2 b  {6 l7 {/ C然而我们发现无论参数怎么变化，最终RGV给CNC上下料的过程始终是一个周期性过程。当然这个似乎很“显然”的事实却是相当难以通过数学严格证明的（参数已知的情况下一般比较容易证明，但是所有的参数都是未知的情况下是很难严格说明的）。我赛后也仔细的思考过，但是也没有得出很漂亮的证明。我最终仅仅是针对给定的三组数据使用了遗传算法对RGV前17次上下料（17次是考虑从初始状态开始循环两圈的最短路径）的最优路径进行了搜索，并且利用穷举证明了这是前17步最优的上下料次序。之后基本上就是不断地循环。
6 r: m3 V0 F, r. Q# a7 D* P
% e! @' o' p: S  F& Y0 a这里的模拟退火遗传算法比较鸡肋，所以我不详细说明，在第三种情况我会详细说明模拟退火遗传算法的原理。

以下给出第一种情况的模拟退火遗传算法算法以及对应的穷举最优证明 ↓↓↓
+ R% c3 @  [1 r' E5 t# -*- coding:UTF-8 -*-
& [% ^' v# K; u4 s6 b& S"""
5 I/ `# d7 V5 c+ a) K( R        作者：囚生CY
        平台：CSDN
8 P- B$ Q: z# p% q* j* v. P: l- s        时间：2018/10/09
5 q4 Y! ^* e# ~        转载请注明原作者
        创作不易，仅供分享
+ c  V. g" Y( A9 F"""
" i0 _6 e# S% w. Y# `
1 I- ~  f% s" h; a+ @- Simport math
import random
0 f; z6 w7 J+ j8 W6 Y3 Wimport itertools

""" 选取一组数据 """
T = 580
3 n* K4 u7 s. B0 `d1 = 23
$ y. t4 H( S5 q8 S7 O, nd2 = 41
$ U6 W$ A( T( ?d3 = 59
Te = 35
To = 30
Tc = 30

9 z6 K3 c0 V8 w/ y' p* C6 }; Y' TCNCT = [To,Te,To,Te,To,Te,To,Te]                                                                                 # CNC上下料时间

7 ?' J: h  h3 g! qN = 50
0 E! f% u* ?1 Z5 t9 S  }L = 17

varP = 0.1
croP = 0.6
+ X' o& `1 N+ y2 _* Q  ^
1 ]3 e# }8 [1 I3 v% b( ZcroL = 4
* E- r+ {. a9 E$ n% we = 0.99
8 l/ T2 G9 {5 G9 |1 i2 h/ z
tm = [
        [0,0,d1,d1,d2,d2,d3,d3],
        [0,0,d1,d1,d2,d2,d3,d3],
. N5 Y8 E% N+ e: W: t        [d1,d1,0,0,d1,d1,d2,d2],
# F9 `; q2 J/ i3 U/ _        [d1,d1,0,0,d1,d1,d2,d2],
3 U; U( L. U0 A' l        [d2,d2,d1,d1,0,0,d1,d1],
/ v4 z9 y) p+ E# ^- t( W7 D3 Z0 W        [d2,d2,d1,d1,0,0,d1,d1],
        [d3,d3,d2,d2,d1,d1,0,0],
; Q2 F6 u3 P5 s        [d3,d3,d2,d2,d1,d1,0,0],
]
" A* {2 O$ @! w6 `& p2 D
def update_state(state,t):
. K% L7 [7 [9 m0 L- a, r        length = len(state)
        for i in range(length):
7 b/ I, X4 P5 R0 M                if state < t:
                        state = 0
                else:
                        state -= t
        return state
! C3 _( m9 _2 s& ?
* f8 r. I8 p# ?& ^' R- ]def time_calc(seq):
        state = [0 for i in range(8)]                                                                                   # 记录CNC状态
6 ]* v: \* V4 j        isEmpty = [1 for i in range(8)]                                                                                 # CNC是否为空？
        currP = 0
        total = 0
        length = len(seq)
: @2 X0 P  }8 \3 L) h8 m- l        for No in seq:
2 |' X4 k7 C# x+ K( @1 \9 `                nextP = No
- P( x& d# y# X/ P" d8 \                t = tm[currP][nextP]
4 j! |4 [9 Z1 o* O7 f2 b" P- ]                total += t                                                                                                                 # rgv移动
& o: k% j" k1 H  X                state = update_state(state,t)                                                                         # 更新state
/ _6 `3 o- X, S) }  S                if state[No]==0:                                                                                                 # 表明CNC等待
                        if isEmpty[No]:                                                                                                 # 当前CNC空
                                t = CNCT[No]
                                isEmpty[No] = 0
9 d! b3 D( ]) y, |1 v3 F                        else:
4 j( N5 ]2 U* Y9 J3 V2 a7 e                                t = CNCT[No]+Tc
                        total += t
, z. B$ S# P/ X8 W                        state = update_state(state,t)
                        state[No] = T
                else:                                                                                                                         # 当前CNC忙
* ]- X( M* [4 x7 i2 L' o                        total += state[No]                                                                                         # 先等当前CNC结束
                        state = update_state(state,state[No])                                                 
                        t = CNCT[No]+Tc
                        total += t
3 P" m& C) R, A7 }- B2 b                        state = update_state(state,t)
/ I- }. G6 j' F& B' q                        state[No] = T
                currP = No
: |* z$ Q! {9 ]        total += tm[currP][0]
        return total

def init_prob(sample):
  B( I2 m+ `# H  `* C2 T. A6 N: v        prob = []
        for seq in sample:
                prob.append(time_calc(seq))
5 j+ Z( i$ _1 p; E4 S8 I2 Z        maxi = max(prob)
. ^3 Z  \( p! q: z8 v3 M        prob = [maxi-prob+1 for i in range(N)]
        temp = 0
        for p in prob:
% \4 L5 c4 I; @/ d. k5 u                temp += p
/ v1 ^6 u) |$ S5 m4 m% h9 c        prob = [prob/temp for i in range(N)]
        for i in range(1,len(prob)):
                prob += prob[i-1]
        prob[-1] = 1                                                                                                                 # 精度有时候很出问题
6 P$ I6 P' p# T- E. K        return prob

& U2 [1 ]2 E6 e) n9 `) o6 l. Edef minT_calc(sample):
3 t% D/ X* F  y# `8 k  G        minT = time_calc(sample[0])
; |+ T- m& ^* k" [2 t        index = 0
# l3 d: ?8 G, h$ p- _        for i in range(1,len(sample)):
2 k9 i3 ^  P% I8 E5 L4 B# B                t = time_calc(sample)
( x9 L5 m2 a9 f                if t < minT:
( I- N  U! b8 Y4 o+ B                        index = i
; ?( g7 s% G: B6 R; Y                        minT = t
        return minT,index
       
4 ?* y: S/ ~$ t8 D1 Udef init():
        sample = []
        for i in range(N):
% O  V8 E  ~0 f7 G/ c$ ^% r+ j! J                sample.append([])
                for j in range(L):
                        sample[-1].append(random.randint(0,7))
$ |5 x; N( N7 g  y4 i" H        return sample

def select(sample,prob):                                                                                                 # 选择
) f& {# `; V* \- a; f" g        sampleEX = []
2 [6 C- Z5 U% K& [2 B, Y        for i in range(N):                                                                                                         # 取出N个样本
                rand = random.random()
                for j in range(len(prob)):
) M2 E2 b  N3 y. V' j7 I                        if rand<=prob[j]:
                                sampleEX.append(sample[j])
2 W, v- @9 s9 z+ E' R! k* D, M$ s                                break
1 x8 m  b1 a! D& C        return sampleEX
, O- A& h0 l. \( z" s5 b& J
def cross(sample,i):                                                                                                         # 交叉
        for i in range(len(sample)-1):
                for j in range(i,len(sample)):
                        rand = random.random()
                        if rand<=croP*(e**i):                                                                                 # 执行交叉
# ^. Z* `/ F( a0 u                                loc = random.randint(0,L-croL-1)
, E! z# p) g2 ?5 ^                                temp1 = sample[loc:loc+croL]
                                temp2 = sample[j][loc:loc+croL]
3 s+ Q  d5 F( F  q- M( R! \                                for k in range(loc,loc+croL):
) x" S& l0 @0 s8 n. ]& r+ |                                        sample[k] = temp2[k-loc]
                                        sample[j][k] = temp1[k-loc]
        return sample
* Q4 M" s: A$ `0 e; J! k, V               
def variance(sample,i):                                                                                                         # 变异算子                                                                                 
, m" x7 P# h7 S8 _        for i in range(len(sample)):
                rand = random.random()
7 o3 ~7 \( U2 r2 S& |                if rand<varP*(e**i):
5 a0 H4 C; `  G% a, H4 m                        rand1 = random.randint(0,L-1)
* M5 k9 g$ W: D0 q1 C5 t) L0 X! R                        rand2 = random.randint(0,L-1)
! |3 I& v  K& e" L! J! l0 _% P* ~. @                        temp = sample[rand1]
                        sample[rand1] = sample[rand2]
                        sample[rand2] = temp
9 F. J+ J" Z+ F* ^" o/ Z        return sample
       
& X- j5 K# A' x0 s  d( M$ fdef main():
$ |" \7 v+ a2 m% }* N        sample = init()
0 r! }# H* _& J/ r        mini,index = minT_calc(sample)
        best = sample[index][:]
        print(best)
8 r, i& E& c2 L  P" n1 x- R        for i in range(10000):
                print(i,'\t',minT_calc(sample),end="\t")
                prob = init_prob(sample)
3 u0 T' b" J  ?* \2 O                sample = select(sample,prob)
                sample = cross(sample,i)
' @) W- Z# ?9 r) a# W                sample = variance(sample,i)
6 U3 h4 o, |) C' }% V7 L5 y5 u9 p; ]                mi,index = minT_calc(sample)
                if mi>mini and random.random()<e**i:                                                         # 精英保留策略
6 `6 P5 }! c) p                        rand = random.randint(0,N-1)
* F& e. p9 A0 h  M! y7 O                        sample[rand] = best[:]
- J3 I  U: A4 x0 v. x- i: }                mini,index = minT_calc(sample)
                best = sample[index][:]
                print(best)
( x- w, m  q; @% p' M        print(sample)
. y2 L( E- O+ W+ ]. T! m
* U, P1 D/ {" O% E: [if __name__ == "__main__":
        main1()
* H! _# W) u. X6 X" B( w2 E        """ 穷举搜索验证 """
        a = list(itertools.permutations([1,2,3,4,5,6,7],7))
        ts = []
        first = [0,1,2,3,4,5,6,7,0]
        for i in a:
3 M- x8 m8 g2 N( k                temp = first+list(i)
                temp.append(0)
" F) g' F8 J6 v8 v" w/ v                t = time_calc(temp)
                ts.append(t)
        print(min(ts))       
' F( f( _  {9 K# t        print(time_calc([0,1,2,3,4,5,6,7,0,1,2,3,4,5,6,7,0]))
       

一道工序有故障
7 N1 L$ ^$ y8 ~( W4 v& {. C+ o8 p
+ L4 S/ F1 e9 s+ t- y9 B4 d这部分是学姐做的，学姐用了偏数学的思考方式，仍然从循环的角度去考虑，主要考虑故障发生是否会影响当前循环，是否需要建立新的循环。因此就没有写代码处理问题了。具体的思路我确实不是很能讲清楚。但是这里面有一个非常大的问题，就是如果出现多台CNC同时发生故障怎么办。关于多台机器同时发生故障的概率，我们通过估算认为以给定的三组数据8小时内会出现这种特殊情况的可能性大约为30%。这个问题是我无法很好严格处理的（当然如果用贪心算法也就没这么多事了）。

" E4 A7 @$ F+ |; r5 ]: G* ?两道工序无故障 & 两道工序有故障

- V( H$ q% s4 d0 `; ?) P这两个部分都是我来处理的，因为使用的方法大致相同，就并在一起说了。

两道工序与一道工序最大的区别在于三点：

1、开始要处理CNC任务分配：分配给第一道工序几台CNC，分配给第二道工序几台CNC？具体怎么布局？
) G& M1 t& M1 e/ P  g" B
$ ]8 g* K# v' R& r& a2、加工过程可能仍然是一个循环，但是这个循环将可能会非常的庞大以至于不可能直观的看出来。

3、两道工序的分配已经是一个严格的NP难问题了（即理论上无法在多项式时间内求得最优解）。
! w+ l  w. y: M, Z' ~
第一点我的想法很单纯——穷举，没错，就是穷举，除了显然不合适的分配方案外，其他方案都试一遍（虽然真的很蠢，但是我真的想不出到底能怎么办了）
5 `0 x8 E  S8 P! P; W
, l2 Y$ c* L# i  y- H5 L, m* L第二点因为不存在循环则使用遗传算法需要设定一个相当长的染色体长度（我们设定的染色体是RGV为各台CNC上下料的次序，如果要考虑全过程的模拟退火遗传算法，则染色体长度大约在300~400左右）。事实上我也尝试了这个方法，结果从我写完这个算法我开始跑，一直跑到比赛结束算法依旧没有收敛[捂脸]。这里给出代码仅供参考（各位朋友要是有好意见也可以提出） ↓↓↓

$ y) P+ D" N1 x# -*- coding:UTF-8 -*-
/ z2 e9 W' ~9 }"""
  r; W; e1 F- H) E5 B. E/ f, Y        作者：囚生CY
        平台：CSDN
        时间：2018/10/09
        转载请注明原作者
        创作不易，仅供分享
"""
8 h: [/ M! B! B6 s) i. f5 Rimport random

: S' J2 g+ D- h) t* }( v' I# 第1组
7 c- @/ v. s6 ^- J2 M"""
d1 = 20
d2 = 33
$ v& ~# g9 ]; d9 R. t( j. V9 ]$ ud3 = 46
6 f. ^6 D( H" LT1 = 400
- s6 ~4 H2 V6 L" _& ~T2 = 378
To = 28
. }1 H. I6 w( z5 CTe = 31
" A7 l6 Z* E! y  W+ H& R: CTc = 25
"""

5 d: \) ]3 l6 {+ }* Y" @' @3 h# 第2组
8 A1 \" ]9 {3 J7 l"""
; v% x0 O" E: A' ~d1 = 23
9 g3 v1 C& |* S4 B+ c/ vd2 = 41
0 r: `8 b/ E( q; V) t9 @7 M8 M  yd3 = 59
T1 = 280
T2 = 500
To = 30
7 l$ f% a& o/ x, q! O% E5 uTe = 35
Tc = 30
0 Q  ~5 X% _3 G9 ?; }: P"""

9 v. [& d% i5 `8 v( {; G3 O  {# 第3组
d1 = 18
d2 = 32
d3 = 46
T1 = 455
T2 = 182
To = 27
Te = 32
% e0 U9 I) ~! s  y9 R& y6 w6 {Tc = 25

" S" p( y& W9 R# s  N. X# \cncT = [To,Te,To,Te,To,Te,To,Te]
tm = [
        [0,0,d1,d1,d2,d2,d3,d3],
- q+ p9 ~; l( O1 O8 k* c7 P0 q8 f        [0,0,d1,d1,d2,d2,d3,d3],
        [d1,d1,0,0,d1,d1,d2,d2],
        [d1,d1,0,0,d1,d1,d2,d2],
        [d2,d2,d1,d1,0,0,d1,d1],
        [d2,d2,d1,d1,0,0,d1,d1],
. [' ]" L4 g' T: c) j4 ^- ^) k' U        [d3,d3,d2,d2,d1,d1,0,0],
        [d3,d3,d2,d2,d1,d1,0,0],
* J/ ?/ M' R( a7 F$ C% |# D]
Type = [1,0,1,0,1,0,0,0]                                                                                                 # CNC刀具分类
, Q# O8 Z" Q0 m! X- x  \
2 ~* M6 F4 c: A  [N = 64
2 n) x' K9 N* H9 q  L# F9 yL = 100
" ]+ L$ `0 X3 f0 ]% Q5 Y% LvarP = 0.1
croP = 0.6
7 I$ S2 v4 T2 \+ BcroL = 2
e = 0.99

* R$ d) G8 g# xdef init_first_round():                                                                                                         # 第一圈初始化（默认把所有第一道CNC按顺序加满再回到当前位置全部加满）
        state = [0 for i in range(8)]                                                                                   # 记录CNC状态（还剩多少秒结束，0表示空闲）
        isEmpty = [1 for i in range(8)]                                                                                 # CNC是否为空
        rgv = 0                                                                                                                                 # rgv状态（0表示空车，1表示载着半成品）
. t* ^  d/ ~5 ]6 T; z& x6 Q        currP = 0
        total = 0
        seq = []
" l6 [9 d# J# }2 m        flag = False
        for i in range(len(Type)):
  d2 n* m) @, A; d: d                if Type==0:
                        seq.append(i)
6 \! r9 M. ]6 V, w! [% b                        flag = True
        currP = seq[0]
        seq.append(currP)
/ m2 j+ d6 m( @+ k        rgv,currP,total = time_calc(seq,state,isEmpty,rgv,currP,total)
: X" Y" n: @# H7 a* x        return state,isEmpty,rgv,currP,total,seq

* Q1 r% A2 O7 `$ R. _5 ldef update(state,t):
; H& f0 o' w' @& K9 b/ |7 l        for i in range(len(state)):
% u, y) W; z6 c) G: l, e                if state < t:
% R4 G2 D0 I' p% i8 r2 N9 x, N                        state = 0
                else:
                        state -= t
$ P: Z9 j( c6 h2 g- O  Q: o* b
def time_calc(seq,state,isEmpty,rgv,currP,total):                                                 # 事实上sequence可能是无效的，所以可能需要
        index = 0
' x; F9 ?! B+ A! Q        temp = 0
        while index<len(seq):
                """ 先移动到下一个位置 """
& A0 K* u3 S: G% l7 y                nextP = seq[index]
' @7 y$ J2 ?+ ~/ z0 p/ r! q& A                t = tm[currP][nextP]
, W) j% C0 A7 q% \                total += t
                update(state,t)
& a8 r3 {6 E0 O; S/ k                if Type[nextP]==0:                                                                                                 # 如果下一个位置是第一道工作点
                        if rgv==1:                                                                                                         # 然而载着半成品
" }3 U' ^1 a1 E" k4 c& X  x. Q                                seq.pop(index)                                                                                         # 去掉这个元素并中止当次循环进入下一个循环
* b3 }9 P' Q  L$ A* C                                continue                               
9 e: ~3 U$ @! W6 m                        if isEmpty[nextP]:                                                                                         # 如果下一个位置是空的
                                t = cncT[nextP]
. s$ M, U& |7 p. A0 e                                total += t
# v- {: F5 L5 T                                update(state,t)
                                state[nextP] = T1                                                                                 # 更新当前的CNC状态
                                isEmpty[nextP] = 0                                                                                 # 就不空闲了
4 F- p, f9 J  z8 T  ^                        else:                                                                                                                 # 如果没有空闲
. M3 ?. F# A$ O  Q9 b1 Q7 m                                if state[nextP] > 0:                                                                         # 如果还在工作就等待结束
; S+ r$ y; k- R7 ^- `                                        t = state[nextP]
: p+ q- ^2 ]2 ^1 [7 v8 @6 X0 D" W                                        total += t
- }* v# L1 p) L( m" Z                                        update(state,t)
                                t = cncT[nextP]                                                                                         # 完成一次上下料
                                total += t
/ A- Y# g5 D: a: t/ v                                update(state,t)
, K0 l- ]) {7 }# w- K                                state[nextP] = T1
* ~1 O/ o- n0 u# r                                rgv = 1
                else:                                                                                                                         # 如果下一个位置是第二道工作点
                        if rgv==0:                                                                                                         # 如果是个空车
. g* T7 ~, w4 v$ h: G                                seq.pop(index)                                                                                         # 删除当前节点
- X0 N- D6 j1 v4 ~9 t8 I/ w0 p                                continue
% r' S: x7 B2 P9 I! l* K                        if isEmpty[nextP]:                                                                                         # 如果下一个位置是空的
                                t = cncT[nextP]
                                total += t
                                update(state,t)
                                state[nextP] = T2
                                isEmpty[nextP] = 0       
                        else:                                                                                                                 # 如果没有空闲
! m$ M1 m9 \3 Y                                if state[nextP] > 0:                                                                         # 如果还在工作就等待结束
                                        t = state[nextP]
) v$ c0 c! Y) I# [# ?+ f# k                                        total += t
" f4 x% a  [2 ]% l                                        update(state,t)
                                t = cncT[nextP]+Tc
5 K5 }' A) v) V# b; w) ^; `7 K                                total += t
                                update(state,t)
- V% _" n. T! e/ n3 U3 F- z+ Z                                state[nextP] = T2
9 R. ]" v! q. q7 Z2 Q                        rgv = 0
                currP = nextP
                temp = total
                index += 1       
        total += tm[currP][Type.index(0)]                                                                         # 最后归零
7 z' \2 D& Z$ V( U8 a# j        return rgv,currP,total
* o* d1 }1 [! L
def init_prob(sample,state,isEmpty,rgv,currP,total):                                         # 计算所有sample的
9 z; L2 y4 }& f2 C4 t- A; ?        prob = []
5 i5 d2 u9 l& \% y5 N# b        for seq in sample:
                t = time_calc(seq,state[:],isEmpty[:],rgv,currP,total)[-1]
                prob.append(t)
  n/ B9 u+ D, @, w        maxi = max(prob)
( T( `$ W; ]* @& J; F% m+ e+ S) ^        prob = [maxi-prob+1 for i in range(N)]
        temp = 0
2 M* w- ?) O3 E- k6 h% Y7 S        for p in prob:
2 E/ I$ X) S' A3 p: |' {, K% M9 W0 i5 F                temp += p
        prob = [prob/temp for i in range(N)]
        for i in range(1,len(prob)):
                prob += prob[i-1]
. w9 z- E" l& d/ w1 C        prob[-1] = 1                                                                                                                 # 精度有时候很出问题
        return prob

9 q4 t& Y% Z; R8 P! K% adef minT_calc(sample,state,isEmpty,rgv,currP,total):
$ `9 a( C* g0 X        minT = time_calc(sample[0],state[:],isEmpty[:],rgv,currP,total)[-1]
        index = 0
        for i in range(1,len(sample)):
                t = time_calc(sample,state[:],isEmpty[:],rgv,currP,total)[-1]
2 d5 E7 e7 ^  B2 m- C$ u                if t < minT:
  g  X4 a* u- H3 [. H: o; l                        index = i
$ X, d' f% g" V+ q2 {2 x6 b                        minT = t
        return minT,index
& K) G/ ]: H( [* t9 q6 x2 X       
def init():                                                                                                                                 # 初始化种群（按照第二道工序，第一道工序，第二道工序，第一道工序顺序排列即可）
4 g0 G0 ?& I( A0 z+ x! C- ?        sample = []
        refer0 = []
- p& |" n0 ]2 U" M        refer1 = []
        for i in range(8):
; j0 c3 W% L( a                if Type==0:
4 |- n* C. h& Z                        refer0.append(i)
                else:
5 d5 d" I) J# z9 Y- v. ]2 {; b                        refer1.append(i)
; Q4 k! c" `& {8 T        for i in range(N):
8 A1 p) q, k4 a                sample.append([])
; f! G8 k' r0 [* b$ ]# w8 @- Y                for j in range(L):
                        if j%2==0:
9 m" X) ^  Z: L# p9 @                                sample[-1].append(refer1[random.randint(0,len(refer1)-1)])
                        else:
                                sample[-1].append(refer0[random.randint(0,len(refer0)-1)])
# a( o2 g4 \8 J& N        return sample

. t$ A$ [) }( S3 ]0 `; N9 U- Mdef select(sample,prob):                                                                                                 # 选择算子
- ?7 g( Q% f! O8 @+ f* N        sampleEX = []
        for i in range(N):                                                                                                         # 取出N个样本
                rand = random.random()
  A. e, N7 V0 J8 ^: B% x                for j in range(len(prob)):
9 o2 G3 U  [+ U" W6 N! Q. W  K                        if rand<=prob[j]:
7 q$ l1 S: B6 {7 J: Y                                sampleEX.append(sample[j])
6 \; d3 T$ p4 h, N& i* O                                break
* U2 ?) C- o9 H        return sampleEX

; N; u  r, Q0 j2 w3 Bdef cross(sample,i):                                                                                                         # 交叉算子
/ r; k8 h7 I  e) i( m+ g) T. _' `        for i in range(len(sample)-1):
                for j in range(i,len(sample)):
$ m& V% k! g" _" P                        rand = random.random()
' k9 {& x) ^* v: T  U                        if rand<=croP*(e**i):                                                                                 # 执行交叉
                                loc = random.randint(0,L-croL-1)
                                temp1 = sample[loc:loc+croL]
  c2 \, T4 q2 b5 P                                temp2 = sample[j][loc:loc+croL]
7 H2 C4 z- W$ J- q, D                                for k in range(loc,loc+croL):
( \' M, p5 G& u2 N) l                                        sample[k] = temp2[k-loc]
$ j- ~: B3 \8 f6 A                                        sample[j][k] = temp1[k-loc]
        return sample
               
def variance(sample,i):                                                                                                         # 变异算子                                                                                 
        for i in range(len(sample)):
                rand = random.random()
                if rand<varP*(e**i):
" T) f2 j$ q0 T                        rand1 = random.randint(0,L-1)
% A& R$ T! G3 }/ _9 T                        randTemp = random.randint(0,int(L/2)-1)
% V% o, A( N" m+ F: ?5 z                        rand2 = 2*randTemp if rand1%2==0 else 2*randTemp+1
                        temp = sample[rand1]
, T$ `: Y% O1 E& Q* A                        sample[rand1] = sample[rand2]
                        sample[rand2] = temp
        return sample
7 g) `6 J- J# b. ]5 l6 L8 _1 a
3 J4 o- l3 ^' x' K8 w7 Hif __name__ == "__main__":
4 D' g5 N/ h7 p1 D        state,isEmpty,rgv,currP,total,seq = init_first_round()
% y6 w( E1 a8 ^3 R- u        print(state,isEmpty,rgv,currP,total)
        sample = init()
  R( N7 A  K! m- Y; f' ~0 _        mini,index = minT_calc(sample,state[:],isEmpty[:],rgv,currP,total)       
        best = sample[index][:]
* h5 Q" e& j2 ?* f1 ]& b        for i in range(100000):
' }/ w2 H5 X9 j                f = open("GA.txt","a")
                tmin = minT_calc(sample,state[:],isEmpty[:],rgv,currP,total)[0]
0 p5 V5 |% Q3 ]" ^                f.write("{}\t{}\n".format(i,tmin))
                print(i,"\t",tmin,end="\t")
& y. V3 Z6 U- q/ {8 @2 f$ B                prob = init_prob(sample,state[:],isEmpty[:],rgv,currP,total)
                sample = select(sample,prob)
# w; ]+ s" B. ?7 z; F" n4 ~                sample = cross(sample,i)
                sample = variance(sample,i)
. R( P% u, L: S6 _                mi,index = minT_calc(sample,state[:],isEmpty[:],rgv,currP,total)
) ~6 k0 m) _2 B" M, D                if mi>mini and random.random()<e**i:                                                         # 精英保留策略
                        rand = random.randint(0,N-1)
                        sample[rand] = best[:]
+ c* s% Q$ G# m                mini,index = minT_calc(sample,state[:],isEmpty[:],rgv,currP,total)
                best = sample[index][:]
$ S# U& b, X) H1 B) k! [                print(best)
" w. w* u0 \1 Y                f.close()
& |) g: H6 o$ G. J' @4 f7 `        print(sample)
遗传算法这条路被堵死后我一度陷入俗套，用最直接的贪心搞了一阵子，觉得用贪心算法（即考虑下一步的最优策略）实在是对不起这种比赛。然后我就变得——更贪心一点了。
4 f: u5 b, n$ e3 V, `6 k
+ X5 J+ K* a# p0 d/ R我试图去寻找接下来K步最优的策略，然后走一步。K=1时算法退化为贪心算法，最终我们设置为K=4（当K>=8时算法速度已经相当缓慢，而4~7的结果大致相同，且K=4的速度基本可以做到2秒内得到结果）。

值得注意的是我假设RGV在两道工序下只能由第一道工序的CNC到第二道工序的CNC（忽略清洗时间情况下），然后回到第一道工序的CNC，这样往复移动（这里我不说明为什么一定要这样，但是我认为确实应该是这样）。在这个规律的引导下我大大减缩了代码量以及计算复杂度。
2 n+ L4 \' Y; v5 ]
% D8 p+ h8 ~! ~  H! G. U; y# |$ _9 k然后到第四种情况我们已经没有多余时间了，只能延续使用情况三的算法，进行了随机模拟的修改，完成了第四种情况的填表。
4 B3 r  H4 Q$ O. ~
/ P( |) B/ t  ~+ y0 ^以下是第三种情况的代码（第四种类似就不上传了）↓↓↓

8 k& t; ~" `0 P) k#coding=gbk
4 ]3 J' ?& S. a! a# aimport random
) Q& L) n1 J9 m+ x# -*- coding:UTF-8 -*-
"""
8 g  ]7 r: u$ S+ l0 ]% v        作者：囚生CY
( [' s$ v1 r7 m        平台：CSDN
        时间：2018/10/09
        转载请注明原作者
) U9 |$ y1 h8 T- _( z/ R6 X1 M. C        创作不易，仅供分享
- ?! |* o1 t& V' s% I& z"""
: y3 X+ r% X" f4 Q4 r# rfrom tranToXls import *
# G5 `/ W* t9 l" _5 O
4 m( x& p9 _( Z) p# 第1组
"""
d1 = 20
d2 = 33
+ V! M! J% J! s4 O, od3 = 46
T1 = 400
T2 = 378
) K& g: W2 e& b+ v! jTo = 28
Te = 31
' b/ @% @7 Q3 M( b  E7 ATc = 25
5 \+ x" U( B- O& [: g. H7 `"""
6 m: R% V* C% U/ X- @  ?- a# 第2组

1 J+ o7 R+ B1 td1 = 23
d2 = 41
d3 = 59
( W4 B+ e+ f! S0 b( kT1 = 280
+ H$ i) O  v$ g- o, f% x) eT2 = 500
2 q; {& p. ^* R  T) b6 HTo = 30
* a3 {* ]+ h. N! _  ~Te = 35
Tc = 30

2 b7 T: f/ ~. J+ Z/ D- L
( Z" f: u" D' P8 p1 e# 第3组
; S9 h6 A9 m: a. u* o0 n
"""
+ [, m0 ?0 Y: pd1 = 18
* z& V. Y3 t& ~; n6 e2 Vd2 = 32
$ c9 {( T( i  J. H$ r1 zd3 = 46
T1 = 455
% N+ \1 E: m+ oT2 = 182
; S) d, R, b+ N* K) [To = 27
% O/ z% @8 q  J( I8 FTe = 32
9 V8 V5 r( s& V7 @" K4 v! LTc = 25
"""

' U' ?/ q( `# wcncT = [To,Te,To,Te,To,Te,To,Te]
7 l1 |% n8 L* z5 o3 i4 j9 Btm = [
        [0,0,d1,d1,d2,d2,d3,d3],
- U; i/ D( S( B" X* b4 }        [0,0,d1,d1,d2,d2,d3,d3],
        [d1,d1,0,0,d1,d1,d2,d2],
        [d1,d1,0,0,d1,d1,d2,d2],
        [d2,d2,d1,d1,0,0,d1,d1],
6 U# H8 @" ?3 G5 }- w: `        [d2,d2,d1,d1,0,0,d1,d1],
1 V& q* Q- x( @; \  Q, t: S5 R: [& m        [d3,d3,d2,d2,d1,d1,0,0],
6 {* v+ u# A8 m0 ~. l9 J        [d3,d3,d2,d2,d1,d1,0,0],
3 b) K, O6 Q+ t, s% M; _]
Type = [0,1,0,1,1,1,0,1]                                                                                                 # CNC刀具分类
- F7 K4 k' Z7 [- O5 i  n
A = []                                                                                                                                         # 储存第一道工序的CNC编号
B = []                                                                                                                                         # 储存第二道工序的CNC编号
% m4 o+ b' s" Q7 o: Kfor i in range(len(Type)):
" a, Z, T2 ~9 o+ @$ S# r$ N: N  M        if Type:
                B.append(i)
        else:
- @8 N" C" B* t- L6 t                A.append(i)
2 X8 \$ u' u& e0 T7 l
$ n. \. e6 n' B" }& A2 D* P, ]. R7 Sdef init_first_round():                                                                                                         # 第一圈初始化（默认把所有第一道CNC按顺序加满再回到当前位置全部加满）
5 G$ ~6 L1 }1 k* n, |        state = [0 for i in range(8)]                                                                                   # 记录CNC状态（还剩多少秒结束，0表示空闲）
0 J. l" X9 f2 H$ c) o0 T& g        isEmpty = [1 for i in range(8)]                                                                                 # CNC是否为空
        log = [0 for i in range(8)]                                                                                         # 记录每台CNC正在加工第几件物料
6 T) H* K: Q" Q: Y& t0 H        count1 = 0
, B' @6 N% \- s3 r0 R        rgv = 0                                                                                                                                 # rgv状态（0表示空车，1表示载着半成品）
2 v( G3 I' t, [) J. @        currP = 0
( N- ~0 {: g+ [  f  `        total = 0
        seq = []
7 S9 i. E2 p& D( ~6 u        flag = False
  {. }; \" \  q$ m& ~& P. {        for i in range(len(Type)):
4 b2 C0 ~4 c* k                if Type==0:
                        seq.append(i)
                        flag = True
        currP = seq[0]
        seq.append(currP)
( c' g8 A  y8 \' ^* N        count1,rgv,currP,total = simulate(seq,state,isEmpty,log,count1,rgv,currP,total)
        return state,isEmpty,log,count1,rgv,currP,total,seq

def update(state,t):
        for i in range(len(state)):
7 M+ `5 N' n: ?& n! g, B; d                if state < t:
2 t3 e: e% B2 j4 z9 [5 @" B6 r, U                        state = 0
                else:
% G6 j2 }' L+ _                        state -= t

def simulate(seq,state,isEmpty,log,count1,rgv,currP,total,fpath="log.txt"):        # 给定了一个序列模拟它的过程以及返回结果（主要用于模拟并记录）
        index = 0
# n" M! S9 Z4 B0 L$ E. v/ ?        temp = 0
        pro1 = {}                                                                                                                         # 第一道工序的上下料开始时间
        pro2 = {}                                                                                                                         # 第二道工序的上下料开始时间
        f = open(fpath,"a")
3 Q8 h. k/ u: }* `- P, v  N        while index<len(seq):
0 J9 g/ _. s9 H8 w' m; D% q                print(isEmpty)
$ d" d5 g' V8 W                nextP = seq[index]
9 H# d- [, h, e7 E8 Y                t = tm[currP][nextP]
3 G" }9 i9 j' |6 [( J                total += t
                update(state,t)
                if Type[nextP]==0:                                                                                                 # 如果下一个位置是第一道工作点
                        count1 += 1
: ?  o3 Z& M# v                        if isEmpty[nextP]:                                                                                         # 如果下一个位置是空的
8 k4 E* @2 V# j                                f.write("第{}个物料的工序一上料开始时间为{}\tCNC编号为{}号\n".format(count1,total,nextP+1))
                                t = cncT[nextP]
6 i6 k' k2 y7 B. M  [0 o2 M4 e) [                                total += t
                                update(state,t)
" o( y; f% Z; {' m" K3 M7 `! d1 u* C, \                                state[nextP] = T1                                                                                 # 更新当前的CNC状态
                                isEmpty[nextP] = 0                                                                                 # 就不空闲了
                        else:                                                                                                                 # 如果没有空闲
$ h/ P4 E! ^) V5 u                                if state[nextP] > 0:                                                                         # 如果还在工作就等待结束
# e  ~' ?4 ]2 `$ ?* w* X7 ]" p/ Q                                        t = state[nextP]
, [) O; y, g' ]! Q1 n6 c                                        total += t
                                        update(state,t)
$ z% ~; Y1 a& `                                f.write("第{}个物料的工序一下料开始时间为{}\tCNC编号为{}号\n".format(log[nextP],total,nextP+1))
) F- O% `: \& b3 z, Z                                f.write("第{}个物料的工序一上料开始时间为{}\tCNC编号为{}号\n".format(count1,total,nextP+1))
                                t = cncT[nextP]                                                                                         # 完成一次上下料
                                total += t
1 j9 r$ ~9 \; }" S2 F                                update(state,t)
                                state[nextP] = T1
                                rgv = log[nextP]
                        log[nextP] = count1
                else:                                                                                                                         # 如果下一个位置是第二道工作点
: C- }+ o2 G# G. U                        if isEmpty[nextP]:                                                                                         # 如果下一个位置是空的
                                f.write("第{}个物料的工序二上料开始时间为{}\tCNC编号为{}号\n".format(rgv,total,nextP+1))
                                t = cncT[nextP]
0 M, q/ i* g% [- f: {! V: j                                total += t
5 ?! s* e1 y. E: ~. R3 k                                update(state,t)
                                state[nextP] = T2
                                isEmpty[nextP] = 0       
                        else:                                                                                                                 # 如果没有空闲
! ~& k4 z$ {5 ]( Z3 W                                f.write("第{}个物料的工序二下料开始时间为{}\tCNC编号为{}号\n".format(log[nextP],total,nextP+1))
: f+ I% t' O9 J8 o/ ]; D  ?                                f.write("第{}个物料的工序二上料开始时间为{}\tCNC编号为{}号\n".format(rgv,total,nextP+1))
& X6 k. }2 U* v* t2 G3 N! w                                if state[nextP] > 0:                                                                         # 如果还在工作就等待结束
                                        t = state[nextP]
; y! P  J% l( k( V                                        total += t
( Y6 a) s) v* q: M, ~                                        update(state,t)
' q9 ~" ^) D$ }                                t = cncT[nextP]+Tc
                                total += t
                                update(state,t)
                                state[nextP] = T2
                        log[nextP] = rgv
                        rgv = 0
1 u! j6 P3 l" |$ Q8 R1 {5 M" P8 c                currP = nextP
8 X" M: b/ v4 F+ }) ^( N+ c                temp = total
, S# ~' q- x0 V/ q                index += 1       
1 E; x' V# D- ~) M5 F4 }        f.close()
        total += tm[currP][Type.index(0)]                                                                         # 最后归到起始点
5 T7 J$ l. u$ h, v0 w- A4 ?2 n1 k" x        return count1,rgv,currP,total
2 ]1 S6 R  f" v4 G1 V
) Z3 Y( H$ p& v/ Y# Fdef time_calc(seq,state,isEmpty,rgv,currP,total):                                                 # 主要用于记录时间
# Y* _1 G. x% e7 ^; e+ K/ u        index = 0
: s) a2 H6 Z6 T4 E, k, h        temp = 0
7 m- I$ G. |( V' @/ W  R& ]        while index<len(seq):
                nextP = seq[index]
                t = tm[currP][nextP]
                total += t
" V. |3 m" o; d- x0 B' P                update(state,t)
* z. `- c* I# X, Z                if Type[nextP]==0:                                                                                                 # 如果下一个位置是第一道工作点
$ {( \# }, H; W4 P- z                        if rgv==1:                                                                                                         # 然而载着半成品
1 s9 s: k$ t, \$ T( o0 q                                seq.pop(index)                                                                                         # 去掉这个元素并中止当次循环进入下一个循环
                                continue                               
                        if isEmpty[nextP]:                                                                                         # 如果下一个位置是空的
) \( p% Y. }! b6 j4 G                                t = cncT[nextP]
                                total += t
                                update(state,t)
1 `5 c- Q, I1 t& t# S7 D                                state[nextP] = T1                                                                                 # 更新当前的CNC状态
+ n, Y* v+ G3 p9 A8 y$ J6 T: Y& V                                isEmpty[nextP] = 0                                                                                 # 就不空闲了
                        else:                                                                                                                 # 如果没有空闲
                                if state[nextP] > 0:                                                                         # 如果还在工作就等待结束
                                        t = state[nextP]
, }3 d" _; M8 m: `, V$ x                                        total += t
  Q- Q' P& n  p7 P                                        update(state,t)
( N/ e+ y  X) w! v, x                                t = cncT[nextP]                                                                                         # 完成一次上下料
                                total += t
                                update(state,t)
- r  v  \1 N$ D: i" c6 D                                state[nextP] = T1
                                rgv = 1
                else:                                                                                                                         # 如果下一个位置是第二道工作点
! G! S5 c1 P7 J  `+ h5 \) V                        if rgv==0:                                                                                                         # 如果是个空车
                                seq.pop(index)                                                                                         # 删除当前节点
& A& o" v4 J! w                                continue
& r/ F2 \2 H. x( `6 l1 @% Z7 z                        if isEmpty[nextP]:                                                                                         # 如果下一个位置是空的
                                t = cncT[nextP]
                                total += t
5 n- h8 Q/ V1 s7 x3 \" z$ Z! U                                update(state,t)
                                state[nextP] = T2
: c1 ?( x  q3 c+ v                                isEmpty[nextP] = 0       
                        else:                                                                                                                 # 如果没有空闲
                                if state[nextP] > 0:                                                                         # 如果还在工作就等待结束
4 e/ q* a: t% p3 L$ R( C                                        t = state[nextP]
1 ?9 m/ n' ~. X9 C                                        total += t
                                        update(state,t)
                                t = cncT[nextP]+Tc
8 ?. i" N1 q' H  l' j( t                                total += t
                                update(state,t)
                                state[nextP] = T2
                        rgv = 0
                currP = nextP
                temp = total
3 h1 f7 _) d, ^  p( N! e# c/ t                index += 1       
7 b0 U* T9 v! K( k6 n0 G+ J- F        return rgv,currP,total

def forward1(state,isEmpty,currP):                                                                                 # 一步最优
        lists = []
        if currP in A:
% N5 e9 S+ q8 Q$ W: I                rgv = 1
5 D0 s$ Q! q" F$ z% n8 W' {                for e1 in B:
; J! u9 `# `& D9 g1 E! y                        lists.append([e1])
       
" c+ i; f% D) s* e        else:
                rgv = 0
% y; y1 @. |. V7 V/ `3 q' B" n: D                for e1 in A:
# t" K% M0 J# \1 ^6 N                        lists.append([e1])
% s! Z' w6 K: g7 H       
0 q% ?& z5 K& G8 d        minV = 28800
        for i in range(len(lists)):
- C& w6 O4 z3 B1 q, n1 f: r                t = time_calc(lists,state[:],isEmpty[:],rgv,currP,0)[-1]
. ?' u' u- u; p2 N$ }$ E# o3 K                if t<minV:
                        minV = t
5 d2 i8 @5 T2 |$ z2 E                        index = i
  y( g4 n6 ~0 s4 x$ J        return lists[index][0]
- T5 j) ?4 O0 G% i: p+ F. }
def forward4(state,isEmpty,currP):                                                                                 # 四步最优
& M& X  g. d& j. S  S9 @        lists = []
7 B- m9 X+ E: v4 }* C* i        """ 遍历所有的可能性 """
0 r# l9 J7 D$ v. J        if currP in A:                                                                                                                 # 如果当前在第二道工序CNC的位置
: W* z2 E5 v) L                rgv = 1
) g& Z+ ^8 I+ A: j9 v2 S8 n                for e1 in B:
                        for e2 in A:
                                for e3 in B:
5 L9 U: Y+ @/ K+ {- D- T7 H  |                                        for e4 in A:
, |1 _6 E7 V9 s" T. C7 F# E. S5 K                                                lists.append([e1,e2,e3,e4])
        else:
                rgv = 0
7 x% J* p9 r( O. m" q8 t0 q* L# m                for e1 in A:
( V- o# \2 j) R                        for e2 in B:
                                for e3 in A:
7 y, v: ]6 a6 c. ?                                        for e4 in B:
                                                lists.append([e1,e2,e3,e4])
        minV = 28800
. D1 A3 G3 r  j8 Y$ W, ]        for i in range(len(lists)):
                t = time_calc(lists,state[:],isEmpty[:],rgv,currP,0)[-1]
1 ]. m% r% C% {0 F- o( X! X                if t<minV:
* {. F; |9 A) G5 g4 ~                        minV = t
                        index = i
        return lists[index][0]                                                                                                 # 给定下一步的4步计算最优

def forward5(state,isEmpty,currP):                                                                                 # 五步最优
        lists = []
; @8 [+ C* P% d5 B        """ 遍历所有的可能性 """
        if currP in A:                                                                                                                 # 如果当前在第二道工序CNC的位置
9 Z1 w! S. }% Y/ \6 l3 M0 L                rgv = 1
" @2 y2 ^% R* k* l( V) t                for e1 in B:
                        for e2 in A:
                                for e3 in B:
; F( z& x& ^1 _* |) F  N+ ?4 Q                                        for e4 in A:
* q; O+ ]* M2 _+ d. _. Y5 S5 t                                                for e5 in B:
9 c7 s3 B9 m# v; B$ D7 i                                                        lists.append([e1,e2,e3,e4,e5])
        else:
                rgv = 0
. p; h/ J$ ]* ~- f                for e1 in A:
3 ^4 f( t5 T! g' }: U# A1 {                        for e2 in B:
* U4 e! k' I9 A" b" I. S9 g; Z                                for e3 in A:
                                        for e4 in B:
3 z/ D) v% X+ w, `) W' u                                                for e5 in A:
" ^& w! J+ @9 H                                                        lists.append([e1,e2,e3,e4,e5])
        minV = 28800
        for i in range(len(lists)):
                t = time_calc(lists,state[:],isEmpty[:],rgv,currP,0)[-1]
5 _7 X0 R- I3 `6 S9 J& W  [* z                if t<minV:
  V2 q5 z7 I$ H                        minV = t
                        index = i
( n" z. D' A  c3 D        return lists[index][0]                                                                                                 # 给定下一步的5步计算最优
4 d  p# s" Q* [" ]2 `% o
# r) n# M/ [9 |0 I) _def forward6(state,isEmpty,currP):                                                                                 # 六步最优
        lists = []
        """ 遍历所有的可能性 """
        if currP in A:                                                                                                                 # 如果当前在第二道工序CNC的位置
                rgv = 1
                for e1 in B:
( m, C% C* a1 O) T                        for e2 in A:
) K5 j+ q9 e; o3 F; d                                for e3 in B:
                                        for e4 in A:
9 D; S3 ^7 \- J! m                                                for e5 in B:
                                                        for e6 in A:
* F9 G; y) R+ g& B                                                                lists.append([e1,e2,e3,e4,e5,e6])
        else:
, ]4 `3 T  X/ a, C7 N  R                rgv = 0
3 Q2 e2 x6 v& @6 i4 t$ a, R1 E; `                for e1 in A:
                        for e2 in B:
9 {7 a/ [: a/ X$ g                                for e3 in A:
                                        for e4 in B:
                                                for e5 in A:
                                                        for e6 in B:
2 R  d; h. \- K8 m0 w5 m0 q                                                                lists.append([e1,e2,e3,e4,e5,e6])
/ e1 }  D. W2 e2 i1 J        minV = 28800
0 k% `- G" h4 a$ d& n1 |: z) b        for i in range(len(lists)):
" d* L$ D2 I; }8 Z3 D                t = time_calc(lists,state[:],isEmpty[:],rgv,currP,0)[-1]
2 r3 B) u8 t! k& n                if t<minV:
3 L, [( I! H- l% Y) Z" i: H; J6 m                        minV = t
                        index = i
3 J9 I/ t( k7 H5 K5 D2 W/ k: |5 k        return lists[index][0]                                                                                                 # 给定下一步的6步计算最优
/ ]3 ^8 t% e  I# O
def forward7(state,isEmpty,currP):                                                                                 # 七步最优
7 e, w" z) I3 U: W4 ?) W. u        lists = []
        """ 遍历所有的可能性 """
: \3 n5 ]2 F# m$ ^3 W+ w9 o        if currP in A:                                                                                                                 # 如果当前在第二道工序CNC的位置
+ f8 d9 I1 T& L9 o% s                rgv = 1
                for e1 in B:
                        for e2 in A:
8 A7 S/ V! c0 @                                for e3 in B:
                                        for e4 in A:
                                                for e5 in B:
/ g5 T; C6 j  A# j                                                        for e6 in A:
! b6 M9 m6 v1 N" H                                                                for e7 in B:
2 d2 o" m9 x, X* A- K                                                                        lists.append([e1,e2,e3,e4,e5,e6,e7])
        else:
                rgv = 0
                for e1 in A:
                        for e2 in B:
; M1 c  R8 @" F                                for e3 in A:
                                        for e4 in B:
                                                for e5 in A:
. F4 @8 f8 v% i: H7 W8 ?- x5 ?                                                        for e6 in B:
                                                                for e7 in A:
" v, ?" W! ?4 b" B: r                                                                        lists.append([e1,e2,e3,e4,e5,e6,e7])
+ j! |5 |& B& i        minV = 28800
        for i in range(len(lists)):
                t = time_calc(lists,state[:],isEmpty[:],rgv,currP,0)[-1]
' r0 l# k! @  i: d5 C5 z$ k1 \+ Q                if t<minV:
                        minV = t
                        index = i
! }. [2 h; w. U        return lists[index][0]                                                                                                 # 给定下一步的7步计算最优
5 e- n, o. Q4 U) ^! U  m: V  g- {
" |9 @; Z8 p1 w3 ~' _  v; V+ R/ |def forward8(state,isEmpty,currP):                                                                                 # 八步最优
        lists = []
( h5 Y9 P( {2 U* v5 o; q        """ 遍历所有的可能性 """
$ p) x% O5 O) G6 {7 |# r" H7 h' `        if currP in A:                                                                                                                 # 如果当前在第二道工序CNC的位置
                rgv = 1
                for e1 in B:
                        for e2 in A:
                                for e3 in B:
3 t' @4 K; d. X% h                                        for e4 in A:
                                                for e5 in B:
, N) V0 w1 D8 Y6 A                                                        for e6 in A:
                                                                for e7 in B:
5 O- [0 N( |% ~# G                                                                        for e8 in A:
7 h7 @4 f) z" s( k3 L( @7 I1 R5 c                                                                                lists.append([e1,e2,e3,e4,e5,e6,e7,e8])
        else:
                rgv = 0
1 f4 ^  [  g: n& k                for e1 in A:
                        for e2 in B:
                                for e3 in A:
                                        for e4 in B:
                                                for e5 in A:
& p+ w3 Z* ^' N8 F; D                                                        for e6 in B:
$ q6 F6 z6 M' h; m& I  v. ?                                                                for e7 in A:
                                                                        for e8 in B:
& N+ j0 g8 L8 U% w( c$ j                                                                                lists.append([e1,e2,e3,e4,e5,e6,e7,e8])
! ^. W8 V9 `2 K" I        minV = 28800
        for i in range(len(lists)):
                t = time_calc(lists,state[:],isEmpty[:],rgv,currP,0)[-1]
                if t<minV:
                        minV = t
                        index = i
* i5 w* g+ T$ I        return lists[index][0]                                                                                                 # 给定下一步的8步计算最优

+ x( f  o+ i+ T2 s0 ~def greedy(state,isEmpty,rgv,currP,total):                                                                 # 贪婪算法
  y$ b  K2 g, U* n6 b1 E+ p* j        line = []
# k; N$ @, S1 d5 q, X* ^        count = 0
) s) ?$ o0 k4 n/ }) P9 A        while True:
; m8 J# \  u7 E                #nextP = forward4(state[:],isEmpty[:],currP)               
                nextP = forward5(state[:],isEmpty[:],currP)               
6 ]" r$ S/ s$ F* a, }: E3 e' s                line.append(nextP)
                rgv,currP,t = time_calc([nextP],state,isEmpty,rgv,currP,0)
) h- S+ z9 m2 P* ?* w4 t                total += t
                count += 1
                if total>=28800:
                        break
2 X* T- o3 S$ P" @: |        return line
; w& ?2 m" f$ k
+ @5 A' v9 X' S* _. G1 R- Aif __name__ == "__main__":
        state,isEmpty,log,count1,rgv,currP,total,seq = init_first_round()
        print(state,isEmpty,log,count1,rgv,currP,total,seq)
5 q# n" B6 T/ R0 l9 q2 w; d$ `; d        line = greedy(state[:],isEmpty[:],rgv,currP,total)
        simulate(line,state,isEmpty,log,count1,rgv,currP,total)
       
        write_xlsx()
! I) C# O- a6 C后记

这次博客有点赶，所以质量有点差，很多点没有具体说清楚。主要最近事情比较多。本来也没想写这篇博客，但是觉得人还是要善始善终，虽然没有人来阅读，但是学习的路上还是要多做小结，另外也是万一有需要的朋友也可以给一些参考。虽然我的水平很差劲，但是我希望能够通过交流学习提高更多人包括我自己的水平。不喜勿喷！
4 }8 U2 {7 p  o' G. p& o---------------------

+ ~- o( C- f( k7 G

+ i9 {( n  w8 [5 P$ X  M
" \5 c4 u* e& j: @5 o& J. P
% }& I6 `3 W2 E/ [  a, ?8 Y
/ ^0 g: `4 T- o

; k* O$ O5 C% Z' `/ T+ g
  }: Y9 H$ h. G# @" |* f