最大二分匹配问题在现实生活中比较普遍,常常出现在任务分配上。例如,有5个员工,4个不同的任务,而不同员工能够完成不同或相同的任务。也就是说,有的员工只会做这个任务,有的员工会做那个任务,有的员工会做一些任务。图解如下:左边代表员工,右边代表任务,连线代表有能力完成。
我们的问题是合理安排员工,尽可能地完成最多的任务数。上图中阴影部分为一种最好的分配方式。前一篇文章中,我们介绍了最大流问题,在这里我们可以将最大二分匹配问题转变成最大流问题。具体如下图所示:
其中每条边的最大流量限制为1,因此要求能完成的最大任务数,相当于求转变后的网络的最大流,而最大流问题在前面已经提及。
具体的程序实现如下:
#include<stdio.h>
#define N 11 //顶点数
/********************************************************\
函数功能:从残留网络中找到从源点s到汇点t的增广路径
输入:残留网络矩阵、记录前一顶点的矩阵
输出:0表示未找到路径,1表示找到了路径
\********************************************************/
int search(int dist1[N][N],int vertex[N])
{
int queue[20]={-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1};//初始化,-1作为标记
int flag[N]={0};//标记顶点是否被访问过
int front=0;//队列头指针
int rear=0;//队列尾指针
int temp=0;
int i;
queue[rear]=0;//将源点入队列
flag[0]=1;
rear++;
while(queue[front]!=-1)//队列不为空
{
temp=queue[front];//取队头元素
for(i=0;i<N;i++)
if(dist1[temp][i]!=0&&flag[i]==0)//广度搜索法
{
queue[rear++]=i;//入队
flag[i]=1;
vertex[i]=temp;//标记当前节点的上一个节点
if(i==10)//找到汇点后就不用再寻找了
return 1;
}
front++;
if(front==N)//找完所有顶点后停止寻找
break;
}
if(queue[rear-1]!=10)//没有找到路径到汇点
return 0;
return 1;
}
/**************************************************************************\
函数功能:修改残余网络矩阵和原始流网络矩阵
输入:原始流网络矩阵、残留网络矩阵、记录前一顶点的矩阵、源点和汇点,最大流值
输出:0表示未找到路径,1表示找到了路径
\***************************************************************************/
int modify(int dist[N][N],int dist1[N][N],int vertex[N],int s,int t,int flow)
{
int i,j;
int min=10000;//记录找到的路径所能通过的最大流
i=vertex[t];
j=t;
while(j!=s)//寻找路径所含边的最大流量值中的最小值
{
if(dist1[i][j]<min)
min=dist1[i][j];
j=i;
i=vertex[i];
}
i=vertex[t];
j=t;
flow=min;//记录最大流量
while(j!=s)
{
if(dist1[i][j]>0)//更改残余图
dist1[j][i]=dist1[i][j];
dist1[i][j]=dist1[i][j]-min;
dist[i][j]=dist[i][j]+min-dist[j][i];//更改原始流网路
if(dist[i][j]<0)
{
dist[j][i]=-dist[i][j];
dist[i][j]=0;
}
j=i;
i=vertex[i];
}
printf("原始流网络矩阵:\n");
for(int i=0;i<N;i++)
{
for(int j=0;j<N;j++)
printf("%d ",dist[i][j]);
printf("\n");
}
printf("残留网络矩阵:\n");
for(int i=0;i<N;i++)
{
for(int j=0;j<N;j++)
printf("%d ",dist1[i][j]);
printf("\n");
}
return flow;
}
void FordFulkerson(int dist[N][N],int dist1[N][N])
{
int vertex[N]={0};//初始化
int flow=0;
while(search(dist1,vertex)==1)//当能找到增广路径时
{
flow=flow+modify(dist,dist1,vertex,0,10,0);
printf("\n");
}
printf("最大流为%d \n",flow);
}
void main()
{
int dist1[N][N]={{0,1,1,1,1,1,0,0,0,0,0},
{0,0,0,0,0,0,1,0,0,0,0},
{0,0,0,0,0,0,1,0,1,0,0},
{0,0,0,0,0,0,0,1,1,1,0},
{0,0,0,0,0,0,0,0,1,0,0},
{0,0,0,0,0,0,0,0,1,0,0},
{0,0,0,0,0,0,0,0,0,0,1},
{0,0,0,0,0,0,0,0,0,0,1},
{0,0,0,0,0,0,0,0,0,0,1},
{0,0,0,0,0,0,0,0,0,0,1},
{0,0,0,0,0,0,0,0,0,0,0}};
int dist[N][N]={0};//初始的流网络为0
FordFulkerson(dist,dist1);
}
程序结果为最大流为3,与前面的图解一致,但是我们发现结果与图解的不同,这说明最大二分匹配可以有不同的解即有不同的分配方案。
注:如果程序出错,可能是使用的开发平台版本不同,请点击如下链接: 解释说明
原文:http://blog.csdn.net/tengweitw/article/details/17848773
作者:nineheadedbird