poj2931 Ombrophobic Bovines dinic

坑点：答案很大会爆int  细节点很多，二分判断最长路径的范围 ，再把每个点拆点，一个记现有的牛，一个记可以住的牛。

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <queue>
#include <cmath>
using namespace std;
typedef long long sint;
#define INF  0x7fffffff
#define fuck 0x7fffffffffffffLL
#define maxn 3000
struct node
{
    int v,next;
    int c;
}e[222222];
int st,ed,en;
int q[maxn*100];
int first[maxn],dis[maxn];
int n,m,f,d,sum;
sint save[800][800],aa[800],bb[800];
void add(int a,int b,sint c)
{
    e[en].c=c;
    e[en].v=b;
    e[en].next=first[a];
    first[a]=en;
    en++;
    e[en].c=0;
    e[en].v=a;
    e[en].next=first[b];
    first[b]=en;
    en++;
}
bool bfs()
{
    int head=1;
    int tail=1;
    memset(dis,-1,sizeof(dis));
    dis[st]=0;
    tail++;
    q[tail]=st;
    int u,v;
    while(head<tail)
    {
        head++;
        u=q[head];
        for(int i=first[u];i!=-1;i=e[i].next)
        {
            v=e[i].v;
            if(dis[v]==-1&&e[i].c)
            {
                dis[v]=dis[u]+1;
                tail++;
                q[tail]=v;
                if(v==ed)
                {
                    return true;
                }
            }
        }
    }
    return false;
}
void init()
{
    en=0;
    st=0;
    ed=2*n+1;
    memset(first,-1,sizeof(first));
}
int dfs(int x,int mx)
{
    int i,a;
    if(x==ed||mx==0) return mx ;
    int ret=0,f,flow=0;
    for(int k=first[x];k!=-1&&ret<mx;k=e[k].next)
    {
        if(e[k].c&&dis[e[k].v]==dis[x]+1&&(f=dfs(e[k].v,min(e[k].c,mx))))
        {
            e[k].c-=f;
            e[k^1].c+=f;
            mx-=f;
            ret+=f;
        }
    }
    if(!ret) dis[x]=-1;
    return ret;
}
int dinic()
{
    int tmp=0;
    int maxflow=0;
    while(bfs())
    {
        while(tmp=dfs(st,INF)) maxflow+=tmp;
    }
    return maxflow;
}
void build(long long maxs)
{
	init();
	for(int i=1;i<=n;i++)
	{
		add(st,i,aa[i]);
		add(i+n,ed,bb[i]);
		add(i,i+n,INF);
	}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			if(save[i][j]<=maxs)
			{
				add(i,j+n,INF);
			}
}
void solve()
{
    sum=0;
    sint ll,rr,ans=-1,mid;
    for(int i=1;i<=n;i++)
    {
        scanf("%d%d",&aa[i],&bb[i]);
        sum+=aa[i];
    }
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
            {
                save[i][j]=fuck+1;
            }
        }
        int a,b;
        sint c;
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d%lld",&a,&b,&c);
            save[a][b]=save[b][a]=min(save[a][b],c);
        }
    for(int k=1;k<=n;k++)
    {
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
            {
                if(i==j||i==k||j==k) continue;
                save[i][j]=min(save[i][j],save[i][k]+save[k][j]);
            }
        }
    }
    ll=0;
    rr=fuck;
    while(ll<=rr)
    {
        mid=(ll+rr)>>1;
        build(mid);
        if(dinic()==sum)
        {
            ans=mid;
            rr=mid-1;
        }
        else
        {
            ll=mid+1;
        }
    }
    cout<<ans<<endl;
}
int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        st=n*2+1;
        ed=n*2+2;
        solve();
    }
    return 0;
}