！HDU 4173 到点的距离不超过2.5，找最多能被满足的点的个数-简单几何

来源：互联网发布：我的淘宝是什么意思编辑：程序博客网时间：2024/06/05 21:55

题意：有n个人参加派对，但是条件是派对地点到他家的距离不能超过2.5，现在要你找一个最佳的派对地点让最多人参加派对。

分析：

题目看起来挺难的，怎么求范围然后包含点？其实一个圆心+半径不就代表一个圆了吗。

枚举求两个点的半径为2.5的圆的圆心，保存起来，然后用这些圆心求所有点到圆心的距离，记录距离小于等于2.5的点的个数，更新结果。200*199/2*200不会超时，计算就是高中基本的数学知识。仔细一些耐心一点就是了。注意计算得到的浮点数判定要用eps。

代码：

#include<iostream>#include<cmath>#include<algorithm>#define eps 1e-8using namespace std;int n,ans,cnt;double x[300],y[300];struct node{double x,y;}center[40000];double dis(double a,double b,double c,double d){return sqrt((a-c)*(a-c)+(b-d)*(b-d));}void getcenter(int i,int j){double midx=(x[i]+x[j])/2.0;double midy=(y[i]+y[j])/2.0;double k,b;if(x[i]-x[j]==0){double x1=sqrt(6.25-(midy-y[i])*(midy-y[i]))+x[i];double x2=x[i]-sqrt(6.25-(midy-y[i])*(midy-y[i]));center[cnt].x=x1,center[cnt++].y=midy;center[cnt].x=x2,center[cnt++].y=midy;}else if(y[i]-y[j]==0){double y1=sqrt(6.25-(midx-x[i])*(midx-x[i]))+y[i];double y2=y[i]-sqrt(6.25-(midx-x[i])*(midx-x[i]));center[cnt].x=midx,center[cnt++].y=y1;center[cnt].x=midx,center[cnt++].y=y2;}else{double k1=(y[i]-y[j])/(x[i]-x[j]);k=-1.0/k1;    b=midy-k*midx;double tmp1=(x[i]+k*y[i]-k*b);double tmp2=(x[i]+k*y[i]-k*b)*(x[i]+k*y[i]-k*b)-(k*k+1)*(x[i]*x[i]+y[i]*y[i]+b*b-2*b*y[i]-6.25);double x1=(tmp1+sqrt(tmp2))/(k*k+1);double x2=(tmp1-sqrt(tmp2))/(k*k+1);double y1=k*x1+b;double y2=k*x2+b;center[cnt].x=x1,center[cnt++].y=y1;center[cnt].x=x2,center[cnt++].y=y2;}}int main(){     while(cin>>n){     cnt=0,ans=1;     for(int i=0;i<n;i++) cin>>x[i]>>y[i];     for(int i=0;i<n;i++){     for(int j=i+1;j<n;j++){     if(dis(x[i],y[i],x[j],y[j])<5.0+eps) getcenter(i,j);     }     }     for(int i=0;i<cnt;i++){     int tot=0;     for(int j=0;j<n;j++){     if(dis(x[j],y[j],center[i].x,center[i].y)<2.5+eps) tot++;     }     ans=max(ans,tot);     }     cout<<ans<<endl;     }}

0 0