【HDU4082】相似三角形计数

1题目链接。题面很长，但是题意比较简单，大意就是给你n个点，找出这些点组成的三角形中相似的又多少个。

2.emmm，首先我们分析一下，n个点理论上来说最多会组成C(n,3)个三角形，但是题目中给的数据很小，所以我们暴力枚举即可。那么我们如何判断两个三角形相似呢？最简单的写法肯定是把三个角求出来，然后比较这些角度即可。暴力枚举虽然是没有错，但是我们在比较相似的时候还是会出现很多的问题，一个简单的处理就是我们把角度排序，那么对比就比较简单了，然后这样我们还是需要比较三次，比较麻烦。所以一个想法就是我们把三角形HASH一下（这个思路好），我们把三角形通过角度hash成一个志，然后用map存下来，这样最后就找到这个map中最大的数据即可。emmm，很好。码之，WA了。我很镇定，因为A了才是让我震惊的。计算几何不踩坑1A对我来说痴人说梦，然后DEBUG2h（codeing 10分钟，找bug 2h）最后我发现点可能重复（omfk），改了之后A了。真好。

#include<bits/stdc++.h>
#include<map>
#include<iostream>
#include<set>
using namespace std;
typedef long long ll;
#pragma warning(disable:4996)
const double eps = 1e-9;
map<ll, int>mp;
struct Point
{
	int x, y;
}p[500];
double ACOS(double x)
{
	if (abs(x) > 1)
		x = 1;
	return acos(x);
}
double dis(Point a, Point b)
{
	return sqrt(pow(a.x - b.x, 2) + pow(a.y - b.y, 2));
}
void solve(Point  a, Point b, Point c)
{
	double la = dis(b, c);
	double lb = dis(a, c);
	double lc = dis(a, b);
	if (la + lb <= lc + eps)return;
	if (lb + lc <= la + eps)return;
	if (la + lc <= lb + eps)return;
	double A = acos((lb*lb + lc * lc - la * la) / (2 * lb*lc));
	double B = acos((la*la + lc * lc - lb * lb) / (2 * la*lc));
	double C = acos((la*la + lb * lb - lc * lc) / (2 * la*lb));
	if (A < eps || B < eps || C < eps)return;
	int t1 = (int)(A * 10000);
	int t2 = (int)(B * 10000);
	int t3 = (int)(C * 10000);
	if (t1 > t2)swap(t1, t2);
	if (t1 > t3)swap(t1, t3);
	if (t2 > t3)swap(t2, t3);
	if (t1 == 0)return;
	long long t = t1 * 1000000 * 1000000 + t2 * 1000000 + t3;
	mp[t]++;
}
const int N = 1000;
int vis[N][N];
int main()
{
	int n;
	while (~scanf("%d", &n) && n)
	{
		memset(vis, 0, sizeof(vis));
		int t = 0;
		for (int i = 0; i < n; i++)
		{
			scanf("%d%d", &p[t].x, &p[t].y);
			if (vis[p[t].x + 100][p[t].y + 100] == 0)
			{
				vis[p[t].x + 100][p[t].y + 100] = 1;
				t++;
			}
		}
		n = t;
		mp.clear();
		for (int i = 0; i < n; i++)
		{
			for (int j = i + 1; j < n; j++)
			{
				for (int k = j + 1; k < n; k++)
				{
					solve(p[i], p[j], p[k]);
				}
			}
		}
		map<ll, int>::iterator iter;
		int ans = 0;
		for (iter = mp.begin(); iter != mp.end(); iter++)
		{
			ans = max(ans, iter->second);
		}
		cout << ans << endl;
	}
	return 0;
}