KM算法模板 hdu 2255

KM算法是在匹配是完备的情况下寻找最优匹配。

首先，先将范围定为最大的情况，如果最大的情况无法满足，就下降一个维度继续匹配。

直到匹配成功。

 #include<cstdio>
 #include<string.h>
 #include<algorithm>
 using namespace std;
 const int maxn=3e2+;
 const int inf=0x3f3f3f3f;
 int G[maxn][maxn],n;
 int lx[maxn],ly[maxn];
 int match[maxn];
 int visx[maxn],visy[maxn];
 void init()
 {
     memset(match,,sizeof(match));
 }
 int dfs(int k)
 {
     visx[k]=;
     for(int i=;i<=n;i++){
         if(!visy[i]&&G[k][i]==lx[k]+ly[i]){
             visy[i]=;
             if(!match[i]||dfs(match[i])){
                 match[i]=k;
                 return ;
             }
         }
     }
     return ;
 }
 int KM()
 {
     for(int i=;i<=n;i++){
         lx[i]=-inf,ly[i]=;
         for(int j=;j<=n;j++)
             lx[i]=max(lx[i],G[i][j]);
     }
     for(int k=;k<=n;k++){
         while(){
             memset(visx,,sizeof(visx));
             memset(visy,,sizeof(visy));
             if(dfs(k)) break;
             int mn=inf;

             for(int i=;i<=n;i++) if(visx[i])
             for(int j=;j<=n;j++) if(!visy[j])
             mn=min(mn,lx[i]+ly[j]-G[i][j]);
             if(mn==inf) return -;
             for(int i=;i<=n;i++) if(visx[i]) lx[i]-=mn;
             for(int i=;i<=n;i++) if(visy[i]) ly[i]+=mn;
         }
     }
     int ans=;
     for(int i=;i<=n;i++) if(match[i])
     ans+=G[match[i]][i];
     return ans;
 }
 int main()
 {
     while(scanf("%d",&n)!=EOF){
         init();
         for(int i=;i<=n;i++)
         for(int j=;j<=n;j++)
         scanf("%d",&G[i][j]);
         printf("%d\n",KM());
     }
     return ;
 }