[AHOI2014]支线剧情 bzoj3876

最新推荐文章于 2018-12-01 17:00:37 发布

原创最新推荐文章于 2018-12-01 17:00:37 发布 · 515 阅读

0 ·

CC 4.0 BY-SA版权

本文介绍了一种基于费用流算法的实现方法，通过使用裸的上下界费用流解决特定类型的网络流问题。该方法利用SPFA算法进行寻路，并通过不断调整路径上的流量来优化总费用。

裸的上下界费用流

#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=50000,inf=1<<30;
int head[maxn],dist[maxn],s,t,prev[maxn],pree[maxn],ans,n,que[maxn];
bool inq[maxn];
struct edge{
	int to,next,c,f;
}g[maxn];
int num=1;
void adde(int u,int v,int f,int c){
	g[++num].to=v;g[num].f=f;g[num].c=c; g[num].next=head[u];head[u]=num;
	g[++num].to=u;g[num].f=0;g[num].c=-c;g[num].next=head[v];head[v]=num;
}
bool spfa(){
	for(int i=0;i<=n+1;i++) dist[i]=inf;
	dist[s]=0;
	que[1]=s;
	inq[s]=1;
	int he=1,ta=1;
	while(he<=ta){
		int x=que[he++];
		for(int i=head[x];i;i=g[i].next){
			int node=g[i].to;
			if(g[i].f&&dist[node]>dist[x]+g[i].c){
				dist[node]=dist[x]+g[i].c;
				prev[node]=x;pree[node]=i;
				if(!inq[node]){
					inq[node]=1;
					que[++ta]=node;
				}
			}
		}
		inq[x]=0;
	}
	return dist[t]<inf;
}
void augment(){
	int node=t,flow=inf;
	while(node!=s){
		flow=min(flow,g[pree[node]].f);
		node=prev[node];
	}
	node=t;
	while(node!=s){
		g[pree[node]].f-=flow;
		g[pree[node]^1].f+=flow;
		node=prev[node];
	}
	int ret=dist[t]*flow;
	ans+=ret;
}
int main(){
	scanf("%d",&n);
	s=0,t=n+1;
	for(int i=1;i<=n;i++){
		int m;
		scanf("%d",&m);
		for(int j=1;j<=m;j++){
			int u,c;
			scanf("%d%d",&u,&c);
			adde(i,u,inf,c);
			adde(s,u,1,c);
		}
		if(i!=1) adde(i,1,inf,0);
		adde(i,t,m,0);
	}
	/*for(int i=0;i<=n+1;i++) {
		for(int j=head[i];j;j=g[j].next){
			printf("to-%d:f-%d:c-%d   ",g[j].to,g[j].f,g[j].c);
		}
		printf("\n");
	}*/
	//spfa();
	//printf("%d",dist[t]);
	//int node=t;
	//while(node!=s) printf("%d",node),node=prev[node];
	while(spfa()) augment();
	//spfa();augment();
	printf("%d",ans);
	return 0;
}