poj2135 - Farm Tour

                               想看更多的解题报告:http://blog.csdn.net/wangjian8006/article/details/7870410
                                转载请注明出处：http://blog.csdn.net/wangjian8006

 

题目大意：FJ带朋友参观自己的农场，从自己的房子出发到农场，再从农场返回自己的房子，要求去回不走同一条路。房子的点数为1，农场为n，在1到n之间有很多点，给出n个顶点，m条边，然后m行每行有三个数，a，b，c代表a到c的路径长度为c，并且a到b是无向边，现在要求从1点到n点在从n点返回1点的最短路

解题思路：其实很好建图，建一个源点到1加一条有向边，路径长度为0，容量为2，其他边建无向边路径长度为c，容量为1，再从n点到汇点加一条有向边，路径为0，容量为2，然后从源点到汇点求最小费用最大流，当然最小费用在图中表示就是最小路径

/*
Memory 968K
Time    0MS
*/
#include <iostream>
#include <queue>
using namespace std;
#define INF INT_MAX
#define min(a,b) (a>b?b:a)
#define MAXV 1100
#define MAXM 40100

typedef struct{
	int s,t,next,w,r;
}Edge;

Edge edge[MAXM];
int source,sink,n,m,mincost,edgesum;
int head[MAXV],d[MAXV],parent[MAXV];

void addedge(int a,int b,int c,int r){
	edge[edgesum].s=a;
	edge[edgesum].t=b;
	edge[edgesum].r=r;
	edge[edgesum].w=c;
	edge[edgesum].next=head[a];
	head[a]=edgesum++;

	edge[edgesum].s=b;
	edge[edgesum].t=a;
	edge[edgesum].r=0;
	edge[edgesum].w=-c;
	edge[edgesum].next=head[b];
	head[b]=edgesum++;
}

int spfa(){
	queue <int>q;
	int v,i,tmp;
	bool vis[MAXV];
	
	for(i=0;i<=sink;i++) d[i]=INF;
	memset(vis,false,sizeof(vis));
	memset(parent,-1,sizeof(parent));

	q.push(source);
	vis[source]=true;
	d[source]=0;
	while(!q.empty()){
		v=q.front();q.pop();
		vis[v]=false;

		for(i=head[v];i!=-1;i=edge[i].next){
			tmp=edge[i].t;
			if(edge[i].r && edge[i].w+d[v]<d[tmp]){
				d[tmp]=edge[i].w+d[v];
				parent[tmp]=i;
				if(!vis[tmp]){
					q.push(tmp);
					vis[tmp]=true;
				}
			}
		}
	}
	return 0;
}

void MCMF(){
	int u;
	mincost=0;
	while(1){
		spfa();
		if(parent[sink]==-1) break;

		u=parent[sink];		//这里不用求出增光路径上的最小流量
        while(u!=-1){		//因为最小流量一定为1
			edge[u].r--;
			edge[u^1].r++;
            u=parent[edge[u].s];
        }
        mincost+=d[sink];
	}
}

int main(){
	int i,a,b,c;
	while(~scanf("%d%d",&n,&m)){
		memset(head,-1,sizeof(head));
		edgesum=0,source=0,sink=n+1;
		for(i=1;i<=m;i++){
			scanf("%d%d%d",&a,&b,&c);
			addedge(a,b,c,1);
			addedge(b,a,c,1);
		}
		addedge(source,1,0,2);
		addedge(n,sink,0,2);

		MCMF();
		printf("%d\n",mincost);
	}
	return 0;
}