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

一開始理所當然的認爲橫排和縱排應該是橫座標和縱座標和的平均值,但是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;
}


發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章