HDU 3488 最小费用流

时间:2023-02-03 07:41:59

点击打开链接

题意:在给出的有向图中,找到环,环的定义为头和尾出现两次,中间其他/它点出现一次,然后每个点属于一个环,问最小的费用

思路:因为肯定有一个方案满足条件,那么也就是每个点都会出发并且回到了自己,对于一个完备匹配的二分图来说,随意找一个起点,那么最后一定会回到这个起点对应的那个右部的点,所以我们跑最小费用流时流量设为点的个数就可以了

#include <queue>
#include <vector>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
#include <functional>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3fll;
const int maxn=410;
typedef pair<int,int> P;
struct edge{
    int to,cap,rev,cost;
    edge();
    edge(int a,int b,int c,int d){to=a,cap=b,cost=c,rev=d;};
};
vector<edge>G[maxn];
int h[maxn],dis[maxn];
int prevv[maxn],preve[maxn];
void add_edge(int st,int en,int cap,int cost){
    G[st].push_back(edge(en,cap,cost,G[en].size()));
    G[en].push_back(edge(st,0,-cost,G[st].size()-1));
}
int min_cost_flow(int st,int en,int f){
    int ans=0;
    memset(h,0,sizeof(h));
    while(f>0){
        priority_queue<P,vector<P>,greater<P> >que;
        for(int i=0;i<maxn;i++) dis[i]=inf;
        dis[st]=0;que.push(P(0,st));
        while(!que.empty()){
            P p=que.top();que.pop();
            int v=p.second;
            if(dis[v]<p.first) continue;
            for(unsigned int i=0;i<G[v].size();i++){
                edge &e=G[v][i];
                if(e.cap>0&&dis[e.to]>dis[v]+e.cost+h[v]-h[e.to]){
                    dis[e.to]=dis[v]+e.cost+h[v]-h[e.to];
                    prevv[e.to]=v;
                    preve[e.to]=i;
                    que.push(P(dis[e.to],e.to));
                }
            }
        }
        if(dis[en]==inf) return -1;
        for(int i=0;i<maxn;i++) h[i]+=dis[i];
        int d=f;
        for(int i=en;i!=st;i=prevv[i]){
            d=min(d,G[prevv[i]][preve[i]].cap);
        }
        f-=d;
        ans+=d*h[en];
        for(int i=en;i!=st;i=prevv[i]){
            edge &e=G[prevv[i]][preve[i]];
            e.cap-=d;
            G[i][e.rev].cap+=d;
        }
    }
    return ans;
}
int main(){
    int T,n,m,u,v,cost;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        for(int i=0;i<maxn;i++) G[i].clear();
        for(int i=1;i<=n;i++) add_edge(0,i,1,0);
        for(int i=1+n;i<=2*n;i++) add_edge(i,2*n+1,1,0);
        for(int i=0;i<m;i++){
            scanf("%d%d%d",&u,&v,&cost);
            add_edge(u,v+n,1,cost);
        }
        int ans=min_cost_flow(0,2*n+1,n);
        printf("%d\n",ans);
    }
    return 0;
}