#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

int n,m;
ll p,c,res;
ll a[50010],t[50010];

ll Pow(ll x,ll a){
    if (!a) return 1;
    if (a==1) return x;
    if (a%2) return Pow(x,a/2)*Pow(x,a/2)*x;
    return Pow(x,a/2)*Pow(x,a/2);
}

int lowbit(int x){
    return x&-x;
}

void add(int x,ll a){
    while (x<=n){
        t[x]+=a,t[x]%=p,x+=lowbit(x);
    }
    return;
}

ll query(int x){
    ll cnt=0;
    while (x){
        cnt+=t[x],cnt%=p,x-=lowbit(x);
    }
    return cnt;
}

int main(){
    freopen("verbinden.in","r",stdin);
    freopen("verbinden.out","w",stdout);
    
    int k,l,r;
    scanf("%d %d %lld %lld",&n,&m,&p,&c);
    for (int i=1;i<=n;i++){
        scanf("%lld",&a[i]);
        add(i,a[i]);
    }
    if (n<=100 && m<=100){
        for (int x=1;x<=m;x++){
            scanf("%d %d %d",&k,&l,&r);
            if (!k){
                for (int i=l;i<=r;i++) a[i]=Pow(c,a[i]);
            }
            else{
                res=0;
                for (int i=l;i<=r;i++) res+=a[i],res%=p;
                printf("%lld\n",res);
            }
        }
    }
    else{
        for (int x=1;x<=m;x++){
            scanf("%d %d %d",&k,&l,&r);
            if (!k){
                ll w;
                for (int i=l;i<=r;i++){
                    w=Pow(c,a[i])-a[i];
                    add(i,w);
                }
            }
            else{
                res=query(r)-query(l-1);
                printf("%lld\n",res);
            }
        }
    }
    
    return 0;
}