Loj 2028 随机序列

Loj 2028 随机序列

• 连续的乘号会将序列分成若干个块,块与块之间用加减号连接:

[(a_1*a_2*...a_i)pm(a_{i+1}*a_{i+2}*...a_j)pm... ]

• 除去所有都是乘号的一种,对于任意一个序列,总有与之对应唯一的另一个序列,它们所有的加减号都相反.
• 把这两个序列的和相加,就只剩下了 (2*(a_1*a_2*...a_i)),其中 (i) 为第一个加减号之前的位置.
• 那么平均下来,相当于每个序列的贡献为 ((a_1*a_2*...a_i)).
• 考虑对每个 (i) ,前面填入了 (i-1) 个乘号,相邻的必须是加减号,后面 (n-i-1) 个符号随便填,共有 (2*3^{n-i-1}) 种.
• 那么用线段树对于 ([1,n-1]) 维护 (v_i=3^{n-i-1}*2*prod_{j=1}^i a_i) 的和,答案即为

[(prod_{i=1}^{n}a_i)+(sum_{i=1}^{n-1}v_i). ]

• 将位置 (i) 的数从 (a) 改为 (b) ,相当于对区间 ([i,n-1]) 乘了 (frac b a) .显然容易实现.
• 注意修改位置 (n) 上的数时不用区间修改.否则 (L>R) ,会炸.

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define mp make_pair
#define pii pair<int,int>
inline int read()
{
	int x=0;
	bool pos=1;
	char ch=getchar();
	for(;!isdigit(ch);ch=getchar())
		if(ch=='-')
			pos=0;
	for(;isdigit(ch);ch=getchar())
		x=x*10+ch-'0';
	return pos?x:-x;
}
const int P=1e9+7;
const int MAXN=1e5+10;
const int MAXV=1e4;
int n,m;
inline int add(int a,int b)
{
	return (a + b) % P;
}
inline int mul(int a,int b)
{
	return 1LL * a * b % P;
}
int fpow(int a,int b)
{
	int res=1;
	while(b)
		{
			if(b&1)
				res=mul(res,a);
			a=mul(a,a);
			b>>=1;
		}
	return res;
}
int inv[MAXV+10],Pow[MAXN];
int a[MAXN],prod[MAXN];
struct node{
	int l,r;
	int sum,tag;
	node()
		{
			tag=1;
		}
}Tree[MAXN<<2];
#define root Tree[o]
#define lson Tree[o<<1]
#define rson Tree[o<<1|1]
void pushup(int o)
{
	root.sum=add(lson.sum,rson.sum);
}
void BuildTree(int o,int l,int r)
{
	root.l=l,root.r=r;
	if(l==r)
		{
			//fpow(3,n-l-1) Pow[n-l-1]
			int tmp=mul(Pow[n-l-1],2);
			tmp=mul(tmp,prod[l]);
			root.sum=tmp;
			return;
		}
	int mid=(l+r)>>1;
	BuildTree(o<<1,l,mid);
	BuildTree(o<<1|1,mid+1,r);
	pushup(o);
}
void Modifiy(int o,int c)
{
	root.sum=mul(root.sum,c);
	root.tag=mul(root.tag,c);
}
void pushdown(int o)
{
	if(root.tag!=1)
		{
			Modifiy(o<<1,root.tag);
			Modifiy(o<<1|1,root.tag);
			root.tag=1;
		}
}
void update(int o,int L,int R,int c)
{
	int l=root.l,r=root.r;
	if(l>R || r<l)
		return;
	if(L<=l && r<=R)
		{
			Modifiy(o,c);
			return;
		}
	pushdown(o);
	int mid=(l+r)>>1;
	if(L<=mid)
		update(o<<1,L,R,c);
	if(R>mid)
		update(o<<1|1,L,R,c);
	pushup(o);
}
void init()
{
	inv[1]=1;
	for(int i=2;i<=MAXV;++i)
		inv[i]=mul(P-P/i,inv[P%i]);
	Pow[0]=1;
	for(int i=1;i<=n;++i)
		Pow[i]=mul(Pow[i-1],3);
}
int main()
{
	n=read(),m=read();
	init();
	prod[0]=1;
	for(int i=1;i<=n;++i)
		prod[i]=mul(prod[i-1],a[i]=read());
	BuildTree(1,1,n-1);
	while(m--)
		{
			int i=read(),v=read();
			int c=mul(v,inv[a[i]]);
			a[i]=v;
			prod[n]=mul(prod[n],c);
			if(i<=n-1)
				update(1,i,n-1,c);
			int ans=add(Tree[1].sum,prod[n]);
			printf("%d
",ans);
		}
	return 0;
}