poj2002Squares(点集组成正方形数）

链接

可以枚举两个点，因为是正方形两外两点可以由已知求出，据说可以根据三角形全等求出下列式子，数学渣不会证。。。

已知： (x1,y1)  (x2,y2)

则：   x3=x1+(y1-y2)   y3= y1-(x1-x2)

x4=x2+(y1-y2)   y4= y2-(x1-x2)

或

x3=x1-(y1-y2)   y3= y1+(x1-x2)

x4=x2-(y1-y2)   y4= y2+(x1-x2)

然后就可以hash或者二分做了，这里只用hash做的

应该算是简单的hash解决冲突的应用，放一个邻接表里。

两点需正反枚举两次，才能保证两种位置的正方形都被枚举到。

最后的结果需要除4，因为重复枚举了。

 1 #include <iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<algorithm>
 5 #include<stdlib.h>
 6 #include<vector>
 7 #include<cmath>
 8 #include<queue>
 9 #include<set>
10 using namespace std;
11 #define N 1010
12 #define mod 99991
13 #define LL long long
14 #define INF 0xfffffff
15 const double eps = 1e-8;
16 const double pi = acos(-1.0);
17 const double inf = ~0u>>2;
18 struct point
19 {
20     int x,y;
21     point(int x=0,int y=0):x(x),y(y){}
22 }p[N],o[N];
23 int next[N],head[mod],t;
24 void insert(int i)
25 {
26     int key = (p[i].x*p[i].x+p[i].y*p[i].y)%mod;
27     next[t] = head[key];
28     o[t].x = p[i].x;
29     o[t].y = p[i].y;
30     head[key] = t++;
31 }
32 int find(point a)
33 {
34     int  key = (a.x*a.x+a.y*a.y)%mod;
35     int i;
36     for(i = head[key] ; i!= -1 ; i = next[i])
37     {
38         if(o[i].x==a.x&&o[i].y == a.y) return 1;
39     }
40     return 0;
41 }
42 bool cmp(point a,point b)
43 {
44     if(a.x==b.x)
45     return a.y<b.y;
46     return a.x<b.x;
47 }
48 int main()
49 {
50     int n,i,j;
51     while(scanf("%d",&n)&&n)
52     {
53         memset(head,-1,sizeof(head));
54         t = 0;
55         for(i = 1; i <= n; i++)
56         {
57             scanf("%d%d",&p[i].x,&p[i].y);
58             insert(i);
59         }
60         //sort(p+1,p+n+1,cmp);
61         int ans = 0;
62         for(i = 1; i <= n ;i++)
63             for(j = 1 ; j <= n; j++)
64             {
65                 if(i==j) continue;
66                 point p1,p2;
67                 p1.x = p[i].x+(p[i].y-p[j].y);
68                 p1.y = p[i].y-(p[i].x-p[j].x);
69                 p2.x = p[j].x+(p[i].y-p[j].y);
70                 p2.y = p[j].y-(p[i].x-p[j].x);
71                 if(!find(p1)) continue;
72                 if(!find(p2)) continue;
73                 ans++;
74             }
75         printf("%d\n",ans/4);
76     }
77     return 0;

View Code