poj2002

题意：最多1000个整点，问可以组成多少个正方形。

枚举在一条边上的其中两个点，计算出另外两个点的位置，然后判断是否存在。最后结果要除以4.

其中计算两个点的坐标利用全等三角形。

顺便学习了一下hash和multimap，参考自http://blog.csdn.net/qq574857122/article/details/9736535， multimap（多重映射器）可参考http://blog.csdn.net/acmmmm/article/details/9735379。

代码中的用了平方同余法， mod = 99991。

#include <iostream>
#include <cmath>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <cstdlib>
#include <iomanip>
using namespace std;

#define rep(i, a, b) for( i = (a); i <= (b); i++)
#define reps(i, a, b) for( i = (a); i < (b); i++)
#define pb push_back
#define ps push
#define mp make_pair
#define CLR(x,t) memset(x,t,sizeof x)
#define LEN(X) strlen(X)
#define F first
#define S second
#define Debug(x) cout<<#x<<"="<<x<<endl;


const double euler_r = 0.57721566490153286060651209;
const double pi = 3.141592653589793238462643383279;
const int inf=~0U>>1;
const int MOD = int(1e9) + 7;
const double EPS=1e-6;

//poj2002
//multimap or hash

const int mod = 99991;

typedef long long LL;

struct P
{
	int x, y;
}p[1010];

bool hash[100000];
multimap <int ,int> m;//multimap的创建
multimap<int,int>::iterator it;

bool have(P g)
{

	int k = g.x * g.x + g.y * g.y;
	k %= mod;
	if(!hash[k])return false;
	else
	{
		int num=m.count(k); //hash值为k的个数
		it=m.find(k); //迭代器到k的位置
		while(num--) //依次判断是不是和a结点值相同
		{
			int s=(*it).S; // s表示mutilmap第二个值，即为结点的下标
			if(p[s].x==g.x && p[s].y==g.y)return true;
			it++;
		}
		return false;
	}

}
int find(int i,int j){

	int ans=0;
	P a=p[i], b=p[j], c, d;

	//cal
	c.x=a.x+(a.y-b.y),c.y=a.y-(a.x-b.x);
	d.x=b.x+(a.y-b.y),d.y=b.y-(a.x-b.x);
    if(have(c)&& have(d)) ans++;

	//cal
	c.x=a.x-(a.y-b.y),c.y=a.y+(a.x-b.x);
	d.x=b.x-(a.y-b.y),d.y=b.y+(a.x-b.x);
	if(have(c)&& have(d)) ans++;

	return ans;
}
void insert(int num)
{
	int k = p[num].x * p[num].x + p[num].y * p[num].y;//平方同余法
	k %= mod;
	hash[k]=true;
	m.insert(pair<int,int>(k,num));//插入
}

int main()
{
	int n, i, j;
	while(scanf("%d",&n),n)
    {
		CLR(hash, 0);
		m.clear();
		rep(i, 1, n){scanf("%d%d",&p[i].x,&p[i].y);insert(i);}

		int ans=0;
		rep(i, 1, n) rep(j, i + 1, n) ans += find(i, j);
		printf("%d\n", ans >> 2);
	}
	return 0;
}