bzoj1877（费用流）

Elaxia每天从寝室出发 跑到学校，保证寝室编号为1，学校编号为N。 Elaxia的晨跑计划是按周期（包含若干天）进行的，由于他不喜欢走重复的路线，所以 在一个周期内，每天的晨跑路线都不会相交（在十字路口处），寝室和学校不算十字路 口。Elaxia耐力不太好，他希望在一个周期内跑的路程尽量短，但是又希望训练周期包含的天 数尽量长。 

题目说的还是挺委婉

1：不走重复路线，一条边只能走一次

2：希望包含的天长，流量最大

3：尽量跑得路程短，每一条边路程作为费用，路程少就是让费用少

4：在十字路口处不能相交，其实就是说点只能走一次，所以要分点

裸的费用流

就是中间的点要拆点，连边让cap为1，限制只能流过去一次

原点和汇点不用拆点

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<cmath>
#include<queue>
#include<iostream>
using namespace std;
const int inf=0x3f3f3f3f;
const int N=1009;
int n,m,id1[209],id2[209],ansf,ansc,head[N],tot,s,t;
struct aa
{
	int pre,to,cap,flow,w;
}edge[100009];
void addedge(int x,int y,int z,int w)
{
	edge[++tot].to=y;edge[tot].cap=z;edge[tot].w=w;edge[tot].pre=head[x];head[x]=tot;
	edge[++tot].to=x;edge[tot].cap=0;edge[tot].w=-w;edge[tot].pre=head[y];head[y]=tot;
}
int pre[N],now[N],dis[N];
bool b[N];
bool spfa()
{
	memset(dis,inf,sizeof(dis));
	memset(b,false,sizeof(b));
	memset(pre,0,sizeof(pre));
	memset(now,0,sizeof(now));
	dis[s]=0;
	queue<int> q;q.push(s);
	while (!q.empty())
	{
		int u=q.front();q.pop();
		b[u]=false;
		for (int i=head[u];i;i=edge[i].pre)
		if (edge[i].cap>edge[i].flow&&dis[edge[i].to]>dis[u]+edge[i].w)
		{
			dis[edge[i].to]=dis[u]+edge[i].w;
			pre[edge[i].to]=u;now[edge[i].to]=i;
			if (!b[edge[i].to]) {b[edge[i].to]=true;q.push(edge[i].to);}
		}
	}
	if (dis[t]!=inf) return true;
	return false;
}
void work()
{
	while (spfa())
	{
		int flow=inf;
		for (int i=t;i!=s;i=pre[i])
			flow=min(flow,edge[now[i]].cap-edge[now[i]].flow);
		ansf+=flow;
		ansc+=flow*dis[t];
		for (int i=t;i!=s;i=pre[i]) 
		edge[now[i]].flow+=flow,edge[((now[i]-1)^1)+1].flow-=flow;
	}	
}
int main()
{
	scanf("%d%d",&n,&m);
	int bj=0,total=0;
	for (int i=1;i<=n;i++) 
	id1[i]=++bj,id2[i]=++bj,addedge(id1[i],id2[i],1,0);
	int x,y,z;
	for (int i=1;i<=m;i++)
	{
		scanf("%d%d%d",&x,&y,&z);
		addedge(id2[x],id1[y],1,z);
	}
	s=id2[1];t=id1[n];
	work();
	printf("%d %d",ansf,ansc);
	return 0;
}

总结

1：题目中很多，看似普通的话语就是题目数学建模的约束，要非常注意