题目描述
考虑维护一个这样的问题:
(1) 给出一个数组A,标号为1~n
(2) 修改数组中的一个位置。
(3) 询问区间[l,r]中所有子集的位运算and之和mod(10 9+7)。
位运算and即为“pascal中的and”和“C/C++中的&”
我们定义集合S={ l , l+1 , ... , r-1 , r}
若集合T,T ∩ S = T,则称T为S的子集
设f(T)=A T1 and A T2 and ... and A Tk (设k为T集大小,若k=0则f(T)=0)
所有子集的位运算and之和即为∑f(T)
那么,现在问题来了。
(1) 给出一个数组A,标号为1~n
(2) 修改数组中的一个位置。
(3) 询问区间[l,r]中所有子集的位运算and之和mod(10 9+7)。
位运算and即为“pascal中的and”和“C/C++中的&”
我们定义集合S={ l , l+1 , ... , r-1 , r}
若集合T,T ∩ S = T,则称T为S的子集
设f(T)=A T1 and A T2 and ... and A Tk (设k为T集大小,若k=0则f(T)=0)
所有子集的位运算and之和即为∑f(T)
那么,现在问题来了。
输入描述:
第一行,一个正整数N 第二行,N个非负整数,为数组A 第三行,一个正整数M,为操作次数 接下来M行格式如下 修改操作: 1 x y,将Ax修改为y 询问操作: 2 l r,区间[l,r]中所有子集的位运算and之和 mod(109+7)
输出描述:
对于每次询问输出一行,为该次询问的答案mod(109+7)。
long long 请使用lld
示例1
输入
3 1 2 3 6 2 1 3 1 1 2 2 1 3 2 2 3 1 2 5 2 1 3
输出
9 15 7 13
说明
第一次询问: Answer =1+2+3+(1 and 2)+(1 and 3)+(2 and 3)+(1 and 2 and 3) =1+2+3+0+1+2+0 =9 第二次询问: Answer =2+2+3+(2 and 2)+(2 and 3)+(2 and 3)+(2 and 2 and 3) =2+2+3+2+2+2+2 =15 第三次询问: Answer =2+3+(2 and 3) =2+3+2 =7 第四次询问: Answer =2+5+3+(2 and 5)+(2 and 3)+(3 and 5)+(2 and 5 and 3) =2+5+3+0+2+1+0 =13
备注:
题解:对于区间与,很多时候考虑枚举每一位上的贡献,对于该题
我们考虑将每个数拆分,对于每一位建一个树状数组(不会超过30)
然后我们按照树状数组进行修改即可,查询的话,对于一个区间上该位为1对该区间所有子集的贡献为
(2^x-1)*2^p (x:该区间这一位为一的数的个数,p:第p位) PS: 不知道为啥卡long long 。。。
#include<stdio.h>
#include<algorithm>
using namespace std;
#define maxn 100005
#define mod 1000000007
#define lowbit(x) x&(-x)
typedef long long ll;
int n,m,a[maxn],p[maxn],c[35][maxn];
void add(int k,int x,int val)
{
while(x<=n)
{
c[k][x]+=val;
x+=lowbit(x);
}
}
int sum(int k,int x)
{
int res=0;
while(x)
{
res+=c[k][x];
x-=lowbit(x);
}
return res;
}
int main(void)
{
p[0]=1;
int i,j,t,x,y;
scanf("%d",&n);
for(i=1;i<=n;i++)
p[i]=1ll*2*p[i-1]%mod;
for(i=1;i<=n;i++)
{
scanf("%d",&a[i]);
for(j=0;j<=30;j++)
if(a[i]&(1<<j))
add(j,i,1);
}
scanf("%d",&m);
while(m--)
{
scanf("%d%d%d",&t,&x,&y);
if(t==1)
{
for(i=0;i<=30;i++)
if(a[x]&(1<<i))
add(i,x,-1);
a[x]=y;
for(i=0;i<=30;i++)
if(a[x]&(1<<i))
add(i,x,1);
}
else
{
ll ans=0;
for(i=0;i<=30;i++)
{
int tmp=sum(i,y)-sum(i,x-1);
ans=(ans+1ll*(p[tmp]-1)*p[i]%mod)%mod;
}
printf("%d\n",ans);
}
}
return 0;
}