#include<bits/stdc++.h>
using namespace std;
const int MAXN=1e5+5,mod=1e9+7;
int n,m,a[MAXN],type,l,r,x,cf[MAXN],ans;
int qpow(int x,int y) {
	int ret=1;
	while(y) {
		if(y&1) {
			ret=ret*x%mod;
		}
		y>>=1;
		x=x*x%mod;
	}
	return ret;
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        
        cf[i]=a[i]-a[i-1];
        
    }
    
    while(m--){
        scanf("%d%d%d%d",&type,&l,&r,&x);
        if(type==1){
            cf[l]+=x;
            cf[r+1]-=x;
        }else if(type==2){
            int k=0,j=0;
            for(int i=1;i<=l-1;i++){
                k+=cf[i];
            }
            j+=k;
            for(int i=l;i<=r+1;i++){
                j+=cf[i];
            }
            for(int i=l+1;i<=r;i++){
                cf[i]=0;
            }
            cf[l]=x-k;
            cf[r+1]=j-x;
        }else{
            int k=0;ans=0;
            for(int i=1;i<=l-1;i++){
                k+=cf[i];
            }
            for(int i=l;i<=r;i++){
                k+=cf[i];
                ans+=qpow(k,x);
                ans%=mod;
            }
            printf("%d",ans);
        }
    }

    return 0;
}