BZOJ4561: [JLoi2016]圆的异或并

题目大意：给你一堆圆，只存在相离和包含的关系，求他们的异或面积并

如果考试的时候卷面像我这么标注，我怎么能看错题呢？！

用一个扫描线，从左到右扫一遍，因为只有这两种关系，所以扫描线上的这些圆相对位置是固定的

我们可以把每个圆想象成两个括号，括号之间不会相交，每当扫描线扫到一个新的圆时，相当于在某处添加一对括号，而当扫描线扫完一个圆时，相当于把一个中间不包含任何东西的括号删除

用set维护一下相对位置就行了

然后加完这个括号之后，假如他后面紧跟着一个左括号，说明这两个圆是并列的，即Fi=Fj

否则这个圆被那个圆直接包含，即Fi=-Fj

目前好像是BZOJ的最短代码，ljss和xuruifan写的都是treap...

#include<iostream>
#include<cstdio>
#include<set>
#include<algorithm>
#include<cmath>
#define N 200010
using namespace std;
long long T;
struct circle{long long x,y,r;};
struct ppp{long long num,x,o;}b[N<<1];
bool cmp(ppp x,ppp y){return x.x<y.x;}
circle c[N];
long long sqr(long long x){return x*x;}
bool operator <(ppp x,ppp y)
{
	double xx=c[x.num].y+x.o*sqrt(sqr(c[x.num].r)-sqr(T-c[x.num].x));
	double yy=c[y.num].y+y.o*sqrt(sqr(c[y.num].r)-sqr(T-c[y.num].x));
	if(xx!=yy) return xx<yy;
	return x.o<y.o;
}
set<ppp>S;
long long f[N];
int main()
{
	long long n;
	scanf("%lld",&n);
	long long i,j;
	for(i=1;i<=n;i++)
	{
		scanf("%lld%lld%lld",&c[i].x,&c[i].y,&c[i].r);
		b[2*i-1]=(ppp){i,c[i].x-c[i].r,1};
		b[2*i]=(ppp){i,c[i].x+c[i].r,-1};
	}
	sort(b+1,b+2*n+1,cmp);
	set<ppp>::iterator it;
	for(i=1;i<=2*n;i++)
	{
		T=b[i].x;
		if(b[i].o==1)
		{
			it=S.upper_bound((ppp){b[i].num,0,-1});
			if(it==S.end()) f[b[i].num]=1;
			else
			{
				if((*it).o==1) f[b[i].num]=-f[(*it).num];
				else f[b[i].num]=f[(*it).num];
			}
			S.insert((ppp){b[i].num,0,-1});
			S.insert((ppp){b[i].num,0,1});
		}
		else
		{
			S.erase((ppp){b[i].num,0,-1});
			S.erase((ppp){b[i].num,0,1});
		}
	}
	long long ans=0;
	for(i=1;i<=n;i++)
	ans+=sqr(c[i].r)*f[i];
	printf("%lld",ans);
}