我不会啊
老老实实去学 曼哈顿距离最小生成树
以一个点为原点建立直角坐标系,在每45度内只会向距离该点最近的一个点连边
这样只会有
O(n)
条边
怎么求最近点 是个二维偏序 排完序后 树状数组维护即可
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
typedef pair<int,int> abcd;
inline char nc(){
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline void read(int &x){
char c=nc(),b=1;
for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}
const int N=10005;
struct edge{
int u,v,w;
edge(int u=0,int v=0,int w=0):u(u),v(v),w(w) { }
bool operator < (const edge &B) const{ return w<B.w; }
}ed[N<<2]; int tot;
int fat[N];
inline void Init(int n){
for (int i=1;i<=n;i++) fat[i]=i;
}
inline int Fat(int u){
return u==fat[u]?u:fat[u]=Fat(fat[u]);
}
inline int Merge(int x,int y){
x=Fat(x); y=Fat(y); if (x==y) return 0;
fat[x]=y; return 1;
}
inline int Krus(int n,int K){
Init(n); sort(ed+1,ed+tot+1);
for (int i=1;i<=tot;i++)
if (Merge(ed[i].u,ed[i].v))
if (!(--K))
return ed[i].w;
}
namespace BIT{
abcd c[N]; int maxn;
inline void init(int n){
maxn=n; for (int i=1;i<=n;i++) c[i]=abcd(1<<30,0);
}
inline void modify(int x,abcd r){
for (int i=x;i;i-=i&-i)
c[i]=min(c[i],r);
}
inline int query(int x){
abcd ret(1<<30,0);
for (int i=x;i<=maxn;i+=i&-i)
ret=min(ret,c[i]);
return ret.second;
}
}
struct P{
int x,y,idx;
bool operator < (const P &B)const { return x==B.x?y<B.y:x<B.x; }
friend int dist(P a,P b){ return abs(a.x-b.x)+abs(a.y-b.y); }
}p[N];
int n,K;
int sx[N],icnt;
inline int Bin(int x){
return lower_bound(sx+1,sx+icnt+1,x)-sx;
}
inline void Build(){
using namespace BIT;
for (int d=0;d<4;d++){
if (d&1) for (int i=1;i<=n;i++) swap(p[i].x,p[i].y);
if (d==2) for (int i=1;i<=n;i++) p[i].x=-p[i].x;
sort(p+1,p+n+1);
icnt=0; for (int i=1;i<=n;i++) sx[++icnt]=p[i].y-p[i].x;
sort(sx+1,sx+icnt+1); icnt=unique(sx+1,sx+icnt+1)-sx-1;
init(icnt);
for (int i=n;i;i--){
int pos=Bin(p[i].y-p[i].x),v=query(pos);
if (v)
ed[++tot]=edge(p[i].idx,p[v].idx,dist(p[i],p[v]));
modify(pos,abcd(p[i].x+p[i].y,i));
}
}
}
int main(){
freopen("t.in","r",stdin);
freopen("t.out","w",stdout);
read(n); read(K);
for (int i=1;i<=n;i++) read(p[i].x),read(p[i].y),p[i].idx=i;
Build();
printf("%d\n",Krus(n,n-K));
return 0;
}