【bzoj 1913】signaling 信号覆盖（计算几何）

传送门biu~
考虑贡献。我们可以发现，对于能组成凸四边形的四个点，它们对答案的贡献是2；对于能组成凹四边形的四个点，它们对答案的贡献是1。问题就转化成了计算凹/凸四边形的数量。
因为凹四边形和凸四边形的数量和是确定的$C_n^4$，而凹四边形还有着“有中心点”这个性质，所以我们选择计算凹四边形。
枚举凹四边形的中心点，将其它点按辐角排序。如果其余三个点到中心点的辐角差不超过$\pi$，那么就说明这不是一个凹四边形。将这部分减去就能求出凹四边形数目了。
代码：

#include<bits/stdc++.h>
#define N 1505
using namespace std;
const double pi=acos(-1);
int n;long long t,C[N][5]={1};
struct point{int x,y;}p[N];
double a[N<<1],ans;
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;++i){
        C[i][0]=1;
        for(int j=1;j<=4;++j){
            C[i][j]=C[i-1][j]+C[i-1][j-1];
        }
    }
    for(int i=1;i<=n;++i) scanf("%d%d",&p[i].x,&p[i].y);
    for(int cnt=0,i=1;i<=n;++i,cnt=0){
        for(int j=1;j<=n;++j){
            if(i==j) continue;
            a[++cnt]=atan2(p[j].y-p[i].y,p[j].x-p[i].x);
            if(a[cnt]<0) a[cnt]+=2*pi;
        }
        sort(a+1,a+cnt+1);
        for(int j=1;j<=cnt;++j) a[cnt+j]=a[j]+2*pi;
        long long x(0);int now(1);
        for(int j=1;j<=cnt;++j){
            while(now<2*cnt && a[now+1]-a[j]<pi) ++now;
            if(now-j>1) x+=C[now-j][2];
        }
        t+=C[cnt][3]-x;
    }
    ans=(double)(t+(C[n][4]-t)*2)/C[n][3]+3;
    printf("%.6lf\n",ans);
    return 0;
}