[bzoj 4561--JLoi2016]圆的异或并

在平面直角坐标系中给定N个圆。已知这些圆两两没有交点，即两圆的关系只存在相离和包含。求这些圆的异或面积并。异或面积并为：当一片区域在奇数个圆内则计算其面积，当一片区域在偶数个圆内则不考虑。

看到这个求面积就应该想到扫描线啊，还是太菜了。
我们发现只要知道每个圆覆盖它的圆的个数的奇偶性就可以计算答案了，并且题目中有个重要的性质：要不包含，要不相离，所以我们如果知道离每个圆的最近的圆弧是上还是下，就可以递推计算了。
解决这个可以用扫描线+set来做，所以这题就做完了。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#include<set>
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0' || ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0' && ch<='9')x=x*10+ch-'0',ch=getchar();
    return x*f;
}
inline void write(int x)
{
    if(x<0)putchar('-'),x=-x;
    if(x>9)write(x/10);
    putchar(x%10+'0');
}
inline void pr1(int x){write(x),putchar(' ');}
inline void pr2(int x){write(x),puts("");}
struct circle
{
    int x,y,r;
}c[200010];
struct node
{
    int p,x,id;
}l[400010],uf;
inline long long pf(int x){return 1LL*x*x;}
int P,ans[200010];
bool cmp(node a,node b){return a.x<b.x;}
bool operator <(node a,node b)
{
    double da=c[a.p].y+a.id*sqrt(pf(c[a.p].r)-pf(P-c[a.p].x)),db=c[b.p].y+b.id*sqrt(pf(c[b.p].r)-pf(P-c[b.p].x));
    if(da==db)return a.id<b.id;
    else return da<db;
}
set<node>s;
typedef set<node>::iterator It;
int main()
{
    //freopen("a.in","r",stdin);
    //freopen("a.out","w",stdout);
    int n=read(),top=0;
    for(int i=1;i<=n;i++)
    {
        c[i].x=read(),c[i].y=read(),c[i].r=read();
        l[++top].p=i,l[top].x=c[i].x-c[i].r,l[top].id=1;
        l[++top].p=i,l[top].x=c[i].x+c[i].r,l[top].id=-1;
    }
    sort(l+1,l+top+1,cmp);
    uf.id=1,uf.x=0;
    for(int i=1;i<=top;i++)
    {
        P=l[i].x;uf.p=l[i].p;
        if(l[i].id==1)
        {
            It wy=s.upper_bound(uf);
            if(wy==s.end())ans[uf.p]=1;
            else
            {
                if((*wy).id==1)ans[uf.p]=-ans[(*wy).p];
                else ans[uf.p]=ans[(*wy).p];
            }
            uf.id=1,s.insert(uf);uf.id=-1,s.insert(uf);
        }
        else
        {
            uf.id=1,s.erase(uf);uf.id=-1,s.erase(uf);
        }
    }
    long long sum=0;
    for(int i=1;i<=n;i++)sum+=1LL*ans[i]*c[i].r*c[i].r;
    printf("%lld\n",sum);
    return 0;
}