#include <bits/stdc++.h>
using namespace std;
int m,n,a[100010],b,c,d,e,o[12][12],jc[12];
long long ans=0;
const int mod=1e9+7;
inline int read()
{
    int x=0,f=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')
            f=-1;
        ch=getchar();
    }
    while(ch>='0' && ch<='9')
        x=x*10+ch-'0',ch=getchar();
    return x*f;
}
struct hh
{
	int lazy1;
	int lazy2=-1;
	long long val;
}t[11][100010<<2];
long long qp(int x,int y)
{
	long long ans=1,temp=x;
	while(y)
	{
		if(y&1)
		{
			ans=(ans*temp)%mod;
		}
		y>>=1;
		temp=(temp*temp)%mod;
	}
	return ans;
}
void build(int l,int r,int u)
{
	if(l==r)
	{
		for(int i=0;i<=10;++i)
		{
			t[i][l].val=qp(a[l],i);
		}
		return;
	}
	int mid=(l+r)/2;
	build(l,mid,u*2);
	build(mid+1,r,u*2+1);
}
void work(int kind,int u,int x)
{
	for(int i=1;i<=kind;++i)
	{
		t[kind][u].val=(t[kind][u].val+((o[kind][i]*qp(x,i))%mod*t[kind-i][u].val)%mod)%mod;
	}
}
void pushdown(int u,int l,int r)
{
	if(t[1][u].lazy2==-1&&!t[1][u].lazy1)
	{
		return;
	}
	for(int i=10;i>=0;--i)
	{
		if(t[i][u].lazy2!=-1)
		{
			t[i][u*2].lazy2=t[i][u].lazy2;
			t[i][u*2+1].lazy2=t[i][u].lazy2;
			t[i][u*2].lazy1=t[i][u].lazy1;
			t[i][u*2+1].lazy1=t[i][u].lazy1;
			t[i][u*2].val=qp(t[i][u].lazy2,i)*(r-l+1)%mod;
			t[i][u*2+1].val=qp(t[i][u].lazy2,i)*(r-l+1)%mod;
			work(i,u*2,t[i][u].lazy1);
			work(i,u*2+1,t[i][u].lazy1);
		}
		else
		{
			t[i][u*2].lazy1+=t[i][u].lazy1;
			t[i][u*2+1].lazy1+=t[i][u].lazy1;
			work(i,u*2,t[i][u].lazy1);
			work(i,u*2+1,t[i][u].lazy1);
		}
		t[i][u].lazy2=-1;
		t[i][u].lazy1=0;
	}
}
void change1(int u,int L,int R)
{
	if(L>=c&&R<=d)
	{
		for(int i=10;i>=0;--i)
		{
			work(i,u,e);
			t[i][u].lazy1+=e;
		}
		return;
	}
	int mid=(L+R)/2;
	if(c<=mid)
	{
		pushdown(u,L,R);
		change1(u*2,L,mid);
	}
	if(d>=mid+1)
	{
		pushdown(u,L,R);
		change1(u*2+1,mid+1,R);
	}
}
void change2(int u,int L,int R)
{
	if(L>=c&&R<=d)
	{
		for(int i=0;i<=10;++i)
		{
			t[i][u].val=qp(e,i)*(R-L+1)%mod;
			t[i][u].lazy2=e;
			t[i][u].lazy1=0;
		}
		return;
	}
	int mid=(L+R)/2;
	if(c<=mid)
	{
		pushdown(u,L,R);
		change2(u*2,L,mid);
	}
	if(d>=mid+1)
	{
		pushdown(u,L,R);
		change2(u*2+1,mid+1,R);
	}
}
void query(int u,int L,int R)
{
	if(L>=c&&R<=d)
	{
		ans=(ans+t[e][u].val)%mod;
		return;
	}
	int mid=(L+R)/2;
	if(c<=mid)
	{
		pushdown(u,L,R);
		query(u*2,L,mid);
	}
	if(d>=mid+1)
	{
		pushdown(u,L,R);
		query(u*2+1,mid+1,R);
	}
}
int main()
{
	m=read();
	n=read();
	for(int i=1;i<=m;++i)
	{
		a[i]=read();
	}
	jc[0]=o[0][0]=1;
	for(int i=1;i<=10;++i)
	{
		jc[i]=jc[i-1]*i;
	}
	for(int i=1;i<=10;++i)
	{
		for(int j=i;j>=0;--j)
		{
			o[i][j]=jc[i]/jc[j]/jc[i-j];
		}
	}
	build(1,m,1);
	for(int i=1;i<=n;++i)
	{
		b=read();
		c=read();
		d=read();
		e=read();
		if(b==1)
		{
			pushdown(1,1,m);
			change1(1,1,m);
		}
		else if(b==2)
		{
			pushdown(1,1,m);
			change2(1,1,m);
		}
		else
		{
			ans=0;
			pushdown(1,1,m);
			query(1,1,m);
			printf("%ld\n",ans);
		}
	}
}