[BZOJ 4418][Shoi2013]扇形面积并(树状数组+二分)
Posted Zars19
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[BZOJ 4418][Shoi2013]扇形面积并(树状数组+二分)相关的知识,希望对你有一定的参考价值。
Description
给定N个同心的扇形,求有多少面积,被至少K个扇形所覆盖。
Solution
打开发现是计算几何还以为是看错题号了QwQ
其实就是遇到一条开始的边+1,遇到一条结束的边-1,a1>a2的话就再加上一个完整的圆
计算每一部分只需要找到第cnt-k+1大的半径,在树状数组上二分就好了
#include<iostream> #include<cstdio> #include<cstring> #include<cstdlib> #include<algorithm> #define MAXN 100005 using namespace std; typedef long long LL; int n,m,k,tot=0,c[MAXN]; int read() { int x=0,f=1;char c=getchar(); while(c<‘0‘||c>‘9‘){ if(c==‘-‘)f=-1;c=getchar(); } while(c>=‘0‘&&c<=‘9‘){ x=x*10+c-‘0‘;c=getchar(); } return x*f; } struct Node { int r,a,f; Node(int r=0,int a=0,int f=0):r(r),a(a),f(f){} bool operator < (const Node& x) const {return a<x.a;} }data[400005]; int lowbit(int x){return x&-x;} void add(int pos,int x) { while(pos<=MAXN) { c[pos]+=x; pos+=lowbit(pos); } } int solve(int sum) { int x=0,now=0; for(int i=(1<<20);i;i>>=1) if(x+i<=MAXN&&now+c[x+i]<sum) now+=c[x+i],x+=i; return x+1; } int main() { n=read(),m=read(),k=read(); for(int i=1;i<=n;i++) { int r=read(),a1=read(),a2=read(); data[++tot]=Node(r,a1,1); data[++tot]=Node(r,a2,-1); if(a1>a2) { data[++tot]=Node(r,-m,1); data[++tot]=Node(r,m,-1); } } sort(data+1,data+1+tot); LL ans=0;int cnt=0; for(int i=1;i<tot;i++) { add(data[i].r,data[i].f); cnt+=data[i].f; if(cnt-k+1<=0)continue; int t=solve(cnt-k+1); ans+=1LL*t*t*(data[i+1].a-data[i].a); } printf("%lld\n",ans); return 0; }
以上是关于[BZOJ 4418][Shoi2013]扇形面积并(树状数组+二分)的主要内容,如果未能解决你的问题,请参考以下文章
bzoj4418[Shoi2013]扇形面积并 扫描线+线段树
[BZOJ 4418][Shoi2013]扇形面积并(树状数组+二分)
4418: [Shoi2013]扇形面积并|二分答案|树状数组