Description
在平面直角坐标系中给定N个圆。已知这些圆两两没有交点，即两圆的关系只存在相离和包含。求这些圆的异或面积并。异或面积并为：当一片区域在奇数个圆内则计算其面积，当一片区域在偶数个圆内则不考虑。
Input
第一行包含一个正整数N，代表圆的个数。接下来N行，每行3个非负整数x,y,r，表示一个圆心在(x,y)，半径为r的圆。保证|x|,|y|,≤10^8，r>0,N<=200000
Output
仅一行一个整数，表示所有圆的异或面积并除以圆周率Pi的结果。
Sample Input
2
0 0 1
0 0 2
Sample Output
3
由于圆只存在包含与相离的关系，所以我们可以用扫描线，把圆看成一对对括号
然后判断一段区域的奇偶性，统计答案即可
至于对括号的维护，我们可以用set维护其相对位置即可
/*program from Wolfycz*/
#include<set>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define Fi first
#define Se second
#define inf 0x7f7f7f7f
#define sqr(x) ((x)*(x))
using namespace std;
typedef long long ll;
typedef long double ld;
typedef unsigned int ui;
typedef unsigned long long ull;
inline char gc(){
static char buf[1000000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++;
}
inline int frd(){
int x=0,f=1; char ch=gc();
for (;ch<'0'||ch>'9';ch=gc()) if (ch=='-') f=-1;
for (;ch>='0'&&ch<='9';ch=gc()) x=(x<<3)+(x<<1)+ch-'0';
return x*f;
}
inline int read(){
int x=0,f=1; char ch=getchar();
for (;ch<'0'||ch>'9';ch=getchar()) if (ch=='-') f=-1;
for (;ch>='0'&&ch<='9';ch=getchar()) x=(x<<3)+(x<<1)+ch-'0';
return x*f;
}
inline void print(int x){
if (x<0) putchar('-'),x=-x;
if (x>9) print(x/10);
putchar(x%10+'0');
}
const int N=2e5;
const double eps=1e-8;
double A[N+10],B[N+10],R[N+10],Nx;
struct node{
int type,ID;
node(int _type=0,int _ID=0){type=_type,ID=_ID;}
void insert(int _type,int _ID){type=_type,ID=_ID;}
double T(){return B[ID]+(double)type*(sqrt(sqr(R[ID])-sqr(Nx-A[ID]))+eps);}
};
bool operator <(node a,node b){return a.T()-b.T()<-eps;}
struct S1{
int type,v,ID;
void insert(int _type,int _v,int _ID){type=_type,v=_v,ID=_ID;}
bool operator <(const S1 &tis)const{return v<tis.v;}
}opr[(N<<1)+10];//operator
bool can[N+10];
set<node>ST;
int main(){
int n=read(),cnt=0; ll Ans=0;
for (int i=1;i<=n;i++){
scanf("%lf%lf%lf",A+i,B+i,R+i);
opr[++cnt].insert( 1,A[i]-R[i],i);
opr[++cnt].insert(-1,A[i]+R[i],i);
}
sort(opr+1,opr+1+cnt);
for (int i=1;i<=cnt;i++){
Nx=opr[i].v;
if (opr[i].type>0){
set<node>::iterator it=ST.insert(node(1,opr[i].ID)).Fi;
if (it==ST.begin()) can[opr[i].ID]=1;
else{
it--;
can[opr[i].ID]=can[(*it).ID]^((*it).type<0);
}
Ans+=(can[opr[i].ID]?1:-1)*sqr(R[opr[i].ID]);
ST.insert(node(-1,opr[i].ID));
}else{
ST.erase(node( 1,opr[i].ID));
ST.erase(node(-1,opr[i].ID));
}
}
printf("%lld\n",Ans);
return 0;
}