poj2391

链接：http://poj.org/problem?id=2391

题解：

二分答案，变成判定性问题，只需把能经过的边在当前图中联通

另外对牛和牛棚要拆点

代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <vector>
#include <queue>
using namespace std;
#define ll long long
const int maxn=1e6;
const int maxn2=;
#define INF 1e16
#define INF2 1e9
int n,m,s,t,l,sum,ax;
int head[maxn2],head2[maxn2],d[maxn2],cow[maxn2],cap[maxn2];
ll dis[maxn2][maxn2];
bool vis[maxn2];
struct re{
    int a,b,c,flow;
}a[maxn];
void arr(int x,int y,int z)
{
    a[++l].a=head[x];
    ax=max(ax,l);
  a[l].b=y;
    a[l].c=z;
    a[l].flow=;
    head[x]=l;
}
queue<int> q;
bool bfs(){
    memset(vis,,sizeof(vis));
    q.push(s);
    d[s]=; vis[s]=;
    while (!q.empty())
    {
        int x=q.front();q.pop();
        int u=head[x];
        while (u)
        {
            int v=a[u].b;
            if (!vis[v]&&a[u].c>a[u].flow)
            {
                vis[v]=;
                d[v]=d[x]+;
                q.push(v);
            }
            u=a[u].a;
        }
    }
    return(vis[t]);
}
int dfs(int x,int y)
{
    if (x==t||y==) return y;
    int flow=,f,tmp;
    int u=head[x];
    while (u)
    {
        int v=a[u].b;
        if (d[x]+==d[v]&&(f=dfs(v,min(y,a[u].c-a[u].flow)))>)
        {
            a[u].flow+=f;
            if (u%) tmp=u+; else tmp=u-;
            a[tmp].flow-=f;
            flow+=f;
            y-=f;
            if (y==) break;
        }
        u=a[u].a;
    }
    head[x]=u;
    return(flow);
}
int maxflow()
{
    int flow=;
    while (bfs())
    {
        flow+=dfs(s,INF2);
        memcpy(head,head2,sizeof(head2));
    }
    return(flow);
}
void floyd()
{
    for (int k=;k<=n;k++)
      for (int i=;i<=n;i++)
        for (int j=;j<=n;j++)
          if (i!=j)
          dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
}
void buildedge(ll x)
{
  memset(head,,sizeof(head));l=;
    for (int i=;i<=n;i++)
    {
        arr(,i,cow[i]); arr(i,,);
        arr(i+n,n*+,cap[i]); arr(n*+,i+n,);
        arr(i,i+n,INF2); arr(i+n,i,);
    }
    for (int i=;i<=n;i++)
      for (int j=;j<=n;j++)
        if (dis[i][j]<=x)
        {
          arr(i,j+n,INF2); arr(j+n,i,);
        }
}
ll solve()
{
    ll ans=-,h=,t=INF-;
    while (h<=t)
    {
        ll mid=(h+t)/;
        buildedge(mid); memcpy(head2,head,sizeof(head));
        if (maxflow()>=sum)
        {
            ans=mid; t=mid-;
        } else h=mid+;
    }
    return ans;
}
int main()
{
    freopen("noip.in","r",stdin);
    freopen("noip.out","w",stdout);
    std::ios::sync_with_stdio(false);
    while (cin>>n>>m)
    {
        sum=; s=; t=*n+;
        for (int i=;i<=n;i++)
        {
            cin>>cow[i]>>cap[i];
            sum+=cow[i];
        }
        for (int i=;i<=n;i++)
          for (int j=;j<=n;j++)
            dis[i][j]=INF;
        ll a,b,c;
        for (int i=;i<=m;i++)
        {
            cin>>a>>b>>c;
            dis[a][b]=min(dis[a][b],c);
            dis[b][a]=min(dis[b][a],c);
        }
        floyd();
        cout<<solve()<<endl;
    }
    return ;
}