UVaLive 4043 Ants（二分图最佳完美匹配、KM）

题意：

有n个蚂蚁n棵树，蚂蚁与树要配对，在配对成功的一对之间连一条线段，要求所有线段不能相交。按顺序输出蚂蚁所匹配的树。

分析：

二分图完美匹配问题。但是题目中有个线段不相交，怎么办？其实这个最佳完美匹配就是答案了。最佳完美匹配是权值和最大，那么我们就把两两点线段的权值搞成他们距离的负数即可。这样就不可能有相交的了。为什么？因为假设有相交，a1-b2，a2-b1，而dist(a1,b1)+dist(a2,b2) 肯定比前面交叉的小，可利用三角形三边关系证明，那么负数就是大了，也就是说交叉的在我们设计的负权那里是小的，所以就是最佳，也就是不可能有交叉的。保证最小的权值,所连的边一定是可以不相交的.
这样分析清楚了之后，就只要直接套用KM就OK了！

代码：

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

const int maxn = 100 + 10;

const double INF = 1e30;

int n;
double w[maxn][maxn];
double Lx[maxn], Ly[maxn];   // 顶标
int left[maxn];          // left[i]为右边第i个点的匹配点编号
bool S[maxn], T[maxn];   // S[i]和T[i]为左/右第i个点是否已标记
//KM模板 O(n^4)
bool match(int i) {
    S[i]=1;
    for(int j=1; j<=n; j++)
        if(abs(Lx[i]+Ly[j]-w[i][j])<1e-9&&!T[j]) {
            T[j]=1;
            if(!left[j]||match(left[j])) {
                left[j]=i;
                return true;
            }
        }
    return false;
}
void update() {
    double t=INF;
    for(int i=1; i<=n; i++)if(S[i])
            for(int j=1; j<=n; j++)if(!T[j])
                    t=min(t,Lx[i]+Ly[j]-w[i][j]);
    for(int i=1; i<=n; i++) {
        if(S[i])Lx[i]-=t;
        if(T[i])Ly[i]+=t;
    }
}

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++) {
        while(true) {
            for(int j=1; j<=n; j++)S[j]=T[j]=0;
            if(match(i))break;
            update();
        }
    }
}

int main() {
    int kase = 0;
    while(scanf("%d", &n) == 1) {
        if(++kase > 1) printf("\n");

        int x1[maxn], y1[maxn], x2[maxn], y2[maxn];
        for(int i = 1; i <= n; i++) scanf("%d%d", &x1[i], &y1[i]);
        for(int i = 1; i <= n; i++) scanf("%d%d", &x2[i], &y2[i]);
        for(int i = 1; i <= n; i++) // ant colony
            for(int j = 1; j <= n; j++) // apple tree
                w[j][i] = -sqrt((double)(x1[i]-x2[j])*(x1[i]-x2[j]) + (double)(y1[i]-y2[j])*(y1[i]-y2[j]));
        KM(); // 最大权匹配
        for(int i = 1; i <= n; i++) printf("%d\n", left[i]);
    }
    return 0;
}