#include<bits/stdc++.h>
#define maxn 100010
#define int long long
using namespace std;
struct node
{
	int left,right;
	int sum,mx;
}t[maxn*4];
int a[maxn],n,m;
void update(int id)
{
	t[id].sum=t[id<<1].sum+t[id<<1|1].sum;
	t[id].mx=max(t[id<<1].mx,t[id<<1|1].mx);
}
void build(int id,int l,int r)
{
	t[id].left=l,t[id].right=r;
	if(l==r)
	{
		t[id].sum=t[id].mx=a[l];
		return;
	}
	int mid=(l+r)/2;
	build(id<<1,l,mid);
	build(id<<1|1,mid+1,r);
	update(id);
}
void change(int id,int l,int r,int x,int d)
{
	if(l==r)
	{
		t[id].sum=t[id].mx=d;
		return;
	}
	int mid=(l+r)/2;
	if(x<=mid) change(id<<1,l,mid,x,d);
	else change(id<<1|1,mid+1,r,x,d);
	update(id);
}
void change2(int id,int l,int r,int x,int y,int d)
{
	if(t[id].mx<d) return;
	if(l==r)
	{
		t[id].sum%=d;
		t[id].mx%=d;
		return;
	}
	int mid=(l+r)/2;
	if(x<=mid) change2(id<<1,l,mid,x,y,d);
	if(y>mid) change2(id<<1|1,mid+1,r,x,y,d);
	update(id);
}
int query(int id,int l,int r)
{
	if(t[id].left==l&&t[id].right==r)
	{
		return t[id].sum;
	}
	if(r<=t[id<<1].right) return query(id<<1,l,r);
	else if(l>=t[id<<1|1].left) return query(id<<1|1,l,r);
	else 
	{
		return query(id<<1,l,t[id<<1].right)+query(id<<1|1,t[id<<1|1].left,r);
	}
	update(id);
}
signed main()
{
	scanf("%lld%lld",&n,&m);
	for(int i=1;i<=n;i++) scanf("%lld",&a[i]);
	build(1,1,n);
	for(int i=1;i<=m;i++)
	{
		int opt,x,y,z;
		scanf("%lld%lld%lld",&opt,&x,&y);
		if(opt==1)
		{
			printf("%d\n",query(1,x,y));
		}
		else if(opt==2)
		{
			scanf("%lld",&z);
			change2(1,1,n,x,y,z);
		}
		else 
		{
			change(1,1,n,x,y);
		}
	}
	return 0;
}

感觉写的没问题