WA了好久,发现没开long long。
n个房子选3个有n*(n-1)*(n-2)/6种方案。考虑枚举3个点,超时。
考虑对于4个点,有2种位置关系:可构成凸四边形或凹四边形。对于每个凸四边形,过3点可以构成4种圆,有两种可以包围另一个点,对于每个凹四边形,有一种可以包围另一个点。答案就等于(凸四边形个数*2+凹四边形个数)/总方案数。
显然凸四边形个数+凹四边形个数=四边形个数=n*(n-1)*(n-2)*(n-3)/24。所以只要求出凹四边形个数,就能得到答案。
对于每个点i,设为凹四边形中间的点,枚举其他3个点不能包围它的方案数。以点i为原点,将其他点顺时针排序,枚举点j,用维护多少点在直线ij同侧,计算出以j为这三个点顺时针第一个点,构成三角形不能包围i的方案数,即可求出凹四边形个数,得到答案。
#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#define ll long long
using namespace std;
struct dot{
ll x,y;
double a;
}a[1505],t[1505];
ll n,Ans,x;
bool operator<(dot a,dot b)
{
return a.a<b.a;
}
ll work(ll x)
{
t[0]=a[x];
for (ll i=1,j=0;i<=n;i++)
if (i!=x) t[++j]=a[i];
for (ll i=1;i<n;i++)
t[i].a=atan2( (t[i].y-t[0].y),(t[i].x-t[0].x) );
sort(t+1,t+n);
ll R=2,p=0,sum=(n-1)*(n-2)*(n-3)/6;
for(ll i=1;i<n;i++)
{
while( (t[i].x-t[0].x)*(t[R].y-t[0].y)-(t[i].y-t[0].y)*(t[R].x-t[0].x)>=0 )
{
R=R%(n-1)+1;p++;
if(R==i)break;
}
sum-=p*(p-1)/2;
p--;
}
return sum;
}
int main()
{
scanf("%lld",&n);
if (n==3){puts("3.00");return 0;}
for (ll i=1;i<=n;i++)
scanf("%lld%lld",&a[i].x,&a[i].y);
for (ll i=1;i<=n;i++)
Ans+=work(i);
x=n*(n-1)*(n-2)*(n-3)/24-Ans;
double k=x*2+Ans;
k/=n*(n-1)*(n-2)/6;
printf("%.6f",k+3);
return 0;
}