哈希表：正方形

摘一下题目：

正方形

总时间限制: 3500ms 内存限制: 65536kB
描述
给定直角坐标系中的若干整点，请寻找可以由这些点组成的正方形，并统计它们的个数。
输入
包括多组数据，每组数据的第一行是整点的个数n(1<=n<=1000)，其后n行每行由两个整数组成，表示一个点的x、y坐标。输入保证一组数据中不会出现相同的点，且坐标的绝对值小于等于20000。输入以一组n=0的数据结尾。
输出
对于每组输入数据，输出一个数，表示这组数据中的点可以组成的正方形的数量。
样例输入
4
1 0
0 1
1 1
0 0
9
0 0
1 0
2 0
0 2
1 2
2 2
0 1
1 1
2 1
4
-2 5
3 7
0 0
5 2
0
样例输出
1
6
1
本来是个大水题，但是首先我把初高中数学知识全给忘了，其次数据量非常大，必须自己手写哈希表。
我最开始的思路是以两个点的中点作为哈希值把这些所有点放进各个桶里，再在每个桶的所有点里面检测。但是这样TLE了……因为n^2级别的数据放进n级别的桶里……再两两比较……代价……emmm
参考了别人的思路，是两两取出点，再计算另外两个点坐标，看是否有这两个点。
我看到的一个哈希表的代码是
最开始的辣鸡代码：

#include <iostream>
#include <map>
#include <vector>
#include <cmath>
#include <cstring>
using namespace std;
const int maxn = 1005;
struct point
{
   
    double x, y;
    point() {
   }
    point(double _x, double _y) : x(_x), y(_y) {
   }
} pts[maxn];
struct points
{
   
    point p1, p2;
    points() {
   }
    points(point &_p1, point &_p2) : p1(_p1), p2(_p2) {
   }
};
double dot(points& pp1, points& pp2)
{
   
    return ((pp1.p1.x - pp1.p2.x) * (pp2.p1.x - pp2.p2.x) + (pp1.p1.y - pp1.p2.y) * (pp2.p1.y - pp2.p2.y));
}
double dis(points& p)
{
   
    return sqrt((p.p1.x - p.p2.x) * (p.p1.x - p.p2.x) + (p.p1.y - p.p2.y) * (p.p1.y - p.p2.y));
}
bool issquare(points& pp1, points& pp2)
{
   
    if ((pp1.p1.x + pp1.p2.x) != (pp2.p1.x + pp2.p2.x) || ((pp1.p1.y + pp1.p2.y) != (pp2.p1.y + pp2.p2.y)))
        return 0;
    if (dot(pp1, pp2) != 0)
        return 0;
    if (dis(pp1) != dis(pp2)) return 0;
    return 1;
}
int hashtable(point a)
{
   
    return (long)((a.x * a.x + a.y * a.y)) %