题目链接:https://vjudge.net/problem/HDU-2255
题意:裸的最大权匹配。
解法:模板,KM算法,复杂度O(nx*nx*ny)
#include <bits/stdc++.h> using namespace std; const int maxn = 310; const int inf = 0x3f3f3f3f; int nx,ny; int g[maxn][maxn]; int linker[maxn], lx[maxn], ly[maxn]; int slack[maxn]; bool visx[maxn],visy[maxn]; bool DFS(int x){ visx[x]=true; for(int y=0; y<ny; y++){ if(visy[y]) continue; int tmp = lx[x]+ly[y]-g[x][y]; if(tmp == 0){ visy[y]=true; if(linker[y]==-1||DFS(linker[y])){ linker[y]=x; return true; } } else if(slack[y]>tmp) slack[y]=tmp; } return false; } int KM(){ memset(linker,-1,sizeof(linker)); memset(ly,0,sizeof(ly)); for(int i=0; i<nx; i++){ lx[i]=-inf; for(int j=0; j<ny; j++) if(g[i][j]>lx[i]) lx[i] = g[i][j]; } for(int x=0; x<nx; x++){ for(int i=0; i<ny; i++) slack[i]=inf; while(1){ memset(visx,false,sizeof(visx)); memset(visy,false,sizeof(visy)); if(DFS(x)) break; int d = inf; for(int i=0; i<ny; i++) if(!visy[i]&&d>slack[i]) d = slack[i]; for(int i=0; i<nx; i++) if(visx[i]) lx[i] -= d; for(int i=0; i<ny; i++){ if(visy[i]) ly[i]+=d; else slack[i]-=d; } } } int res = 0; for(int i=0; i<ny; i++) if(linker[i]!=-1) res += g[linker[i]][i]; return res; } int main() { int n; while(scanf("%d",&n)!=EOF) { for(int i=0; i<n; i++) for(int j=0; j<n; j++) scanf("%d", &g[i][j]); nx = ny = n; printf("%d ", KM()); } return 0; }