#include<bits/stdc++.h>
#define int long long
using namespace std;
int n, m, a[100005], c[15], x[15][15], jie[15], inv[15], op;
const int mod = 1000000007;
inline int ksm(int p, int q) {
	if (!q) return 1;
	if (q == 1) return p;
	if (q & 1) return p * ksm(p * p % mod, q / 2) % mod;
	return ksm(p * p % mod, q / 2);
}
inline int C(int p, int q) {
	return ((jie[p] * inv[q] % mod) * inv[p - q] % mod);
}
struct node {
	int l, r, num[11], addtag, changetag;
} b[400005];
inline void update(int p) {
	for (int i = 1; i <= 10; i++) {
		b[p].num[i] = b[2 * p].num[i] + b[2 * p + 1].num[i];
	}
}
inline void push_down(int p) {
	if (b[p].l == b[p].r) return;
	if (b[p].changetag != 1000000008) {
		for (int i = 1; i <= 10; i++) {
			b[2 * p].num[i] = b[2 * p].num[i - 1] * b[p].changetag % mod;
			b[2 * p + 1].num[i] = b[2 * p + 1].num[i - 1] * b[p].changetag % mod;
		}
		b[2 * p].addtag = b[2 * p + 1].addtag = 0;
		b[2 * p].changetag = b[2 * p + 1].changetag = b[p].changetag;
		b[p].changetag = 1000000008;
	}
	if (b[p].addtag != 0) {
		c[0] = 1;
		c[1] = b[p].addtag;
		for (int i = 2; i <= 10; i++) {
			c[i] = c[i - 1] * b[p].addtag % mod;
		}
		for (int i = 10; i >= 1; i--) {
			op = 0;
			for (int j = i; j >= 0; j--) {
				op = op + ((x[i][j] * b[2 * p].num[j] % mod) * c[i - j]) % mod;
				op %= mod;
			}
			b[2 * p].num[i] = op;
		}
		for (int i = 10; i >= 1; i--) {
			op = 0;
			for (int j = i; j >= 0; j--) {
				op = op + ((x[i][j] * b[2 * p + 1].num[j]) % mod * c[i - j]) % mod;
				op %= mod;
			}
			b[2 * p + 1].num[i] = op;
		}
		b[2 * p].addtag += b[p].addtag;
		b[2 * p + 1].addtag += b[p].addtag;
		b[p].addtag = 0;
	}
}
inline void build(int p, int l, int r) {
	b[p].l = l;
	b[p].r = r;
	b[p].changetag = 1000000008;
	b[p].addtag;
	b[p].num[0] = r - l + 1;
	if (l == r) {
		b[p].num[1] = a[l];
		for (int i = 2; i <= 10; i++) {
			b[p].num[i] = b[p].num[i - 1] * a[l] % mod;
		}
		return;
	}
	build(2 * p, l, (l + r) >> 1);
	build(2 * p + 1, ((l + r) >> 1) + 1, r);
	update(p);
}
inline void add(int p, int l, int r, int w) {
	//	cout << p << " " << b[p].l << " " << b[p].r << ' ' << l << ' ' << r << ' ' << w << endl;
	if (b[p].l >= l && b[p].r <= r) {
		c[0] = 1;
		c[1] = w;
		for (int i = 2; i <= 10; i++) {
			c[i] = c[i - 1] * w % mod;
		}
		for (int i = 10; i >= 1; i--) {
			op = 0;
			for (int j = i; j >= 0; j--) {
				op = op + ((x[i][j] * b[p].num[j] % mod) * c[i - j]) % mod;
				op %= mod;
			}
			b[p].num[i] = op;
		}
		b[p].addtag += w;
		return;
	}
	push_down(p);
	int mid = (b[p].l + b[p].r) >> 1;
	if (r <= mid) add(2 * p, l, r, w);
	else if (l > mid) add(2 * p + 1, l, r, w);
	else {
		add(2 * p, l, mid, w);
		add(2 * p + 1, mid + 1, r, w);
	}
	update(p);
}
inline void change(int p, int l, int r, int w) {
	//	cout << p << ' ' << l << ' ' << r << "%%^" << b[p].l << ' ' << b[p].r << "-" << w << " " << b[p].num[0] << endl;
	if (b[p].l >= l && b[p].r <= r) {
		for (int i = 1; i <= 10; i++) {
			b[p].num[i] = b[p].num[i - 1] * w % mod;
		}
		b[p].changetag = w;
		b[p].addtag = 0;
		//		cout << b[p].num[0] << "*" << b[p].num[1] << "$" << b[p].num[2] << endl;
		return;
	}
	push_down(p);
	int mid = (b[p].l + b[p].r) >> 1;
	if (r <= mid) change(2 * p, l, r, w);
	else if (l > mid) change(2 * p + 1, l, r, w);
	else {
		change(2 * p, l, mid, w);
		change(2 * p + 1, mid + 1, r, w);
	}
	update(p);
}
inline int query(int p, int l, int r, int w) {
	//	cout << l << "+" << r << "=" << w << " " << b[p].l << " " << b[p].r << ' ' << b[p].num[w] << endl;
	if (b[p].l >= l && b[p].r <= r) {
		return b[p].num[w];
	}
	push_down(p);
	int mid = (b[p].l + b[p].r) >> 1;
	if (r <= mid) return query(2 * p, l, r, w);
	else if (l > mid) return query(2 * p + 1, l, r, w);
	else {
		return (query(2 * p, l, mid, w) + query(2 * p + 1, mid + 1, r, w)) % mod;
	}
}
signed main() {
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);
	cin >> n >> m;
	for (int i = 1; i <= n; i++) {
		cin >> a[i];
	}
	build(1, 1, n);
	jie[0] = 1;
	inv[0] = 1;
	for (int i = 1; i <= 10; i++) {
		jie[i] = jie[i - 1] * i;
		inv[i] = ksm(jie[i], mod - 2);
	}
	for (int i = 1; i <= 10; i++) {
		for (int j = 0; j <= i; j++) {
			x[i][j] = C(i, j);
		}
	}
	int op, p, q, w;
	for (int i = 1; i <= m; i++) {
		cin >> op >> p >> q >> w;
		if (op == 1) {
			add(1, p, q, w);
		} else if (op == 2) {
			change(1, p, q, w);
		} else {
			cout << query(1, p, q, w) << "\n";
		}
	}
}
//5 9
//1 1 2 4 6
//3 1 3 1
//2 1 3 1
//3 1 3 2
//2 1 2 3
//3 1 2 2
//2 3 5 1
//3 3 5 3
//2 5 5 3
//3 1 2 2