BZOJ 3876(支线剧情) 有下界费用流

Description

【故事背景】

宅男JYY非常喜欢玩RPG游戏，比如仙剑，轩辕剑等等。不过JYY喜欢的并不是战斗场景，而是类似电视剧一般的充满恩怨情仇的剧情。这些游戏往往

都有很多的支线剧情，现在JYY想花费最少的时间看完所有的支线剧情。

【问题描述】

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看完所有支线剧情所需要的最少时间。

Sample Input

6
2 2 1 3 2
2 4 3 5 4
2 5 5 6 6
0
0
0

Sample Output

24

HINT

 JYY需要重新开始3次游戏，加上一开始的一次游戏，4次游戏的进程是

1->2->4，1->2->5，1->3->5和1->3->6。

对于100%的数据满足N<=300,0<=Ki<=50,1<=Tij<=300,Sigma(Ki)<=5000

这是一道有下限的费用流。因为有下限，建图就有所改变。
首先，增加一个源点S和一个汇点T。对于当前点x，将它与原来的起点相连，费用为0，容量为inf，这样就修改了源点， 将它与当前的汇点T相连，费用为0，容量为x的出度(保证跑完所有支路)。对于x到达的每个点，将x与它相连，费用为读入的费用，容量为inf(可以重复跑) ，再将现在的源点S与它相连，费用为读入费用，容量为1(反向弧保证下界)。这样图就建完了。 

#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
using namespace std;

#define maxn (300+10)
int N,S,T,inf=1000000;
int head[maxn],dep[maxn],p[maxn],vis[maxn];

struct node{
	int u,v,w,c,flow,next;
}e[100000];

int k=0;
void adde(int u,int v,int w,int c){
	e[k].u=u; e[k].v=v; e[k].w=w; e[k].c=c;
	e[k].flow=0; e[k].next=head[u]; head[u]=k++;
	
	e[k].u=v; e[k].v=u; e[k].w=-w; e[k].c=0;
	e[k].flow=0; e[k].next=head[v]; head[v]=k++;
}

void input(){
	memset(head,-1,sizeof(head));
	scanf("%d",&N);
	S=0; T=N+1;
	int x,a,b;
	for(int i=1;i<=N;i++){
		scanf("%d",&x);
		for(int j=1;j<=x;j++){
			scanf("%d%d",&a,&b);
			adde(S,a,b,1);
			adde(i,a,b,inf);
		}
		adde(i,T,0,x);
		adde(i,1,0,inf);
	}
	
}

void solve(){
	queue<int>q;
	int ans=0;
	while(1){
		for(int i=0;i<=T;i++)dep[i]=inf;
		q.push(S);
		dep[S]=0;
		while(!q.empty()){
			int u=q.front();
			q.pop(); vis[u]=0;
			
			for(int i=head[u];i!=-1;i=e[i].next){
				int v=e[i].v;
				if(dep[v]>dep[u]+e[i].w&&e[i].c>e[i].flow){
					dep[v]=dep[u]+e[i].w;
					p[v]=i;
					if(vis[v]==0){
						vis[v]=1; q.push(v);
					}
				}
			}
		}
		
		if(dep[T]==inf)break;
		int tt=T,a=inf;
		while(tt!=0){
			a=min(a,e[p[tt]].c-e[p[tt]].flow);
			tt=e[p[tt]].u;
		}
		ans+=dep[T]*a;
		tt=T;
		while(tt!=0){
			e[p[tt]].flow+=a;
			e[p[tt]^1].flow-=a;
			tt=e[p[tt]].u;
		}
	}
	printf("%d",ans);
}

int main(){
	input();
	solve();
	return 0;
}