HDU 2255 奔小康赚大钱(最佳二分图匹配)

来源:互联网 发布:软件服务商合作模式 编辑:程序博客网 时间:2024/04/28 10:56
#include<iostream>using namespace std;const int maxn=300+5;int w[maxn][maxn],n;int lx[maxn],ly[maxn];int Left[maxn];bool s[maxn],t[maxn];bool match(int i){    s[i]=true;    for(int j=1;j<=n;j++) if(lx[i]+ly[j]==w[i][j]&&!t[j]){        t[j]=true;        if(!Left[j]||match(Left[j])){            Left[j]=i;            return true;        }    }    return false;}void update(){    int a=(1<<30);    for(int i=1;i<=n;i++) if(s[i])        for(int j=1;j<=n;j++) if(!t[j])            a=min(a,lx[i]+ly[j]-w[i][j]);    for(int i=1;i<=n;i++){        if(s[i]) lx[i]-=a;        if(t[i]) ly[i]+=a;    }}void KM(){    for(int i=1;i<=n;i++){        Left[i]=lx[i]=ly[i]=0;        for(int j=1;j<=n;j++)            lx[i]=max(lx[i],w[i][j]);    }    for(int i=1;i<=n;i++){        for(;;)        {            for(int j=1;j<=n;j++) s[j]=t[j]=0;            if(match(i)) break;else update();        }    }}int main(){    while(~scanf("%d",&n))    {           for(int i=1;i<=n;i++)            for(int j=1;j<=n;j++) scanf("%d",&w[i][j]);        KM();        int ans=0;        for(int i=1;i<=n;i++) ans+=w[Left[i]][i];        printf("%d\n",ans);     }    return 0;}

0 0
原创粉丝点击