bzoj4642: 泡泡

Description

"OI真的像是一条奇趣横生的路啊,也许它是绕过了高考的大山,也许确实有通往大学的捷径。但我,真的,真的只在
乎那路上美丽的泡泡。"
--TB

TB喜欢所有自然的事物。比如说松爷的仙人掌，Picks的多项式导论，当然，还有OI路上美丽的泡泡。这些泡泡可
以视作某一平面上的一些圆。由于泡泡的特殊性质，当两个泡泡在这一平面上相切的时候，TB认为这对泡泡是自然
的，然而如果它们相交或者包含的话，泡泡就会破裂而无法继续存在（即不会存在相交或包含的情况）。TB想知道
有多少对泡泡是自然的。

Input

输入文件的第一行包含一个正整数n，表示泡泡的个数。
接下来n行，每行三个整数x,y,r，表示一个泡泡的圆心和半径。
对于所有数据，|x|，|y|，r<=10^9, n<=500000
数据保证所有的泡泡都是存在的，既不会出现相交或者包含的关系。

Output

一行，表示有多少对自然的泡泡。
首先旋转坐标系避免一些特判

用扫描线处理，由于圆互不包含/相切，扫描线上用平衡树维护线上圆的y坐标

当线上新增加一个圆时，检查y坐标在其上/下的圆是否与其相切

当线上删除一个圆时，检查y坐标在其上、下的两圆是否相切

以上两个判定可能重复，最后要去重

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<set>
typedef long double ld;
const ld _0=1e-8,_s=std::sin(0.43415),_c=std::cos(0.43415);
int ans=0;
struct cir{
ld x,y,r;
}cs[500010];
struct event{
ld x;
int id,type;
}es[1000010];
std::pair<int,int>as[1500050];
int ap=0;
bool operator<(event a,event b){
return a.x<b.x;
}
struct cmp{bool operator()(int x,int y){return cs[x].y<cs[y].y;}};
std::set<int,cmp>line;
void chk(int a,int b){
ld x=cs[a].x-cs[b].x,y=cs[a].y-cs[b].y;
if(std::fabs(std::sqrt(x*x+y*y)-cs[a].r-cs[b].r)<_0){
if(a>b){int c=a;a=b;b=c;}
as[ap++]=std::make_pair(a,b);
}
}
int main(){
int n=_int();
for(int i=0;i<n;i++){
int x,y,r;
scanf("%d%d%d",&x,&y,&r);
cs[i].x=_c*x+_s*y;
cs[i].y=-_s*x+_c*y;
cs[i].r=r;
es[i<<1]=(event){cs[i].x-cs[i].r,i,1};
es[i<<1^1]=(event){cs[i].x+cs[i].r,i,0};
}
n<<=1;
std::sort(es,es+n);
for(int i=0;i<n;i++){
event w=es[i];
if(w.type){
std::set<int,cmp>::iterator it=line.upper_bound(w.id);
if(it!=line.end())chk(*it,w.id);
if(it!=line.begin())--it,chk(*it,w.id);
line.insert(w.id);
}else{
std::set<int,cmp>::iterator it=line.find(w.id);
if(it!=line.begin()){
int a=*--it;++it;++it;
if(it!=line.end())chk(*it,a);
--it;
}
line.erase(it);
}
}
std::sort(as,as+ap);
if(ap)ans=1;
for(int i=1;i<ap;i++)if(as[i]!=as[i-1])++ans;
printf("%d\n",ans);
return 0;
}