#include <bits/stdc++.h>
using namespace std;
const long long mod=1000000007;
int n,l,r,nq;
long long a,b,v[300010];
void rd(int &x)
{
	x=0;
	char c=getchar();
	for(;c>'9'||c<'0';c=getchar());
	for(;c<='9'&&c>='0';c=getchar())
		x=(x<<3)+(x<<1)+c-'0';
}
void rd(long long &x)
{
	x=0ll;
	char c=getchar();
	for(;c>'9'||c<'0';c=getchar());
	for(;c<='9'&&c>='0';c=getchar())
		x=(x<<3ll)+(x<<1ll)+c-'0';
}
void pt(long long x)
{
	if(x>=10ll)
		pt(x/10ll);
	putchar(x%10ll+'0');
}
struct matrix
{
	int n,m;
	long long a[4][4];
	void init(int k)
	{
		for(int i=1;i<=n;i++)
			for(int j=1;j<=m;j++)
				if(i!=j)
					a[i][j]=0;
				else
					a[i][j]=k;
	}
	friend matrix operator*(matrix a,matrix b)
	{
		matrix c;
		c.n=a.n,c.m=b.m;
		c.init(0);
		for(int i=1;i<=c.n;i++)
			for(int k=1;k<=b.n;k++)
				for(int j=1;j<=c.m;j++)
					(c.a[i][j]+=a.a[i][k]*b.a[k][j]%mod)%=mod;
		return c;
	}
	friend matrix operator+(matrix a,matrix b)
	{
		matrix c;
		c.n=a.n,c.m=a.m;
		c.init(0);
		for(int i=1;i<=c.n;i++)
			for(int j=1;j<=c.m;j++)
				c.a[i][j]=(a.a[i][j]+b.a[i][j])%mod;
		return c;
	}
	friend matrix operator!(matrix a)
	{
		matrix c;
		c.n=a.m,c.m=a.n;
		for(int i=1;i<=c.n;i++)
			for(int j=1;j<=c.m;j++)
				c.a[i][j]=a.a[j][i];
		return c;
	}
	friend matrix operator ^(matrix a,long long b)
	{
		matrix c;
		c.n=c.m=a.n,c.init(1);
		for(;b;b>>=1ll)
		{
			if(b&1ll)
				c=c*a;
			a=a*a;
		}
		return c;
	}
}ept,pl,mi,tpl,tmi,f,tf,dm;
long long ni(long long x)
{
	long long res=1ll,b=mod-2;
	for(;b;b>>=1ll)
	{
		if(b&1ll)
			res=res*x%mod;
		x=x*x%mod;
	}
	return res;
}
void make(matrix &a)
{
	a.n=3,a.m=3;
}
string op;
struct nd
{
	int l,r;
	matrix s,tl,tr;
}T[1200010];
void pushup(int x)
{
	T[x].s=T[x<<1].s+T[x<<1|1].s;
}
inline void pushdown(int x)
{
	T[x<<1].s=T[x].tl*T[x<<1].s*T[x].tr;
	T[x<<1|1].s=T[x].tl*T[x<<1|1].s*T[x].tr;
	T[x<<1].tl=T[x].tl*T[x<<1].tl;
	T[x<<1|1].tl=T[x].tl*T[x<<1|1].tl;
	T[x<<1].tr=T[x<<1].tr*T[x].tr;
	T[x<<1|1].tr=T[x<<1|1].tr*T[x].tr;
	T[x].tl.init(1),T[x].tr.init(1);
}
void build(int p,int l,int r)
{
	T[p].l=l,T[p].r=r;
	make(T[p].tl),make(T[p].tr);
	T[p].tl.init(1),T[p].tr.init(1);
	if(l==r)
	{
		T[p].s=(pl^(v[l-1]-1))*f;
		T[p].s=T[p].s*tf*(tpl^(v[l+1]-3));
		return;
	}
	int mid=l+r>>1;
	build(p<<1,l,mid);
	build(p<<1|1,mid+1,r);
	pushup(p);
}
void change(int p,int L,int R,matrix lv,matrix rv)
{
	int l=T[p].l,r=T[p].r;
	if(l>R||r<L)
		return;
	if(L<=l&&r<=R)
	{
		T[p].s=lv*T[p].s*rv;
		T[p].tl=lv*T[p].tl;
		T[p].tr=T[p].tr*rv;
		return;
	}
	pushdown(p);
	int mid=l+r>>1;
	if(L<=mid)
		change(p<<1,L,R,lv,rv);
	if(R>mid)
		change(p<<1|1,L,R,lv,rv);
	pushup(p);
}
matrix query(int p,int L,int R)
{
	int l=T[p].l,r=T[p].r;
	if(L<=l&&r<=R)
		return T[p].s;
	pushdown(p);
	int mid=l+r>>1;
	matrix res;
	make(res);
	res.init(0);
	if(R>mid)
		res=query(p<<1|1,L,R);
	if(L<=mid)
		res=res+query(p<<1,L,R);
	return res;
}
int main()
{
	rd(n),rd(nq),rd(a),rd(b);
	for(int i=1;i<=n;i++)
		rd(v[i]);
	make(pl),make(f);
	make(mi),make(ept),make(dm);
	ept.init(0),dm.init(1);
	pl.a[1][1]=pl.a[2][1]=pl.a[3][3]=1;
	pl.a[1][2]=a;
	pl.a[1][3]=b;
	f.a[1][1]=2;
	f.a[2][1]=f.a[3][1]=1;
	tpl=!pl;
	tf=!f;
	if(a)
	{
		mi.a[1][2]=mi.a[3][3]=1;
		mi.a[2][1]=ni(a);
		mi.a[2][2]=mod-ni(a);
		mi.a[2][3]=mod-b*ni(a)%mod;
	}
	else
	{
		mi.a[1][2]=mi.a[2][2]=mi.a[3][3]=1;
		mi.a[2][3]=mod-b;
	}
	tmi=!mi;
	build(1,2,n-1);
	while(nq--)
	{
		cin>>op>>l>>r;
		if(op=="plus")
			change(1,l+1,r+1,pl,dm),change(1,l-1,r-1,dm,tpl);
		else if(op=="minus")
			change(1,l+1,r+1,mi,dm),change(1,l-1,r-1,dm,tmi); 
		else if(op=="query")
			pt(query(1,l+1,r-1).a[1][1]),putchar('\n');
	}
	return 0;
}