poj 2135 Farm Tour
来源:互联网 发布:软件著作权 发明专利 编辑:程序博客网 时间:2024/05/16 07:02
无向图最小费用流:因为每条边只能经过一次,所以每条边的流量为1,逆向边的流量为0,在用SPFA松弛路径的时候检查流量是否大于0,再作进一步的判断是否更新源点到当前点的距离,同时因为该图为无向图,建边要注意正反都要建,假设在1~n的过程中走过a-b这一条费用为正的边,那么从n~1的过程中一定不会再走b-a这条费用为正的边,因为还存在一条a-b的逆向边b-a的费用为负值,更新的的时候它一定会取这条负值边而不会取b-a的正值边,如果他取了这条负值边,就相当于没有流量流过这条边,综上所诉,可以保证每条边只经过了一次。#include<queue>#include<stdio.h>#include<string.h>#include<iostream>#define N 1010#define M 10010#define inf 1<<30using namespace std;struct node{int u,v,cost,flag,next;}edge[4*M];int pre[M],vis[N],head[M],dis[N];int n,m,u,v,c,e,ans,s,t;void add(int u,int v,int c,int f){ edge[e].u=u; edge[e].v=v; edge[e].cost=c; edge[e].flag=f; edge[e].next=head[u]; head[u]=e++; edge[e].u=v; edge[e].v=u; edge[e].cost=-c; edge[e].flag=0; edge[e].next=head[v]; head[v]=e++;}void spfa(){ memset(pre,-1,sizeof(pre)); memset(vis,0,sizeof(vis)); for(int i=0;i<=t;i++) dis[i]=inf; dis[0]=0; queue<int>q; q.push(0); vis[0]=1; while(!q.empty()) { u=q.front(); q.pop(); vis[u]=0; for(int i=head[u];i!=-1;i=edge[i].next) { int v=edge[i].v; if(edge[i].flag&&dis[v]>dis[u]+edge[i].cost) { dis[v]=dis[u]+edge[i].cost; pre[v]=i; if(!vis[v]) { q.push(v); vis[v]=1; } } } }}int work(){ for(int i=0;i<2;i++) { spfa(); ans+=dis[t]; for(int j=pre[t];j!=-1;j=pre[edge[j].u]) { edge[j].flag-=1; edge[j^1].flag+=1; } } return ans;}int main(){ while(scanf("%d%d",&n,&m)!=EOF) { s=0,t=n+1,e=0,ans=0; memset(head,-1,sizeof(head)); while(m--) { scanf("%d%d%d",&u,&v,&c); add(u,v,c,1); add(v,u,c,1); } add(0,1,0,2); add(n,t,0,2); printf("%d\n",work()); } return 0;}