蒟蒻线段树30分求调
查看原帖
蒟蒻线段树30分求调
889845
Wunsch楼主2023/2/1 20:29

没用结构体代码有点丑望别嫌弃

个人感觉就是乘那部分的问题

但蒟蒻调不出来

#include<iostream>
#include<cstdio>
#include<cmath>
#include<string>
#include<cstring>
#define ll long long
#define LS o<<1
#define RS LS|1

using namespace std;

ll n,m,t[500001],a[500001],laz[500001],lazc[500001],x,y,k,p,q;

void build(ll o,ll l,ll r){
	if(l==r){cin>>t[o];return;}
	ll mid=l+r>>1;
	build(LS,l,mid);
	build(RS,mid+1,r);
	t[o]=t[LS]%p+t[RS]%p;
}

void ddown(ll o,ll l,ll r){//乘下放 
    if(l==r||lazc[o]==1)return;
    ll mid=l+r>>1;
	lazc[LS]=((lazc[LS]%p)*lazc[o]%p)%p;
	lazc[RS]=((lazc[RS]%p)*lazc[o]%p)%p;
    t[LS]=((t[LS]%p)*lazc[o]%p)%p;t[LS]%=p;
    t[RS]=((t[RS]%p)*lazc[o]%p)%p;t[RS]%=p;
    lazc[o]=1;
}

void down(ll o,ll l,ll r){//加下放 
	if(l==r||!laz[o])return;
	ll mid=l+r>>1;
	laz[LS]+=laz[o];
	laz[RS]+=laz[o];
	t[LS]+=laz[o]*(mid-l+1)%p;t[LS]%=p;
	t[RS]+=laz[o]*(r-mid)%p;t[RS]%=p;
	laz[o]=0;
}

void cchenge(ll o,ll l,ll r,ll x,ll y,ll c){//乘 
    down(o,l,r);
    ddown(o,l,r);
    if(x<=l&&r<=y){lazc[o]=(lazc[o]*c)%p;t[o]%=p;t[o]=((t[o]%p)*(c%p))%p;return;}
    down(o,l,r);//多下放一次总错不了虽然说没什么用
    ddown(o,l,r);
    ll mid=l+r>>1;
    if(x<=mid)cchenge(LS,l,mid,x,y,c);
    if(y>mid)cchenge(RS,mid+1,r,x,y,c);
    t[o]=t[LS]+t[RS];
}

void change(ll o,ll l,ll r,ll x,ll y,ll sum){//加 
    down(o,l,r);
    ddown(o,l,r);
	if(x<=l&&r<=y){laz[o]+=sum;t[o]+=((r-l+1)*sum)%p;return;}
    down(o,l,r);
    ddown(o,l,r);
	ll mid=l+r>>1;
	if(x<=mid)change(LS,l,mid,x,y,sum);
	if(y>mid)change(RS,mid+1,r,x,y,sum);
	t[o]=(t[LS]%p+t[RS]%p)%p;
}

ll qury(ll o,ll l,ll r,ll x,ll y){
    down(o,l,r);
    ddown(o,l,r);
	if(x<=l&&r<=y)return t[o];
	ll ret=0;
	ll mid=l+r>>1;
	if(x<=mid)ret+=qury(LS,l,mid,x,y);
	if(y>mid)ret+=qury(RS,mid+1,r,x,y);
	return ret;
}

int main(){
	cin>>n>>m>>p;
    for(int i=1;i<=500000;i++){
        lazc[i]=1;
    }
	build(1,1,n);// 
	for(int i=1;i<=m;i++){
		cin>>q;
		if(q==1){
			cin>>x>>y>>k;cchenge(1,1,n,x,y,k);
    	}
		else if(q==2){
			cin>>x>>y>>k;change(1,1,n,x,y,k);
		}
        else{
            cin>>x>>y;cout<<qury(1,1,n,x,y)%p<<endl;
        }
	}
	return 0;
}
2023/2/1 20:29
加载中...