HDU 2255 （KM算法）

最新推荐文章于 2021-09-28 11:10:55 发布

凉秋girl

最新推荐文章于 2021-09-28 11:10:55 发布

阅读量262

点赞数

分类专栏：算法

本文链接：https://blog.csdn.net/qq_36124802/article/details/55669505

版权

算法专栏收录该内容

35 篇文章 0 订阅

订阅专栏

the reason of failure：1、超时，while(scanf(%d,&n))是错误的写法，而里面改为cin则为正确的了，但是while(scanf("%d,&n)==1)是正确的。

thinking：

输入例子:

4
21 8 9 13
17 15 34 19
16 22 51 43
32 17 18 22

就是先把lx初始化为每行的最大值，lx[1]=21,lx[2]=34,lx[3]=51,lx[4]=32

然后linker[1]=1,linker[3]=2,但是在第3行找到最大的时候，linker[3]在之前已经存在了，那么就dfs深度搜索，本来slack[1]=51-16,slack[2]=51-22.slack[4]=51-43,

然后在第2行与第3行之间找到最小的slack[],然后linker，这样是求二分图最大权值的方法

输出为117

代码:

#include <iostream>
#include <string.h>
#include <math.h>
using namespace std;
int n;
int cost[500][500];
int lx[500],ly[500];
int visx[500],visy[500];
int linker[500];
int slack[500];
int bfs(int x){
    visx[x]=1;
    int tmp,i;
    for(i=1;i<=n;i++){
        if(visy[i])continue;
        tmp=lx[x]+ly[i]-cost[x][i];
        if(tmp==0){
            visy[i]=1;
            if(linker[i]==0||bfs(linker[i])){
                linker[i]=x;
                return 1;
            }
        }else{
            if(tmp<slack[i])slack[i]=tmp;
        }
    }
    return 0;
}
long long maxcost(){
    memset(lx,0,sizeof(lx));
    memset(ly,0,sizeof(ly));
    int i,j;
    for(i=1;i<=n;i++)
        for(j=1;j<=n;j++)
        if(cost[i][j]>lx[i])lx[i]=cost[i][j];
    for(int k=1;k<=n;k++){
    while(1){
            memset(visx,0,sizeof(visx));
            memset(visy,0,sizeof(visy));
            for(i=1;i<=n;i++)
                slack[i]=99999;
        if(bfs(k))break;
        else{
            int d=99999;
            for(i=1;i<=n;i++)
                if(d>slack[i])d=slack[i];
            for(i=1;i<=n;i++){
                if(visx[i])lx[i]-=d;
                if(visy[i])ly[i]+=d;}
        }
    }
    }
    long long sum=0;
    for(i=1;i<=n;i++)
        sum+=cost[linker[i]][i];
        return sum;
}
int main(){
    freopen("in.txt","r",stdin);
    int i,j,k,l,f1,f2,f3,t1,t2,t3;
    while(scanf("%d",&n)==1){
            memset(cost,0,sizeof(cost));
            memset(linker,0,sizeof(linker));
        for(i=1;i<=n;i++)
        for(j=1;j<=n;j++){
            scanf("%d",&t1);
            cost[i][j]=t1;
    }
    printf("%lld\n",maxcost());
}
    return 0;
}