UVALIVE 4970 最小权匹配

首先贴一下这道题的BNU地址,UVA地址自己找吧。

http://acm.bnu.edu.cn/bnuoj/problem_show.php?pid=11852

题意:这道题的意思就是,给你N个棋子的坐标,这些棋子的走法是象棋中的马的走法。然后再给你N个坐标终点。

问所有的棋子走到其中一个坐标上,不能有重复,最少的步数是多少。

思路:这道题直接搜显然爆,因为他的坐标范围都是int ,所以我们就要考虑怎么处理出两个点之间的最短步数。

我直接把问题扔给队友了,然后他推出了一条公式,太神不能多说。

问题简化成给你两个点(x , y ) ,(x1 , y1),求这两个点之间最少需要多少步。

然后他经过一系列的数学推导,就搞定了。

过程我也不会,等下看代码。

然后接下来每个棋子到每个终点坐标的距离都处理出来了。

那么根据N = 15 。那么很显然状压DP就可以搞定,求出最小值。

但是考虑到每个棋子只有一个位置,那么很显然棋子和终点位置是二分图,那么很显然我们可以用最小权匹配来搞这个问题。

KM的姿势比状压优美许多。

 

#include <set>

#include <map>

#include <stack>

#include <cmath>

#include <queue>

#include <cstdio>

#include <string>

#include <vector>

#include <iomanip>

#include <cstring>

#include <iostream>

#include <algorithm>

#define Max 2505

#define FI first

#define SE second

#define ll long long

#define PI acos(-1.0)

#define inf 1ll << 60ll

#define LL(x) ( x << 1 )

#define bug puts("here")

#define PII pair<int,int>

#define RR(x) ( x << 1 | 1 )

#define mp(a,b) make_pair(a,b)

#define mem(a,b) memset(a,b,sizeof(a))

#define REP(i,s,t) for( int i = ( s ) ; i <= ( t ) ; ++ i )



using namespace std;



ll f(ll x,ll y) {//推出来的,太神不多说

    x=abs(x);

    y=abs(y);

    if(x>y) {

        swap(x,y);

    }

    if(x==0&&y==1) {

        return 3;

    }

    if(x==2&&y==2) {

        return 4;

    }

    ll z=max((y+1)/2ll,(x+y+2)/3ll);

    if((z-x-y)&1) {

        z++;

    }

    return z;

}

int n ;

ll x[22] , y[22] ;

ll xx[22] , yy[22] ;

ll Map[22][22] ;

ll lx[22] ,ly[22] ,linkx[22] , linky[22] ;

bool visx[22] , visy[22] ;

int find(int now){

    visx[now] = 1 ;

    for (int i = 0 ; i < n ; i ++ ){

        if(!visy[i] && Map[now][i] - lx[now] - ly[i] == 0){

            visy[i] = 1 ;

            if(linky[i] == -1 || find(linky[i])){

                linkx[now] = i ;

                linky[i] = now ;

                return 1 ;

            }

        }

    }

    return 0 ;

}

ll KM(){

    mem(linkx , -1) ;

    mem(linky , -1) ;

    mem(ly ,0) ;

    for (int i = 0 ; i < n ; i ++ ){

        while(1){

            mem(visx ,0) ;mem(visy , 0) ;

            if(find(i))break ;

            ll d = inf ;

            for (int j = 0 ; j < n ; j ++ ){

                if(visx[j])

                    for (int k = 0 ; k < n ; k ++ ){

                        if(!visy[k])

                            d = min(d , Map[j][k] - lx[j] - ly[k]) ;

                    }

            }

            for (int j = 0 ; j < n ; j ++ ){

                if(visx[j])lx[j] += d ;

                if(visy[j])ly[j] -= d ;

            }

        }

    }

    ll ans = 0 ;

    for (int i = 0 ; i < n ; i ++ ){

        if(linky[i] != -1 )ans += Map[linky[i]][i] ;

    }

    return ans ;

}

int main() {

    int ca = 0 ;

    while(cin >> n , n ){

        for (int i = 0 ; i < n ; i ++ )cin >> x[i] >> y[i] ;

        for (int i = 0 ; i < n ; i ++ )cin >> xx[i] >> yy[i] ;

        for (int i = 0 ; i < n ; i ++ )lx[i] = inf ;

        for (int i = 0 ; i < n ; i ++ ){

            for (int j = 0 ; j < n ; j ++ ){

                Map[i][j] = f(x[i] - xx[j] , y[i] - yy[j]) ;

                lx[i] = min(lx[i] , Map[i][j]) ;

            }

        }

        cout << ++ ca << ". " ;

        cout << KM() << endl;

    }

    return 0 ;

}


 

 

你可能感兴趣的:(live)