分块30pts求助
查看原帖
分块30pts求助
704275
无名之雾楼主2024/6/15 17:06
#include<bits/stdc++.h>
using namespace std;
#define int long long
inline int read(){
   int s=0,w=1;
   char ch=getchar();
   while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
   while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();
   return s*w;
}
inline void write(int x){
    if(x==0){putchar('0');return;}
	int len=0,k1=x,c[10005];
	if(k1<0)k1=-k1,putchar('-');
	while(k1)c[len++]=k1%10+'0',k1/=10;
	while(len--)putchar(c[len]);
}
const int N=1e6+5;
int pos[N],st[N],ed[N],sum[N],add[N],mul[N],a[N],mod;
void upd1(int l,int r,int k){
	int p=pos[l],q=pos[r];
	if(p==q){
		if(add[p]!=0){
			for(int i=st[p];i<=ed[q];i++)a[i]=(a[i]+add[p])%mod;
			add[p]=0;
		}
		for(int i=l;i<=r;i++){
			sum[p]+=(a[i]*(k-1))%mod;
			sum[p]%=mod;
			a[i]=(a[i]*k)%mod;
		}
	}	
	else{
		for(int i=p+1;i<=q-1;i++){
			if(add[i]!=0)add[i]=(add[i]*k)%mod;
			mul[i]=(mul[i]*k)%mod;
		}
		if(add[p]!=0){
			for(int i=st[p];i<=ed[p];i++)a[i]=(a[i]+add[p])%mod;
			sum[p]+=add[p]*(ed[p]-st[p]+1)%mod;
			sum[p]%=mod;
			add[p]=0;
		}
		for(int i=l;i<=ed[p];i++){
			sum[p]+=(a[i]*(k-1))%mod;
			sum[p]%=mod;
			a[i]=(a[i]*k)%mod;
		}
		if(add[q]!=0){
			for(int i=st[q];i<=ed[q];i++)a[i]=(a[i]+add[q])%mod;
			sum[q]+=add[q]*(ed[q]-st[q]+1)%mod;
			sum[q]%=mod;
			add[q]=0;
		}
		for(int i=st[q];i<=r;i++){
			sum[q]+=(a[i]*(k-1))%mod;
			sum[q]%=mod;
			a[i]=(a[i]*k)%mod;
		}
	}
}
void upd2(int l,int r,int k){
	int p=pos[l],q=pos[r];
	if(p==q){
		if(mul[p]!=1){
			for(int i=st[p];i<=ed[p];i++)a[i]=(a[i]*mul[p])%mod;
			sum[p]=(sum[p]*mul[p])%mod;
			mul[p]=1;
		}
		for(int i=l;i<=r;i++)a[i]=(a[i]+k)%mod;
		sum[p]=(sum[p]+k*(r-l+1)%mod)%mod; 
	}
	else{
		for(int i=p+1;i<=q-1;i++)add[i]=(add[i]+k)%mod;
		if(mul[p]!=1){
			for(int i=st[p];i<=ed[p];i++)a[i]=(a[i]*mul[p])%mod;
			sum[p]=(sum[p]*mul[p])%mod;
			mul[p]=1;
		}
		for(int i=l;i<=ed[p];i++)a[i]=(a[i]+k)%mod;
		sum[p]=(sum[p]+k*(ed[p]-l+1))%mod;
		if(mul[q]!=1){
			for(int i=st[q];i<=ed[q];i++)a[i]=(a[i]*mul[q])%mod;
			sum[q]=(sum[q]*mul[q])%mod;
			mul[q]=1;
		}
		for(int i=st[q];i<=r;i++)a[i]=(a[i]+k)%mod;
		sum[q]=(sum[q]+k*(r-st[q]+1))%mod;
	}
}
void ask(int l,int r){
	int p=pos[l],q=pos[r];
	int ans=0;
	if(p==q){
		for(int i=l;i<=r;i++)ans+=a[i];
		ans=(ans%mod*mul[p])%mod;
		ans=(ans+add[p]*(r-l+1))%mod;
	}
	else{
		for(int i=p+1;i<=q-1;i++)ans=(ans+(sum[i]*mul[i])%mod+add[i]*(ed[i]-st[i]+1))%mod;
		for(int i=l;i<=ed[p];i++)ans=(ans+a[i]*mul[p])%mod;
		ans=(ans+add[p]*(ed[p]-l+1))%mod;
		for(int i=st[q];i<=r;i++)ans=(ans+a[i]*mul[q])%mod;
		ans=(ans+add[q]*(r-st[q]+1))%mod;
	}
	cout<<ans%mod<<"\n";
} 
signed main(){
	int n=read(),m=read();
	mod=read();
	for(int i=1;i<=n;i++)a[i]=read(); 
	int block=sqrt(n);
	int t=n/block;
	if(n%block)t++;
	for(int i=1;i<=t;i++){
		st[i]=(i-1)*block+1;
		ed[i]=i*block; 
	}	
	ed[t]=n;
	for(int i=1;i<=t;i++) mul[i]=1; 
	for(int i=1;i<=n;i++)pos[i]=(i-1)/block+1;
	for(int i=1;i<=t;i++){
		for(int j=st[i];j<=ed[i];j++){
			sum[i]+=a[j];
			sum[i]%=mod;
		}
	}
	for(int i=1;i<=m;i++){
		int op=read();
		if(op==1){
			int x=read(),y=read(),k=read();
			upd1(x,y,k);
		}
		if(op==2){
			int x=read(),y=read(),k=read();
			upd2(x,y,k);
		}
		if(op==3){
			int x=read(),y=read();
			ask(x,y);
		}
//		cout<<"___\n";
//		for(int i=1;i<=n;i++) cout<<"i"<<i<<endl<<(a[i]*mul[pos[i]]+add[pos[i]])%mod<<"\n"<<sum[pos[i]]<<" "<<a[i]<<" "<<mul[pos[i]]<<" "<<add[pos[i]]<<endl;
//		cout<<"\n___"<<endl;
	}
	return 0;
}
2024/6/15 17:06
加载中...