1927: [Sdoi2010]星际竞速

题目链接

题解：

每个点拆成两个点，

st->每个点的入点建一条流量为1，费用为0的边

st->每个点的出点建一条流量为1，费用为瞬间移动的思考时间

对于每一条路径，从编号小的点的入点向编号大的点的出点连一条流量为1，费用为花费时间的边

每个点的出点->ed建一条流量为1，费用为0的边

考虑由于求最大流，所以所有点的出点向汇点的连边都会流满。每个点的出点的到达方式即为此点的到达方式：

1.由源点流来：指在某一时刻瞬移到该星球。

2.由入点流来：指由其他星球沿航路到该星球。

至于为什么第2种情况合法，自己画一个图应该就懂了。

代码：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<queue>
using namespace std;
#define chu(i) i+n
#define INF 1e9
const int N=1605;
const int M=31000;
int n,m,st,ed;
struct node{
	int x,y,z,cost,next,other;
}sa[M*2];int len=0,first[N];
void ins(int x,int y,int z,int cost)
{
	len++;
	sa[len].x=x;
	sa[len].y=y;
	sa[len].z=z;
	sa[len].cost=cost;
	sa[len].next=first[x];
	first[x]=len;
	sa[len].other=len+1;
	
	len++;
	sa[len].x=y;
	sa[len].y=x;
	sa[len].z=0;
	sa[len].cost=-cost;
	sa[len].next=first[y];
	first[y]=len;
	sa[len].other=len-1;
}
int pre[N],frpe[N],dis[N];
bool vis[N];
queue<int>q;
bool spfa()
{
	memset(pre,-1,sizeof(pre));pre[st]=0;
	memset(vis,0,sizeof(vis));
	for(int i=0;i<=ed;i++) dis[i]=INF;
	dis[st]=0;
	vis[st]=1;
	q.push(st);
	while(!q.empty())
	{
		int x=q.front();q.pop();
		for(int i=first[x];i!=-1;i=sa[i].next)
		{
			int y=sa[i].y;
			if(sa[i].z>0&&dis[y]>dis[x]+sa[i].cost)
			{
				dis[y]=dis[x]+sa[i].cost;
				if(!vis[y])
				{
					q.push(y);
					vis[y]=1;
				}
				pre[y]=x;
				frpe[y]=i;
			}
		}
		vis[x]=0;
	}
	if(dis[ed]>=INF) return false;
	return true;
}
void MCMF()
{
	int ans=0,liu=0;
	while(spfa())
	{
		int minn=INF;
		for(int i=ed;i!=st;i=pre[i])
		{
			if(minn>sa[frpe[i]].z)
			minn=sa[frpe[i]].z;
		}
		ans+=minn*dis[ed];
		liu+=minn;
		for(int i=ed;i!=st;i=pre[i])
		{
			sa[frpe[i]].z-=minn;
			sa[sa[frpe[i]].other].z+=minn;
		}
	}
	printf("%d\n",ans);
}
int main()
{
	scanf("%d%d",&n,&m);
	st=0;ed=n*2+1;
	int a;
	memset(first,-1,sizeof(first));
	for(int i=1;i<=n;i++)
	{
		scanf("%d",&a);
		ins(st,i,1,0);
		ins(st,chu(i),1,a);
		ins(chu(i),ed,1,0);
	}
	
	int x,y,z;
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d%d",&x,&y,&z);
		if(x>y) swap(x,y);
		ins(x,chu(y),1,z);
	}
	MCMF();
}