题目要求是给你n个点,求n个点形成的所有三角形中包含的点的总数的平均值。
题目可以理解为求每个点能在多少三角形内。我们枚举一个点x,然后将其他点按照极角排序,那么每一个点i到与点i成pi的范围内的所有点组成的三角形都不包含点x,然后我们就可以求出所有的解。
里面有一些小的需要注意的问题,就是有的点形成的极角会在(0,-pi),我们求解的时候会有一点小小的麻烦,我们可以将其+2pi来解决。
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
typedef long long LL;
using namespace std;
int N;
const double pi = acos(-1.0);
const int maxn = 1500*2;
double X[maxn],Y[maxn],K[maxn];
LL solve(int n)
{
int cnt = 0;
for(int i=0; i<N; i++)
{
if(i==n)continue;
K[cnt++] = atan2(Y[i]-Y[n],X[i]-X[n]);
K[cnt++] = atan2(Y[i]-Y[n],X[i]-X[n])+2*pi;
}
sort(K,K+cnt);
// int _count = 0;
int pos = 0;
LL num =1;
LL res = 0;
for(int i=0; i<N-1; i++)
{
num--;
while(K[pos]<K[i]+pi)
{
pos++;
num++;
}
res+=(num-1)*(num-2)/2;
}
LL sum = (LL)(N-1)*(LL)(N-2)*(LL)(N-3)/(LL)6;
return sum - res;
}
int main()
{
int kas = 1;
while(~scanf("%d",&N)&&N)
{
LL res = 0;
for(int i=0; i<N; i++)
{
scanf("%lf%lf",&X[i],&Y[i]);
}
for(int i=0; i<N; i++)
{
res += solve(i);
}
LL sum = (LL)(N)*(LL)(N-1)*(LL)(N-2)/(LL)6;
printf("City %d: %.2lf\n",kas++,(double)res/(double)sum);
}
return 0;
}