HDU 3395 Special Fish(KM算法)

来源:互联网 发布:毛利寿三郎 越知月光 编辑:程序博客网 时间:2024/05/21 21:47
#include<iostream>#include<cstring>using namespace std;const int maxn=100+5;int w[maxn][maxn],n;int lx[maxn],ly[maxn];int Left[maxn];bool s[maxn],t[maxn];int val[maxn];char map[maxn][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)&&n)    {        memset(w,0,sizeof(w));        for(int i=1;i<=n;i++) scanf("%d",&val[i]);        for(int i=1;i<=n;i++) scanf("%s",map[i]+1);        for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) if(map[i][j]=='1') w[i][j]=(val[i]^val[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