分析:
扇形面积
=
p
2
m
×
π
r
2
=\frac{p}{2m}\times \pi r^2
=2mp×πr2
p
p
p为扇形所占的份数 最后答案乘上
2
m
p
\frac{2m}{p}
p2m 这样就化掉了
要找的就是第
k
k
k大半径 用
a
r
a_r
ar记录半径为
r
r
r的数量 从
[
−
m
,
m
]
[-m,m]
[−m,m]扫过去 遇到起始点就将
a
r
+
1
a_r+1
ar+1 终点就
a
r
−
1
a_r-1
ar−1
然后线段树维护
CODE:
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long ll;
const int N=4e5+5;
struct node{
int st,val,r;
}a[N<<1];
int n,m,qwq,tot,head[N],next[N],tree[N<<1];
void up(int x){tree[x]=tree[x<<1]+tree[x<<1|1];}
void update(int x,int l,int r,int p,int k)
{
if(l==r)
{
tree[x]+=k;
return;
}
int mid=(l+r)>>1;
if(p<=mid) update(x<<1,l,mid,p,k);
else update(x<<1|1,mid+1,r,p,k);
up(x);
}
ll query(int x,int l,int r,int p)
{
if(l==r)
{
if(tree[x]>=p) return l;
else return 0;
}
int mid=(l+r)>>1;
if(tree[x<<1|1]>=p) return query(x<<1|1,mid+1,r,p);
else return query(x<<1,l,mid,p-tree[x<<1|1]);
}
int main(){
scanf("%d%d%d",&n,&m,&qwq);
for(int i=1,R,S,T;i<=n;i++)
{
scanf("%d%d%d",&R,&S,&T);
if(S==-m||S==m)
{
a[++tot]=(node){-m,1,R};
a[++tot]=(node){T,-1,R};
}
else
{
if(T==-m) T=m;
if(S<=T)
{
a[++tot]=(node){S,1,R};
a[++tot]=(node){T,-1,R};
}
else
{
a[++tot]=(node){S,1,R};
a[++tot]=(node){m,-1,R};
a[++tot]=(node){-m,1,R};
a[++tot]=(node){T,-1,R};
}
}
}
for(int i=1;i<=tot;i++)
{
int k=a[i].st+m;
next[i]=head[k];
head[k]=i;
}
ll ans=0;
for(int i=-m;i<m;i++)
{
for(int k=head[i+m];k;k=next[k])
update(1,1,N,a[k].r,a[k].val);
ll res=query(1,1,N,qwq);
ans+=res*res;
}
printf("%lld",ans);
return 0;
}