曼哈顿距离最小生成树

将n个点分成K部分，求最小的X，使得每部分中的任意点，都至少能找到一个相邻点 曼哈顿距离小于等于X

转换成求最小生成树上的第K大边。

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#define N 505050
using namespace std;
int n,K,tot=0,edgenum,fa[N],tree[N],treep[N];
struct point{int x,y,id,w;}p[N],a[N];
struct edge{int u,v,w;}b[N];
void add(int u,int v,int w){
    tot++;b[tot].u=u;b[tot].v=v;b[tot].w=w;
}
int dist(point a,point b){return abs(a.x-b.x)+abs(a.y-b.y);}
int find(int x){
    if(x==fa[x])return fa[x];
    fa[x]=find(fa[x]);
    return fa[x];
}
int cmp1(point a,point b){
    if(a.x==b.x)return a.y<b.y;
    return a.x<b.x;
}
int cmp2(point a,point b){
    return a.x<b.x;
}
int cmp3(edge a,edge b){return a.w<b.w;}
int query(int x){
    int p=-1,tmp=2000000000;
    for(int i=x;i<=n;i+=i&(-i)){
        if(tree[i]<tmp){
            tmp=tree[i];
            p=treep[i];
        }
    }
    return p;
}
void modify(int x,int val,int p){

    for(int i=x;i>0;i-=i&(-i)){
        if(val<tree[i])tree[i]=val,treep[i]=p;
    }
    return;
}
int main(){
    scanf("%d%d",&n,&K);
    for(int i=1;i<=n;i++){
        scanf("%d%d",&p[i].x,&p[i].y);
        p[i].id=i;
    }
    for(int dir=1;dir<=4;dir++){
        if(dir==2 || dir==4){
            for(int i=1;i<=n;i++){
                swap(p[i].x,p[i].y);
            }
        }
        if(dir==3)for(int i=1;i<=n;i++)p[i].x=-p[i].x;
        sort(p+1,p+n+1,cmp1);
        for(int i=1;i<=n;i++)a[i].x=p[i].y-p[i].x,a[i].id=i;
        sort(a+1,a+n+1,cmp2);
        int nowc=0;
        for(int i=1;i<=n;i++)tree[i]=2000000000,treep[i]=-1;
        for(int i=1;i<=n;i++){
            if(i==1 || a[i].x!=a[i-1].x)nowc++;
            p[a[i].id].w=nowc;
        }
        for(int i=n;i>=1;i--){
            int x=query(p[i].w);
            if(x>0)add(p[i].id,p[x].id,dist(p[i],p[x]));
            modify(p[i].w,p[i].x+p[i].y,i);
        }
    }
    int ans;
    sort(b+1,b+tot+1,cmp3);int now=0;
    for(int i=1;i<=n;i++)fa[i]=i;
    for(int i=1;i<=tot;i++){
        int u=b[i].u,v=b[i].v;
        int x=find(u),y=find(v);
        if(x!=y){
            fa[x]=y;
            now++;
            if(now==n-K){ans=b[i].w;break;}
        }
    }
    printf("%d\n",ans);
    return 0;
}

一件很奇怪的事情：

关于欧几里得距离最小生成树 网上到处是王栋老师WC论文中的Delaunay三角剖分法

而jiangshibiao大哥证明了60度这个结论，并且小于六十度这个分法也是正确的（据说）。

那其实完全可以用这个更好写的曼哈顿生成树写法去做欧几里得生成树。

想验证一下正确性，但是没有找到裸题。只有一个写不动的bzoj4219.。。。。。。。。 

曼哈顿这个模版应该错不了了，51nod和poj数据不水的话）

---

upd:

听了专题才知道。原来jsb这个四十五度做六十度的推广是错的。><