#include<bits/stdc++.h>
using namespace std;
int c[15][15];
int sum[400005][15];
const int mod=1000000007;
int a[100005][15],lz1[4000005],lz2[4000005];
int d[15];
void pushup(int x){
	for(int i=0;i<=10;i++)sum[x][i]=(sum[x<<1][i]+sum[x<<1|1][i])%mod;
}
void build(int l,int r,int x){
	if(l==r){
		for(int i=0;i<=10;i++)sum[x][i]=a[l][i];
		return;
	}
	int mid=(l+r)/2;
	build(l,mid,x<<1);
	build(mid+1,r,x<<1|1);
	pushup(x);
}
void pushdown(int x){
	if(lz2[x]<0x3f3f3f3f){
		int v=lz2[x];
		for(int i=1;i<=10;i++)sum[x<<1][i]=1LL*sum[x<<1][i-1]*v%mod,sum[x<<1|1][i]=1LL*sum[x<<1|1][i-1]*v%mod;
		lz2[x<<1]=lz2[x<<1|1]=lz2[x];
		lz1[x<<1]=lz1[x<<1|1]=0;
	}
	else if(lz1[x]!=0){
		int v=lz1[x];
		d[0]=1;
		for(int i=1;i<=10;i++)d[i]=1LL*d[i-1]*v%mod;
		for(int i=10;i>=0;i--)for(int j=0;j<i;j++)sum[x<<1][i]=(sum[x<<1][i]+1LL*sum[x<<1][j]*c[i][j]%mod*d[i-j])%mod,sum[x<<1|1][i]=(sum[x<<1|1][i]+1LL*sum[x<<1|1][j]*c[i][j]%mod*d[i-j])%mod;;
		if(lz2[x<<1]==0x3f3f3f3f)lz1[x<<1]=(lz1[x<<1]+v)%mod;
		else lz2[x<<1]=(lz2[x<<1]+v)%mod;
		if(lz2[x<<1|1]==0x3f3f3f3f)lz1[x<<1|1]=(lz1[x<<1|1]+v)%mod;
		else lz2[x<<1|1]=(lz2[x<<1|1]+v)%mod;
	}
	lz1[x]=0;
	lz2[x]=0x3f3f3f3f;
}
void update(int l,int r,int x,int L,int R,int op,int v){
	if(L<=l && r<=R){
		if(op==1){
			d[0]=1;
			for(int i=1;i<=10;i++)d[i]=1LL*d[i-1]*v%mod;
			for(int i=10;i>=0;i--)for(int j=0;j<i;j++)sum[x][i]=(sum[x][i]+1LL*sum[x][j]*c[i][j]%mod*d[i-j])%mod;
			if(lz2[x]==0x3f3f3f3f)lz1[x]=(lz1[x]+v)%mod;
			else lz2[x]=(lz2[x]+v)%mod;
		}
		else{
			for(int i=1;i<=10;i++)sum[x][i]=1LL*sum[x][i-1]*v%mod;
			lz2[x]=v;
			lz1[x]=0;
		}
		return;
	}
	int mid=(l+r)/2;
	pushdown(x);
	if(L<=mid)update(l,mid,x<<1,L,R,op,v);
	if(R>mid)update(mid+1,r,x<<1|1,L,R,op,v);
	pushup(x);
}
int query(int l,int r,int x,int L,int R,int v){
	if(L<=l && r<=R)return sum[x][v];
	int mid=(l+r)/2;
	int ans=0;
	pushdown(x);
	if(L<=mid)ans=(ans+query(l,mid,x<<1,L,R,v))%mod;
	if(R>mid)ans=(ans+query(mid+1,r,x<<1|1,L,R,v))%mod;
	pushup(x);
	return ans;
}
int main(){
//	freopen("t3.in","r",stdin);
	int n,m;
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)scanf("%d",&a[i][1]);
	for(int i=1;i<=n;i++){
		a[i][0]=1;
		for(int j=2;j<=10;j++)a[i][j]=1LL*a[i][j-1]*a[i][1]%mod;
	}
	for(int i=0;i<=10;i++){
		c[i][0]=1;
		for(int j=1;j<=i;j++)c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;
	}
	memset(lz2,0x3f,sizeof(lz2));
	build(1,n,1);
	for(int i=1;i<=m;i++){
		int op,l,r,x;
		scanf("%d%d%d%d",&op,&l,&r,&x);
		if(op<=2)update(1,n,1,l,r,op,x);
		else printf("%d\n",query(1,n,1,l,r,x));
	}
	return 0;
}
/*
5 12
1 1 2 4 6
1 1 5 2
3 1 5 0
3 1 5 1
3 1 5 2
3 1 5 3
3 1 5 4
3 1 5 5
3 1 5 6
3 1 5 7
3 1 5 8
3 1 5 9
3 1 5 10
*/