pku2112 最大流，弗洛伊德，二分答案

第一次写二分答案的最大流。其实二分答案是一个很神奇的东西，很多问题本来没法解决然后二分答案变成了判定性问题就可以很简单的解决了。膜拜cai0715，信偲牛，得永生！

题目大意：

  n个挤奶机，m个奶牛，挤奶机和挤奶机之间，挤奶机和奶牛之间，奶牛和奶牛之间都有距离，求满足让所有奶牛都喝到牛奶的情况下，所走最大距离的最小值。

简略题解：

  其实一想就是最大流，但是这个最大流怎么写呢？我想了很长时间终于没想出来，然后看了偲哥的blog，顿悟了，原来二分答案是如此的神奇。

  先floyd求出任意两点间的最短距离。

  然后从源点向奶牛连一条流量为1费用为0的边，从挤奶机向汇点连一条流量为k费用为0的边。

  我们二分枚举答案。

  然后构图，如过奶牛和挤奶机的距离不超过答案，那么连边，流量为1。

  然后求最大流，如果最大流=m那么r=mid-1 else l=mid+1

我发现最近写程序随便一写就100+行了，额，果然是一坨noi

附程序：

program pku_2112; const inf=1000000001; var a,c,tc:array[0..250,0..250] of longint; d,num:array[0..250] of longint; n,m,max,st,ed:longint; {---------------------------} procedure floyd; var i,j,k:longint; begin for k:=st to ed do for i:=st to ed do for j:=st to ed do if (i<>j) and (j<>k) and (i<>k) then if a[i,k]+a[k,j]<a[i,j] then a[i,j]:=a[i,k]+a[k,j]; end; {---------------------------} function min(x,y:longint):longint; begin if x<y then exit(x); exit(y); end; {---------------------------} function dfs(u,flow:longint):longint; var i,t:longint; begin if u=ed then exit(flow); dfs:=0; for i:=st to ed do if (c[u,i]>0) and (d[u]=d[i]+1) then begin t:=dfs(i,min(flow-dfs,c[u,i])); dec(c[u,i],t); inc(c[i,u],t); inc(dfs,t); if dfs=flow then exit(dfs); end; if d[st]=n+m+2 then exit; dec(num[d[u]]); if num[d[u]]=0 then d[st]:=n+m+2; inc(d[u]); inc(num[d[u]]); end; {---------------------------} function sap(limit:longint):boolean; var i,j,ans,flow:longint; begin fillchar(d,sizeof(d),0); fillchar(num,sizeof(num),0); num[0]:=n+m+2; c:=tc; ans:=0; flow:=0; for i:=1 to m do for j:=1 to n do if a[i+n,j]<=limit then c[i+n,j]:=1 else c[i+n,j]:=0; while d[st]<n+m+2 do begin flow:=dfs(st,inf); inc(ans,flow); end; if ans=m then exit(true) else exit(false); end; {---------------------------} procedure init; var i,j,limit:longint; begin fillchar(c,sizeof(c),0); fillchar(a,sizeof(a),50); max:=0; readln(n,m,limit); for i:=1 to n+m do for j:=1 to n+m do begin read(a[i,j]); if a[i,j]>max then max:=a[i,j]; if a[i,j]=0 then a[i,j]:=inf; end; st:=0;ed:=n+m+1; for i:=1 to m do begin c[st,i+n]:=1;a[st,i+n]:=0;end; for i:=1 to n do begin c[i,ed]:=limit;a[i,ed]:=0;end; tc:=c; end; {---------------------------} procedure main; var l,r,mid:longint; begin l:=1;r:=(n+m)*max+1; repeat mid:=(l+r) div 2; if sap(mid) then r:=mid-1 else l:=mid+1; until l>r; writeln(l); end; {---------------------------} begin init; floyd; main; end.