洛谷P3747 [六省联考2017]相逢是问候

传送门

题解

扩展欧拉定理。

线段树维护，已经全改到底了的节点就不管，不然暴力修改下去。

//Achen
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<queue>
#include<cmath>
const int N=+;
#define For(i,a,b) for(int i=(a);i<=(b);i++)
#define Rep(i,a,b) for(int i=(a);i>=(b);i--)
typedef long long LL;
using namespace std;
LL n,m,p,pp,c,tp,cnt,phi[],pr[][N],pr2[][N],mini[];

template<typename T>void read(T &x)  {
    char ch=getchar(); x=; T f=;
    while(ch!='-'&&(ch<''||ch>'')) ch=getchar();
    if(ch=='-') f=-,ch=getchar();
    for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; x*=f;
}

LL a[N],sg[N<<],tot[N<<];

int get_phi(int x) {
    int up=sqrt(x),rs=x;
    For(i,,up) {
        if(!(x%i)) {
            rs=rs-rs/i;
            while(!(x%i)) x/=i;
        }
    }
    if(x!=) rs=rs-rs/x;
    return rs;
}

LL get(int i,int a) { //c^a%phi[i]
    return pr2[i][a/pp]*pr[i][a%pp]%phi[i];
}

#define lc x<<1
#define rc x<<1|1
#define mid ((l+r)>>1)
void build(int x,int l,int r) {
    if(l==r) { sg[x]=a[l]; return; }
    build(lc,l,mid); build(rc,mid+,r);
    sg[x]=(sg[lc]+sg[rc])%p;
    tot[x]=min(tot[lc],tot[rc]);
}

int lz[N<<];
void down(int x,int l_len,int r_len) {
    if(!lz[x]) return;
    if(l_len) { sg[lc]=l_len*c%p; tot[lc]=cnt; lz[lc]=; }
    if(r_len) { sg[rc]=r_len*c%p; tot[rc]=cnt; lz[rc]=; }
    lz[x]=;
}

void update(int x,int l,int r,int ql,int qr) {
    if(tot[x]>=cnt) return;
    if(l>=ql&&r<=qr&&l==r) {
        tot[x]++;
        sg[x]=a[l]>=phi[tot[x]+]?a[l]%phi[tot[x]+]+phi[tot[x]+]:a[l];
        Rep(i,tot[x],) {
            int now=sg[x]>=mini[i];
            sg[x]=get(i,sg[x]);
            if(now&&i!=) sg[x]+=phi[i];
        }
        return ;
    }
    if(ql<=mid) update(lc,l,mid,ql,qr);
    if(qr>mid) update(rc,mid+,r,ql,qr);
    sg[x]=(sg[lc]+sg[rc])%p;
    tot[x]=min(tot[lc],tot[rc]);
}

LL qry(int x,int l,int r,int ql,int qr) {
    if(l>=ql&&r<=qr) return sg[x];
    down(x,mid-l+,r-mid);
    if(qr<=mid) return qry(lc,l,mid,ql,qr);
    if(ql>mid) return qry(rc,mid+,r,ql,qr);
    return (qry(lc,l,mid,ql,qr)+qry(rc,mid+,r,ql,qr))%p;
}

//#define DEBUG
int main() {
#ifdef DEBUG
    freopen("bzoj4869相逢是问候.in","r",stdin);
    freopen("bzoj4869相逢是问候.out","w",stdout);
#endif
    read(n); read(m); read(p); read(c);
    For(i,,n) read(a[i]);
    build(,,n);
    tp=p; phi[cnt=]=tp;
    while(tp!=) {
        phi[++cnt]=get_phi(tp);
        tp=phi[cnt];
    } phi[++cnt]=; phi[cnt+]=;
    pp=1e4;
    For(i,,cnt) {
        tp=; int tc=;
        if(c!=) {
            for(;;) {
                tp*=c; tc++;
                if(tp>=phi[i]) {
                    mini[i]=tc;
                    break;
                }
            }
        } else mini[i]=1e8;
        pr[i][]=pr2[i][]=;
        For(j,,pp)
            pr[i][j]=pr[i][j-]*c%phi[i];
        For(j,,pp)
            pr2[i][j]=pr2[i][j-]*pr[i][pp]%phi[i];
    } if(c!=) { mini[cnt]=; mini[cnt-]=; }
    while(m--) {
        int o,l,r;
        read(o); read(l); read(r);
        if(!o) update(,,n,l,r);
        else printf("%lld\n",qry(,,n,l,r)%p);
    }
    return ;
}