分块求条
查看原帖
分块求条
704275
无名之雾楼主2024/6/15 14:31
#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;
			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;
			add[p]=0;
		}
		for(int i=l;i<=ed[p];i++){
			sum[p]+=(a[i]*(k-1))%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;
			add[q]=0;
		}
		for(int i=st[q];i<=r;i++){
			sum[p]+=(a[i]*(k-1))%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){
		for(int i=l;i<=r;i++)a[i]=(a[i]+k)%mod;
		sum[p]=(sum[p]+k*(r-l+1))%mod; 
	}
	else{
		for(int i=p+1;i<=q-1;i++)add[i]=(add[i]+k)%mod;
		for(int i=l;i<=ed[p];i++)a[i]=(a[i]+k)%mod;
		sum[p]=(sum[p]+k*(ed[p]-l+1))%mod;
		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*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])%mod;
		ans=(ans*mul[p])%mod;
		ans=(ans+add[p]*(r-l+1))%mod;
		for(int i=st[q];i<=r;i++)ans=(ans+a[i])%mod;
		ans=(ans*mul[q])%mod;
		ans=(ans+add[q]*(r-l+1))%mod;
	}
	cout<<ans<<"\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<=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];
		}
	}
	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);
		}
	}
	return 0;
}
2024/6/15 14:31
加载中...