poj-2002-Squares--二分

确定两点(a,b),(c,b)..可确定另外两点，使之组成正方形

大概如上图所示，可得另外两点坐标

x1=a+b-d;

y1=b+c-a;

x2=c+b-d;

y2=d+c-a;

还有一种情况，在下面。。

x1=a-b+d;

y1=b-c+a;

x2=c-b+d;

y2=d-c+a;

然后二分检查下，是否存在这两个点。。

每个情况算了4次，，结果在除以4.。

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
int n;
struct node
{
	int x,y;
}p[1010],p1,p2;
int cmp(node a,node b)
{
	if(a.x==b.x)
		return a.y<b.y;
	return a.x<b.x;
}
int check(node a,node b)
{
	int l=0,r=n-1,f=0;
	while(l<=r)
	{
		int mid=(l+r)/2;
		if(p[mid].x==a.x){
			if(p[mid].y==a.y){
				f=1;
				break;
			}
			else if(p[mid].y>a.y)
				r=mid-1;
			else
				l=mid+1;
		}
		else if(p[mid].x>a.x)
			r=mid-1;
		else
			l=mid+1;
	}
	if(f==0)
		return 0;
	l=0,r=n-1,f=0;
	while(l<=r)
	{
		int mid=(l+r)/2;
		if(p[mid].x==b.x){
			if(p[mid].y==b.y){
				f=1;
				break;
			}
			else if(p[mid].y>b.y)
				r=mid-1;
			else
				l=mid+1;
		}
		else if(p[mid].x>b.x)
			r=mid-1;
		else
			l=mid+1;
	}
	if(f==0)
		return 0;
	else
		return 1;
}
int main()
{
	int i,j,k;
	while(scanf("%d",&n),n)
	{
		for(i=0;i<n;i++)
			scanf("%d%d",&p[i].x,&p[i].y);
		sort(p,p+n,cmp);
		int sum=0;
		for(i=0;i<n;i++)
		{
			for(j=i+1;j<n;j++)
			{
				p1.x=p[j].x+p[i].y-p[j].y;
				p1.y=p[j].y+p[j].x-p[i].x;
				p2.x=p[i].x+p[i].y-p[j].y;
				p2.y=p[i].y+p[j].x-p[i].x;
				if(check(p1,p2))
					sum++;
				p1.x=p[j].x-p[i].y+p[j].y;
				p1.y=p[j].y-p[j].x+p[i].x;
				p2.x=p[i].x-p[i].y+p[j].y;
				p2.y=p[i].y-p[j].x+p[i].x;
				if(check(p1,p2))
					sum++;
			}
		}
		printf("%d\n",sum/4);
	}
}