CF 1466F-Line Distance 几何关系转化+二分+树状数组

注释有详细的推导过程，一些结论可以多画一画

	
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;

/*
    CF1446F
        题意：给定一些点，计算这些点两两之间连成的直线与原点的距离第K大的是多少
        思路：直接暴力必炸->不去枚举点对，而去二分距离mid，若能在线性复杂度以下得到有多少点对构成的直线与原点的距离<mid，就可以以此调整二分区间
                    ->若以原点为圆心，mid为半径，做一个圆，则所有与圆相交的直线到原点的距离都小于mid
                    ->若一条直线与圆相交，则直线上两个在园外的点向圆做两条切线，若记这两点为A, B，切于A1, A2, B1, B2，则A1A2与B1B2两条弦不相交
                        反之，若直线与圆相离，则A1A2与B1B2两条弦相交（！！重要结论，可以把直线与圆的关系问题转化为弦相交问题）
                    ->若两条弦相交，记A1, A2, B1, B2的极角为a1, a2, b1, b2，一定有a1<b1<a2<b2 (即由极角构成的区间一定有重叠部分)
                        ！！重要结论，将弦相交转化为区间交
                    ->现在问题变成了给定一些区间，求有多少对区间有重叠部分，但又不是包含
                    ->将区间按的左端点从小到大排序，对于第i个区间，前1~i-1个区间的左端点都小于当前的左端点，那么只要统计有多少个右端点落在L[i]~R[i]之间即可
                    ->用树状数组可以解决这个问题
*/

const int maxn = 1e5 + 5;
const double eps=1e-8;
int N; LL K; double pi;
int X[maxn], Y[maxn];
struct Seg {
    double radl, radr;
    void correct() {
        if(radl>2*pi) radl-=2*pi;
        if(radr>2*pi) radr-=2*pi;
        if(radl<0) radl+=2*pi;
        if(radr<0) radr+=2*pi;
        if(radl>radr) swap(radl, radr);
    }
    bool operator<(const Seg& rhs) const {
        return radl<rhs.radl;
    }
} ses[maxn]; double rads[maxn*2];
unordered_map<double, int> mp;

int c[maxn*2];
void add(int pos, int val) {
    while(pos<=N*2) {
        c[pos]+=val;
        pos+=(pos&(-pos));
    }
}

int sum(int pos) {
    int res = 0;
    while(pos) {
        res+=c[pos];
        pos-=(pos&(-pos));
    }
    return res;
}

int main() {
    scanf("%d%lld", &N, &K);
    LL ALL=1ll*N*(N-1)/2; pi=acos(-1);
    for(int i=1; i<=N; i++) scanf("%d%d", &X[i], &Y[i]);
    double mxr=1e5, mir=0;//一定注意至少最大半径至少是1e4*sqrt(2)以上
    while(mxr-mir>eps) {
        double mid=(mxr+mir)/2.0;
        int rcnt = 0, rscnt = 0; mp.clear(); memset(c, 0, sizeof(c));
        for(int i=1; i<=N; i++) {
            double dis=sqrt(X[i]*X[i]+Y[i]*Y[i]);
            if(dis<=mid) continue;
            double ang=atan2(Y[i], X[i])+pi;//atan2的极角范围为[-pi, pi]，加上pi后就可以转变到[0, 2pi]
            double dan=acos(mid/dis);
            ses[++rcnt]=Seg{ang+dan, ang-dan};//计算两个切点的极角
            ses[rcnt].correct();//将极角限制在[0, 2pi]以内
            rads[++rscnt]=ses[rcnt].radl;
            rads[++rscnt]=ses[rcnt].radr;
        }
        sort(ses+1, ses+rcnt+1);
        sort(rads+1, rads+rscnt+1);
        for(int i=1; i<=rscnt; i++) mp[rads[i]]=i;//利用map对浮点数的区间离散化
        LL cnt = 0;
        for(int i=1; i<=rcnt; i++) {
            cnt+=sum(mp[ses[i].radr]-1)-sum(mp[ses[i].radl]);//计算与当前区间相交区间个数
            add(mp[ses[i].radr], 1);
        }
        if(ALL-cnt>=K) mxr=mid;
        else mir=mid;
    }
    printf("%.8lf\n", mir);
    return 0;
}