poj2135

题意：给你一个图，让你求一条最小的路径从1到N再回到1。边不能重复走。 

思路：从1到N的最小路径可以直接最短路算法，那么从N到1也可以最短路算法，只要把边取反就好了，网络流中我们在每两个点对应连两条边(正向的，反向的，容量都设为1，表示只能走一次，无向图嘛)，在源点到1连一条边，费用为0，容量为2，表示要走两次。

裸的费用流。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<queue>
#define INF 99999999
using namespace std;
const int N=2003;
struct node{
	int x,y,z,next,other,cost;
}sa[N*N/2];int len=0,first[N];
int n,m,st,ed;
void ins(int x,int y,int z,int v)
{
	len++;
	sa[len].x=x;
	sa[len].y=y;
	sa[len].z=z;
	sa[len].cost=v;
	sa[len].next=first[x];
	sa[len].other=len+1;
	first[x]=len;
	len++;
	sa[len].x=y;
	sa[len].y=x;
	sa[len].z=0;
	sa[len].cost=-v;
	sa[len].next=first[y];
	sa[len].other=len-1;
	first[y]=len;
}
int vis[N],pre[N],dis[N],frpe[N];
bool spfa(int s,int t)
{
	memset(vis,0,sizeof(vis));
	memset(pre,-1,sizeof(pre));
	for(int i=0;i<=t;i++)
	dis[i]=INF;
	pre[s]=0;
	dis[s]=0;
	vis[s]=1;
	queue<int>q;
	q.push(s);
	while(!q.empty())
	{
		
		int x=q.front();
		q.pop();	vis[x]=0;
		for(int i=first[x];i!=-1;i=sa[i].next)
		{
			int y=sa[i].y;
			if(sa[i].z && 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;
			}
			
		}
	
	}//printf("!");
	if(dis[t]>=INF) return false;
	return true;
}
void MCMF(int s,int t)
{
	int ans=0;
	//
	while(spfa(st,ed))
	{//printf("!");
		int minl=INF;
		for(int i=t;i!=0;i=pre[i])
		{
		//	printf("%d\n",i);
		//	system("pause");
			if(minl>sa[frpe[i]].z)
			minl=sa[frpe[i]].z;
		}
		ans+=minl*dis[t];
		for(int i=t;i!=0;i=pre[i])
		{
			sa[frpe[i]].z-=minl;
			sa[sa[frpe[i]].other].z+=minl;
		}
	}
	printf("%d\n",ans);
}
int main()
{
	while(scanf("%d%d",&n,&m)!=EOF)
	{
		memset(first,-1,sizeof(first));
		len=0;
		for(int i=1;i<=m;i++)
		{
			int x,y,z; 
			scanf("%d%d%d",&x,&y,&z);
			ins(x,y,1,z);		
			ins(y,x,1,z);
		}
		st=0;ed=n+1;
		ins(st,1,2,0);
		ins(n,ed,2,0);
		//printf("!");
		MCMF(st,ed);
	}
}