[校内模拟][bitset][计算几何]ball

传送门

我们可以发现最优解一定在一个三角形上，且这个三角形最短边就是答案
那么我们按边的长度排序，用bitset查询什么时候会形成一个三角形
然后就做完了

Code：

#include<bits/stdc++.h>
#define pb push_back
using namespace std;
inline int read(){
    int res=0,f=1;char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') f=-f;ch=getchar();}
    while(isdigit(ch)) {res=(res<<1)+(res<<3)+(ch^48);ch=getchar();}
    return res*f;
}
const int N=4005;
struct point{
    int x,y;
    point(int _x=0,int _y=0):x(_x),y(_y){}
    friend inline point operator +(const point &a,const point &b){return point(a.x+b.x,a.y+b.y);}
    friend inline point operator -(const point &a,const point &b){return point(a.x-b.x,a.y-b.y);}
    friend inline int operator *(const point &a,const point &b){return a.x*b.y-a.y*b.x;}
    inline int dis(){return x*x+y*y;}
}p[N];
bitset<N>v[N];
struct node{
    short fi,se;int d;
    friend inline bool operator <(const node &a,const node &b){return a.d>b.d;}
}vec[N*N/2];
int n,tot;
int main(){
    n=read();
    for(int i=1;i<=n;i++) p[i].x=read(),p[i].y=read();
    for(int i=1;i<=n;i++)
    for(int j=1;j<i;j++) vec[++tot]=node{i,j,(p[i]-p[j]).dis()};
    sort(vec+1,vec+tot+1);
    for(int i=1;i<=tot;i++){
        node x=vec[i];
        if((v[x.fi]&v[x.se]).count()){printf("%.10lf",sqrt((p[x.fi]-p[x.se]).dis())/2);return 0;}
        v[x.fi].set(x.se),v[x.se].set(x.fi);
    }
}