#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=1e5+3,T=200,MOD=1e9+7,INF=2e9;
int n,m;
LL a[N];
int l[N],r[N],belong[N];
int turn[N];
LL tag[N],sum[N][20],C[12][12];
void add(LL &x,LL y){
	x+=y;
	if(x>=MOD) x-=MOD;
	if(x<0) x+=MOD;
}
void pre(int x){
	for(int i=l[x];i<=r[x];i++){
		if(turn[x]!=INF) a[i]=turn[x];
		add(a[i],tag[x]);
	}
	turn[x]=INF,tag[x]=0;
}
void work(int x){
	for(int k=0;k<=10;k++) sum[x][k]=0;
	for(int i=l[x];i<=r[x];i++){
	    int mul=1;
		for(int k=0;k<=10;k++){
			add(sum[x][k],mul);
			mul=1ll*mul*a[i]%MOD; 
		}
	}
}
void modify1(int s,int t,int val){
	int L=belong[s],R=belong[t];
	if(L==R){
		pre(L);
		for(int i=s;i<=t;i++) add(a[i],val);
		work(L);
		return;
	}
	pre(L),pre(R);
	for(int i=s;i<=r[L];i++) add(a[i],val);
	for(int i=l[R];i<=t;i++) add(a[i],val);
	work(L),work(R);
	for(int i=L+1;i<=R-1;i++) add(tag[i],val);
}
void modify2(int s,int t,int val){
	int L=belong[s],R=belong[t];
	pre(L),pre(R);
	if(L==R){
		pre(L);
		for(int i=s;i<=t;i++) a[i]=val;
		work(L);
		return;
	}
	for(int i=s;i<=r[L];i++) a[i]=val;
	for(int i=l[R];i<=t;i++) a[i]=val;
	work(L),work(R);
	for(int i=L+1;i<=R-1;i++) turn[i]=val,tag[i]=0;
}
LL Q(int s,int t,int val){
	int L=belong[s],R=belong[t];
	LL ans=0;
	if(L==R){
		pre(L);work(L);
		for(int i=s;i<=t;i++){
	    	LL mul=1;
	    	for(int j=1;j<=val;j++) mul=mul*a[i]%MOD;
	    	add(ans,mul); 
    	}
		return ans;
	}
	pre(L),pre(R);
	work(L),work(R);
	for(int i=s;i<=r[L];i++){
		LL mul=1;
		for(int j=1;j<=val;j++) mul=mul*a[i]%MOD;
		add(ans,mul); 
	}
	for(int i=l[R];i<=t;i++){
		LL mul=1;
		for(int j=1;j<=val;j++) mul=mul*a[i]%MOD;
		add(ans,mul); 
	}
	for(int i=L+1;i<=R-1;i++){
		LL mul=1;
		LL la=ans;
		for(int k=0;k<=val;k++){
		    add(ans,C[val][k]*mul%MOD*sum[i][val-k]%MOD);
			mul=mul*tag[i]%MOD;
		}
	}
	return ans;
}
int Main()
{
	for(int i=0;i<=10;i++) C[i][0]=1;
	for(int i=1;i<=10;i++) 
	    for(int j=1;j<=10;j++)
	        C[i][j]=C[i-1][j]+C[i-1][j-1];
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]),turn[i]=INF;
    int num=n/T+1;
    for(int i=1;i<=num;i++){
    	l[i]=(i-1)*T+1,r[i]=i*T;
	    for(int j=l[i];j<=r[i];j++)
		    belong[j]=i;
		work(i);
	}
    while(m--){
    	int op,l,r,x;
    	scanf("%d%d%d%d",&op,&l,&r,&x);
    	if(op==1) modify1(l,r,x);
    	if(op==2) modify2(l,r,x);
    	if(op==3) printf("%lld\n",Q(l,r,x));
	}
	return 0;
}
int main()
{
    int Case=1;
    while(Case--) Main();
	return 0;
}