半径R覆盖最多点

struct point {
    double x, y;
};
point p[N];
struct alpha {
    double v;
    bool flag;
    bool friend operator <(const alpha &a, const alpha &b) {
        return a.v < b.v;
    }
} alp[N * ];
double dis(point a, point b) {
    return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));
}
int n,s;
bool judge(double R) {
    int ans = ;
    double theta, phi, dist;
    for(int i = ; i < n; i++) {
        int k = ;
        for(int j = ; j < n; j++) {
            if(i == j) continue;
            dist = dis(p[i], p[j]);
            if(dist > 2.0 * R)
                continue;
            theta = atan2(p[j].y - p[i].y, p[j].x - p[i].x);
            if(theta < )
                theta +=  * pi;
            phi = acos(dist / (2.0 * R));
            alp[k].v = theta - phi +  * pi;
            alp[k].flag = true;
            alp[k + ].v = theta + phi +  * pi;
            alp[k + ].flag = false;
            k += ;
        }
        sort(alp, alp + k);
        int sum = ;
        for(int j = ; j < k; j++) {
            if(alp[j].flag)
                sum ++;
            else
                sum --;
            if(sum > ans)
                ans = sum;
        }
    }
    return ans +  >= s;//ans+1 为最多点
}