BZOJ3110: [Zjoi2013]K大数查询

喜闻乐见的简单树套树= =第一维按权值建树状数组,第二维按下标建动态开点线段树,修改相当于第二维区间加,查询在树状数组上二分,比一般的线段树还短= =可惜并不能跑过整体二分= =另外bzoj上的数据有负数= =其他树套树方法也是可以的爱怎么套怎么套= =

#include<cstdio>
#define J (i+j>>1)
#define I (J+1)
typedef unsigned ll;
const int N=1e5+5;
ll n,m,q,i,j,k,s,t,u,v;
struct node{
ll s,t;
node*i,*j;
}e[N*320];
node*a=e,*r[N];
void vary(node*&o,int s,int t,int i=1,int j=n){
if(!o)o=a++;
o->s+=t-s+1;
if(s==i&&t==j)++o->t;
else if(t<I)
vary(o->i,s,t,i,J);
else if(s>J)
vary(o->j,s,t,I,j);
else{
vary(o->i,s,J,i,J);
vary(o->j,I,t,I,j);
}
}
void ask(node*o,int s,int t,int i=1,int j=n){
if(!o)return;
if(s==i&&t==j)u+=o->s;
else{
u+=o->t*(t-s+1);
if(t<I)
ask(o->i,s,t,i,J);
else if(s>J)
ask(o->j,s,t,I,j);
else{
ask(o->i,s,J,i,J);
ask(o->j,I,t,I,j);
}
}
}
int main(){
scanf("%d%d",&n,&q);
m=n<<1^1;
while(q--){
scanf("%d%d%d%d",&j,&s,&t,&k);
if(j==2){
for(i=65536,v=j=0;i;i>>=1)
if(j+i<=m){
u=0,ask(r[j+i],s,t);
if(v+u<k)
j+=i,v+=u;
}
printf("%d\n",n-j);
}else
for(i=n-k+1;i<=m;i+=i&-i)
vary(r[i],s,t);
}
}
上一篇:perl模块安装


下一篇:如何优雅地使用Markdown (Sublime 3 + MarkdownEditing+OmniMarkupPreviewer)