UVA 1405 The Great Wall Game(最佳二分图匹配)

来源:互联网 发布:中智软件开发有限公司 编辑:程序博客网 时间:2024/06/06 02:40

一开始理所当然的认为横排和纵排应该是横坐标和纵坐标和的平均值,但是WA了好几次。检查了一下,好像也就这么一个可疑点,算算时间复杂度也不会T,就直接暴力了,然后就过了= =。

#include<iostream>using namespace std;const int maxn=300+5;const int inf=(1<<30);int w[maxn][maxn],n;int lx[maxn],ly[maxn];int Left[maxn];bool s[maxn],t[maxn];int id[maxn][maxn];struct point{int x,y;}p[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();        }    }}void clear() {for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) w[i][j]=-inf;}int main(){    int kase=0;    while(scanf("%d",&n)&&n)    {        for(int i=1;i<=n;i++){            scanf("%d%d",&p[i].x,&p[i].y);            // cout<<p[i].x<<" "<<p[i].y<<endl;        }        int ans=inf,tmp=0;        //横排和纵排        for(int i=1;i<=n;i++){            clear();            for(int k=1;k<=n;k++){                int x=p[k].x,y=p[k].y;                for(int j=1;j<=n;j++){                    w[k][j]=-(abs(x-i)+abs(y-j));                }            }            KM();            int tmp=0;            for(int k=1;k<=n;k++) tmp-=w[Left[k]][k];            ans=min(ans,tmp);            clear();            for(int k=1;k<=n;k++){                int x=p[k].x,y=p[k].y;                // cout<<x<<" "<<y<<endl;                for(int j=1;j<=n;j++){                    w[k][j]=-(abs(x-j)+abs(y-i));                    // cout<<k<<" "<<x<<" "<<y<<" "<<i<<" "<<j<<" "<<w[k][j]<<endl;                }            }            KM();            tmp=0;            for(int k=1;k<=n;k++) tmp-=w[Left[k]][k];            ans=min(ans,tmp);        }        //        for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) w[i][j]=-inf;        tmp=0;        for(int i=1;i<=n;i++){            int xx=p[i].x,yy=p[i].y;            for(int j=1;j<=n;j++){                w[i][j]=-(abs(xx-j)+abs(yy-j));            }        }        KM();        for(int i=1;i<=n;i++) tmp-=w[Left[i]][i];        ans=min(ans,tmp);                for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) w[i][j]=-inf;        tmp=0;        for(int i=1;i<=n;i++){            int xx=p[i].x,yy=p[i].y;            for(int j=1;j<=n;j++){                w[i][j]=-(abs(xx-j)+abs(yy-(n-j+1)));            }        }        KM();        for(int i=1;i<=n;i++) tmp-=w[Left[i]][i];        ans=min(ans,tmp);        // if(kase++) printf("\n");        printf("Board %d: %d moves required.\n",++kase,ans);        printf("\n");    }    return 0;}


0 0
原创粉丝点击