Tyvj1461

题目链接

分析：
KDtree最远点对

tip

新学的算法就要多加练习
注意：kdtree中的左右儿子不是*2，*2+1这么简单了
一定要调用t[bh].l,t[bh].r
最远距离的计算还是要注意：

这里写代码片
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long

using namespace std;

const int N=100010;
struct node{
    ll mx[2],mn[2],d[2],l,r;
};
node t[N];
int n,root,cmpd;
ll x,y;
ll an;

int cmp(const node &a,const node &b)
{
    return (a.d[cmpd]<b.d[cmpd])||((a.d[cmpd]==b.d[cmpd])&&(a.d[!cmpd]<b.d[!cmpd]));
}

ll max(ll a,ll b){if (a>b) return a;else return b;}

void update(int bh)
{
    int lc=t[bh].l;
    int rc=t[bh].r;
    if (lc)
    {
        t[bh].mn[0]=min(t[bh].mn[0],t[lc].mn[0]);
        t[bh].mn[1]=min(t[bh].mn[1],t[lc].mn[1]);
        t[bh].mx[0]=max(t[bh].mx[0],t[lc].mx[0]);
        t[bh].mx[1]=max(t[bh].mx[1],t[lc].mx[1]);
    }
    if (rc)
    {
        t[bh].mn[0]=min(t[bh].mn[0],t[rc].mn[0]);
        t[bh].mn[1]=min(t[bh].mn[1],t[rc].mn[1]);
        t[bh].mx[0]=max(t[bh].mx[0],t[rc].mx[0]);
        t[bh].mx[1]=max(t[bh].mx[1],t[rc].mx[1]);
    }
}

int build(int l,int r,int D)
{
    cmpd=D;
    int mid=(l+r)>>1;
    nth_element(t+1+l,t+1+mid,t+1+r,cmp);
    t[mid].mn[0]=t[mid].mx[0]=t[mid].d[0];
    t[mid].mn[1]=t[mid].mx[1]=t[mid].d[1];
    if (l!=mid) t[mid].l=build(l,mid-1,!D);
    if (r!=mid) t[mid].r=build(mid+1,r,!D);
    update(mid);
    return mid; 
}

ll sqr(ll x){return x*x;}
ll dis(int p,ll x,ll y)
{
    ll d=0;
    d=sqr(max(abs(t[p].mn[0]-x),abs(t[p].mx[0]-x)));
    d+=sqr(max(abs(t[p].mn[1]-y),abs(t[p].mx[1]-y)));
    return d;
}

void ask(int now)
{
    ll d0,dl,dr;
    d0=(x-t[now].d[0])*(x-t[now].d[0])+(y-t[now].d[1])*(y-t[now].d[1]);
    if (t[now].l) dl=dis(t[now].l,x,y);
    else dl=0;
    if (t[now].r) dr=dis(t[now].r,x,y);
    else dr=0;
    if (d0>an) an=d0;
    if (dl>dr)
    {
        if (dl>an) ask(t[now].l);
        if (dr>an) ask(t[now].r);
    }
    else
    {
        if (dr>an) ask(t[now].r);
        if (dl>an) ask(t[now].l); 
    }
}

int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;i++)
        scanf("%lld%lld",&t[i].d[0],&t[i].d[1]);
    root=build(1,n,0);
    ll ans=0;
    for (int i=1;i<=n;i++)
    {
        x=t[i].d[0]; y=t[i].d[1];
        an=0;
        ask(root);
        if (an>ans) ans=an;
    }
    printf("%lld",ans);
    return 0;
}