POJ 1106 Transmitters(计算几何)

最新推荐文章于 2022-02-17 21:58:54 发布

转载最新推荐文章于 2022-02-17 21:58:54 发布 · 85 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/naix-x/p/3363630.html

本文深入探讨了计算几何领域的经典问题——如何在圆内找到最多数量的直线，这些直线能够平分圆。通过使用枚举法、距离公式及交叉乘积等数学工具，作者详细阐述了解决该问题的算法步骤，并提供了代码实现，帮助读者理解并应用到实际场景中。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目链接

切计算几何，感觉计算几何的算法还不熟。此题，枚举线段和圆点的直线，平分一个圆

 1 #include <iostream>
 2 #include <cstring>
 3 #include <cstdio>
 4 #include <cstdlib>
 5 #include <cmath>
 6 using namespace std;
 7 #define eps 1e-8
 8 struct point
 9 {
10     double x,y;
11 }p[201];
12 double dis(point p1,point p2)
13 {
14     return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
15 }
16 double xmult(point p1,point p2,point p0)
17 {
18     return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
19 }
20 int main()
21 {
22     point r;
23     double d;
24     int n,m,i,j,ans,temp;
25     while(scanf("%lf%lf%lf",&r.x,&r.y,&d)!=EOF)
26     {
27         if(d < 0) break;
28         scanf("%d",&n);
29         m = 0;
30         for(i = 0;i < n;i ++)
31         {
32             scanf("%lf%lf",&p[i].x,&p[i].y);
33         }
34         for(i = 0;i < n;i ++)
35         {
36             if(dis(p[i],r) > d)
37             continue;
38             p[m++] = p[i];
39         }
40         ans = 0;
41         for(i = 0;i < m;i ++)
42         {
43             temp = 0;
44             for(j = 0;j < m;j ++)
45             {
46                 if(xmult(p[i],p[j],r) >= 0)
47                 temp ++;
48             }
49             ans = max(ans,temp);
50         }
51         printf("%d\n",ans);
52     }
53     return 0;
54 }