POJ 2391 Ombrophobic Bovines【网络流建模汇总】拆点最大流

题意：一个无向图。点i处有Ai头牛，牛棚能够容纳Bi头牛，从i到j有个时间，求一个最短时间T，让所有的牛都能够进入牛棚

这个题呢，思路其实很简单：二分时间就能找到最小值，对于每个时间点T0，判断是不是可行就好

可行的意思：在给定的时间内，跑出来的最大流等于各个点的牛的数量之和

时间影响了什么？影响了从i点到j点能不能去。时间不够的话，就去不了

然后看看题目中的数据，n不超过200：所以用floyd算法得到任意两点之间的最短时间

这样建模就很清晰了：

拆点：

对于任意不同的两点i，j：如果time【i】【j】不超过当前二分的时间T0，说明第i个点的牛在时间内可以到第j个牛棚

所以连边（i1，j2，INF）容量为INF是说，第i个牛棚的所有牛都可以去（当然用Ai也可以）

必须有的边：

源点到所有的点：（s，i1，Ai），这个是牛的数量的限制

所有的点到汇点：（i2，t，Bi），这个是牛棚中的最大值的限制

然后，拆点之后，自己和自己必须连边：因为距离为0，无论时间多少都是可以去到的点

（i1，i2，INF）

其实，这个题的难点不在拆点和连边：那么在哪儿？

在如何确定二分的L，R中的R（也就是最大时间上限）：是题目中的1e9这个值吗？

考虑一组非常特别特殊的数据：

10 9
1 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 1
1 2 1000000000
2 3 1000000000
3 4 1000000000
4 5 1000000000
5 6 1000000000
6 7 1000000000
7 8 1000000000
8 9 1000000000
9 10 1000000000

看到这种特别恶心的数据：我们知道最后的答案是

9000000000

因为从1到10只有一条路径，必须慢慢的走

那么，拓展看来，如果都是这种边，最大的时间可能是1000*1e9=1e12（注意，这个值已经超过int了，需要用long long来搞）

所以把值设为这个应该就是没有问题的（当然二分是比较快的，再设大几个数量级其实更保险，反正多算不了几次）

讨论了这些，注意好细节，应该可以上代码上模板1A了吧：

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<stdio.h>
#include<cstring>
#include<string.h>
using namespace std;

#define LL long long
const int maxm=50000;
const int maxn=1050;
const int INF=0x3f3f3f3f;
const LL LLINF=1e16;
int s,t,tot;
int n,m;
int f[maxn];
int p[maxn];
int all;
LL L,R,ans,mid,dist[maxn][maxn];

struct Edge{
    int to,nxt,cap,flow;
}edge[maxm*10];
int tol,Head[maxm*10];

void init(){
    memset(Head,-1,sizeof(Head));
    tol=2;
}

void addedge(int u,int v,int w,int rw=0){
    edge[tol].to=v;
    edge[tol].cap=w;
    edge[tol].flow=0;
    edge[tol].nxt=Head[u];
    Head[u]=tol++;

    edge[tol].to=u;
    edge[tol].cap=rw;
    edge[tol].flow=0;
    edge[tol].nxt=Head[v];
    Head[v]=tol++;
}

int Q[maxn*10],dep[maxn*10],cur[maxn*10],sta[maxn*10];

bool bfs(int s,int t,int n){
    int Front=0,Tail=0;
    memset(dep,-1,sizeof(dep[0])*(n+1));
    dep[s]=0;
    Q[Tail++]=s;
    while(Front<Tail){
        int u=Q[Front++];
        for(int i=Head[u];i!=-1;i=edge[i].nxt){
            int v=edge[i].to;
            if (edge[i].cap>edge[i].flow&&dep[v]==-1){
                dep[v]=dep[u]+1;
                if (v==t) return true;
                Q[Tail++]=v;
            }
        }
    }
    return false;
}

int dinic(int s,int t,int n){
    int maxflow=0;
    while(bfs(s,t,n)){
        for(int i=0;i<n;i++) cur[i]=Head[i];
        int u=s,tail=0;
        while(cur[s]!=-1){
            if (u==t){
                int tp=INF;
                for(int i=tail-1;i>=0;i--)
                    tp=min(tp,edge[sta[i]].cap-edge[sta[i]].flow);
                maxflow+=tp;
                for(int i=tail-1;i>=0;i--){
                    edge[sta[i]].flow+=tp;
                    edge[sta[i]^1].flow-=tp;
                    if (edge[sta[i]].cap==edge[sta[i]].flow) tail=i;
                }
                u=edge[sta[tail]^1].to;
            }
            else if (cur[u]!=-1&&edge[cur[u]].cap>edge[cur[u]].flow&&dep[u]+1==dep[edge[cur[u]].to]){
                sta[tail++]=cur[u];
                u=edge[cur[u]].to;
            }
            else{
                while(u!=s&&cur[u]==-1) u=edge[sta[--tail]^1].to;
                cur[u]=edge[cur[u]].nxt;
            }
        }
    }
    return maxflow;
}

void add(LL x){
    init();
    for(int i=1;i<=n;i++){
        addedge(s,i,f[i]);
        addedge(i+n,t,p[i]);
        addedge(i,i+n,INF);
    }
    for(int i=1;i<=n;i++)
        for(int j=i+1;j<=n;j++)
        if (dist[i][j]<=x){
            addedge(i,j+n,INF);
            addedge(j,i+n,INF);
        }
}

int main(){
    freopen("input.txt","r",stdin);
    while(scanf("%d%d",&n,&m)!=EOF){
        s=0;t=2*n+1;tot=t+1;all=0;
        for(int i=1;i<=n;i++){
            scanf("%d%d",&f[i],&p[i]);
            all+=f[i];
        }
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                dist[i][j]=LLINF;
        int u,v,w;
        for(int i=1;i<=m;i++){
            scanf("%d%d%d",&u,&v,&w);
            if (w<dist[u][v])
                dist[u][v]=dist[v][u]=w;
        }
        for(int k=1;k<=n;k++)
        for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            if (dist[i][j]>dist[i][k]+dist[k][j])
                dist[i][j]=dist[i][k]+dist[k][j];
        ans=-1;L=0;R=LLINF-1;
        int maxflow;
        while(L<=R){
            mid=(L+R)>>1;
            add(mid);
            maxflow=dinic(s,t,tot);
            if (maxflow==all){
                ans=mid;
                R=mid-1;
            }
            else L=mid+1;
        }
        printf("%lld\n",ans);
    }
    return 0;
}