「CQOI2016」K 远点对

/*
考虑暴力  可以n ^ 2枚举点对 然后用一个容量为2k的小根堆来维护第k大
kd-tree呢就是来将这个暴力优化， 每次先找远的并且最远距离不如堆顶的话就不继续找下去
貌似挺难构造数据卡的
*/

#include<iostream>
#include<queue>
#include<cmath>
#include<cstdio>
#include<algorithm>
#include<cstring>
#define M 100010
#define ll long long
#define inf 1000000000000000ll
#define sqr(x) 1ll * (x)* (x)
using namespace  std;
int read() {
    int nm = , f = ;
    char c = getchar();
    for(; !isdigit(c); c = getchar()) if(c == '-') f = -;
    for(; isdigit(c); c = getchar()) nm = nm *  + c - '';
    return nm * f;
}
int n, k, q, root, lc[M], rc[M];
ll minn[M][], maxx[M][];
struct Note {
    ll d[];
    bool operator < (const Note &b) const {
        return  d[q] != b.d[q] ? d[q] < b.d[q] : d[q^] < b.d[q ^ ];
    }
} tmp, ver[M];

priority_queue<ll, vector<ll>, greater<ll> >que;

void pushup(int x) {
    for(int i = ; i <= ; i++) {
        minn[x][i] = maxx[x][i] = ver[x].d[i];
        minn[x][i] = min(minn[x][i], min(minn[lc[x]][i], minn[rc[x]][i]));
        maxx[x][i] = max(maxx[x][i], max(maxx[lc[x]][i], maxx[rc[x]][i]));
    }
}

int build(int l, int r, int kx) {
    if(l > r) return ;
    int mid = (l + r) >> ;
    q = kx;
    nth_element(ver + l, ver + mid, ver + r + );
    lc[mid] = build(l, mid - , kx ^ );
    rc[mid] = build(mid + , r, kx ^ );
    pushup(mid);
    return mid;
}

ll get(int x) {
    if(!x) return ;
    ll rt = ;
    for(int i = ; i <= ; i++) rt += max(sqr(maxx[x][i] - tmp.d[i]), sqr(minn[x][i] - tmp.d[i]));
    return rt;
}

ll calc(Note a, Note b) {
    return sqr(a.d[] - b.d[]) + sqr(a.d[]-b.d[]);
}

void query(int x) {
    if(!x) return;
    ll dl = get(lc[x]), dr = get(rc[x]), d = calc(tmp, ver[x]);
    if(d > que.top()) que.pop(), que.push(d);
    if(dl > dr) {
        if(dl > que.top()) query(lc[x]);
        if(dr > que.top()) query(rc[x]);
    } else {
        if(dr > que.top()) query(rc[x]);
        if(dl > que.top()) query(lc[x]);
    }
}
int main() {
    n = read(), k = read();
    for(int i = ; i <= n; i++) ver[i].d[] = read(), ver[i].d[] = read();
    for(int i = ; i <= ; i++) minn[][i] = inf, maxx[][i] = -inf;
    root = build(, n, );
    for(int i = ; i <=  * k; i++) que.push();
    for(int i = ; i <= n; i++) tmp = ver[i],query(root);
    cout << que.top();
    return ;
}