【bzoj1914】[Usaco2010 OPen]Triangle Counting 数三角形 计算几何

补集转化
总数-不能过原点的三角形
我们发现先按照极角排序，然后枚举每个点，这个点与原点连线的一侧任选两个点与这个点围成的三角形都不能过原点
那么怎么计算呢？
其实按顺序枚举每个点，每个点只计算一遍，则正好不重不漏的计算了每一种情况
我们发现只要顺序枚举就可以了，按极角序枚举时，每次向后移一个点，会在原来的点的个数上-1
假如枚举到点i，记一个右端点r，在当前节点与原点所在的直线右侧的点的集合一定是[i+1,r]
注意！！！r是会循环回来的，就是会从n循环到1,

好难想呀！！！这都什么神题呀！！！

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<iostream>
#define maxn 100010

using namespace std;

struct yts
{
	long long x,y;
	double k;
}a[maxn];

long long operator*(yts a,yts b)
{
	return a.x*b.y-a.y*b.x;
}

int n,cnt;
long long ans,sum;

bool cmp(yts x,yts y)
{
	return x.k<y.k;
}

int main()
{
	scanf("%d",&n);
	for (int i=1;i<=n;i++) {scanf("%lld%lld",&a[i].x,&a[i].y);a[i].k=atan2(a[i].y,a[i].x);}
	sort(a+1,a+n+1,cmp);
	int r=2,cnt=0;
	for (int i=1;i<=n;i++)
	{
		while (r!=i && a[i]*a[r]>=0) r=r%n+1,cnt++;
		ans+=(long long)cnt*(cnt-1)/2;
		cnt--;
	}
	printf("%lld\n",(long long)n*(n-1)*(n-2)/6-ans);
	return 0;
}