2016 UESTC Training for Data Structures G - 郭大侠与阴阳家 CDOJ 1337 强行map

G - 郭大侠与阴阳家

给n个二维平面上的点，n属于[1,2000]，x和y是[1,1000000000]，唔，这题好像没用什么数据结构，可以用map，但是感觉强行map，不map感觉完全没问题

 

坑点有两个，重点，共线，重点简单，把所有点sort一遍，然后unique一发就可以，

共线，就需要先建立向量然后再处理了，

首先，能 组成平行四边形的条件是这两条向量平行且模长相等，所以我每个向量我就存储他的斜率和模长，还有向量所在直线与y轴的交点，当然这就涉及到一个问题，斜率为无穷时，而且此时向量所在直线与y轴没有焦点或有无穷个交点，这时候我就取一个很大的数1000010000，因为按所给题目的条件，这个数比所以可能的斜率都大，然后我对于这种情况，我存储它与x轴的交点，因为我存交点是为了判断共线，所以我只要这个交点能区分平行向量是否共线就行

然后，把这些向量sort，按斜率升序排序，斜率相同按模长升序排序，

然后就依次扫过去，处理那些平行且模长相等的向量，再来一个二重循环，枚举它们能组成多少个平行四边形

然后因为我们只判断了一组向量，实际中平行四边形的两组对边都是平行且模长相等的，所以我们得到的答案再除2就好了，另：反正我是用long long存的

代码：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
#define ll long long
#define maxn 2005
#define maxk 1000010000
int n;
struct Node
{
	ll x, y;
	bool operator <(Node a)const
	{
		return x < a.x || (x == a.x&&y < a.y);
	}
	bool operator ==(Node a)const
	{
		return (x == a.x) && (y == a.y);
	}
}node[maxn];
struct Edge
{
	ll length;
	double k, b;
	bool operator <(Edge t)const
	{
		if (fabs(k - t.k)>1e-6)
			return k < t.k;
		else
			return length < t.length;
	}
}edge[maxn*maxn];
inline bool equal_v(Edge m, Edge n)
{
	return fabs(m.k - n.k) < 1e-6&&m.length == n.length;
}
inline ll comb_2(ll x)
{
	return x*(x - 1) / 2;
}
int main()
{
	//freopen("input.txt", "r", stdin);
	scanf("%d", &n);
	for (int i = 0; i < n; ++i)
		scanf("%lld%lld", &node[i].x, &node[i].y);
	sort(node, node + n);
	n = unique(node, node + n) - node;
	int cnt = 0;
	ll ans = 0;
	for (int i = 0; i < n; ++i)
	{
		for (int j = i + 1; j < n; ++j)
		{
			if (node[j].x == node[i].x)
			{
				edge[cnt].k = maxk;
				edge[cnt].b = node[i].x;
			}
			else
			{
				edge[cnt].k = (node[j].y - node[i].y)*1.0 / (node[j].x - node[i].x);
				edge[cnt].b = node[i].y*1.0 - edge[cnt].k*node[i].x;
			}
			edge[cnt].length = (node[j].x - node[i].x)*(node[j].x - node[i].x) + (node[j].y - node[i].y)*(node[j].y - node[i].y);
			++cnt;
		}
	}
	sort(edge, edge + cnt);
	/*for (int i = 0; i < cnt; ++i)
	printf("11  %f %f %f\n", edge[i].k, edge[i].a, edge[i].b);
	printf("\n");*/
	/*for (int i = 0; i < cnt; ++i)
	printf("%f %f %f\n", edge[i].k, edge[i].a, edge[i].b);
	printf("\n");*/
	for (int i = 0; i < cnt;)
	{
		int j;
		for (j = i + 1; j < cnt; ++j)
		{
			if (!equal_v(edge[i], edge[j]))
				break;
		}
		for (int k = i; k < j; ++k)
		{
			for (int l = k + 1; l < j; ++l)
			{
				if (fabs(edge[k].b - edge[l].b)>1e-6)
					++ans;
			}
		}
		i = j;
	}
	ans >>= 1;
	printf("%lld\n", ans);
	//system("pause");
	//while (1);
	return 0;
}