POJ 2112 Optimal Milking(Floyd + 二分 + 最大流）

题目连接：http://poj.org/problem?id=2112

想了N久，建图的方法还是没想完全正确，看了别人的建图方法，恍然大悟啊，建图很重要啊，要多多练习。。。

题意：有C头牛，K台加工奶的机器，每台机器每天可以加工M头牛产的奶，机器和牛在牧场的不同的位置分布着，给出他们之间任何两点的距离，现在让你给出一个分配的方案，该方案须满足距离机器最远的牛距离最小。

题解：

我们可以二分枚举所有的可能距离d，用最大流判断否存在一种分配方案。

先用Floyd算法求出任意两点之间的最短路，如果从机器 Ki 到牛 Cj 的距离不大于d,那么就从Ki到Cj建边，权值为1，建立一个超级源点S，从S建立到所有机器的边，权值为M,建立一个超级汇点U，建立所有牛到U的边，权值为1，求S到U的最大流Maxflow，如果Maxflow等于牛的个数C，则该方案符合要求，二分求出最小的d即可。

参考代码：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;

#define NN 405
#define MAX 100000
#define INS 1<<29
#define CLR(arr,v) memset(arr,v,sizeof(arr))

template<int MaxV,int MaxE>
class MaxFlow{
public:
	void Clear(){
		pos = 0;
		CLR(h,-1); CLR(cnt,0); CLR(Flow,0);
		CLR(gap,0);
	}
	void add(int u,int v,int f){
		num[pos] = v;
		cap[pos] = f;
		next[pos] = h[u];
		h[u] = pos++;

		num[pos] = u;
		cap[pos] = 0;
		next[pos] = h[v];
		h[v] = pos++;
	}
	int GetMaxFlow(int s,int t){
		int maxflow = 0,cur = s,minflow = INS;
		memcpy(Cur,h,sizeof(h)); 
		Init(t);
		while(gap[cur] < t){
			int &i = Cur[cur];
			for(;i != -1;i = next[i]){
				if(cap[i] - Flow[i] > 0 && gap[cur] - gap[ num[i] ] == 1){
					minflow = min(minflow,cap[i] - Flow[i]);
					pre_v[ num[i] ] = cur;
					pre_e[ num[i] ] = i;
					cur = num[i];
					if(cur == t){
						maxflow += minflow;
						while(cur != s){
							Flow[ pre_e[cur] ] += minflow;
							Flow[ pre_e[cur]^1 ] -= minflow;
							cur = pre_v[cur];
						}
						minflow = INS;
					}
					break;
				}	
			}	
			if(i == -1){
				if(--cnt[ gap[cur] ] == 0) return maxflow;
				i = h[cur];
				gap[cur] = INS;
				for(int j = h[cur];j != -1;j = next[j])
					if(gap[ num[j] ] + 1 < gap[cur] && cap[j] > Flow[j])
						gap[cur] = gap[ num[j] ] + 1;
				if(gap[cur] < INS) ++cnt[ gap[cur] ];
				cur = pre_v[cur];
			}
		}
		return maxflow;
	}
private:
	int h[MaxV],Cur[MaxV],num[MaxE],cap[MaxE],next[MaxE],gap[MaxV],cnt[MaxV],pre_v[MaxV],pre_e[MaxE],Flow[MaxE],que[MaxV];
	int pos;
	void Init(int n){
		int head = 0,total = 0;
		cnt[0] = 1;
		que[total++] = n;
		while(head < total){
			int p = que[head++];
			for(int i = h[p]; i != -1;i = next[i]){
				if((i&1) && !gap[ num[i] ]){
					gap[ num[i] ] = gap[p] + 1;
					cnt[ gap[ num[i] ] ]++;
					que[total++] = num[i];	
				}
			}
		}
	}
};

int map[NN][NN];

void Floyd(int s,int t){
	for(int k = s;k <= t;++k)
		for(int i = s;i <= t;++i)
			for(int j = s;j <= t;++j)
				if(map[i][k] + map[k][j] < map[i][j])
					map[i][j] = map[i][k] + map[k][j];
}

MaxFlow<NN,200000> g;

int main()
{
	int K,C,M;
	while(~scanf("%d%d%d",&K,&C,&M))
	{
		int n = K + C,sum = 0;
		for(int i = 1;i <= n;++i)
			for(int j = 1;j <= n;++j)
			{
				scanf("%d",&map[i][j]);
				if(!map[i][j]) map[i][j] = INS;
			}
		Floyd(1,n);
		for(int i = 1;i <= n;++i)
			for(int j = 1;j <= n;++j)
				if(map[i][j] != INS) sum += map[i][j];
		int L = 0,R = sum,s = 0,t = K+C+1;
		while(L < R)
		{
			int mid = (L + R) / 2;
			g.Clear();
			for(int i = 1;i <= K;++i)
				g.add(0,i,M);
			for(int i = 1;i <= C;++i)
				g.add(K+i,t,1);
			for(int i = 1;i <= K;++i)
				for(int j = 1;j <= C;++j)
					if(map[i][K+j] <= mid) g.add(i,K+j,1);
			if(g.GetMaxFlow(s,t) == C) R = mid;
			else L = mid + 1;
		}
		printf("%d\n",R);
	}
	return 0;
}