求函数

题目链接：求函数

样例
输入

2 3
1 1
1 0
1 2 114514 1919810
2 1 2
2 1 1

输出

2148838
2

解题思路

维护一个线段树，前n项积，这里不能用树状数组，因为是单点更新所以这里不能用树状数组，emmm或许也可以我太菜了；维护的公式看代码就可以了

代码

#include<bits/stdc++.h>
using namespace std;
int n,m;
#define ll long long
const ll mod=1e9+7;
ll k[200020],b[200020];
struct node{
	ll k,b;
}tre[1000000];
ll ans=0;
void build(int x,int l,int r)
{
	if(l==r){
		tre[x].k=k[l];
		tre[x].b=b[l];
		return ;
	}
	int mid=(l+r)>>1;
	build(x<<1,l,mid);
	build(x<<1|1,mid+1,r);
	tre[x].b=(tre[x<<1].b*tre[x<<1|1].k%mod+tre[x<<1|1].b)%mod;//区间b值
	tre[x].k=tre[x<<1|1].k*tre[x<<1].k%mod;//区间k值
}
void up(int op,int l,int r,int x,int k,int b)
{
	if(l==r&&r==x)
	{
		tre[op].k=k;
		tre[op].b=b;
		return ;
	}
	int mid=(l+r)>>1;
	if(mid>=x) up(op<<1,l,mid,x,k,b);
	else up(op<<1|1,mid+1,r,x,k,b);
	tre[op].b=(tre[op<<1].b*tre[op<<1|1].k%mod+tre[op<<1|1].b)%mod;//更新b值
	tre[op].k=tre[op<<1|1].k*tre[op<<1].k%mod;//更新k值
	return;
}
void query(int x,int l,int r,int ql,int qr)
{
	if(ql<=l&&r<=qr)
	{
		ans=(tre[x].b+tre[x].k*ans%mod)%mod;
		return ;
	}
	int mid=(l+r)>>1;
	if(mid>=ql) query(x<<1,l,mid,ql,qr);
	if(qr>mid) query(x<<1|1,mid+1,r,ql,qr);
	return ;
}
int main()
{
	ios::sync_with_stdio(0);
	cin>>n>>m;
	for(int i=1;i<=n;i++)cin>>k[i];
	for(int i=1;i<=n;i++)cin>>b[i];
	build(1,1,n);
	
	int op,pos,K,B,l,r;
	while(m--)
	{
		cin>>op;
		if(op&1){
			cin>>pos>>K>>B;
			up(1,1,n,pos,K,B);
		}
		else{
			cin>>l>>r;
			ans=1;
			query(1,1,n,l,r);
			cout<<ans<<"\n";
		}
	}
	return 0;
}