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[n];
    head[n]=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;
}