hdu 1853 Cyclic Tour

最小費用流問題：

題意：找出若干個環覆蓋所有的點，使得總的花費最小
因爲每個點只能經過一次，又要使花費最少，所以很快就可想到拆點求最小費用流
因爲要形成環，所以每個點的入度和出度都不能爲0，每一個點的出度都要對應一個點的入度。
把每個點分成兩部分，一部分控制入度，一部分控制出度。
建立超級源點和超級匯點。
讓超級源點和每一個點的出度相連，容量爲1，費用爲0，這保證了，一個點的出度只能對應一個點的入度。
同理，讓超級匯點和每一個點的入度相連。
建好圖後，運行最小費用流的算法即可。




#include<queue>   
#include<stdio.h>   
#include<iostream>   
using namespace std;  
#define N 300
#define M 50000   
#define inf 10000000   
struct node{    
int u,v,cost,flag,next;    
}edge[M];    
int pre[M],vis[N],head[N],dis[N];    
int n,m,e,ans,u,v,c,s,t;  
  
void add_edge(int u,int v,int f,int c)    
{    
     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=s;i<=t;i++)    
        dis[i]=inf;    
     dis[s]=0;    
     queue<int>q;    
     q.push(s);    
     vis[s]=1;      
     while(!q.empty())    
     {      
        int 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<n;i++)    
   {    
      spfa(); 
      if(dis[t]==inf) return 0; 
      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)
    {
         ans=0; s=0; t=2*n+1;e=0; 
         memset(head,-1,sizeof(head));
         for(int i=0;i<m;i++)
         {
            scanf("%d%d%d",&u,&v,&c);
            add_edge(u,v+n,1,c);
         }
         for(int i=1;i<=n;i++)
         {
            add_edge(0,i,1,0);
            add_edge(i+n,2*n+1,1,0);
         }
         if(work())
            printf("%d\n",ans);
         else
            printf("-1\n");
    }
    return 0;
}