poj2391 二分+floyed+网络流Dinic（附数据网址）

《POJ 2391 Ombrophobic Bovines》
【题目大意】
给定一个无向图，点i处有Ai头牛，点i处的牛棚能容纳Bi头牛，求一个最短时间T使得在T时间内所有的牛都能进到某一牛棚里去。（1 <= N <= 200, 1 <= M <= 1500, 0 <= Ai <= 1000, 0 <= Bi <= 1000, 1 <= Dij <= 1,000,000,000）
题目戳这里：https://vjudge.net/problem/POJ-2391
数据戳这里：http://contest.usaco.org/MAR05_4.htm
顺便吐槽一下：我最开始一直没过，看了数据才发现有个极限大数据，一条链，所有距离都是最大，所以dis数组一定要开long long
坑死了。
可以很清楚的知道，我们先用floyd处理两两牛棚之间的最短距离，N的三次方N为200完全大丈夫，萌大奶。
弄完floyed后，然后我们一看根据人生的经验就知道这十分像一道网络流的题，从源点流出多少头牛，流进汇点就是牛棚的容量。
但还不够，因为牛走到牛棚有距离，所以我们需要二分这个距离，把所有距离小于等于二分距离的边连起来。
但是这里会有一个问题，比如1到2可以走，2到3可以走而1到3不能走，但按照这样建图的话，1是可以流到3的，这是不严谨的，并且会wa，所以我们把点拆成两个，1连到2”，2连到3”就可以完美解决这个问题了。
代码如下：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
#define del(a,b) memset(a,b,sizeof(a))
#define N 1005
#define M 100505
#define inf 0x7fffffff
using namespace std;
int head[N],num=1,n;
int m,A[N],B[N],a,b,sum;
long long dis[N][N],c;
struct Edge{
    int v,next,flow,cap;
};
Edge e[N*N*10];
void adde(int i,int j,int w){
    e[++num].v=j;
    e[num].next=head[i];
    e[num].cap=w;
    e[num].flow=0;
    head[i]=num;
}
void add(int u,int v,int w){
    adde(u,v,w);
    adde(v,u,0);
}
struct Dinic{
    int s,t,d[N];
    bool bfs(long long mid){
        queue<int>q;
        while(!q.empty())q.pop();
        del(d,0);q.push(s);
        d[s]=1;
        while(!q.empty()){
            int u=q.front();q.pop();
            for(int i=head[u];i;i=e[i].next){
                int v=e[i].v,U=u,V=v;
                if(v>n)V-=n;if(u>n)U-=n;
                if(e[i].cap>e[i].flow&&!d[v]&&dis[U][V]<=mid){
                    d[v]=d[u]+1;
                    q.push(v);
                }
            }
        }
        if(d[t])return true;
        else return false;
    }
    int dfs(int u,int a,long long mid){
        if(u==t||a==0)return a;
        int flow=0,f;
        for(int i=head[u];i;i=e[i].next){
            int v=e[i].v,U=u,V=v;
            if(v>n)V-=n;if(u>n)U-=n;
            if(dis[U][V]<=mid&&d[v]==d[u]+1&&(f=dfs(v,min(a,e[i].cap-e[i].flow),mid))>0){
                e[i].flow+=f;
                e[i^1].flow-=f;
                flow+=f;
                a-=f;
                if(a==0)return flow;
            }
        }
        if(!flow)d[u]=-1;
        return flow;
    }
    int maxflow(long long mid){
        int flow=0;
        for(int i=2;i<=num;i++)e[i].flow=0;
        while(bfs(mid))flow+=dfs(s,inf,mid);
        return flow;
    }
};
Dinic date;
int main(){
    long long maxw=0;
    memset(dis,63,sizeof dis);
    scanf("%d%d",&n,&m);
    date.s=0;date.t=2*n+1;
    for(int i=1;i<=n;i++){
        scanf("%d%d",&A[i],&B[i]);
        sum+=A[i];
        add(date.s,i,A[i]);
        add(i+n,date.t,B[i]);
    }
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
    add(i,j+n,inf);
    for(int i=1;i<=m;i++){
        scanf("%d%d%lld",&a,&b,&c);
        maxw+=c;
        dis[a][b]=min(dis[a][b],c);
        dis[b][a]=min(dis[b][a],c);
    }
    for(int k=1;k<=n;k++)
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
    if(dis[i][j]>dis[i][k]+dis[k][j])
    dis[i][j]=dis[i][k]+dis[j][k];
    for(int i=1;i<=n;i++)
    dis[i][i]=dis[0][i]=dis[i][n+1]=dis[i][0]=dis[n+1][i]=0;
    long long l=0,r=maxw,ans=-1;
    while(l<=r){
        long long mid=(l+r)>>1;
    //  cout<<mid<<endl;
        if(date.maxflow(mid)==sum)ans=mid,r=mid-1;
        else l=mid+1;
    }
    cout<<ans<<endl;
    return 0;
}