poj-2002 Squares

题意很明确，有n个点，输入n个点的坐标x,y得出这些点能构成几个正方形。

思路：由于坐标很大直接存二维数组是肯定不行的，所以较为通用的方法是用哈希表，该题给的点的数量有点少，只有1000个点，有可能开一个1000的数组存点能过，不过还是推荐哈希表过。本人表头存法是用的坐标的平方和，表头开了1w+10，稍微大些来减少表中在查找的次数。并且对任意两点求另两点中，两点的选取作为对角点，求另两个对角点，减少重复的查找。

对于坐标的计算，设对角线上两点为A(x1,y1)，B(x2,y2)，未知点为C(x3,y3)，D(x4,y4)。AB中点O为((x1+x2)/2,(y1+y2)/2),向量OB=((x2-x1)/2,(y2-y1)/2)。向量OC=((y1-y2)-2,(x2-x1)/2)，C坐标即OC+O，可得x3=(y1-y2+x1+x2)/2,y3=(x2-x1+y1+y2)/2,同理D为x4=(y2-y1+x1+x2)/2,y4=(x1-x2+y1+y2)/2.如果C，D点中有非整数点，直接跳过查找。

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#define INF 0x3f3f3f3f
#define mod 1000000007

using namespace std;

struct node
{
int x,y;
int next;
}d[1010];
int head[10010];
int top,flag;
void creat(int x,int y)
{
int n=(x*x+y*y)%10007;
d[top].x=x;
d[top].y=y;
d[top].next=head
;
head
=top++;
}
void Find(int x,int y)
{
int i=head[(x*x+y*y)%10007];
while(i)
{
if(d[i].x==x&&d[i].y==y)
{flag++;break;}
i=d[i].next;
}
}
int main()
{
int n,x,y,i,j;
while(scanf("%d",&n),n)
{
top=1;
for (i=0;i<n;i++)
{
scanf("%d%d",&x,&y);
creat(x,y);
}
int x1,x2,y1,y2,x3,y3,x4,y4;
int ans=0;
for (i=1;i<top-1;i++)
{
for (j=i+1;j<top;j++)
{
flag=0;
x1=d[i].x;y1=d[i].y;
x2=d[j].x;y2=d[j].y;
x3=(y1-y2+x1+x2);
y3=(x2-x1+y1+y2);
x4=(y2-y1+x1+x2);
y4=(x1-x2+y1+y2);
if(x3%2!=0||y3%2!=0||x4%2!=0||y4%2!=0) continue;
x3/=2,y3/=2;
x4/=2,y4/=2;
Find(x3,y3);
Find(x4,y4);
if(flag==2) ans++;
}
}
printf("%d\n",ans/2);
}
return 0;
}