生活随笔
收集整理的這篇文章主要介紹了
poj 1981(单位圆覆盖点)
小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.
題意:一個單位圓最多能覆蓋平面上多少點。
解題思路:一個覆蓋最多點的圓,必然至少有兩個點在圓上。枚舉兩個點,求過這兩個點的單位圓,判斷有多少個點在圓中,枚舉N^2,判斷N
參考博客:http://www.cnblogs.com/-sunshine/archive/2012/10/11/2719859.html
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;const int maxn = 305;
const double eps = 1e-8;
struct Point
{double x,y;Point(){}Point(double _x,double _y){x = _x, y = _y;}
}p[maxn];
int n;double dis(Point p1,Point p2)
{return sqrt((p1.x - p2.x) * (p1.x - p2.x) + (p1.y - p2.y) * (p1.y - p2.y));
}Point get_circle(Point p1,Point p2)
{Point mid = Point((p1.x+p2.x)/2,(p1.y+p2.y)/2);double angle = atan2(p1.x-p2.x,p2.y-p1.y);double d = sqrt(1-dis(p1,mid)*dis(p1,mid));return Point(mid.x+d*cos(angle),mid.y+d*sin(angle));
}int main()
{while(scanf("%d",&n),n){for(int i = 0; i < n; i++)scanf("%lf%lf",&p[i].x,&p[i].y);int ans = 1;for(int i = 0;i < n; i++){for(int j = i + 1; j < n; j++){if(dis(p[i],p[j]) > 2.0) continue;Point central = get_circle(p[i],p[j]);int cnt = 0;for(int k = 0;k < n; k++)if(dis(central,p[k]) < 1.0 + eps)cnt++;ans=max(ans,cnt);}}printf("%d\n",ans);}return 0;
}
總結(jié)
以上是生活随笔為你收集整理的poj 1981(单位圆覆盖点)的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
如果覺得生活随笔網(wǎng)站內(nèi)容還不錯,歡迎將生活随笔推薦給好友。