poj 1981

最新推荐文章于 2020-10-09 19:15:53 发布

sky_zdk

最新推荐文章于 2020-10-09 19:15:53 发布

阅读量400

点赞数

CC 4.0 BY-SA版权

分类专栏： ACM poj 文章标签：计算几何

本文链接：https://blog.youkuaiyun.com/sky_zdk/article/details/68962684

ACM poj 专栏收录该内容

201 篇文章

订阅专栏

本文介绍了一种解决点集覆盖问题的方法，通过计算不同点之间的距离和角度来确定最大覆盖范围。使用C++实现算法，并考虑了点间距离限制及角度偏移，最终输出覆盖最多点的数量。

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

在判断的时候不要用等于号！！！！

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#define Mn 300
using namespace std;
#define sqr(x) ((x) * (x))
int n;
struct point
{
	double x,y;
}p[300+16];
double dis(const point a,const point b)
{
	return sqrt(sqr(a.x-b.x)+sqr(a.y-b.y));
}
struct alpha
{
    double v;
    bool flag;
    bool friend operator <(const alpha &a,const alpha &b)
    {
       return a.v < b.v;
    }
}alp[Mn+15];
void solve()
{
    int i,j;
    int MAX =1;
    for( i = 0;i < n;i++)
    {
       int t = 0;
       for( j = 0;j < n;j++)
       {
           if(i == j)
              continue;
           double theta,phi,D;
           D = dis(p[i],p[j]);
           if(D > 2.0 )
              continue;
           theta = atan2(p[j].y - p[i].y,p[j].x - p[i].x);
           phi = acos(D / (2.0 ));
           alp[t].v = theta - phi;
           alp[t++].flag = true;
           alp[t].v = theta + phi;
           alp[t++].flag = false;
       }
       sort(alp,alp + t);
       int sum = 1;
       for( j = 0;j < t;j++)
       {
           if(alp[j].flag)
              sum ++;
           else
              sum --;
           if(sum > MAX)
              MAX = sum;
       }
    }
    printf("%d\n",MAX);
}

int main()
{
    int i;
    while(~scanf("%d",&n)&&n)

    {
       for(i=0;i<n;i++)
		   scanf("%lf%lf",&p[i].x,&p[i].y);
       solve();
    }
}