bzoj 4597: [Shoi2016]随机序列

来源：互联网发布：windows7旗舰版优化编辑：程序博客网时间：2024/05/01 09:25

对于一个表达式序列 ****+??????? 一定有一个表达式序列 ****-??????? 与它相加为0，所以对答案有贡献的是第一个非*之前的序列乘上这个序列出现的次数。对于前i个数组成的序列，出现次数=2*3^(n-i-1)，前n个数组成的序列出现次数为1.所以用线段树维护这个式子的值就好了，具体实现看代码

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
 
#define ll long long
#define inf 1e9
#define eps 1e-8
#define md 1000000007
#define N 100010
using namespace std;
ll ji[4*N],sum[4*N],k[N],a[N];
void update(int i)
{
ji[i]=ji[i<<1]*ji[i<<1|1]%md;
sum[i]=(sum[i<<1]+sum[i<<1|1]*ji[i<<1])%md;
}
 
void build(int i,int l,int r)
{
if (l==r) { ji[i]=a[l]; sum[i]=ji[i]*k[l]%md; return;}
int mid=(l+r)>>1;
build(i<<1,l,mid); build(i<<1|1,mid+1,r);
update(i);
}
 
void modify(int i,int l,int r,int x,int d)
{
if (l==r) { ji[i]=d; sum[i]=ji[i]*k[l]%md; return;}
int mid=(l+r)>>1;
if (x<=mid) modify(i<<1,l,mid,x,d);
  else modify(i<<1|1,mid+1,r,x,d);
update(i);
}
 
int main()
{
int n,Q;
scanf("%d%d",&n,&Q);
for (int i=1;i<=n;i++) scanf("%lld",&a[i]);
k[n]=1; k[n-1]=2; for (int i=n-2;i;i--) k[i]=k[i+1]*3%md;
build(1,1,n);
for (int i=1;i<=Q;i++)
{
int x,d;
scanf("%d%d",&x,&d);
modify(1,1,n,x,d);
printf("%lld\n",sum[1]);
}
return 0;
}

0 0