#include <bits/stdc++.h>
using namespace std;
const long long maxn = 100005;
long long n, m;
long long c[maxn];
struct sg_tr {
  long long l, r, sum, mx;
};
sg_tr a[maxn << 2];
inline long long grd() {
  long long x = 0, f = 1;
  char ch = getchar();
  while (ch < '0' || ch > '9') {
    if (ch == '-') f = -1;
    ch = getchar();
  }
  while (ch >= '0' && ch <= '9') {
    x = (x << 1) + (x << 3) + (ch ^ 48);
    ch = getchar();
  }
  return x * f;
}
inline void upd(long long p) {
  a[p].sum = a[(p << 1)].sum + a[((p << 1) | 1)].sum;
  a[p].mx = std::max(a[(p << 1)].mx, a[((p << 1) | 1)].mx);
}
void bld(long long p, long long l, long long r) {
  a[p].l = l, a[p].r = r;
  if (l == r) {
    a[p].sum = a[p].mx = c[l];
    return;
  }
  long long mid = (l + r) >> 1;
  bld((p << 1), l, mid);
  bld(((p << 1) | 1), mid + 1, r);
  upd(p);
}
void mdf(long long p, long long ql, long long qr, long long mod) {
  if (a[p].mx < mod) return;
  if (a[p].l == a[p].r) {
    a[p].sum %= mod;
    a[p].mx %= mod;
    return;
  }
  long long mid = (a[p].l + a[p].r) >> 1;
  if (qr <= mid)
    mdf((p << 1), ql, qr, mod);
  else if (ql > mid)
    mdf(((p << 1) | 1), ql, qr, mod);
  else
    mdf((p << 1), ql, mid, mod), mdf(((p << 1) | 1), mid + 1, qr, mod);
  upd(p);
}
void replc(long long p, long long pos, long long val) {
  if (a[p].l == a[p].r) {
    a[p].sum = a[p].mx = val;
    return;
  }
  long long mid = (a[p].l + a[p].r) >> 1;
  if (pos <= mid)
    replc((p << 1), pos, val);
  else
    replc(((p << 1) | 1), pos, val);
  upd(p);
}
long long qry(long long p, long long ql, long long qr) {
  if (a[p].l == ql && a[p].r == qr) return a[p].sum;
  long long mid = (a[p].l + a[p].r) >> 1;
  if (qr <= mid)
    return qry((p << 1), ql, qr);
  else if (ql > mid)
    return qry(((p << 1) | 1), ql, qr);
  return qry((p << 1), ql, mid) + qry(((p << 1) | 1), mid + 1, qr);
}
int main() {
  n = grd(), m = grd();
  for (register long long i = 1; i <= n; i++) c[i] = grd();
  bld(1, 1, n);
  while (m--) {
    long long opt = grd(), x = grd(), y = grd();
    switch (opt) {
      case 1:
        printf("%lld\n", qry(1, x, y));
        break;
      case 3:
        replc(1, x, y);
        break;
      case 2:
        long long z = grd();
        mdf(1, x, y, z);
        break;
    }
  }
  return 0;
}
