题意:在给出的有向图中,找到环,环的定义为头和尾出现两次,中间其他/它点出现一次,然后每个点属于一个环,问最小的费用
思路:因为肯定有一个方案满足条件,那么也就是每个点都会出发并且回到了自己,对于一个完备匹配的二分图来说,随意找一个起点,那么最后一定会回到这个起点对应的那个右部的点,所以我们跑最小费用流时流量设为点的个数就可以了
#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; }