BZOJ4561: [JLoi2016]圆的异或并

来源:互联网 发布:软件出现乱码怎么办 编辑:程序博客网 时间:2024/05/29 02:31

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

如果考试的时候卷面像我这么标注,我怎么能看错题呢?!

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

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

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

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

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


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

#include<iostream>#include<cstdio>#include<set>#include<algorithm>#include<cmath>#define N 200010using 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);}

0 0
原创粉丝点击