BZOJ 1798: [Ahoi2009]Seq 维护序列seq|线段树

这道题和2631对标记的处理方法是一样的..以下直接复制

将乘法和加法合并起来算做一个标记，$mul$数组和$plus$数组两个数组算做一个标记，设权值为$x$那么标记的意思就是把所有子节点的的$x$都变为：

$$mul*x+plus$$ 这样只要多维护一个 $size$数组就可以维护，下传标记也很好写
附上链接 http://blog.csdn.net/ws_yzy/article/details/50734951

由于开$longlong$据说会很慢，诶..然后我就想到了用直接$int$直接快速乘不就好了吗，然后发现$int$的快速乘直接TLE！！！30s的限制都过不去..是不是应该改名叫慢速乘
$longlong$实测只需要4s左右..smg
突然发现没有线段树分类.......

#include<cstdio>  
#include<cstdlib>  
#include<cstring>  
#include<cmath>  
#include<queue>  
#include<set>  
#include<map>  
#include<vector>  
#include<algorithm>  
#include<iostream>  
#define N 503333  
using namespace std;  
int sc()  
{  
    int i=0,f=1; char c=getchar();  
    while(c>'9'||c<'0'){if(c=='-')f=-1; c=getchar();}  
    while(c>='0'&&c<='9') i=i*10+c-'0',c=getchar();  
    return i*f;  
}  
struct W{int l,r;long long sum,plu,mul;}tr[N];
int v[N],n,m,P;
void build(int x,int l,int r)
{
    tr[x].l=l,tr[x].r=r,tr[x].plu=0,tr[x].mul=1;
    if(l==r)
    {
        tr[x].sum=v[l]%P;
        return;
    }
    int mid=l+r>>1;
    build(x<<1,l,mid),build(x<<1|1,mid+1,r);
    tr[x].sum=(tr[x<<1].sum+tr[x<<1|1].sum)%P;
    //cout << l <<" "<< r <<" "<< tr[x].sum<< endl;
}
void cal(int x,long long MUL,long long PLU)
{
    tr[x].sum=(MUL*tr[x].sum+PLU*(tr[x].r-tr[x].l+1))%P;
    tr[x].plu=(MUL*tr[x].plu+PLU)%P;
    tr[x].mul=MUL*tr[x].mul%P;
}
void push_down(int x)
{
    if(tr[x].l==tr[x].r)return;
    if(tr[x].plu!=0||tr[x].mul!=1)
        cal(x<<1,tr[x].mul,tr[x].plu),
        cal(x<<1|1,tr[x].mul,tr[x].plu);
    tr[x].mul=1,tr[x].plu=0;
}
void change(int x,int l,int r,int MUL,int PLU)
{
    push_down(x);
    int L=tr[x].l,R=tr[x].r,mid=L+R>>1;
    if(L==l&&R==r)
    {
        cal(x,MUL,PLU);
        return;
    }
    if(r<=mid) change(x<<1,l,r,MUL,PLU);
    else if(l>mid) change(x<<1|1,l,r,MUL,PLU);
    else change(x<<1,l,mid,MUL,PLU),change(x<<1|1,mid+1,r,MUL,PLU);
    tr[x].sum=(tr[x<<1].sum+tr[x<<1|1].sum)%P;
}   
int query(int x,int l,int r)
{
    push_down(x);
    int L=tr[x].l,R=tr[x].r,mid=L+R>>1;
    if(L==l&&R==r)return tr[x].sum;
    if(r<=mid)return query(x<<1,l,r);
    else if(l>mid) return query(x<<1|1,l,r);
    else return (query(x<<1,l,mid)+query(x<<1|1,mid+1,r))%P;
}   
int main()
{
    n=sc(),P=sc();
    for(int i=1;i<=n;i++)v[i]=sc();
    build(1,1,n);
    m=sc();
    while(m--)
    {
        int opt=sc();
        if(opt==1)
        {
            int l=sc(),r=sc(),c=sc();
            change(1,l,r,c,0);
        }
        else if(opt==2)
        {
            int l=sc(),r=sc(),c=sc();
            change(1,l,r,1,c);
        }
        else 
        {
            int l=sc(),r=sc();
            printf("%d\n",query(1,l,r));
        }
    }
    return 0;
}