BZOJ 4080 Wf2014 Sensor Network 随机化

最新推荐文章于 2021-02-25 20:25:37 发布

原创最新推荐文章于 2021-02-25 20:25:37 发布 · 1.9k 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#BZOJ #BZOJ4080 #随机化

BZOJ 专栏收录该内容

667 篇文章

订阅专栏

本文介绍了一种使用随机化贪心算法解决平面内最大点集问题的方法，该方法通过随机序列和局部贪心策略找到满足两点间距离不超过给定阈值的最大点集。代码实现包括点类定义、距离计算、随机化排序和贪心选择过程。

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

题目大意：给定平面上的 $n$ 个点，求一个最大的点集，使得两两之间距离不超过 $d$

爆搜T到死，加什么剪枝都没用……

随机化大法好

每次随机一个序列，依次贪心加入，然后更新答案

据说很靠谱？反正写完直接过了

#include <bitset>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 110
using namespace std;
typedef bitset<M> abcd;
struct Point{
    int x,y;
    friend istream& operator >> (istream &_,Point &p)
    {
        return scanf("%d%d",&p.x,&p.y),_;
    }
    friend bool operator < (const Point &p1,const Point &p2)
    {
        return p1.x < p2.x ;
    }
    friend int Distance(const Point &p1,const Point &p2)
    {
        return (p1.x-p2.x)*(p1.x-p2.x) + (p1.y-p2.y)*(p1.y-p2.y) ;
    }
}points[M];
int n,d;
abcd a[M];
abcd now,able,ans;
int main()
{
    int i,j;
    cin>>n>>d;
    for(i=1;i<=n;i++)
        cin>>points[i];
    for(i=1;i<=n;i++)
        for(j=i+1;j<=n;j++)
            if(Distance(points[i],points[j])<=d*d)
                a[i][j]=a[j][i]=true;
    static int order[M];
    for(i=1;i<=n;i++)
    {
        able[i]=true;
        order[i]=i;
    }
    for(j=1;j<=1000;j++)
    {
        abcd able=::able;
        now.reset();
        for(i=1;i<=n;i++)
            if(able[order[i]])
            {
                now[order[i]]=true;
                able&=a[order[i]];
            }
        if(now.count()>ans.count())
            ans=now;
        random_shuffle(order+1,order+n+1);
    }
    cout<<ans.count()<<endl;
    for(i=1;i<=n;i++)
        if(ans[i])
            printf("%d ",i);
    return 0;
}