易错点：
- 如果是double类型的数据，判断是否为零需要if(fabs(gap)<1e-10).
- 既然是double就不能使用常见的delta了，因为能更新delta不代表是最终结果.
- 如果要开double类型的INF可以直接用1e12.
- 中间值转换成double类型可以使用value*1.0.
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
const int N=120;
const double INF=1e12;
double w[N][N],la[N],lb[N],delta;
int match[N],n;
bool va[N],vb[N];
bool dfs(int x){
va[x]=1;
for(int y=1;y<=n;y++){
if(vb[y])continue;
double gap=la[x]+lb[y]-w[x][y];
if(fabs(gap)<1e-10){
vb[y]=1;
if(!match[y]||dfs(match[y])){
match[y]=x;
return 1;
}
}
}
return 0;
}
int ans[N];
void KM(){
for(int x=1;x<=n;x++){
la[x]=-INF;
for(int y=1;y<=n;y++){
la[x]=max(la[x],w[x][y]);
}
}
for(int x=1;x<=n;x++){
while(1){
memset(va,0,sizeof(va));
memset(vb,0,sizeof(vb));
if(dfs(x))break;
delta=INF;
for(int i=1;i<=n;i++)
if(va[i])
for(int j=1;j<=n;j++){
if(!vb[j])
delta=min(delta,la[i]+lb[j]-w[i][j]);
}
for(int i=1;i<=n;i++){
if(va[i])la[i]-=delta;
if(vb[i])lb[i]+=delta;
}
}
}
for(int y=1;y<=n;y++){
ans[match[y]]=y;
}
return;
}
struct point{
int x,y;
}black[N],white[N];
double dis(point a,point b){
return sqrt((b.x-a.x)*(b.x-a.x)*1.0+(b.y-a.y)*(b.y-a.y)*1.0);
}
int main(){
scanf("%d",&n);
for(int i=1;i<=n;i++){
int x,y;
scanf("%d%d",&x,&y);
black[i].x=x,black[i].y=y;
}
for(int i=1;i<=n;i++){
int x,y;
scanf("%d%d",&x,&y);
white[i].x=x,white[i].y=y;
}
for(int x=1;x<=n;x++)
for(int y=1;y<=n;y++){
w[x][y]=-dis(black[x],white[y]);
}
KM();
for(int x=1;x<=n;x++)
printf("%d
",ans[x]);
return 0;
}