cf #308 D. Vanya and Triangles (几何)

题目：http://codeforces.com/contest/554/problem/D

题意：给你n个点的坐标，n<2000,问这n个点能形成多少个三角形。

分析：简单的几何题。枚举所有直线，求出k和b，将k和b映射成hash值，存起来。然后枚举直线，求出所有不能做成三角形的情况，再用所有的可能情况减掉这些不可能的情况即可。

重点不在上面。在于精度问题。已经在代码里面标出来了。

代码：

#include <bits/stdc++.h>
using namespace std;

typedef long long LL;
typedef unsigned long long ULL;
const LL INF = 1E3+9;

const int maxn = 2007;

inline int MyAbs(int a)
{
	return a>0?a:-a;
}

inline LL cal(int x)
{
	return x*(x-1ll)*(x-2ll)/6ll;
}

inline LL getHash(double d)
{
	return LL(d*1E14);
}

typedef pair <LL,LL> pll;

map <pll ,set <int> > Map;
int x[maxn],y[maxn];
int main()
{
//	freopen("1.txt","r",stdin);
	int n;
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		scanf("%d%d",&x[i],&y[i]);
	for(int i=1;i<=n;i++)
	{
		for(int j=i+1;j<=n;j++)
		{
			int dx=x[i]-x[j];
			int dy=y[i]-y[j];
			static pll key;
			if(0==dx)
				key=make_pair(getHash(INF),getHash(x[i]));
			else
			{
				double k=double(dy)/dx;
				double b=(y[i]*x[j]-y[j]*x[i])/double(x[i]-x[j]);
				//double b=y[i]-k*x[i]  这样写会有比较严重的精度损失  
				key = make_pair(getHash(k),getHash(b));
			}
			Map[key].insert(i);
			Map[key].insert(j);
		}
	}
	LL ans=cal(n);
	for(const auto &it:Map)
		ans-=cal(it.second.size());
	cout<<ans;
	return 0;
}