【POJ2135】Farm Tour 最小费用最大流

题目大意：给你n个点（n <= 1000）和m条边（m <= 10000），现在从点1出发经过一些边后到点n，然后再经过一些边回到点1，每一条边只能走一次，每条边花费的时间不同。求最小路径。

解题思路：最小费用最大流。加一个源点S和一个汇点T，源点S与点1连一条容量为2费用为0点边，点n与汇点T也连一条容量为2费用为0的边，对于题目给定的那些边。每条边的容量为1，费用为每条边花费的时间。注意这是双向边。然后求最小费用。

需要注意的地方：

一、由于具有不同费用的平行边，所以，必须使用邻接表来存，而不能使用邻接矩阵

二、图是双向的，因为此题是无向图，所以建边的时候如果建两条费用都是正的边的话，退流时无法修正费用

所以应该建4条边：

第一对：

a->b cost 1

b->a -cost 0

第二对：

b->a cost 1

a->b -cost 0

 

 

#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<climits> #include<queue> #define CLR(arr, val) memset(arr,val,sizeof(arr)) using namespace std; const int INF=100000000; template<int MaxV,int MaxE> class MinCostMaxFlow { public: MinCostMaxFlow():top(0) { CLR(head,-1); CLR(next,-1); } void Add1(int u,int v,int tcap,int tcost) { cap[top]=tcap; tonum[top]=v; next[top]=head[u]; head[u]=top; cost[top]=tcost; ++top; cap[top]=0; tonum[top]=u; next[top]=head[v]; head[v]=top; cost[top]=-tcost; ++top; } void Add2(int u,int v,int tcap,int tcost) { Add1(u,v,tcap,tcost); Add1(v,u,tcap,tcost); } int GetMaxFlow(int s,int t,int n) { int max_flow=0; while(SPFA(s, t, n)) { int cur=t; int minc=INF; while(cur!=s) { minc=min(minc,cap[pre_edge[cur]]-flow[pre_edge[cur]]); cur=pre[cur]; } max_flow+=minc; cur=t; while(cur!=s) { flow[pre_edge[cur]]+=minc; flow[rev(pre_edge[cur])]-=minc; cur=pre[cur]; } } return max_flow; } int CurFlowMinCost()const { int mcost=0; for(int i=0;i!=top;i++) if(flow[i]>0) mcost+=flow[i]*cost[i]; return mcost; } private: inline int rev(int n) { return n+1 - 2* (n&1); } bool SPFA(int s,int t,int n) { CLR(pre,-1); fill(len,len+n,INF); queue<int> Q; Q.push(s); len[s]=0; int cur; while(!Q.empty()) { cur=Q.front(); Q.pop(); for(int i=head[cur];i!=-1;i=next[i]) { if(len[cur]+cost[i]<len[tonum[i]] && cap[i]-flow[i]>0) { pre[tonum[i]]=cur; pre_edge[tonum[i]]=i; len[tonum[i]]=len[cur]+cost[i]; Q.push(tonum[i]); } } } return pre[t]!=-1; } int pre[MaxV]; int pre_edge[MaxV]; int q[MaxV]; int len[MaxV]; int top; int head[MaxV]; int next[MaxE]; int tonum[MaxE]; int cap[MaxE]; int flow[MaxE]; int cost[MaxE]; }; MinCostMaxFlow<1010,40010> g; int main() { //freopen("in.txt","r",stdin); int v,e,a,b,c; scanf("%d%d",&v,&e); int s=0,t=v+1; for(int i=0;i!=e;i++) { scanf("%d%d%d",&a,&b,&c); g.Add2(a,b,1,c); } g.Add1(s,1,2,0); g.Add1(v,t,2,0); g.GetMaxFlow(s,t,t+1); cout<<g.CurFlowMinCost()<<endl; }