poj2002--Squares(n个点求正方形个数)

标签：des   style   blog   io   ar   os   sp   for   strong   

Description

Input

Output

Sample Input

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

Sample Output

1
6
1

给出n个点，求出有可以组成多少个正方形？

枚举对角的两个点，然后求解出其他的两个点，将这两个点带入到n个点中查找，可以hash 或 二分。

已知对角线的点 (x1,y1) (x3,y3)  求出中点( (x1+x3)/2 , (y1+y3)/2 ) ->（X,Y） 用对角线的一个点减去中点得到(x,y)，那么其他的两个点就是( x1-Y,y1+X ) （x1+Y,y1-X）

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <vector>
using namespace std ;
#define eqs 1e-9
struct node{
    double x , y ;
}p[1100] ;
bool cmp(node a,node b)
{
    return ( a.x < b.x || ( a.x == b.x && a.y < b.y ) ) ;
}
bool judge(double x,double y,int n)
{
    int low = 0 , mid , high = n-1 ;
    while( low <= high )
    {
        mid = (low + high) / 2 ;
        if( fabs(p[mid].x-x) < eqs && fabs(p[mid].y-y) < eqs )
            return true ;
        else if( p[mid].x-x > eqs || ( fabs(p[mid].x-x) < eqs && p[mid].y-y > eqs ) )
            high = mid - 1 ;
        else
            low = mid + 1 ;
    }
    return false ;
}
int main()
{
    int n , i , j , num ;
    double x , y , xx , yy ;
    while( scanf("%d", &n) && n )
    {
        num = 0 ;
        for(i = 0 ; i < n ; i++)
        {
            scanf("%lf %lf", &p[i].x, &p[i].y) ;
        }
        sort(p,p+n,cmp) ;
        for(i = 0 ; i < n ; i++)
        {
            for(j = i+1 ; j < n ; j++)
            {
                if( i == j ) continue ;
                x = (p[i].x+p[j].x)/2 ;
                y = (p[i].y+p[j].y)/2 ;
                xx = p[i].x - x ;
                yy = p[i].y - y ;
                if( judge(x+yy,y-xx,n) && judge(x-yy,y+xx,n) )
                {
                    //printf("(%.1lf,%.1lf) (%.1lf,%.1lf) (%.1lf,%.1lf) (%.1lf,%.1lf)\n", p[i].x, p[i].y , p[j].x, p[j].y,x+yy,y-xx,x-yy,y+xx ) ;
                    num++ ;
                }
            }
        }
        printf("%d\n", num/2) ;
    }
    return 0;
}

poj2002--Squares(n个点求正方形个数)

标签：des   style   blog   io   ar   os   sp   for   strong   

原文地址：http://blog.csdn.net/winddreams/article/details/41853127