HDU 1077 枚举圆心， 数学

从所有点中依次取两点，然后判断到两点m, n距离为半径1的点centre（易知center在m, n连线的中垂线上）。然后判断其余有多少点在以centre为圆心的圆内，记作ans，最后找出最大的ans,就是要求的值。

#include <iostream>
#include <cmath>
const double eps=1e-6;
const int MAXN=305;
struct point
{
double x,y;
};
struct point ps[MAXN];
int N;
double get_dis(int a,int b)
{
return (ps[a].x-ps[b].x)*(ps[a].x-ps[b].x)+(ps[a].y-ps[b].y)*(ps[a].y-ps[b].y);
}
void get_center(int a,int b)
{
double ox,oy,dx,dy,r,temp;
ox=(ps[a].x+ps[b].x)/2;
oy=(ps[a].y+ps[b].y)/2;
dx=ps[b].x-ps[a].x;
dy=ps[b].y-ps[a].y;
ps
.x=ox;ps
.y=oy;
temp=get_dis(N,b);
r=sqrt(1.0-temp);
if(fabs(dx)<eps)//竖直的情况
ps
.x-=r;//+或-都可以
else
{
double ang=atan(-dy/dx);
ps
.x-=r*sin(ang);//+或-都可以
ps
.y-=r*cos(ang);
}
}
int main()
{
int T,i,j,k,max;
scanf("%d",&T);
while(T--)
{
scanf("%d",&N);
for(i=0;i<N;i++)
scanf("%lf%lf",&ps[i].x,&ps[i].y);
int ans=1;
for(i=0;i<N;i++)//枚举O(n^2)
{
for(j=i+1;j<N;j++)
{
if(get_dis(i,j)>=4) continue;
get_center(i,j);
for(k=0,max=0;k<N;k++)
{
if(N-k+max<=ans) break;
double temp=sqrt(get_dis(k,N));
if(temp<=1.000001)
max++;
}
if(max>ans)
ans=max;
}
}
printf("%d\n",ans);
}
return 0;
}