修仙录 3.22

今天的题有点难呀
不过鲍丽芬非常良心有30！！！

---

jzoj 6078 魔法阵

https://jzoj.net/senior/#contest/show/2677/2

其实用线段树就好啦
注意如果在修改中出现负数，要处理一下

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int MAXN=1e5+5;
const int mod=1e9+7;

int n,Q,opt,x,y;
int A[MAXN],B[MAXN],C[MAXN];
long long ans=1;

long long Pow(long long a,int b){
	long long ans=1;
	while(b){
		if(b&1) ans=ans*a%mod;
		a=a*a%mod;
		b>>=1;
	}
	return ans;
}

struct line_tree{
	struct node{
		int B,C,mn,b,pb,pm,num,laz;
		long long prd;
	}tr[MAXN*4];
	void up(int p){
		int l=p<<1,r=p<<1|1;
		if(tr[l].b<tr[r].b) tr[p].b=tr[l].b,tr[p].pb=tr[l].pb;
		else tr[p].b=tr[r].b,tr[p].pb=tr[r].pb;
		if(tr[l].mn<tr[r].mn) tr[p].mn=tr[l].mn,tr[p].pm=tr[l].pm;
		else tr[p].mn=tr[r].mn,tr[p].pm=tr[r].pm;
		tr[p].C=tr[l].C,tr[p].num=tr[l].num+tr[r].num;
		tr[p].prd=tr[l].prd*tr[r].prd%mod;
	}
	void solve(int p,int l){
		if(tr[p].B>tr[p].C){
			tr[p].mn=tr[p].B-tr[p].C;
			tr[p].pm=tr[p].pb=l;
			tr[p].num=1,tr[p].b=tr[p].B,tr[p].prd=tr[p].C;
			ans=ans*tr[p].C%mod;
		}
		else{
			tr[p].mn=1e9;
			tr[p].pm=tr[p].pb=0;
			tr[p].num=0,tr[p].b=1e9,tr[p].prd=1;
			ans=ans*tr[p].B%mod;
		}
	}
	void build(int p,int l,int r){
		if(l==r){
			tr[p].B=B[l],tr[p].C=C[l];
			solve(p,l);
			return ;
		}
		int mid=l+r>>1;
		build(p<<1,l,mid);
		build(p<<1|1,mid+1,r);
		up(p);
	}
	void turn(int p,int v){
		tr[p].C=v;
		if(tr[p].b!=1e9) tr[p].mn=tr[p].b-v,tr[p].pm=tr[p].pb;
		else tr[p].mn=1e9,tr[p].pm=0;
		tr[p].prd=Pow(v,tr[p].num);
		tr[p].laz=1;
	}
	void down(int p){
		if(!tr[p].laz) return ;
		turn(p<<1,tr[p].C);
		turn(p<<1|1,tr[p].C);
		tr[p].laz=0;
	}
	int find(int p,int l,int r,int x){
		if(tr[p].C>=x) return -1;
		if(l==r) return l;
		down(p);
		int mid=l+r>>1;
		int ans=find(p<<1|1,mid+1,r,x);
		if(ans==-1) ans=find(p<<1,l,mid,x);
		return ans;
	}
	void change(int p,int l,int r,int L,int R,int x){
		if(L<=l&&r<=R){
			ans=ans*Pow(tr[p].prd,mod-2)%mod;
			turn(p,x);
			ans=ans*tr[p].prd%mod;
			return ;
		}
		down(p);
		int mid=l+r>>1;
		if(L<=mid) change(p<<1,l,mid,L,R,x);
		if(mid<R) change(p<<1|1,mid+1,r,L,R,x);
		up(p);
	}
	void change_node(int p,int l,int r,int x,int y){
		if(l==r){
			if(tr[p].num) ans=ans*Pow(tr[p].prd,mod-2)%mod;
			else ans=ans*Pow(tr[p].B,mod-2)%mod;
			tr[p].B=y;
			solve(p,l);
			return ;
		}
		down(p);
		int mid=l+r>>1;
		if(x<=mid) change_node(p<<1,l,mid,x,y);
		else change_node(p<<1|1,mid+1,r,x,y);
		up(p);
	}
	void keep(int p,int l,int r,int x){
		if(l==r){
			ans=ans*Pow(tr[p].prd,mod-2)%mod;
			solve(p,l);
			return ;
		}
		down(p);
		int mid=l+r>>1;
		if(x<=mid) keep(p<<1,l,mid,x);
		else keep(p<<1|1,mid+1,r,x);
		up(p);
	}
	void make(){
		while(tr[1].mn<0) keep(1,1,n,tr[1].pm);
	}
}st;

int main(){
	freopen("magic.in","r",stdin);
	freopen("magic.out","w",stdout);
	cin>>n>>Q;
	for(int i=1;i<=n;i++) scanf("%d",&A[i]),C[i]=max(C[i-1],A[i]);
	for(int i=1;i<=n;i++) scanf("%d",&B[i]);
	st.build(1,1,n);
	while(Q--){
		scanf("%d%d%d",&opt,&x,&y);
		if(!opt){
			int t=st.find(1,1,n,y);
			if(t>=x) st.change(1,1,n,x,t,y);
		}
		else st.change_node(1,1,n,x,y);
		st.make();
		printf("%lld\n",ans);
	}
	return 0;
}