我们发现,对于一个有向图,他若是形成了一个换。那么我们把他拆点,成二分图。
会发现这个二分图一定会形成一个完备匹配。
那么这题就明了了。我们建立二分图,其中边的权值取负数,然后跑一遍KM就可以了。
我邻接表用多了,没怎么用邻接矩阵。这里刚开始忘给矩阵付初始的值-inf了。。。。。。
代码如下
#include<iostream> #include<algorithm> #include<queue> #include<cstring> using namespace std; const int max_n = 250; const int inf = 1e9; int nx,ny; int linker[max_n],lx[max_n],ly[max_n],slack[max_n]; int visx[max_n],visy[max_n],w[max_n][max_n]; int DFS(int x){ visx[x]=1; for(int y=1;y<=ny;y++){ if(visy[y]) continue; int tmp=lx[x]+ly[y]-w[x][y]; if(tmp==0){ visy[y]=1; if(linker[y]==-1 || DFS(linker[y])){ linker[y]=x; return 1; } }else if(slack[y]>tmp){ slack[y]=tmp; } } return 0; } int KM(){ int i,j; memset(linker,-1,sizeof(linker)); memset(ly,0,sizeof(ly)); for(i=1;i<=nx;i++) for(j=1,lx[i]=-inf;j<=ny;j++) if(w[i][j]>lx[i]) lx[i]=w[i][j]; for(int x=1;x<=nx;x++){ for(i=1;i<=ny;i++) slack[i]=inf; while(1){ memset(visx,0,sizeof(visx)); memset(visy,0,sizeof(visy)); if(DFS(x)) break; int d=inf; for(i=1;i<=ny;i++) if(!visy[i] && d>slack[i]) d=slack[i]; for(i=1;i<=nx;i++) if(visx[i]) lx[i]-=d; for(i=1;i<=ny;i++) if(visy[i]) ly[i]+=d; else slack[i]-=d; } } int res=0; for(i=1;i<=ny;i++) if(linker[i]!=-1) res+=w[linker[i]][i]; return res; } int n,m; int main(){ int T;scanf("%d",&T); while (T--){ scanf("%d %d",&n,&m); for (int i=1;i<=n;++i)for (int j=1;j<=n;++j)w[i][j]=-inf; for (int i=1;i<=m;++i){ int u,v,c; scanf("%d %d %d",&u,&v,&c); w[u][v]=max(-c,w[u][v]); }nx=ny=n; printf("%d\n",-KM()); } }