CodeVS2133星际竞速

题目

　　http://codevs.cn/problem/2133/

题解

　　为何SDOI的题总是让人感觉如此舒爽！

　　言归正传，这道题的确让我获益匪浅。把每个点拆成两个，假如是Xi和Yi，设能力爆发的费用是Bi，任意边为(u,v,w)，那么连边：(S,Xi,1,0)(S,Yi,1,Bi)(Xu,YV,1,w)最后跑一遍最小费用最大流即可。

　　解释：题目规定，每个点只能访问一次，那么Yi的作用就是这个，一个点i被访问的充要条件是由一个大小为1的流流经Yi并流进了T。那么一个点被访问有两种方式：一种是能力爆发，即边(S,Yi,1,Bi)；另一种是通过边来访问，即边(Xu,Yv,1,w)，这里的c=1是没必要的，只是为了统一。其实最开始的想法应该是把所有的边建在左边然后连向对应的终点，但是为了节省节点数，就在左边拆了一排Xi这样的点，这道题的拆点和平常所说的拆点限流没有丝毫关系。

代码

//费用流 
#include <cstdio>
#include <algorithm>
#include <queue>
#define maxn 1000000
#define inf 0x3f3f3f3f
using namespace std;
int N, M, d[maxn], in[maxn], head[maxn], next[maxn], c[maxn], w[maxn], to[maxn], S, T, pre[maxn],
	cost, tot=1;
queue<int> q;
void adde(int a, int b, int cc, int ww)
{to[++tot]=b;c[tot]=cc;w[tot]=ww;next[tot]=head[a];head[a]=tot;}
void adde2(int a, int b, int cc, int ww){adde(a,b,cc,ww);adde(b,a,0,-ww);}
bool spfa()
{
	int p, x;
	for(int i=1;i<=T;i++)d[i]=inf;
	in[S]=1;d[S]=0;q.push(S);
	while(!q.empty())
	{
		in[x=q.front()]=0;q.pop();
		for(p=head[x];p;p=next[p])
			if(c[p] and d[to[p]]>d[x]+w[p])
			{
				d[to[p]]=d[x]+w[p];
				pre[to[p]]=p;
				if(!in[to[p]])in[to[p]]=1,q.push(to[p]);
			}
	}
	return d[T]!=inf;
}
int augment()
{
	int f=inf, x;
	for(x=T;x!=S;x=to[pre[x] xor 1])f=min(f,c[pre[x]]);
	for(x=T;x!=S;x=to[pre[x] xor 1])c[pre[x]]-=f, c[pre[x] xor 1]+=f;
	cost+=f*d[T];
	return f;
}
void build()
{
	int i, x, y, z;
	scanf("%d%d",&N,&M);
	S=N*2+1, T=S+1;
	for(i=1;i<=N;i++)scanf("%d",&x),adde2(S,N+i,1,x),adde2(S,i,1,0),adde2(N+i,T,1,0);
	for(i=1;i<=M;i++)scanf("%d%d%d",&x,&y,&z),adde2(x<y?x:y,(x<y?y:x)+N,1,z);
}
int main()
{
	build();
	while(spfa())augment();
	printf("%d",cost);
	return 0;
}