POJ2002

 

题目大意：给出散点，求构成正方形的个数有多少？

思路：给出的是散点，通过2个点来枚举，但是通过公式求得的点   如果是通过线性查找来验证肯定是超时的，因为执行10^9次（O（N^3）），所以要考虑哈希表和二分查找

//Memory Time
//652K  1438MS

#include<iostream>
using namespace std;

const int prime=1999;  //长度为2n区间的最大素数 (本题n=1000)

//其他prime可取值:
// 1n  区间: 997   1704ms
// 2n  区间: 1999  1438ms
// 8n  区间: 7993  1110ms
// 10n 区间: 9973  1063ms
// 30n 区间: 29989 1000ms
// 50n 区间: 49999 1016ms
// 100n区间: 99991 1000ms

//为了尽量达到key与地址的一一映射，hash[]至少为1n，
//当为1n时，空间利用率最高，但地址冲突也相对较多，由于经常要为解决冲突开放寻址，使得寻找key值耗时O(1)的情况较少
//当n太大时，空间利用率很低，但由于key分布很离散，地址冲突也相对较少，使得寻找键值耗时基本为O(1)的情况

typedef class
{
public:
  int x,y;
}Node;

typedef class HashTable
{
public:
  int x,y;   //标记key值对应的x,y
  HashTable* next;  //当出现地址冲突时，开放寻址

  HashTable()  //Initial
  {
   next=0;
  }
}Hashtable;

Node pos[1001];
Hashtable* hash[prime];   //hash[]是指针数组，存放地址

void insert_vist(int k)
{
int key=((pos[k].x * pos[k].x)+(pos[k].y * pos[k].y))%prime +1;   //+1是避免==0
                                                                   //使key从[0~1998]后移到[1~1999]
if(!hash[key])
{
  Hashtable* temp=new Hashtable;
  temp->x=pos[k].x;
  temp->y=pos[k].y;
  hash[key]=temp;
}
else   //hash[key]已存地址，地址冲突
{
  Hashtable* temp=hash[key];
 
  while(temp->next)     //开放寻址，直至next为空
   temp=temp->next;

  temp->next=new HashTable;   //申请新结点，用next指向，记录x、y
  temp->next->x=pos[k].x;
  temp->next->y=pos[k].y;
}
return;
}

bool find(int x,int y)
{
int key=((x * x)+(y * y))%prime +1;

if(!hash[key])   //key对应的地址不存在
  return false;
else
{
  Hashtable* temp=hash[key];

  while(temp)
  {
   if(temp->x==x && temp->y==y)
    return true;

   temp=temp->next;
  }
}

return false;
}

int main(void)
{
int n;
while(cin>>n)
{
  if(!n)
   break;

  memset(hash,0,sizeof(hash));   //0 <-> NULL

  for(int k=1;k<=n;k++)
  {
   cin>>pos[k].x>>pos[k].y;
   insert_vist(k);   //插入哈希表，标记散点
  }

  int num=0;  //正方形的个数
  for(int i=1;i<=n-1;i++)
   for(int j=i+1;j<=n;j++)
   {
    int a=pos[j].x-pos[i].x;
    int b=pos[j].y-pos[i].y;

    int x3=pos[i].x+b;
    int y3=pos[i].y-a;
    int x4=pos[j].x+b;
    int y4=pos[j].y-a;
   
    if(find(x3,y3) && find(x4,y4))
     num++;

    x3=pos[i].x-b;
    y3=pos[i].y+a;
    x4=pos[j].x-b;
    y4=pos[j].y+a;

    if(find(x3,y3) && find(x4,y4))
     num++;
   }

  cout<<num/4<<endl;  //同一个正方形枚举了4次
}
return 0;
}

另转（二分查找）：http://kenby.iteye.com/blog/949273