POJ 2112 Optimal Milking（最大流+Floyd+二分||二分图匹配）

题目链接：POJ2112

题意：农场主John 将他的K(1≤K≤30)个挤奶器运到牧场，在那里有C(1≤C≤200)头奶牛，在奶牛和挤奶器之间有一组不同长度的路。K个挤奶器的位置用1～K的编号标明，奶牛的位置用K+1～K+C 的编号标明。每台挤奶器每天最多能为M(1≤M≤15)头奶牛挤奶。编写程序，寻找一个方案，安排每头奶牛到某个挤奶器挤奶，并使得C 头奶牛需要走的所有路程中的最大路程最小。每个测试数据中至少有一个安排方案。每条奶牛到挤奶器有多条路。

题解：看了好久才懂“C 头奶牛需要走的所有路程中的最大路程最小”这句话的意思。每头奶牛到挤奶器之间都有一个最短距离的路程，现在要求 所有奶牛的最短路程中的最大值最小。 先用Floyd求出各个实体之间的最短距离，然后二分枚举最短距离，用最大流求出最大能完成的奶牛个数，小于总的奶牛个数，说明距离太小，大于说明距离可以继续缩小。

最大流建图：源点和奶牛相连，容量为1，机器和汇点相连，容量为M，奶牛和机器之间，连接 距离小于当前枚举距离的边，容量为1。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int MAX=200+30+10;
const int INF=0x3f3f3f3f;
int dis[MAX][MAX];
int g[MAX][MAX];
int level[MAX];
int N;
int S,E;
int k,c,m;

int BFS()//构建层次网络
{
	memset(level,0xff,sizeof(level));
	queue<int> que;
	que.push(S);
	level[S]=0;
	while(!que.empty())
	{
		int u=que.front();
		que.pop();
		for(int v=0;v<N;v++)
		{
			if(level[v]<0&&g[u][v]>0)//存在路径而且没有被标记层次
			{
				level[v]=level[u]+1;
				que.push(v);
			}
		}
	}
	if(level[E]>0) return 1;
	return 0;
}
int DFS(int x,int low)
{
	if(x==E) return low;
	for(int i=0;i<N;i++)
	{
		int a;
		if(g[x][i]>0&&level[i]==level[x]+1&&(a=DFS(i,min(low,g[x][i]))))
		{
			g[x][i]-=a;
			g[i][x]+=a;
			return a;
		}
	}
	return 0;
}
int Dinic()
{
	int ans=0;
	while(BFS())
	{
		int tans;
		while(tans=DFS(S,INF)) ans+=tans;
	}
	return ans;
}

void floyd()
{
	for(int m=1;m<=k+c;m++)
	{
		for(int i=1;i<=k+c;i++)
		{
			for(int j=1;j<=k+c;j++)
			{
				dis[i][j]=min(dis[i][j],dis[i][m]+dis[m][j]);
			}
		}
	}
}
void build(int maxdis)
{
	memset(g,0,sizeof(g));
	for(int i=1;i<=k;i++)
	{
		g[i][E]=m;
	}
	for(int i=k+1;i<=k+c;i++)
	{
		g[S][i]=1;
	}
	for(int i=k+1;i<=k+c;i++)
	{
		for(int j=1;j<=k;j++)
		{
			if(dis[i][j]<=maxdis) g[i][j]=1;
		}
	}
}
int main()
{
	//freopen("in.txt","r",stdin);
	//freopen("out.txt","w",stdout);
	while(scanf("%d%d%d",&k,&c,&m)!=EOF)
	{
		for(int i=1;i<=k+c;i++)
		{
			for(int j=1;j<=k+c;j++)
			{
				scanf("%d",&dis[i][j]);
				if((dis[i][j]==0&&i!=j)) dis[i][j]=INF;
			}
		}
		floyd();
		N=k+c+2;S=0;E=N-1;
		int l=0,r=200*MAX;
		while(l<r)
		{
			int mid=(l+r)/2;
			build(mid);
		/*	for(int i=0;i<N;i++)
			{
				for(int j=0;j<N;j++)
				{
					cout<<g[i][j]<<" ";
				}
				cout<<endl;
			}*/
			int ans=Dinic();
			//cout<<ans<<endl;
			if(ans>=c)r=mid;
			else l=mid+1;
		//	cout<<l<<" "<<r<<endl;
		}
		cout<<r<<endl;
	}
	return 0;
}

学了二分图，又用二分图多重匹配写了一遍。依旧要用floyd跑一遍最短距离，然后二分最大距离，然后在牛和机器之间建立边，匹配的时候如果加一个小于距离的限制条件。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int MAX=200+30+10;
const int MAXC=200+5;
const int MAXK=30+5;
const int INF=0x3f3f3f3f;

int map[MAX][MAX];
int g[MAXC][MAXK];
int k,c,m;
int mid;

void floyd(int n)
{
	for(int k=0;k<n;k++)
	{
		for(int i=0;i<n;i++)
		{
			for(int j=0;j<n;j++)
			{
				map[i][j]=min(map[i][j],map[i][k]+map[k][j]);
			}
		}
	}
}
int match[MAXK][MAXC];
int vis[MAXK];
bool dfs(int u)
{
	for(int v=0;v<k;v++)
	{
		if(!vis[v]&&g[u][v]<=mid)
		{
			vis[v]=1;
			if(match[v][0]<m)
			{
				match[v][++match[v][0]]=u;
				return true;
			}
			for(int i=1;i<=m;i++)
			{
				if(dfs(match[v][i]))
				{
					match[v][i]=u;
					return true;
				}
			}
		}
	}
	return false;
}

int hungray()
{
	int res=0;
	memset(match,0,sizeof(match));
	for(int i=0;i<c;i++)
	{
		memset(vis,0,sizeof(vis));
		if(dfs(i)) res++;
	}
	return res;
}

int main()
{
	//freopen("in.txt","r",stdin);
	//freopen("out.txt","w",stdout);
	while(scanf("%d%d%d",&k,&c,&m)!=EOF)
	{
		memset(map,0x3f,sizeof(map));
		for(int i=0;i<k+c;i++)
		{
			for(int j=0;j<k+c;j++)
			{
				scanf("%d",&map[i][j]);
				if(map[i][j]==0&&i!=j) map[i][j]=INF;
			}
		}
		floyd(k+c);
/*		for(int i=0;i<k+c;i++)
		{
			for(int j=0;j<k+c;j++)
				cout<<map[i][j]<<" ";
			cout<<endl;
		}*/
		memset(g,0x3f,sizeof(g));
		for(int i=0;i<c;i++)
		{
			for(int j=0;j<k;j++)
			{
				g[i][j]=map[k+i][j];
			//		cout<<g[k+i][j]<<" ";
			}
		//	cout<<endl;
		}
		int l=0,r=50000;
		while(l<r)
		{
			mid=(l+r)/2;
			if(hungray()>=c) r=mid;
			else l=mid+1;
		}
		printf("%d\n",r);
	}
	return 0;
}