#include<bits/stdc++.h>
using namespace std;
struct tree{
	int l,r;
	long long x;
	long long sum;
}tr[800020];
const long long Mod=1ll<<32;
long long lcm(long long x,long long y){
	if(x==-1||y==-1) return -1;
	long long g=__gcd(x,y);
	if((__int128)(x/g*y)>(__int128)(1e18)) return -1;
	return x/g*y;
}
long long a[200005];
long long X;
int n,m;
void build(int l,int r,int x){
    if(l==r){
		tr[x].l=l;
		tr[x].r=r;
		tr[x].sum=a[l];
		tr[x].x=a[l];
		return;
	}
	int mid=(l+r)/2;
	build(l,mid,x<<1);
	build(mid+1,r,x<<1|1);
	tr[x].l=l;
	tr[x].r=r;
	tr[x].sum=tr[x<<1].sum+tr[x<<1|1].sum;
	tr[x].x=lcm(tr[x<<1].x,tr[x<<1|1].x);
	return;
}
long long fin(int l,int r,int x){
	if(tr[x].l>=l&&tr[x].r<=r) return tr[x].sum;
	int mid=(l+r)/2;
	long long sum=0;
	if(mid>=l) sum+=fin(l,r,x<<1)%Mod;
	if(mid<r) sum+=fin(l,r,x<<1|1)%Mod;
	return sum%Mod;
}
void change(int l,int r,int x){
	if(tr[x].l==tr[x].r){
		tr[x].sum=tr[x].x=__gcd(tr[x].x,X);
		return;
	}
	int mid=(l+r)/2;
	if(mid>=l){
		if(tr[x<<1].x==-1||X%tr[x<<1].x!=0) change(l,r,x<<1);
	}
	if(mid<r){
		if(tr[x<<1|1].x==-1||X%tr[x<<1|1].x!=0) change(l,r,x<<1|1);
	}
	tr[x].sum=tr[x<<1].sum+tr[x<<1|1].sum;
	tr[x].x=lcm(tr[x<<1|1].x,tr[x<<1].x);
	return;
}
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++){
		scanf("%lld",&a[i]);
	}
	build(1,n,1);
	for(int i=1;i<=m;i++){
		int op;
		scanf("%d",&op);
		if(op==1){
			int l,r;
			scanf("%d%d%lld",&l,&r,&X);
			change(l,r,1);
		}
		else{
			int l,r;
			scanf("%d%d",&l,&r);
			printf("%lld\n",fin(l,r,1)%Mod);
		}
	}
	return 0;
}