hdu 1853/3488/3435 最小费用最大流 环覆盖图

N个城市，M条单向路，用若干个环覆盖整个图，求环的最小权值

s向u连一条流为1费用为0的边 u+n向t连一条流为1费用为0的边 u向u+n连一条流为1费用为w的边

进行最小费用最大流运算 如果流等于节点数 说明是可以完全覆盖的 求出来的就是最小代价

因为我们可以这样考虑 如果用环覆盖整个图的话 所有点的入度出度都是2 

然后在上面费用流的计算中 也是这样的 所以完成问题的转换

但是这个方案的速度比较慢 更快的方案是km匹配

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;
int sumFlow;
const int MAXN=502;
const int MAXM=10002;
const int INF=1000000000;
struct Edge
{
    int u,v,cap,cost,next;
}edge[MAXM<<2];
int NE;
int head[MAXN],dist[MAXN],pp[MAXN];
bool vis[MAXN];
void init()
{
    NE=0;
    memset(head,-1,sizeof(head));
}
void ae(int u,int v,int cap,int cost)
{
    edge[NE].u=u;edge[NE].v=v;edge[NE].cap=cap;edge[NE].cost=cost;
    edge[NE].next=head[u];head[u]=NE++;

    edge[NE].u=v;edge[NE].v=u;edge[NE].cap=0;edge[NE].cost=-cost;
    edge[NE].next=head[v];head[v]=NE++;
}
bool SPFA(int s,int t,int n)
{
    int i,u,v;
    queue<int>q;
    memset(vis,0,sizeof(vis));
    memset(pp,-1,sizeof(pp));
    for(i=0;i<=n;i++) dist[i]=INF;
    vis[s]=1;dist[s]=0;
    q.push(s);
    while(!q.empty())
    {
        u=q.front();q.pop();vis[u]=0;
        for(i=head[u];i!=-1;i=edge[i].next)
        {
            v=edge[i].v;
            if(edge[i].cap&&dist[v]>dist[u]+edge[i].cost)
            {
                dist[v]=dist[u]+edge[i].cost;
                pp[v]=i;
                if(!vis[v])
                {
                    q.push(v);
                    vis[v]=1;
                }
            }
        }
    }
    if(dist[t]==INF) return 0;
    return 1;
}
pair<int,int> MCMF(int s,int t,int n)
{
    int flow=0,i,minflow,mincost=0;
    while(SPFA(s,t,n))
    {
        minflow=INF;
        for(i=pp[t];i!=-1;i=pp[edge[i].u])
            minflow=min(edge[i].cap,minflow);
        flow+=minflow;
        for(i=pp[t];i!=-1;i=pp[edge[i].u])
            edge[i].cap-=minflow,edge[i^1].cap+=minflow;
        mincost+=dist[t]*minflow;
    }
    return make_pair(flow,mincost);
}
int main()
{
    int n,m,u,v,c;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        init();
        int vs=0,vt=n+n+1,NV=vt+1;
        for(int i=1;i<=n;i++)
        {
            ae(vs,i,1,0);
            ae(i+n,vt,1,0);
        }
        while(m--)
        {
            scanf("%d%d%d",&u,&v,&c);
            ae(u,v+n,1,c);
        }
        pair<int,int>ans;
        ans=MCMF(vs,vt,NV);
        if(ans.first!=n) printf("-1\n");
        else printf("%d\n",ans.second);
    }
    return 0;
}