bzoj3876 支线剧情【上下界最小费用可行流】

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=3876

解题思路

很明显是求上下界最小费用可行流做。
因为每条边都要走一次，所以将原图中的边下届设为1，费用即为边长；
又因为可以从任何点退出，所以从所有点向汇点连边，无上下界，费用为0；
基本图建立完毕，再建立有超级源，超级汇的图跑最小费用可行流。
记得最后还要加上原图中每条边下届（即是1）乘以其费用，因为网络流中是忽略了下届的费用的。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
#include<vector>
#define ll long long
using namespace std;

int getint() 
{
    int i=0,f=1;char c;
    for(c=getchar();(c<'0'||c>'9')&&c!='-';c=getchar());
    if(c=='-')c=getchar(),f=-1;
    for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
    return i*f;
}

const int N=500,M=20000,INF=0x3f3f3f3f;
int n,ans,s,t,S,T,totflow;
int tot=1,first[N],nxt[M],to[M],cap[M],cost[M];
int in[N],out[N],dis[N];
bool exist[N],visit[N];
queue<int>q;

inline void add(int x,int y,int f,int w)
{
    nxt[++tot]=first[x],first[x]=tot,to[tot]=y,cap[tot]=f,cost[tot]=w;
    nxt[++tot]=first[y],first[y]=tot,to[tot]=x,cap[tot]=0,cost[tot]=-w;
}

bool SPFA()
{
    memset(dis,INF,sizeof(dis));
    memset(visit,false,sizeof(visit));
    dis[S]=0,exist[S]=1;
    q.push(S);
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        exist[u]=false;
        for(int e=first[u];e;e=nxt[e])
        {
            int v=to[e];
            if(dis[v]>dis[u]+cost[e]&&cap[e])
            {
                dis[v]=dis[u]+cost[e];
                if(!exist[v])exist[v]=true,q.push(v);
            }
        }
    }
    return dis[T]!=INF;
}

int dinic(int u,int flow)
{
    if(u==T)
    {
        ans+=dis[T]*flow;
        return flow;
    }
    int res=0;
    visit[u]=true;
    for(int e=first[u];e;e=nxt[e])
    {
        int v=to[e];
        if(dis[v]==dis[u]+cost[e]&&cap[e]&&!visit[v])
        {
            int delta=dinic(v,min(flow-res,cap[e]));
            cap[e]-=delta,cap[e^1]+=delta;
            res+=delta;if(res==flow)break;
        }
    }
    return res;
}

int maxflow()
{
    while(SPFA())dinic(S,INF);
}

int main()
{
    //freopen("lx.in","r",stdin);
    int x,y,z;
    n=getint();
    S=0,s=1,t=n+1,T=n+2;
    for(x=1;x<=n;x++)
    {
        out[x]=getint();
        for(int i=1;i<=out[x];i++)
        {
            y=getint(),z=getint();
            in[y]++;
            add(x,y,INF,z);
            ans+=z;
        }
        add(x,t,INF,0);
    }
    add(t,s,INF,0);
    for(x=1;x<=n;x++)
        if(in[x]>out[x])add(S,x,in[x]-out[x],0);
        else if(in[x]<out[x])add(x,T,out[x]-in[x],0);
    maxflow();
    cout<<ans;
    return 0;
}