poj2135

来源:互联网 发布:淘宝企业店铺条件 编辑:程序博客网 时间:2024/05/19 22:25

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

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

裸的费用流。

#include<cstdio>#include<cstdlib>#include<cstring>#include<queue>#define INF 99999999using 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);}}