题目:
题意:
给出一个二进制串,我们需要维护对其局部的升序或降序排序,并要支持查询局部的二进制数在十进制下的数是多少
分析:
对于排序,因为是二进制数所以只含有
0
0
0或
1
1
1,降序就是把所有
1
1
1移到最前面,升序反之
查询操作我们根据子树的
0
、
1
0、1
0、1来计算当前节点的权值,然后进行维护
代码:
#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<queue>
#include<algorithm>
#include<cmath>
#include<ctime>
#define LL long long
#define mo 1000000007
using namespace std;
inline LL read() {
LL d=0,f=1;char s=getchar();
while(s<'0'||s>'9'){if(s=='-')f=-1; s=getchar();}
while(s>='0'&&s<='9'){d=d*10+s-'0'; s=getchar();}
return d*f;
}
LL ksm(LL x,LL y)
{
LL s=1;
while(y)
{
if(y&1) (s*=x)%=mo;
(x*=x)%=mo;y>>=1;
}
return s%mo;
}
LL n;
char c[200005];
struct Tree{
LL w[800005],s[800005],lazy[800005];
void build(LL k,LL L,LL R)
{
if(L==R) {s[k]=c[L-1]-'0';w[k]=((c[L-1]-'0')*ksm(2,n-L))%mo;return;}
LL mid=(L+R)>>1;
build(k*2,L,mid);build(k*2+1,mid+1,R);
w[k]=(w[k*2]+w[k*2+1])%mo;
s[k]=s[k*2]+s[k*2+1];
return;
}
void pushdown(LL k,LL L,LL R)
{
LL mid=(L+R)>>1;
if(!lazy[k]) return;
if(lazy[k]==1)
{
w[k*2]=(ksm(2,n-L+1)-1-(ksm(2,n-mid)-1)+mo)%mo;
w[k*2+1]=(ksm(2,n-mid)-1-(ksm(2,n-R)-1)+mo)%mo;
s[k*2]=(mid-L+1);s[k*2+1]=(R-mid);
}
else w[k*2]=w[k*2+1]=s[k*2]=s[k*2+1]=0;
lazy[k*2]=lazy[k*2+1]=lazy[k];
lazy[k]=0;
return;
}
LL qq(LL k,LL L,LL R,LL l,LL r)
{
if(l<=L&&R<=r) return s[k];
pushdown(k,L,R);
LL mid=(L+R)>>1;
if(mid<l) return qq(k*2+1,mid+1,R,l,r);
else if(mid>=r) return qq(k*2,L,mid,l,r);
else return qq(k*2,L,mid,l,r)+qq(k*2+1,mid+1,R,l,r);
}
void change(LL k,LL L,LL R,LL l,LL r,LL x)
{
if(l<=L&&R<=r)
{
if(!x) lazy[k]=2,w[k]=0,s[k]=0;
else
{
lazy[k]=1;
s[k]=R-L+1;
w[k]=(ksm(2,n-L+1)-1-(ksm(2,n-R)-1)+mo)%mo;
}
return;
}
pushdown(k,L,R);
LL mid=(L+R)>>1;
if(mid<l) change(k*2+1,mid+1,R,l,r,x);
else if(mid>=r) change(k*2,L,mid,l,r,x);
else change(k*2,L,mid,l,r,x),change(k*2+1,mid+1,R,l,r,x);
w[k]=(w[k*2]+w[k*2+1])%mo;
s[k]=s[k*2]+s[k*2+1];
return;
}
LL query(LL k,LL L,LL R,LL l,LL r)
{
if(l<=L&&R<=r) return w[k]%mo;
pushdown(k,L,R);
LL mid=(L+R)>>1;
if(mid<l) return query(k*2+1,mid+1,R,l,r)%mo;
else if(mid>=r) return query(k*2,L,mid,l,r)%mo;
else return (query(k*2,L,mid,l,r)+query(k*2+1,mid+1,R,l,r))%mo;
}
}t;
int main()
{
scanf("%s",c);
n=strlen(c);
t.build(1,1,n);
LL q=read();
for(LL i=1;i<=q;i++)
{
LL op=read(),l=read(),r=read();
if(op==2)
{
LL s=t.qq(1,1,n,l,r);
if(!s||s==r-l+1) continue;
LL l1=l,r1=r-s,l2=r-s+1,r2=r;
t.change(1,1,n,l1,r1,0);
t.change(1,1,n,l2,r2,1);
}
if(op==1)
{
LL s=t.qq(1,1,n,l,r);
if(!s||s==r-l+1) continue;
LL l1=l,r1=l+s-1,l2=l+s,r2=r;
t.change(1,1,n,l1,r1,1);
t.change(1,1,n,l2,r2,0);
}
if(op==3)
printf("%lld\n",t.query(1,1,n,l,r)*ksm(ksm(2,n-r)%mo,mo-2)%mo);
}
return 0;
}