数学建模指派问题论文

1、 目录一问题重述 . 2二模型假设 . 2三匈牙利法陈述 . 2四问题分析 . 3五问题实现 . 51问题重述 . 52 问题求解 . 52.1由匈牙利法构造目标函数. 52.2模型建立 . 63 模型解析 . 64 程序实现 . 7六结果显示及min求解 . 17七模型深入 . 171 模型建立 . 182 进行求解 . 183程序分析 . 19八模型检验 . 19九整体总结 . 20十参考文献 . 20一问题重述指派问题亦称平衡指派问题仅研究人数与事数相等、一人一事及一事一人的情形。现有的不平衡指派问题将研究范围扩大到人数与事数可以不等、一人一事或一人多事及一事一人的情形。日常活动中也不乏

2、人数与事数可以不等、一人多事及一事多人的情形，这类事务呈现了广义指派问题的实际背景。平衡指派问题是特殊形式的平衡运输问题，可运用匈亚利法、削高排除法和缩阵分析法等特殊方法求解。另一方面，正是平衡指派问题的这种特殊性，使得不平衡指派问题不能按常规技术转化为平衡指派问题。因此，各种不平衡指派问题需要确立相应的有效解法问题的提出及其数学模型广义指派问题并非奇特和抽象的构想，相反，该问题可以从司空见惯的日常事务中引出。现在我们就运用匈牙利法，去实现n个人，n件工作的指派问题。 二模型假设1 假设一共有n个人，n件工作，即人数与工作数相等。2 假设每个人的都能从事某项工作，但是付出的代价不同。3 假设求

3、解代价最小的解。4甲乙丙丁四个人，abcd四项工作，要求每人只能做一项工作，每项工作只由一人完成，问如何指派总时间最短？ 三匈牙利法陈述第一步：找出矩阵每行的最小元素，分别从每行中减去这个最小元素； 第二步：再找去矩阵每列的最小元素，分别从各列减去这个最小元素； 第三步：经过这两步变换后，矩阵的每行每列至少都有了一个零元素，接着根据以下准则进行试指派，找出覆盖上面矩阵中所有零元素至少需要多少条直线；（1）从第一行开始，若该行只有一个零元素打上（）号。对打（）号零元素所在列划一条直线。若该行没有零元素或有两个以上零元素（已划去的不计在内），则转下一行，一直到最后一行为止；（2）从第一列开始，若该

4、列只有一个零元素就对这个零元素打上（）号（同样不考虑已划去的零元素），对打（）号零元素所在行划一条直线。若该列没有零元素或 还有两个以上零元素，则转下一列，并进行到最后一列；（3）重复（1）、（2）两个步骤，可能出现三种情况： 矩阵每行都有一个打（）号零元素，很显然，按照上述步骤得到的打（）的零元素都位于不同行不同列，因此就找到了问题的答案； 有多于两行或两列存在两个以上零元素，即出现了零元素的闭回路，这个时候可顺着闭回路的走向，对每个间隔的零元素打上（）号，然后对所有打（）号零元素或所有列或所在行划一条直线。 矩阵中所有零元素或打上（）号，或被划去，但打（）号零元素个数小于m。 第四步：为了

5、设法使每行都有一个打（）的零元素，就要继续对矩阵进行变换；（1）从矩阵未被直线覆盖的元素找出最小元素k；（2）对矩阵的每行，当该行有直线覆盖时，令ui=0，无直线覆盖的，令ui=k；（3）对矩阵的每列，当该列有直线覆盖时，令vj=-k，无直线覆盖的，令vj=0；（4）得列一个变换后的矩阵，其中每个元素bij=aij-ui-vj。第五步：回到第三步，反复进行，一直到矩阵中每一行都有一个打（）的零元素为止，即找到最优分配方案为止。 四问题分析指派问题的标准形式（以人和事为例）如下。有n个人和n项任务，已知第i个人做第j件事的费用为cij，要求确定人和事之间的一一对应的指派方案，使完成这n项任务的费

6、用最少。一般把目标函数的系数写为矩阵形式，称矩阵c=(cij)nnc11c21.=.n 为系数矩阵（coefficient matrix），也称为效益矩阵或价值矩阵。矩阵的元素cij（i,j=1,2,n）表示分配第i个人去完成第j项任务时的效益。一般地，以xij表示给定的资源分配用于给定活动时的有关效益（时间，费用，价值等），且0,xij=1,不分配第i单位资源用于第分配第i单位资源用于第j项活动j项活动，i,j=1,2,.,n 然后我们求解最小（最大（这里不再讨论）代价和模型，nnijmin（max）z=nci=1j=1ijxiji=1,2,.

7、,n(1)s.t.xj=1ni=1=1,(2)xij=1,j=1,2,.,ni,j=1,2,.,n(3)(4) 当然，作为可行解，矩阵的每列元素中都有且只有一个1，以满足约束条件式（3）。每行元素中也有且只有一个1，以满足约束条件（2）。指派问题n!个可行解。nnijxij=0或1,如果要求解最大值maxz=ci=1j=1)，xij时，我们将构造一个新的矩阵(cij=m-cij，其中m是一个足够大的常数。一般取cij中最大的元素作为m，使cij求解minz=(mi=1j=1nn-cij)xij，所得的解(xij)就是原问题的解。事实上，由nnijijcxi=1j=1=(mi=1nj=1nnn-

8、cij)xijnijn=mxi=1j=1ni=1-i=1cj=1ijxij nij=mn-cj=1xij可的此结论。 五问题实现1问题重述已知问题甲乙丙丁四个人，abcd四项工作，要求每人只能做一项工作，每项工作只由一人完成，问如何指派总时间最短？每个人的对每项工作的代价如下： 2 问题求解开始求解 2.1由匈牙利法构造目标函数 9220ijmin=di=1j=1*xij 引入0-1变量xij，xij =1：第i人做第j项工作xij =0：第i人不做第j项工作 约束条件：xij=10第i个工作从事第j个工作第i个工作不能被从事(i=1,2,92 j=1,2,20)2.2模型建立 即一项任务只由

9、一个人完成 x11+x21+x31+x41=1x12+x22+x32+x42=1x13+x23+x33+x43=1x14+x24+x34+x44=1 一人只能完成一项任务x11+x12+x13+x14=1x21+x22+x23+x24=1x31+x32+x33+x34=1x41+x42+x43+x44=1 求出目标函数minz=3x11+5x12+8x13+4x14+6x21+8x22+5x23+4x24+2x31+5x32+8x33+5x34+9x41+2x42+5x43+2x443 模型解析 根据指派问题的最优性定理，求最优解的问题可以转换为求效益矩阵的 大1元素组的问题。匈牙利法的一般计

10、算步骤为： 步骤1：对效益矩阵进行初等变换，使每行每列都出现0元素。1. 从效益矩阵a中每一行减去该行的最小元素；2. 再在所得矩阵中每一列减去该列的最小元素，得矩阵d；步骤2：将矩阵d中0元素置为1元素，非0元素置为0元素，得矩阵e。 步骤3：确定独立1元素组。1. 在矩阵e中含有1元素的各行中选择1元素最少的行，比较该行中各1元素所在的列中1元素的个数，选择1元素的个数最少的那一列中的1元素；2. 将所选的1元素所在的行和列的元素置为0；3. 重复第2步和第3步，直到没有1元素为止，即得到一个独立1元素组。 步骤4：判断是否为最大独立1元素组。1. 如果所得独立1元素组为原效益矩阵的最大独

11、立1元素组（即1元素的个数等于矩阵的阶数），则已得到最优解，停止计算；2. 如果所得独立1元素组还不是原效益矩阵的最大独立1元素组，那么继续寻找可扩路的方法对其进行扩张，进行下一步； 步骤5：利用寻找可扩路方法确定最大独立1元素组。1. 做最少的直线覆盖矩阵d的所有0元素；2. 在没有被直线覆盖的部分找出最小元素，在没有被直线覆盖的各行减去此最小元素，在没有被直线覆盖的各列上加上此最小元素，得到一个新的矩阵，返回第二步。 4 程序实现这里我们运用c+语言进行编程进行求解 /*/* zhipai2.cppauthor: 路遥date:2012-12-1description:需要在同样的目录下建

12、立一个input.txt的文件夹在里面写入每个人从事不同的工作代价，首行要写入几阶方程，然后是每个人的不同工作代价，用空格隔开。每人一行，见下面数字形式。用匈牙利法，实现n个人，n件工作的指派问题。input:input.txt 格式为代价矩阵维度n，和矩阵内容，例如： 43 5 8 46 8 5 42 5 8 59 2 5 2output:控制台输出指派矩阵，例如0 0 0 10 0 1 01 0 0 00 1 0 0*/*/ #include <iostream>#include <fstream> using namespace std; #define max_

13、n 100intnn;/输入矩阵的维度nn voidprintmatrix(int amax_nmax_n)inti,j;for(i=0;i<nn;i+)for(j=0;j<nn;j+)cout<<aij<<" "cout<<endl; int main()int cmax_nmax_n,bmax_nmax_n;intquanmax_nmax_n,chamax_nmax_n;introwzeromax_n,colzeromax_n;introwcheckmax_n,colcheckmax_n;inti,j,k; ifstrea

14、m input("input.txt");if(!input)cout<<"open input file failed."<<endl;system("pause");exit(1); /读取输入文件 input>>nn; for(i=0;i<nn;i+) for(j=0;j<nn;j+) input>>cij; bij=cij; input.close(); /每行减去该行最小 for(i=0;i<nn;i+) int min=bi0; for(j=1;j<n

15、n;j+) if(bij<min) min=bij; for(j=0;j<nn;j+) bij -= min; /每列减去该列最小 for(j=0;j<nn;j+) int min=b0j; for(i=1;i<nn;i+) if(bij<min) min=bij; for(i=0;i<nn;i+) bij -= min; /开始尝试进行试指派assign: /初始化标记数组和计数数组for(i=0;i<nn;i+)rowzeroi=colzeroi=0;rowchecki=colchecki=0;for(j=0;j<nn;j+)quanij=c

16、haij=0; /计算每行每列0元素个数for(i=0;i<nn;i+)for(j=0;j<nn;j+)if(bij=0)rowzeroi+;colzeroj+; bool flag;/直到尽可能多的0元素都被圈出和划掉为止doflag=false;/寻找只有一个0元素的行，加圈划叉for(i=0;i<nn;i+)if(rowzeroi=1)/该行只有一个0元素/cout<<"rowzero found: "<<i<<endl;flag=true; /找到0元素,加圈划叉for(j=0;j<nn;j+)if(bij

17、=0 &&quanij=0 && chaij=0) quanij=1;rowzeroi-;colzeroj-;for(k=0;k<nn;k+)if(bkj=0 &&quankj=0 &&chakj=0) chakj=1; rowzerok-; colzeroj-; break; /break; /寻找只有一个0元素的列，加圈划叉 for(j=0;j<nn;j+) if(colzeroj=1)/该列只有一个0元素 flag=true; /找到0元素,加圈划叉 for(i=0;i<nn;i+) if(bij=0 &a

18、mp;&quanij=0 && chaij=0) quanij=1; rowzeroi-; colzeroj-; for(k=0;k<nn;k+) if(bik=0 &&quanik=0 && chaik=0) chaik=1; rowzeroi-; colzerok-; break; /break; while (flag); /判断是否还有0元素未被标记intzeronotmarked = 0;for(i=0;i<nn;i+)zeronotmarked += rowzeroi;while(zeronotmarked !=

19、0)/*从剩有0元素最少的行(列)开始，比较这行各0元素所在列中0元素的数目，选择0元素少的那列的这个0元素加圈(表示选择性多的要“礼让”选择性少的)。然后划掉同行同列的其它0元素。可反复进行，直到所有0元素都已圈出和划掉为止。*/ intleastzerorow=0;while(rowzeroleastzerorow=0)leastzerorow+;for(i=leastzerorow+1;i<nn;i+)if(rowzeroi!=0 &&rowzeroi<rowzeroleastzerorow) leastzerorow=i;i=leastzerorow;/得到

20、0元素最少的行标 intleastzerocol=0;while(colzeroleastzerocol=0)leastzerocol+;for(j=0;j<nn;j+)if(bij=0 &&quanij=0 && chaij=0 &&colzeroj<colzeroleastzerocol)leastzerocol=j;j=leastzerocol;/得到0元素最少的列标 /圈出bij,划掉i行和j列其余0元素quanij=1;rowzeroi-;colzeroj-;zeronotmarked-; for(k=0;k<nn;k

21、+)if(bik=0 &&quanik=0 && chaik=0)chaik=1;rowzeroi-;colzerok-;zeronotmarked-; for(k=0;k<nn;k+)if(bkj=0 &&quankj=0 && chakj=0)chakj=1;rowzerok-;colzeroj-;zeronotmarked-; /printmatrix(quan);/若quan元素的数目countquan等于矩阵的阶数nn,则得出最优解，否则画线intcountquan=0;for(i=0;i<nn;i+)for

22、(j=0;j<nn;j+)countquan += quanij;if(countquan<nn)/需要作直线覆盖0元素 /对没有quan的行打勾for(i=0;i<nn;i+)int temp=0;for(k=0;k<nn;k+)temp+=quanik; if(temp=0) rowchecki=1; /* 重复执行(1)(2)，直到打不出新的勾为止： (1)对已打勾的行中含cha的元素的列打勾 (2)对已打勾的列中含quan的元素的行打勾 */ intnewcheck=0; int check; do check=newcheck; /(1) for(i=0;i&

23、lt;nn;i+) if(rowchecki=1) for(j=0;j<nn;j+) if(chaij=1 &&colcheckj=0) colcheckj=1; newcheck+; /(2) for(j=0;j<nn;j+) if(colcheckj=1) for(i=0;i<nn;i+) if(quanij=1 &&rowchecki=0) rowchecki=1; newcheck+; while (check!=newcheck); /(5)对rowcheck=0的行画横线，colcheck=1的列画纵线 /这就得到覆盖所有0元素的最

24、少直线数numlinesintnumlines=0;for(i=0;i<nn;i+)if(rowchecki=0)numlines+;for(j=0;j<nn;j+)if(colcheckj=1)numlines+; if(numlines != countquan)cout<<"直线个数不等于画圈0元素个数，试指派有误"<<endl; /若 numlines = countquan<nn，须再变换当前的系数矩阵，以找到nn个独立的0元素if(numlines=countquan&&numlines<nn)/*

25、变换矩阵bij以增加0元素：设没有被直线覆盖(rowchecki=1 &&colcheckj=0) 的有所有元素中的最小元素为min2，然后打勾各行都减去min2；打勾各列都加上min2，将得到的矩阵重新进行试指派*/i=0;while(rowchecki=0) i+;j=0;while(colcheckj=1) j+; int min2=bij;for(i=0;i<nn;i+)for(j=0;j<nn;j+)if(bij<min2 &&rowchecki=1 &&colcheckj=0) min2=bij; for(i=0;i

26、<nn;i+) if(rowchecki=1) for(j=0;j<nn;j+) bij -= min2; for(j=0;j<nn;j+) if(colcheckj=1) for(i=0;i<nn;i+) bij += min2; /将变换后的矩阵bij重新试指派 goto assign; /此时countquan=nn,输出指派结果 cout<<"指派结果如下："<<endl; printmatrix(quan); system("pause"); return 0;六结果显示及min求解 min（z）

27、=4+5+2+2=13即求出最小代价。完成假设要求 七模型深入将约束条件由整数数据变为小数数据且目标函数由最小值化为最大值问题的求解。假设有4件工作分派给4个人来做，每项工作只能由一人来做，每个人只能做一项工作。下表为各人对各项工作所具有的工作效率。问应该如何安排人选，及发挥个人特长又能使总的效率最大。每个人完成各项任务需要的时间1 模型建立数学模型为maxz=0.6x11+0.2x12+0.3x13+0.1x14+0.7x21+0.4x22+0.3x23+0.2x24+0.8x31+1.0x32+0.7x33+0.3x34+0.7x41+0.7x42+0.5x43+0.4x44x11+x12

28、+x13+x14=1 s.t.x21+x22+x23+x24=1x31+x32+x33+x34=1x41+x42+x43+x44=1 x11+x21+x31+x41=1x12+x22+x32+x42=1x13+x23+x33+x43=1x14+x24+x34+x44=1xij=0或 1，i，j=1，2，3，42 进行求解本次我们使用matlab求解具体程序如下：function y,fval=maxzp(m,c) %m为c中元素的最大值 m,n=size(c);c=m+zeros(n)-c; %将求极小值的目标函数系数矩阵转换成求极大值的系 数矩阵 m=1.0; c=c; f=c(:);aeq

29、=zeros(2*n,n*n); for i=1:naeq(1:n,1+(i-1)*n:i*n)=eye(n,n); endfor i=1:naeq(n+i,1+(i-1)*n:i*n)=ones(1,n); endbeq=ones(2*n,1); lb=zeros(n*n,1); ub=ones(n*n,1);x=linprog(f,aeq,beq,lb,ub); y=reshape(x,n,n); y=y;y=round(y);sol=zeros(n,n); for i=1:n for j=1:nif y(i,j)=1sol(i,j)=c(j,i);endendendfval=sum(so

30、l(:);fval=m*n-fval; %求出极大值的目标函数值其中，c=1.00.70.5 0.4;>> y,fval=maxzp(m,c)输出结果为：optimization terminated.y =0 0 1 01 0 0 00 1 0 00 0 0 1fval =2.40003程序分析显然运用matlab专用软件进行求解简单很多，本程序是我们通过学习很快写出的，不过需要一定matlab运用基础。深刻体会到数学工具的强大。八模型检验在上述假设中我们可以准确求出最合理的安排，这种规划方法一定程度上减少了

31、误差存在，但是我们考虑到当人数上升到100以后，担心c+程序运行会很费时，当然matlab可能不会存在这个问题，所以推荐大家运用matlab进行数学求解。 九整体总结通过运用匈牙利算法和0-1整数规划同时对指派问题求解，我们发现用0-1整数规划的方法来求解可以更简单，也更方便程序的阅读和理解。用0-1整数规划可以给出运筹学中指派问题的最优解，对于求目标函数极大值与极小值的程序有略微的区别，原程序为求目标函数的极小值，而当求极大值时则可改变一下程序。这里我们自己编的c+程序如果要想求得最大值可能需要改变一下矩阵数据，让数据全部取负，然后计算就可以，最后结果再次取负就可以。同样对于matlab程序

32、并不是直接给目标函数的系数矩阵加负号，而是要做一些简单的矩阵变换。当然matlab是其求指派问题的最优解的方法都是简单易行的，方便求解，可以说是matlab程序的一个优点。但是我们的c+程序，也是很锻炼能力的，在编写过程中她完全依靠匈牙利法，这有助于我们更好地掌握匈牙利法，可以说编程的过程就是一个再现的过程，这相对于matlab来说无疑很有助于学习，所以我们坚持用c+写完了整体程序，并顺利测通，可用性很强。 十参考文献1 胡运权运筹学基础及应用m北京：高等教育出版社，2009。2 薛毅、耿美英 运筹学与实验m 北京：北京电子工业出版社, 2008。3 常巍. matlab2007 基础与提高.

33、 北京：电子工业出版社，2007。4 王海英等. 图论算法及matlab实现. 北京：北京航空航天大学出版社，2010。 蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇

34、蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄

35、肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅

36、膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆

37、膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄

38、膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄

39、芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂

40、芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃

41、莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃

42、莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁

43、莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆芈荿螈袅莁薅蚄羄肀莇薀羄膃薃袈羃莅莆袄羂蒇蚁螀羁膇蒄蚆羀艿蚀薂罿莁蒂袁羈肁蚈螇肈膃蒁蚃肇芆蚆蕿肆蒈葿羈肅膈莂袄肄芀薇螀肃莂莀蚆肃肂薆薂膂膄莈袀膁芇薄螆膀荿莇蚂腿聿薂蚈膈芁蒅羇膇莃蚀袃膇蒅蒃蝿膆膅虿蚅螂芇蒁薁袁莀蚇衿袀聿蒀螅袀节蚅螁衿莄薈蚇袈蒆莁羆袇膆薆袂袆蚁螂膁莈螃羇肆莇蒃螀羂莆薅羆莁莅螈螈芇莅袀肄膃莄蕿袇聿莃蚂肂羅莂螄袅芄蒁蒄肀膀蒀薆袃肆葿蚈聿肂葿袁袂莀蒈薀螄芆蒇蚃羀膂蒆螅螃肈蒅蒅羈羄薄薇螁芃薃虿羆腿薃螂蝿膅薂薁肅肁薁蚃袈荿薀螆肃芅蕿袈袆膁薈薈肁肇芅蚀袄羃芄螂肀节芃蒂袂芈节蚄膈膄芁螇羁肀芁衿螄荿芀蕿罿芅艿蚁螂膁莈螃羇肆莇蒃螀羂莆薅羆莁莅螈螈芇莅袀肄膃莄蕿袇聿莃蚂肂羅莂螄袅芄蒁蒄肀膀蒀薆袃肆葿蚈聿肂葿袁袂莀蒈薀螄芆蒇蚃羀膂蒆螅螃肈蒅蒅羈羄薄薇螁芃薃虿羆腿薃螂蝿膅薂薁肅肁薁蚃袈荿薀螆肃芅蕿袈袆膁薈薈肁肇芅蚀袄羃芄螂肀节芃蒂袂芈节蚄膈膄芁螇羁肀芁衿螄荿芀蕿罿芅艿蚁螂膁莈螃羇肆莇蒃螀羂莆薅羆莁莅螈螈芇莅袀肄膃莄蕿袇聿莃蚂肂羅莂螄袅芄蒁蒄肀膀蒀薆袃肆葿蚈聿肂葿