poj 2932 Coneology

最新推荐文章于 2020-10-02 19:30:05 发布

ch2010091

最新推荐文章于 2020-10-02 19:30:05 发布

阅读量1.1k

点赞数 1

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/ch2010091/article/details/7795171

计算几何专栏收录该内容

3 篇文章

订阅专栏

本文介绍了一种解决圆覆盖问题的扫描线算法，通过使用set集合存储不相交的圆，并利用圆心位置进行排序来高效判断圆之间的覆盖关系。

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

依然是圆有关的扫描线.感觉额.上一篇说得自己都看不太懂的样子.

这个题目吧。是要问你有哪个圆没有被其他的圆覆盖住.说白了就是他的被覆盖层数为0.和上一篇的题相比.简单了点。

同样的.没有相交的圆..

啊.直接说做法.加入set的集合的元素是.额.是不相互包含的圆.也就是说两两关系是相离.那么一个圆只要和它上下两圆无关系，自然也不会和这之外的圆有关系了.

然后这里判断一个圆的上下，就没有那么麻烦了.按圆心位置高下来就行.即能和当前圆发生关系的就只有它上下的圆。毕竟加入set的集合的圆都是没有关系的.

额......老是感觉没说清的感觉......

好吧.硬着头皮上我的cuo代码.....set的操作害我RE了好几次...

#include<iostream>
#include<cstdio>
#include<cstring>
#include<set>
#include <cmath>
#include<algorithm>
using namespace std;
const int L=50001;
const double eps=0.0;
struct point
{
    double x,y,r;
    int va,id;
}cir[L];
struct node
{
    double x;
    int id,flag;
}lr[L*2];
set<int> seg;
typedef set<int>::iterator it;
bool cmp(point a,point b)
{
    return a.y<b.y || a.y==b.y && a.x<b.x;
}
bool cmp1(node a,node b)
{
    return a.x<b.x;
}
bool pr(point a,point b)
{
    return a.x==b.x && a.y==b.y;
}
double dis(point a,point b)
{
    a.x-=b.x;a.y-=b.y;
    return sqrt(a.x*a.x+a.y*a.y);
}
int check(point a,point b)
{
    if(pr(a,b))
        if(a.r>b.r-eps) return 1;
    if(a.r<b.r) return 0;
    double d=dis(a,b);
    if(d<a.r-b.r-eps) return 1;
    return 0;
}
void slove(int n)
{
    int i,j;
    for(i=0;i<n;i++)
    {
        if(lr[i].flag==0)
        {
            it idx=seg.insert(lr[i].id).first;
            it be=idx,ed=idx;
            if(seg.size()>=2)
            {
                if(be==seg.begin())
                {
                    j=*++be;
                   if(check(cir[j],cir[lr[i].id]))
                        {   
                            cir[lr[i].id].va=1;
                            seg.erase(lr[i].id);
                        }
                }
                else if(++ed==seg.end())
                {
                    j=*--be;
                    if(check(cir[j],cir[lr[i].id]))
                        {
                            cir[lr[i].id].va=1;
                            seg.erase(lr[i].id);
                        }
                }
                else
                {
                    ed--;
                    if((check(cir[*--be],cir[lr[i].id])) || (check(cir[*++ed],cir[lr[i].id])))
                        {
                            cir[lr[i].id].va=1;
                            seg.erase(lr[i].id);
                        }
                }
            }
        }
        else
            seg.erase(lr[i].id);
    }
}
int ans[L];
int main()
{
    int i,j,n,m;
    while(~scanf("%d",&n))
    {
        j=0;
        for(i=0;i<n;i++)
        {
            scanf("%lf %lf %lf",&cir[i].r,&cir[i].x,&cir[i].y);
            cir[i].id=i;
            cir[i].va=0;
        }
        sort(cir,cir+n,cmp);
        for(i=0;i<n;i++)
        {
            lr[j].id=i;lr[j].x=cir[i].x-cir[i].r;lr[j++].flag=0;
            lr[j].id=i;lr[j].x=cir[i].x+cir[i].r;lr[j++].flag=1;
        }
        sort(lr,lr+j,cmp1);
        slove(j);
        j=0;
        for(i=0;i<n;i++)
            if(cir[i].va==0) ans[j++]=cir[i].id+1;
        printf("%d\n",j);
        sort(ans,ans+j);
        for(i=0;i<j;i++)
            printf("%d%c",ans[i],i==j-1?'\n':' ');
    }
    return 0;
}