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;}

原创粉丝点击