Uva live 4043 Ants( KM+简单的计算几何)

一定要清楚的一个问题,就是:(如图)两条黑色的线的和一定小于蓝色的和,但是平方的和之间就不一定是这个关系了!


代码:

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <cmath>
using namespace std;

const double INF = 0x7fffffff;
const double eps = 1e-9;
const int N = 110;
int n;
double w[N][N], lx[N], ly[N];
int l[N], r[N];
int wx[N], wy[N], bx[N], by[N];
bool S[N], T[N];

bool eq( double x, double y ) {
    return fabs ( x - y ) < eps;
}
bool match( int i )
{
    S[i] = true;
    for ( int j = 1; j <= n; ++j ) if ( eq( lx[i] + ly[j], w[i][j] ) && !T[j] ) {
        T[j] = true;
        if ( !l[j] || match(l[j]) ) {
            l[j] = i;
            //r[i] = j;
            return true;
        }
    }
    return false;
}
void update()
{
    double a = INF;
    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()
{
    //memset(l, 0, sizeof(l));
    //memset(r, 0, sizeof(r));
    for ( int i = 1; i <= n; ++i ) {
        l[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()
{
    bool first = false;
    while ( scanf("%d", &n) != EOF ) {
        if ( first ) printf("\n");
        first = true;
        for ( int i = 1; i <= n; ++i ) scanf("%d%d", &wx[i], &wy[i]);
        for ( int i = 1; i <= n; ++i ) scanf("%d%d", &bx[i], &by[i]);
        for ( int i = 1; i <= n; ++i ) 
            for ( int j = 1; j <= n; ++j ) 
                w[j][i] = -sqrt((double)(wx[i]-bx[j])*(wx[i]-bx[j])+(double)(wy[i]-by[j])*(wy[i]-by[j]));
        KM();
        for ( int i = 1; i <= n; ++i ) printf("%d\n", l[i]);
    }
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值