bzoj 3876[AHOI 2014]支線ストーリー(上下界の料金フローあり)

もともと上下のネットの流れをブラシしたいだけだったが、うっかり上下の費用の流れを作ってしまった.
初めてエッジ-infのバージョンを書いたが、infの値が不思議なようだ.どうせ私はだめだ==
湘江からの分析
分析:毎回ストーリーポイント1から、最小限の時間を費やして、すべてのエッジを歩き終えて、1つのエッジを少なくとも1回歩くので、下界が1の費用の流れです.構図:まずスーパーソース点s(0)を確立し、スーパーポイントt(n+1)t->s費は0、流量はINF対辺u->v、費はc、流量下界は1、上界は正無限に変更:s->v費はc、流量は1(下界が1であることを示す)u->v費はc、流量はINF(何度も歩けることを示す)各点u->tに対して0,流量uの出度(各連1本のu->vの辺に等価,1本のu方向tの辺流量をu->vの下限1の最終効果加算)u->1(原図のソース点)に0,流量INF(このソースポイントを除く)ループでi=1の場合、1->1のエッジを接続しないでください.そうしないと、あるデータが死んでしまいます.

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<bitset>
#define LL long long
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define efo(i,x) for(int i=last[x];i!=0;i=e[i].next)
using namespace std;
inline LL read()
{
	LL d=0,f=1;char s=getchar();
	while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
	while(s>='0'&&s<='9'){d=d*10+s-'0';s=getchar();}
	return d*f;
}
#define S 0
#define T n+1
#define N 305
#define inf 1000000
struct edge
{
	int y,c,f,next;
}e[N*N];
int last[N],ne=1;
int n,m,ans=0;

int q[N*9],pre[N],lo[N],dis[N];
bitset<N>inq;
bool spfa(int s,int tt)
{
	inq.reset();
	fo(i,0,N)dis[i]=inf;
	memset(pre,-1,sizeof(pre));
	memset(lo,-1,sizeof(lo));
	int h=0,t=1;
	q[1]=s;dis[s]=0;inq[s]=1;
	while(h<t)
	{
		int now=q[++h];
		efo(i,now)
		if(e[i].f&&dis[now]+e[i].c<dis[e[i].y])
		{
			dis[e[i].y]=dis[now]+e[i].c;
			pre[e[i].y]=now;
			lo[e[i].y]=i;
			if(inq[e[i].y]==0)
			{
				q[++t]=e[i].y;
				inq[e[i].y]=1;
			}
		}
		inq[now]=0;
	}
	if(dis[tt]==inf)return 0;
	else return 1;
}

void mcf(int s,int t)
{
	while(spfa(s,t))
	{
		int x=t,mi=inf;
		while(pre[x]!=-1)
		{
			mi=min(e[lo[x]].f,mi);
			x=pre[x];
		}
		x=t;
		while(pre[x]!=-1)
		{
			ans+=mi*e[lo[x]].c;
			e[lo[x]].f-=mi;
			e[lo[x]^1].f+=mi;
//			cout<<x<<' ';
			x=pre[x];
		}
//		cout<<endl;
	}
}

void add(int x,int y,int f,int c)
{
	e[++ne].y=y;e[ne].f=f;e[ne].c=c;e[ne].next=last[x];last[x]=ne;
}

void add2(int x,int y,int f,int c)
{
	add(x,y,f,c);add(y,x,0,-c);
}

void init()
{
	n=read();
	fo(i,1,n)
	{
		m=read();
		fo(j,1,m)
		{
			int x=read(),y=read();
			add2(S,x,1,y);
			add2(i,x,inf,y);
		}
		add2(i,T,m,0);
		if(i!=1)add2(i,1,inf,0);
	}
	add2(T,S,inf,0);
}

int main()
{
	init();
	mcf(S,T);
	cout<<ans<<endl;
	return 0;
}