日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當(dāng)前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

poj 1981(单位圆覆盖点)

發(fā)布時間:2025/3/16 编程问答 19 豆豆
生活随笔 收集整理的這篇文章主要介紹了 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)容還不錯,歡迎將生活随笔推薦給好友。