poj2932:Coneology（扫描线）

传送门

题意：
给n个两两不相交的圆，求不被其他圆所包含的圆。(n$\le$40000 )

题解：扫描线
一个圆只能被他y左边的前驱后继所包含（当然圆的左右端点要包含这个圆）。直接套set。

#include<cstdio>
#include<algorithm>
#include<set>
using namespace std;
const int Maxn=4e4+50;
int n,tot,ans[Maxn],cnt;
typedef pair<double,int> pii;
typedef set<pii>::iterator it;
set<pii>S;
struct cir{
    double x,y,r;
}c[Maxn];
struct node{
    double pos;
    int val,id;
    node(int pos=0,int val=0,int id=0):pos(pos),val(val),id(id){}
    friend inline bool operator <(const node &a,const node &b){
        return a.pos<b.pos||(a.pos==b.pos&&c[a.id].r>c[b.id].r);
    }
}p[Maxn*2];
inline bool in(int x,int y){
    return c[x].r>c[y].r&&((c[y].y-c[x].y)*(c[y].y-c[x].y)+(c[y].x-c[x].x)*(c[y].x-c[x].x)<=(c[x].r-c[y].r)*(c[x].r-c[y].r));
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%lf%lf%lf",&c[i].r,&c[i].x,&c[i].y);
        p[++tot]=node(c[i].x-c[i].r,1,i);p[++tot]=node(c[i].x+c[i].r,-1,i);
    }
    sort(p+1,p+tot+1);
    for(int i=1;i<=tot;i++){
        if(p[i].val==-1){
            if(S.count(make_pair(c[p[i].id].y,p[i].id)))S.erase(S.find(make_pair(c[p[i].id].y,p[i].id)));
            continue;
        }
        it pre=S.lower_bound(make_pair(c[p[i].id].y,p[i].id));
        if((pre!=S.end()&&in(pre->second,p[i].id))||(pre!=S.begin()&&in((--pre)->second,p[i].id)))continue;
        S.insert(make_pair(c[p[i].id].y,p[i].id));ans[++cnt]=p[i].id;
    }
    printf("%d\n",cnt);sort(ans+1,ans+cnt+1);
    for(int i=1;i<=cnt;i++)printf("%d ",ans[i]);
}