HDU 2255 奔小康赚大钱 最大权匹配

题目链接: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\n", KM());
    }
    return 0;
}

 

posted @ 2017-09-05 19:32  zxycoder  阅读(128)  评论(0编辑  收藏  举报