BZOJ3876 || 洛谷P4043 [AHOI2014/JSOI2014]支线剧情【有上下界最小费用流】

Time Limit: 10 Sec
Memory Limit: 256 MB

Description

JYY现在所玩的RPG游戏中，一共有N个剧情点，由1到N编号，第i个剧情点可以根据JYY的不同的选择，而经过不同的支线剧情，前往Ki种不同的新的剧情点。当然如果为0，则说明i号剧情点是游戏的一个结局了。
JYY观看一个支线剧情需要一定的时间。JYY一开始处在1号剧情点，也就是游戏的开始。显然任何一个剧情点都是从1号剧情点可达的。此外，随着游戏的进行，剧情是不可逆的。所以游戏保证从任意剧情点出发，都不能再回到这个剧情点。由于JYY过度使用修改器，导致游戏的“存档”和“读档”功能损坏了，
所以JYY要想回到之前的剧情点，唯一的方法就是退出当前游戏，并开始新的游戏，也就是回到1号剧情点。JYY可以在任何时刻退出游戏并重新开始。不断开始新的游戏重复观看已经看过的剧情是很痛苦，JYY希望花费最少的时间，看完所有不同的支线剧情。

Input

输入一行包含一个正整数N。
接下来N行，第i行为i号剧情点的信息；
第一个整数为，接下来个整数对，Bij和Tij，表示从剧情点i可以前往剧
情点，并且观看这段支线剧情需要花费的时间。

Output

输出一行包含一个整数，表示JYY看完所有支线剧情所需要的最少时间

---

题目分析

对于每条剧情点$i$到$j$的支线剧情，连边$i,j$，下界1，上界inf，费用为给定时间
每个剧情点 i 向虚拟汇点连边，下界0，上界inf，费用0
并设1号剧情点为虚拟源点

按有上下界网络流的套路构造网络，求$ss$(超源)到$tt$(超汇)的最小费用流
注意此时求出的费用是 假设每条边已有大小等于下界的流量 后补流的费用
此题中是假设每个支线剧情都已观看了一遍，所以还要加上每个支线剧情所需时间的总和

---

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
using namespace std;
typedef long long lt;

int read()
{
    int f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
    while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
    return f*x;
}

const int inf=1128481603;
const int maxn=100010;
int n,m;
struct node{int v,f,c,nxt;}E[maxn<<1];
int head[maxn],tot=1;
int dis[maxn],incf[maxn];
int vis[maxn],pre[maxn];
int deg[maxn];

void add(int u,int v,int f,int c)
{
    E[++tot].nxt=head[u];
    E[tot].v=v; E[tot].f=f; E[tot].c=c;
    head[u]=tot;
}

bool bfs(int s,int t)
{
    memset(dis,67,sizeof(dis)); dis[s]=0;
    queue<int> q; q.push(s); incf[s]=inf;

    while(!q.empty())
    {
        int u=q.front();
        q.pop(); vis[u]=0;
        for(int i=head[u];i;i=E[i].nxt)
        {
            int v=E[i].v;
            if(E[i].f&&dis[v]>dis[u]+E[i].c)
            {
                dis[v]=dis[u]+E[i].c;
                incf[v]=min(incf[u],E[i].f); pre[v]=i;
                if(!vis[v]) q.push(v),vis[v]=1;
            }
        }
    }
    return dis[t]!=inf;
}

int dfs(int s,int t)
{
    int u=t;
    while(u!=s)
    {
        int i=pre[u];
        E[i].f-=incf[t]; E[i^1].f+=incf[t];
        u=E[i^1].v;
    }
    return incf[t]*dis[t];
}

int dicnic(int s,int t)
{
    int fee=0;
    while(bfs(s,t)) fee+=dfs(s,t);
    return fee;
}

int main()
{
    n=read(); int ans=0;
    for(int u=1;u<=n;++u)
    {
    	int k=read();
    	while(k--)
    	{
    		int v=read(),c=read(); ans+=c;
    		add(u,v,inf,c); add(v,u,0,-c);//l=1,r=inf
    		deg[u]--; deg[v]++;
        }
    }
    
    int s=1,t=n+1;
    add(t,s,inf,0); add(s,t,0,0);
    for(int i=1;i<=n;++i)
    add(i,t,inf,0),add(t,i,0,0);
    
    int ss=n+2,tt=ss+1;
    for(int i=1;i<=n;++i)
    {
        if(deg[i]>0) add(ss,i,deg[i],0),add(i,ss,0,0);
        else if(deg[i]<0) add(i,tt,-deg[i],0),add(tt,i,0,0);
    }
    
    printf("%d",ans+dicnic(ss,tt));
    return 0;
}