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）

一开始想拆点建图，0到x集合为汇，值为各个区域的牛数量， Y到终点连边，值为各个区域的容量，然后就是看怎么连x和y了

我一开始把可以连接的X和Y连起来，把可以互达的点在Y集合点那里连边，这样很麻烦，先跑一遍floyd把点到点的最短路求出来，然后直接X和Y集合可达即相连就行

二分结果，再建图，把在mid以内的X点对Y点连起来，跑最大流 判断结果即可

注意要用long long

一开始还没看清题意，一条路上可以同时走无数的牛，我一开始以为只能走一头，还敲MCMF去了。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
#define LL long long
#define INF 1LL<<60
using namespace std;
int f,p;
const int maxn=500;
struct Edge
{
    int from,to,cap,flow;
};
struct Dinic
{
    vector<Edge>edges;
    vector<int> G[maxn];
    int vis[maxn];
    int cur[maxn];
    int d[maxn];
    void init(int n)
    {
        edges.clear();
        for (int i=0; i<=n; i++)
        {
            G[i].clear();
        }
    }
    void addedge(int from,int to,int cap)
    {
        int m;
        edges.push_back((Edge)
        {
            from,to,cap,0
        });
        edges.push_back((Edge)
        {
            to,from,0,0
        });
        m=edges.size();
        G[from].push_back(m-2);
        G[to].push_back(m-1);
    }
    bool bfs(int s,int t)
    {
        memset(vis,0,sizeof vis);
        queue<int> q;
        q.push(s);
        d[s]=0;
        vis[s]=1;
        while (!q.empty())
        {
            int u=q.front();
            q.pop();
            for (int i=0; i<G[u].size(); i++)
            {
                Edge& e=edges[G[u][i]];
                if (!vis[e.to] && e.cap>e.flow)
                {
                    vis[e.to]=1;
                    d[e.to]=d[u]+1;
                    q.push(e.to);
                }
            }
        }
        return vis[t];
    }
    int dfs(int x,int a,int t)
    {
        if (x==t || a==0) return a;
        int flow=0,f;
        for (int& i=cur[x]; i<G[x].size(); i++)
        {
            Edge& e=edges[G[x][i]];
            if (d[x]+1==d[e.to] && (f=dfs(e.to,min(a,e.cap-e.flow),t))>0)
            {
                e.flow+=f;
                edges[G[x][i]^1].flow-=f;
                flow+=f;
                a-=f;
                if (a==0) break;
            }
        }
        return flow;
    }
    int maxflow(int s,int t)
    {
        int flow=0;
        while (bfs(s,t))
        {

            memset(cur,0,sizeof cur);
            flow+=dfs(s,100000000,t);
        }
              return flow;
    }
} mcmf;
int A[210],B[210];
LL path[210][210];
LL N;
void floyd()
{
    for (int i=1; i<=f; i++)
    {
        for (int j=1; j<=f; j++)
        {
            for (int k=1; k<=f; k++)
            {
                if (j==k) continue;
                path[j][k]=min(path[j][k],path[j][i]+path[i][k]);
                N=max(N,path[j][k]);
            }
        }
    }
}
int main()
{
    //freopen("POJ_2391.in","r",stdin);
    int a,b;
    LL c;
    while (scanf("%d%d",&f,&p)!=EOF)
    {
        int cur=0;
        for (int i=1; i<=f; i++)
        {
            scanf("%d%d",&A[i],&B[i]);
            cur+=A[i];
            for(int j=1; j<=f; j++) path[i][j]=INF;
        }
        for (int i=1; i<=p; i++)
        {
            scanf("%d%d%lld",&a,&b,&c);
            path[a][b]=min(path[a][b],c);
            path[b][a]=min(path[b][a],c);
        }
        N=0;
        floyd();
        LL l,r,mid;
        l=1,r=N;
        //cout<<l<<" "<<r<<endl;
        LL ans=-1;
        while(l<r)
        {
            mcmf.init(2*f+10);
            for (int i=1; i<=f; i++)
            {
                mcmf.addedge(0,i,A[i]);
                mcmf.addedge(i,i+f,1<<30);
            }
            for (int i=1; i<=f; i++)
            {
                mcmf.addedge(i+f,2*f+5,B[i]);
            }
            mid=(r+l)>>1;
            for (int i=1;i<=f;i++){
                for (int j=1;j<=f;j++){
                    if (path[i][j]>mid || i==j) continue;
                    mcmf.addedge(i,f+j,1<<30);
                }
            }
            int res=mcmf.maxflow(0,2*f+5);
            if (res>=cur){
                //cout<<res<<" "<<cur<<endl;
                //cout<<mid<<endl;
                ans=mid;
                r=mid;
            }
            else{
                l=mid+1;
            }

        }
        printf("%lld
",ans);
    }
    return 0;
}