#include<bits/stdc++.h>
#define int long long
using namespace std;
struct node
{
	int l,r;
	int sum;
	int lcm;
}tr[1000010];
int a[1000010];
void push_up(int p)
{
	tr[p].sum=tr[p*2].sum+tr[p*2+1].sum;
	tr[p].lcm=min((__int128)2e18,(__int128)tr[p<<1].lcm*tr[p<<1|1].lcm/__gcd(tr[p<<1].lcm,tr[p<<1|1].lcm));
}
void build(int l,int r,int p)
{
	tr[p].l=l,tr[p].r=r;
	if(l==r)
	{
		tr[p].sum=tr[p].lcm=a[l];
		return;
	}
	int mid=(l+r)/2;
	build(l,mid,p*2);
	build(mid+1,r,p*2+1);
	push_up(p);
}
void update(int l,int r,int k,int p)
{
	if(tr[p].l>=l&&tr[p].r<=r)
	{
		tr[p].sum=tr[p].lcm=__gcd(tr[p].lcm,k);
		return;
	}
	int mid=(tr[p].l+tr[p].r)/2;
	if(mid>=l)
		update(l,r,k,p*2);
	if(mid<r)
		update(l,r,k,p*2+1);
	push_up(p);
}
int query(int l,int r,int p)
{
	if(tr[p].l>=l&&tr[p].r<=r)
		return tr[p].sum;
	int mid=(tr[p].l+tr[p].r)/2,sum=0;
	if(mid>=l)
		sum+=query(l,r,p*2);
	if(mid<r)
		sum+=query(l,r,p*2+1);
	return sum;
}
signed main()
{
	int n,m;
	cin>>n>>m;
	for(int i=1;i<=n;i++)
		cin>>a[i];
	build(1,n,1);
	while(m--)
	{
		int op,x,y;
		cin>>op>>x>>y;
		if(op==1)
		{
			int k;
			cin>>k;
			update(x,y,k,1);
		}
		else
		{
			cout<<query(x,y,1)<<endl;
		}
	}
	return 0;
}