生产调度问题及其优化算法

1、个人资料整理仅限学习使用生产调度问题及其优化算法<采用遗传算法与MATLAB 编程）信息 014孙卓明二零零三年八月十四日-0-/14个人资料整理仅限学习使用生产调度问题及其优化算法背景及摘要这是一个典型的 Job-Shop动态排序问题。目前调度问题的理论研究成果主要集中在以Job-Shop问题为代表的基于最小化完工时间的调度问题上。一个复杂的制造系统不仅可能涉及到成千上万道车间调度工序，而且工序的变更又可能导致相当大的调度规模。解空间容量巨大， N个工件、 M 台机器的问题包含 种排列。由于问题的连环嵌套性，使得用图解方法也变得不切实际。传统的运筹学方法，即便在单目标优化的静态调度问题

2、中也难以有效应用。b5E2RGbCAP本文给出三个模型。首先通过贪婪法手工求得本问题最优解，既而通过编解码程序随机模拟优化方案得出最优解。最后采用现代进化算法中有代表性发展优势的遗传算法。文章有针对性地选取遗传算法关键环节的适宜方法，采用MATLAB软件实现算法模拟，得出优化方案，并与计算机随机模拟结果加以比较显示出遗传算法之优化效果。对车间调度系列问题的有效解决具有一定参考和借鉴价值。 p1EanqFDPw一问题重述某重型机械厂产品都是单件性的，其中有一车间共有 A ， B，C，D四种不同设备，现接受 6件产品的加工任务，每件产品接受的程序在指定的设备上加工，其工序与加工周期如下表： <

3、;S-设备号、 T-周期） DXDiTa9E3d工序12345678产品STSTSTSTSTSTSTST1C8A2B4C24D62A4D 5B 3C 43C3D7A15B20A84B 7C6D21A1D16C35D10B4C 8D 4A12C6D 16A 1B 4A 7C 3D 5A 2C 5A 8(表一 >条件： 1、每件产品必须按规定的工序加工，不得颠倒；2、每台设备在同一时间只能担任一项任务。<每件产品的每个工序为一个任务）问题：做出生产安排，希望在尽可能短的时间里，完成所接受的全部任务。要求：给出每台设备承担任务的时间表。注：在上面，机器 A， B， C，D即为机器 1，

4、2， 3， 4，程序中以数字 1， 2， 3，4表示，说明时则用 A， B， C， D-1-/14个人资料整理仅限学习使用二模型假设1每一时刻，每台机器只能加工一个工件，且每个工件只能被一台机器所加工 ，同时加工过程为不间断；2所有机器均同时开工，且工件从机器I到机器 J 的转移过程时间损耗不计；3各工件必须按工艺路线以指定的次序在机器上加工多次；4操作允许等待，即前一操作未完成，则后面的操作需要等待，可用资源有限。三符号说明及初始数据表达分析- 第i 个工件 <i=1 6）- 机器顺序阵表示 i 工件的第 j 个操作的机器号-第j 台机器<j=14）-工件排列阵表示 i 机器上第

5、 j 次加工的工件号-加工时间阵为i 工件的第 j 个操作的时间周期- 整个任务完成时间整理数据后得到：= CABCD000=8 2424 6000RTCrpUDGiT ADBC0000 453400005PCzVD7HxA CDABA000 3715208000jLBHrnAILg BCDADC00 76 21116300xHAQX74J0X DBCDACD0 1048412610LDAYtRyKfE ABACDACA 14735258Zzz6ZB2Ltk上述二阵直接从题目得出，而则是我们要求的。关于工件的加工时间表 ：( 表二 >产品 / 工件 <i ） :123456总计总净

6、加工时间 <周期）441653544535247加工工序总数 <个）54567835关于机器的加工时间表：(表三>机器 / 设备 (j>:ABCD总计总净加工时间60427075247加工操作次数10610935分析：由于各产品总净加工时间和各机器总净加工时间之中最大值为75 ，而总计为 247，-2-/14个人资料整理仅限学习使用那么 总时间 C 介于 75 ，247 。同时各工件加工繁杂程度不一，各机器的任务量也有轻重之别。合理的调度排序是对于节省时间和资源是必要的。dvzfvkwMI1希望最优化答案是75，这样达到最小值，如果答案是75，那么意味着机器D不间断工作

7、 , 直至全部加工任务完成。rqyn14ZNXI四贪婪法快速求解AB如果按照一定规则排序，当多个工件出现“抢占”同一机器的局面的时候 , 我们可以制定如下的工序安排规则 :1. 优先选择总剩余时间或总剩余操作较多的工件。 <如果出现总剩余加工时间多者总剩余操作数反而较少的情况时, 按照程度具体情况具体分析）。EmxvxOtOco2. 机器方面来说 , 尽量避免等待空闲时间，优先考虑剩余净加工时间或者剩余加工总次数较多的机器 , 尤其是机器 D,即倘若能够使机器 D不间断工作且其他机器完工时间均不多余 75时, 那么就可以得到最优解。 SixE2yXPq5首先按照最优化时间为 75的设想避

8、免 D出现等待，排序后得到升以下具体排列顺序。各机器承担任务表为 ( 其中粗体字为对应工件产品号，括号内为对应时间周期段>：操作 1操作 2操作 3操作 4操作 5操作 6操作 7操作 8操作 9操作 106216345636(1>(2-5>(12-13>(14-20>(21-35>(36>(43-54>(55-56>(57-64>(66-73>465132(1-7>(8-11>(12-15>(16-19>(36-55>(56-58>C3145615624(1-3>(4-11>(1

9、2-17>(18-25>(26-28>(29-52>(55-60>(61-65>(66-69>(70-72>D534562415(1-10>(11-17>(18-38>(39-42>(43-47>(48-52>(53-68>(69-74>(75>(表四>A 机器1 4 27151 12288B 机器744 4203C机器38683246543D 机器10721455166101020304050607080-3-/14个人资料整理仅限学习使用( 图一>上图为加工周期图 <甘特

10、图），标注数字为相应操作的周期，完工时间为第75 周期。五计算机随机模拟 <编程）1编码：随机产生生产的工序操作优先顺序，进行编码，如：K=43566231 45<注：同时作为下文的染色体之用）意思为：工件 4优先被考虑进行第一次操作，然后 3进行其第一步操作，然后 5操作， 6操作，再 6操作其第二步工序，依次进行。如果前后互相不冲突，则可同时在不同机器上操作。6ewMyirQFL通过排列组合得出，总共有类似K 的排列序列 2多种！当然，这其中只对应解 75， 247，意味着有大量排列序列对应同一加工方案，而大量加工方案又对应同一时间解。 kavU42VRUs2解码：即对编码进行

11、翻译，产生具体可操作工序安排方案，这里采用活动化解码算法。例如工件 2第i 步操作 <记为<2， i），且在机器 A 上进行）被安排在工件 3第j 步操作 <记为 JM<3,j ）后面进行，那么如果安排好y6v3ALoS89<3,j）后，只要<2， i）在工件 2已经排序好的操作之后进行，那么操作 <2，i ）可插入到机器 A 处最前可安置的时间段进行。在这里，一个编码序列对应一个加工方案，而一个加工方案可对应一个或多个编码序列，这就是二者之关系。3编程：通过一组随机编码产生一生产加工优先序列，通过解码过程产生相应加工方案及其总耗费时间 C .N次模拟

12、后即可得出解 C的概率密度分布情况以及相对最优解 <N个C的最小值，如 80，77等，甚至出现 75）。 M2ub6vSTnP4计算机模拟所得数据分析a. 进行 100次模拟得出最优解情况：(共运行 10次 >82，83，82， 84，78， 80，81，83， 87，82 平均值 82.2，每回耗时约 3秒 b. 进行 1000次模拟得出最优解情况： (共运行 10次 >-4-/14个人资料整理仅限学习使用80，79，78， 78，79， 79，76，80， 77，78 平均值 78.4 , 每回耗时 25秒 c. 进行 10000次模拟得出最优解情况： (共运行 10次&

13、gt;76，77，77， 75，76， 76，77，76， 76，77 平均值 76.3, 每回耗时 4分钟 d. 模拟 1000000次得到的解 C的概率密度分布情况为： <如图二所示）(图二 ><注： 75处为 17次<概率为 17/1000000=1/58823 ）， 76处为 40次， 77处 167次 ）结论：如果想将 2中排序序列以平均出现一次的可能性进行模拟，即运行 2 次，计算机需运行将近 150万亿年！当然，我们没有这个必要，因为我们只需要运行数万次，就很可能得到最优解 75，<在随机模拟 1000000次后出现 17次75，那么意味着概率大约

14、17/1000000=1/58823 ，另外 76处为 40次， 77处167次）。 0YujCfmUCw六遗传算法模型建立和步骤解法遗传算法 <GeneticAlgorithm ）作为一种优化算法特别适合于对象模型难于建立、搜索空间非常庞大的复杂问题的优化求解。它和模糊控制技术一样，虽然在理论上还没有完善，但是在实践中已经得到了广泛的应用。遗传算法的基本思想是：模仿生物系统 “适者生成 " 的原理，通过选择、复制、交叉、变异等简单操作的多次重复来达到去劣存优的目的，从而获得问题的优化结果。遗传算法的实现由两个部分组成，一是编码与解码，二是遗传操作。其中遗传操作又包括选择、复制

15、、交叉、变异等步骤。eUts8ZQVRd本文根据实际情况采取了1-6 整数编码。数字1，2，3，4，5，6 分别代表6 件待加工产品。本文遗传算法基本流程:通过编码，解码程序随机产生N 个 <有一定数量 , 如 50 或 100）个体构成初始种群a)从初始中群中选取 2 个具有最优染色体 <最有排序方案）的个体作为临时个体 <父代）；-5-/14个人资料整理仅限学习使用b)如果此 2 个体中有一个个体通过解码操作能够实现最优排序<即使总时间为 75 周期），那么结束此算法，得到最优解；sQsAEJkW5Tc)对 2 个临时个体以一定方式 ( 循环交叉 >执行染色体

16、交叉变换和变异选择<小概率 , 互换操作），产生2 个新的个体； GMsIasNXkAd)对父代和子代共4 个个体进行选择，从中选出最佳的2 个个体，做为下一代的父代；e) 重复执行第二步 (b> 操作；f) 如果执行完 M步后仍然未得出答案 75，那么将目前的最优解作为本算法的最优解答案。1编码和解码 < 同上）2交叉变换 (crossover>对 2 个父代临时个体进行染色体交叉变换，采用循环交叉方法<Cycle crossoverCX），如父代染色体为：X：9 2 6 4 7 3 5 8 1和 Y：3 4 5 8 1 6 7 2 9，如果随机选到第二位开始交

17、叉，那么X的 2对应 Y的 4，X的 4对应 Y的 8，X的 8对应 Y的 2 ，这样就确定了以上为不变的染色体，其余位置的染色体互换位置，最后得到:3 2 5 4 1 6 7 8 9,:9 4 6 8 7 3 5 2 1,实现交叉变换。TIrRGchYzg3变异选择 <mutation ）采用 互换操作 <SWAP），即随机交换染色体中两不同基因的位置。如上面的染色体为：:3 2 5 4 1 6 7 8 9。随机产生变换位置号，如产生随机数3 和 5，那么交换数字后得到染色体：3 2 1 4 5 6 7 8 9， 变异概率取0.1。 7EqZcWLZNX4选择操作 <sel

18、ection）对父代 2 个体 f1,f2和子代 2 个体 f3,f4进行选择，通过编码操作确定具有最优解的2 个个体，成为新一代f1 和 f2。 lzq7IGf02E如此，通过多次编码和解码随机产生一定数量的个体，选取2 个最佳个体进行交叉变换操作，产生 2 个新个体，然后对 4 个个体进行选择，产生下一代，如果某时刻通过解码操作得出最优解 <所有解的下限，这里是 75 周期），那么算法结束，否则循环进行，直至预先给定的循环次数达到为止，以最后得到的最优解作为最终最优解。zvpgeqJ1hk七遗传算法模拟 <采用 MATLAB工具编程）主程序如下： <子程序见略）% 本程序

19、为主程序，调用以下各分支程序task= 'Welcome! Wait a moment please! - Writer:孙卓明, 信息 014',NrpoJac3v1f1=zeros(1,35>。 f2=zeros(1,35> 。while f1=f2。%此步避免初始染色体 f1,f2相同，导致以下死循环1nowfTG4KI(N>。 %minminmax1,s1=chushijie种群初始化；基于操作的编码策略；活动化解码算法。chushijie(N>参数 N 为初始种群数 fjnFLDa5Zof1=s1 。 minminmax1,%选取的第一个初始个

20、体-6-/14个人资料整理仅限学习使用minminmax2,s2=chushijie(N>再次种群初始化。%f2=s2。 minminmax2,%选取的第二个初始个体end 。Mfor e=1:% e=1:M进行 M次遗传操作<交叉-变异-选择>。tfnNhnE6e5D=jiaocha(f1,f2>。%交叉变化 <循环交叉操作， cycle crossover CX）, 选取HbmVN777sLf1 。 f2 。“染色体”无需变动部分基因f3,f4=jiaocha_bianyi(f1,f2,D>。%生成交叉后的“染色体”, 并进行变异选择V7l4jRB8Hs

21、f3。 f4 。f1,f2=xuanze(f1,f2,f3,f4>。%选择：对父代f1,f2和子代f3,f4进行解码，得出2 个83lcPA59W9f1 。 f2 。较优个体，成为下一代的父代mZkklkzaaPminmaxf1,a1,b1=tongbujinzhan(f1>。% 求 该 时 刻 个 体f1的 最 优 时 间 ( 因 为f1优 于f2> AVktR43bpwif minmaxf1=75。f1,a1,b1,minminmax1,minminmax2,minminmax_last=minmaxf1,ORjBnOwcEdtask='Finish! Succe

22、ssful!Best answer!Congratulation! ' , return。2MiJTy0dTTend 。end。f1,a1,b1,minminmax1,minminmax2,minminmax_last=minmaxf1,gIiSpiue7Aif minminmax_last>=90。 task='Finish! Action again please!',end。uEh0U1Yfmhif minminmax_last>=80&&minminmax_last<90。 task='Finish!' , en

23、d。IAg9qLsgBXif minminmax_last<80。task='Finish! Successful!', end。WwghWvVhPE八遗传算法模拟结果首先给出最优方案 :在进行某次 n=100,m=200 的操作后得到模拟最优结果75 周期时间：minminmax1 =83minminmax2 =78<二个初始较优个体解）f1 =4 5 6 6 3 1 3 6 4 5 6 1 3 2 5 4 5 3 1 5 2 6 4 5 6 4 6 6 4 3 2 2 5 11 asfpsfpi4k<f1为各工件优先选择顺序排列, 即“染色体”）a1= 5

24、35396400000<a1,b1为四台机器空闲周期段）ooeyYZTjj115240000000175365000000000000000b1=113842650000020350000000185468000000000000000minminmax =75<最终最优解 ）其中机器 A 空闲时间段为 :5-11,35-38,39-42,64-65。 机器 B 则为 :15-20,24-35。 机器 C 为 :17-18,53-54,65-68。机器 D 无空闲。 BkeGuInkxI-7-/14个人资料整理仅限学习使用以下为 : 取不同 N 和 M值情况下数据优化过程以及时间上

25、的比较:(表五>NNM第一次运行第二次运行第三次运行第四次运行第五次运行平均运行时间111104-114- 9898-105- 9392-93-92100-132- 9586-86- 84/1110106-97- 86108-100- 89102-87- 8799-90- 9088-104- 84/1110094-81- 8181-102- 7891-105- 91101-84- 8090-101- 90/111000107-100- 7892-101- 76101-100 -8288-97- 8691-93- 878.5(s>111000088-105- 77103-81- 77

26、89-99- 84107-104- 7893-105- 7880(s>10101090-108- 90104-91- 83104-100- 9395-98- 87105-106- 87/101010098-96- 9693-99- 9088-90-80105-92- 8091-95- 85/10101000101-96- 7891-89- 8096-104- 87105-105- 8488-99- 789.5(s>10101000099-92- 7797-95- 7596-104- 7689-99- 7691-101- 7590(s>10010010095-100- 8698

27、-90- 80104-99- 7893-88 -8192-99- 80/1001001000109-98- 8591-100- 82100-99- 77114-101 -8496-110- 7611(s>1001001000096-101- 78101-86- 76101-97- 8099-110- 7699-111- 77100(s>说明 :以最后一行第一次运行“96-101-78 ”为例 ,96 和 101 分别为 2 次 N 取100 时得到的100*2 个随机解中的最优解 ,78 为经过 M=10000代的遗传( 交叉变异选择 >后得到的最终最优解。PgdO0sRlM

28、o明显地发现 : M的增大所产生的优化效果明显好于N 增大产生的优化效果M 从 1-10-100-1000-10000的变化使最优解有从9*-8*-7*的变化规律，N 的 1-10-100 的变化则不明显。因此相对于随机取解 ,经过遗传算法优化之后效果是显著的。另外对采用遗传算法前后模拟所得数据进行比较:第一次第二次第三次第四次第五次均值平均时间N=1000 ,M=0787981807925(S>794N=1,N=1,M=1000807577768278.08.5(S>(表六 >由上表看来 , 虽然在均值方面相差不显著 , 但是时间上采用遗传算法之后节约了约三分之二的运行时间

29、 , 效率显而易见。 3cdXwckm15九模型优缺点及改进模型优点，采用遗传算法可对一个理论上无法求尽的NP 问题以较有效方法在较短时间内得到较精确解。缺点，采用遗传算法还不全面，只是一个简单模型，未考虑收敛性, 适配值等，对参数设置与最优解关系研究还不够深入。h8c52WOngM模型本身还具有漏动，仍需改进，较好设想为采用和模拟退火算法的混合遗传算法，由于时间紧迫，在此就不再深入给出模型及分析了。v4bdyGious-8-/14个人资料整理仅限学习使用参考文献：1车间调度与遗传算法王凌清华大学出版社2数值计算的算法与分析张可村赵英良科学出版社3Permutation Based GAs a

30、nd Ordered GreedPeter G. Anderson,J0bm4qMpJ94MATLAB6.0与科学计算王沫然电子工业出版社5C程序设计 <第二版）潭浩强清华大学出版社信息014孙卓明2003年8月 14日本文网上下载地址 ( 个人主页 >: 附录：子程序一： ( 种群初始化，得较优个while x<=7。体 >Jm(I,x>=Jm(I,x+1>。functionx=x+1。minminmax,ss=chushijie(n>end。% 种群初始化 , 以下为编码和解码过程，同时Jm(I,8>=0。对 n 次选取最优化个体XVauA9

31、grYPend。Jm=3 1 2 3 4 0 0 0。14230000。34Jm=3 1 2 3 4 0 0 0。14230000。 3 4121000。23414300。423413121000。23414300。42341340 。1213 4131。bR9C6TJscw40 。1213 4131。 DJ8T7nHuGTminminmax=200。T=824246000。45340000。 3 7for d=1:n 。15208000。7621116300。1048412610。147352 5s=0。8。%编 码程 序： 基于 操作 的编 码策略%解码程序：活动化解码算法pN9LBDdt

32、rdQF81D7bvUAk=1 。for i=1:6。for t=1:35。k(i>=1。I=randint(1,1,1,6>。end 。while Jm(I,1>=0。for q=1:4。I=randint(1,1,1,6>。for x=1:9。end。a(q,x>=0。 b(q,x>=0 。flag_(q>=0 。s(k>=I。end 。k=k+1。end 。x=1。for p=1:6。-9-/14个人资料整理仅限学习使用flag(p>=0。flag(s(i>>=flag(s(i>>+t。end 。flag_(q

33、>=flag(s(i>>。for i=1:35。end 。q=Jm(s(i>,k(s(i>>>。t=T(s(i>,k(s(i>>>。end 。z=0。 v=0 。k(s(i>>=k(s(i>>+1。for x=1:9。end。ifminmax=0。max(flag(s(i>>,a(q,x>>+t<=b(q,x>&&z=0。for q=1:4。ifif minmax<flag_(q>。flag(s(i>><=a(q,x>

34、&&a(q,x>+t=b(q,x>。minmax=flag_(q>。flag(s(i>>=b(q,x>。forend 。y=x:8 。a(q,y>=a(q,y+1>。end。b(q,y>=b(q,y+1>。 4B7a9QFw9hifend。minminmax>minmax。z=1 。% 得出初始最优解 Yl4HdOAA61elseifminminmax=minmax。 ss=s。flag(s(i>><=a(q,x>&&a(q,x>+t<b(q,x>。end

35、 。a(q,x>=a(q,x>+t。 flag(s(i>>=a(q,x>。 z=1 。end。ix6iFA8xoXelseifflag(s(i>>>a(q,x>&&a(q,x>+t=b(q,x>。flag(s(i>>=b(q,x>。 b(q,x>=b(q,x>-t。 z=1 。子程序二： <父代交叉变换）wt6qbkCyDEfunction D=jiaocha(f1,f2>%交叉变elseif化 <循环交叉操作，CX） , 选取“染色体”无需变flag(s(i>

36、;>>a(q,x>&&a(q,x>+t<b(q,x>。动部分基因 ch4PJx4BlIfor y=8:x+2。s=randint(1,1,1,35>。a(q,y>=a(q,y-1>。b(q,y>=b(q,y-1>。while f1(s>=f2(s>。end。s=randint(1,1,1,35>。b(q,x+1>=b(q,x>。b(q,x>=flag(s(i>>。end。a(q,x+1>=flag(s(i>>+t。for p=1:34。flag(s

37、(i>>=a(q,x+1>。z=1。Kp5zH46zRkD(p>=0。end。end 。end 。z=0 。 j=1 。D(j>=s 。j=j+1。end 。for i=1:35。if z=0。if f1(i>=f2(s>&&z=0。if flag(s(i>><=flag_(q>。D(j>=i。j=j+1 。z=1 。flag(s(i>>=flag_(q>+t。end。flag_(q>=flag_(q>+t。end 。elseif flag(s(i>>>fla

38、g_(q>。if f2(D(j-1>>=f1(s>。for x=1:9。return。if a(q,x>=0&&v=0。end 。a(q,x>=flag_(q>。for m=1:34。b(q,x>=flag(s(i>>。 v=1。z=0。end。end。for i=1:35。10/14个人资料整理仅限学习使用if f1(i>=f2(D(j-1>> && z=0。end。w=0。 for t=3:j。if (i-D(t-1>>>0|(i-D(t-1>><

39、;0。w=w+1。 end 。end。子程序四：<四者中选取最优二个if w=j-2。体）D(j>=i。 j=j+1 。 z=1。functionend。f1,f2=xuanze(f1,f2,f3,f4>%end。对父代 f1,f2和子代 f3,f4进行解码，得出 2个end 。较优个体，成为下一代的父代E836L11DO5if f2(D(j-1>>=f1(s>&&z=1。min1=0。min2=0。min3=0。min4=0。h=0。 g=0。return。min1=tongbujinzhan(f1>。end。min2=tongbuj

40、inzhan(f2>。end 。min3=tongbujinzhan(f3>。end。min4=tongbujinzhan(f4>。if min1<=min2&&min1<=min3&&min1<=min4。h=f1。ifmin2<=min3&&min2<=min4。g=f2 。子程序三： <变异选择，形成子代）elseif min3<=min2&&min3<=min4。 g=f3 。functionelseif min4<=min2&&min

41、4<=min3。 g=f4 。f3,f4=jiaocha_bianyi(f1,f2,D>%end。生成交叉后的“染色体”, 并进行变异选择elseifqd3YfhxCzomin2<=min1&&min2<=min3&&min2<=min4。g2=f1 。 g1=f2 。h=f2。ifmin1<=min3&&min1<=min4。for i=1:34。g=f1 。if D(i>>0。elseifmin3<=min1&&min3<=min4。g1(D(i>>

42、=f1(D(i>>。g2(D(i>>=f2(D(i>> 。g=f3 。end 。elseifmin4<=min1&&min4<=min3。end。g=f4 。f3=g1 。 f4=g2 。end。c=randint(1,1,1,100>。elseifif c=1。min3<=min1&&min3<=min2&&min3<=min4 。d1=randint(1,1,1,35>。h=f3。ifmin1<=min2&&min1<=min4。d2=r

43、andint(1,1,1,35>。g=f1 。while d1=d2。elseifmin2<=min1&&min2<=min4。d2=randint(1,1,1,35>。g=f2 。end 。elseifmin4<=min1&&min4<=min2。m=f3(d1>。f3(d1>=f3(d2>。 f3(d2>=m 。g=f4 。elseif c=2。end。d1=randint(1,1,1,35>。elseifd2=randint(1,1,1,35>。min4<=min1&&a

44、mp;min4<=min2&&min4<=min3。while d1=d2。h=f4。ifmin1<=min2&&min1<=min3。d2=randint(1,1,1,35>。g=f1 。end 。 m=f4(d1> 。f4(d1>=f4(d2> 。 f4(d2>=m。elseifmin2<=min1&&min2<=min3。11/14个人资料整理仅限学习使用g=f2 。ifelseifmin3<=min1&&min3<=min2。max(flag(s

45、(i>>,a(q,x>>+t<=b(q,x>&&z=0。g=f3 。ifend。flag(s(i>><=a(q,x>&&a(q,x>+t=b(q,x>。end 。flag(s(i>>=b(q,x>。 for y=x:8。end。a(q,y>=a(q,y+1>。b(q,y>=b(q,y+1>。f1=h 。f2=g 。LOZMkIqI0wwhile f1=f2。end。 z=1。minminmax3,s1=chushijie(30>。elseiff2=s1 。flag(s(i>><=a(q,x>&&