bzoj 4561: [JLoi2016]圆的异或并

Description

在平面直角坐标系中给定N个圆。已知这些圆两两没有交点，即两圆的关系只存在相离和包含。求这些圆的异或面

积并。异或面积并为：当一片区域在奇数个圆内则计算其面积，当一片区域在偶数个圆内则不考虑。

Solution

我们可以把异或当作容斥的一个过程....现在要确定每一个圆的系数.

由于不存在相交关系 , 圆与圆之间相对顺序是确定的 , 所以可以用 \(set\) 来维护相对顺序.

我们对 \(x\) 做扫描线来维护一个圆的插入和删除.

我们把一个圆拆成上下两个圆弧 , 新插入一个圆时 , 判断上面的第一个圆弧是上圆弧还是下圆弧.

容易发现 : 如果是上圆弧那么系数相反 , 如果是下圆弧则相同.

求出每个圆的系数 , 最后累加答案就行了.

注意一个细节 , 上下圆弧 \(y\) 相同时 , 要把下圆弧放在下面...

#include<bits/stdc++.h>
#define sqr(x) (1ll*(x)*(x))
using namespace std;
template<class T>void gi(T &x){
int f;char c;
for(f=1,c=getchar();c<'0'||c>'9';c=getchar())if(c=='-')f=-1;
for(x=0;c<='9'&&c>='0';c=getchar())x=x*10+(c&15);x*=f;
}
const int N=200010;
int n,f
,y
,r
,X,x
;
struct data{int id,k;}q[N*2];
inline bool operator <(data p,data q){
return x[p.id]+r[p.id]*p.k<x[q.id]+r[q.id]*q.k;
}
struct node{int id,k;};
inline bool operator <(node p,node q){
double y1=y[p.id]+p.k*sqrt(sqr(r[p.id])-sqr(X-x[p.id]));
double y2=y[q.id]+q.k*sqrt(sqr(r[q.id])-sqr(X-x[q.id]));
if(y1<y2)return 1;if(y1>y2)return 0;
return p.k<q.k;
}
set<node>S;
set<node>::iterator it;
int main(){
freopen("pp.in","r",stdin);
freopen("pp.out","w",stdout);
cin>>n;
int cnt=0;
for(int i=1;i<=n;i++){
gi(x[i]),gi(y[i]),gi(r[i]);
q[++cnt]=(data){i,-1},q[++cnt]=(data){i,1};
}
sort(q+1,q+cnt+1);
for(int i=1;i<=cnt;i++){
X=x[q[i].id]+q[i].k*r[q[i].id];
if(q[i].k==-1){
int x=q[i].id;
node t1=(node){x,1},t2=(node){x,-1};
it=S.lower_bound(t1);
if(it!=S.end()){
if(it->k>0)f[x]=-f[it->id];
else f[x]=f[it->id];
}
else f[x]=1;
S.insert(t1),S.insert(t2);
}
else S.erase((node){q[i].id,1}),S.erase((node){q[i].id,-1});
}
long long ans=0;
for(int i=1;i<=n;i++)ans+=1ll*f[i]*r[i]*r[i];
cout<<ans;
return 0;
}