运筹学指派问题的匈牙利法Word文档下载推荐.docx

运筹学指派问题的匈牙利法Word文档下载推荐.docx

《运筹学指派问题的匈牙利法Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《运筹学指派问题的匈牙利法Word文档下载推荐.docx（18页珍藏版）》请在冰豆网上搜索。

2.算法流程或步骤：

1.将原始效益矩阵C的每行、每列各元素都依次减去该行、该列的最小元素，使每行、每列都至少出现一个0元素，以构成等价的效益矩阵C’。

2.圈0元素。

在C’中未被直线通过的含0元素最少的行（或列）中圈出一个0元素，通过这个0元素作一条竖（或横）线。

重复此步，若这样能圈出不同行不同列的n个0元素，转第四步，否则转第三步。

3.调整效益矩阵。

在C’中未被直线穿过的数集D中，找出最小的数d，D中所有数都减去d，C’中两条直线相交处的数都加的d。

去掉直线，组成新的等价效益矩阵仍叫C’，返回第二步。

4.令被圈0元素对应位置的=1,其余=0，这就是一种最优分配。

最低总耗费是C中使=1的各位置上各元素的和。

算法流程图：

Y

N

3.算法源程序：

#include<

iostream.h>

typedefstructmatrix

{

floatcost[101][101];

intzeroelem[101][101];

floatcostforout[101][101];

intmatrixsize;

intpersonnumber;

intjobnumber;

}matrix;

matrixsb;

intresult[501][2];

voidtwozero（matrix&

sb）;

voidjudge（matrix&

sb,intresult[501][2]）;

voidrefresh（matrix&

voidcirclezero（matrix&

matrixinput（）;

voidoutput（intresult[501][2],matrixsb）;

voidzeroout（matrix&

matrixinput（）

matrixsb;

intm;

intpnumber,jnumber;

inti,j;

floatk;

charw;

cout<

<

"

指派问题的匈牙利解法:

endl;

求最大值，请输入1;

求最小值，请输入0:

cin>

>

m;

while（m!

=1&

&

m!

=0）

{

cout<

请输入1或0:

cin>

}

请输入人数（人数介于1和100之间）:

pnumber;

while（pnumber<

1||pnumber>

100）

{

请输入合法数据:

请输入工作数（介于1和100之间）:

jnumber;

while（jnumber<

1||jnumber>

请输入"

pnumber<

行"

jnumber<

列的矩阵,同一行内以空格间隔,不同行间以回车分隔,以$结束输入:

\n"

;

for（i=1;

i<

=pnumber;

i++）

for（j=1;

j<

=jnumber;

j++）

{

cin>

sb.cost[i][j];

sb.costforout[i][j]=sb.cost[i][j];

}

w;

if（jnumber>

pnumber）

for（i=pnumber+1;

for（j=1;

{

sb.cost[i][j]=0;

sb.costforout[i][j]=0;

}

else

if（pnumber>

jnumber）

for（i=1;

for（j=jnumber+1;

{

sb.cost[i][j]=0;

sb.costforout[i][j]=0;

}

sb.matrixsize=pnumber;

if（pnumber<

sb.matrixsize=jnumber;

sb.personnumber=pnumber;

sb.jobnumber=jnumber;

if（m==1）

k=0;

for（i=1;

=sb.matrixsize;

if（sb.cost[i][j]>

k）

k=sb.cost[i][j];

sb.cost[i][j]=k-sb.cost[i][j];

returnsb;

}

sb）

intp;

for（i=0;

sb.cost[i][0]=0;

for（j=1;

sb.cost[0][j]=0;

if（sb.cost[i][j]==0）

sb.cost[i][0]++;

sb.cost[0][j]++;

sb.cost[0][0]++;

for（j=0;

sb.zeroelem[i][j]=0;

k=sb.cost[0][0]+1;

while（sb.cost[0][0]<

k=sb.cost[0][0];

if（sb.cost[i][0]==1）

for（j=1;

if（sb.cost[i][j]==0&

sb.zeroelem[i][j]==0）

break;

sb.zeroelem[i][j]=1;

sb.cost[i][0]--;

sb.cost[0][j]--;

sb.cost[0][0]--;

if（sb.cost[0][j]>

0）

for（p=1;

p<

p++）

if（sb.cost[p][j]==0&

sb.zeroelem[p][j]==0）

{

sb.zeroelem[p][j]=2;

sb.cost[p][0]--;

sb.cost[0][j]--;

sb.cost[0][0]--;

}

if（sb.cost[0][j]==1）

for（i=1;

if（sb.cost[i][0]>

if（sb.cost[i][p]==0&

sb.zeroelem[i][p]==0）

sb.zeroelem[i][p]=2;

sb.cost[i][0]--;

sb.cost[0][p]--;

if（sb.cost[0][0]>

twozero（sb）;

judge（sb,result）;

intp,q;

intm,n;

matrixst;

if（sb.cost[i][0]>

break;

if（i<

=sb.matrixsize）

st=sb;

if（sb.cost[i][j]==0&

for（q=1;

q<

q++）

if（sb.cost[i][q]==0&

sb.zeroelem[i][q]==0）

{

sb.zeroelem[i][q]=2;

sb.cost[i][0]--;

sb.cost[0][q]--;

sb.cost[0][0]--;

}

for（p=1;

if（sb.cost[p][j]==0&

sb.zeroelem[p][j]=2;

sb.cost[p][0]--;

sb.cost[0][j]--;

k=sb.cost[0][0]+1;

while（sb.cost[0][0]<

k=sb.cost[0][0];

for（p=i+1;