UVA11529 - Strange Tax Calculation

链接

https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2524

题解

这题答案从直观上来看，分子每个三角形内部点数之和，分母是三角形数目
其实分子也等于每个点所在三角形数之和


可以发现
如果一个点在一个三角形的内部，那么我连接这个点和三角形上的任何一个点之后，另外两个点都分居连线的两侧
而当一个点在一个三角形的外部时，我就可以连接这个点和一个三角形的顶点，使得其余两个点在这条连线的一侧
利用这一点就可以解决这道题目
假设我现在要统计$A$点在多少三角形的内部，我以$A$为坐标原点建立坐标系
其余点的坐标变为$(x_i-x_A,y_i-y_A)$
极角排序之后，枚举一个点$j$，然后统计角度在${\theta }_j$到${\theta }_j+\pi$之间的点有多少个，然后搞一搞就好了

代码

//计算几何
#include <bits/stdc++.h>
#define maxn 5050
#define eps 1e-8
#define cl(x) memset(x,0,sizeof(x))
using namespace std;
int x[maxn], y[maxn], N;
long long cnt;
double th[maxn];
const double  pi=acos(-1);
int read(int x=0)
{
	int c, f=1;
	for(c=getchar();!isdigit(c);c=getchar())if(c=='-')f=-1;
	for(;isdigit(c);c=getchar())x=(x<<1)+(x<<3)+c-48;
	return f*x;
}
void init()
{
	int i;
	cnt=0;
	for(i=1;i<=N;i++)x[i]=read(), y[i]=read();
}
void work()
{
	int i, j, p, tot, t;
	for(i=1;i<=N;i++)
	{
		tot=t=0;
		for(j=1;j<=N;j++)if(j!=i)th[++tot]=atan2(y[j]-y[i],x[j]-x[i]);
		sort(th+1,th+tot+1);
		for(j=1;j<=tot;j++)th[tot+j]=th[j]+2*pi;
		for(j=1,p=1;j<=tot;j++)
		{
			for(;th[p+1]-th[j]<pi-eps;p++);
			t+=(p-j)*(p-j-1)/2;
		}
		cnt+=(N-1)*(N-2)*(N-3)/6-t;
	}
	printf("%.2lf\n",double(cnt)/N/(N-1)/(N-2)*6);
}
int main()
{
	int kase=0;
	while(N=read())
	{
		printf("City %d: ",++kase);
		init();
		work();
	}
	return 0;
}