#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e5 + 5;
struct Seg {
  long long sum, ma;
};
Seg seg[maxn << 2];
void solve();
int main() {
  solve();
  return 0;
}
void pushup(int rt) {
  seg[rt].sum = seg[rt << 1].sum + seg[rt << 1 | 1].sum;
  seg[rt].ma = max(seg[rt << 1].ma, seg[rt << 1 | 1].ma);
}
void build(int l, int r, int rt) {
  if (l == r) {
    scanf("%I64d", &seg[rt].ma);
    seg[rt].sum = seg[rt].ma;
    return;
  }
  int m = (l + r) >> 1;
  build(l, m, rt << 1);
  build(m + 1, r, rt << 1 | 1);
  pushup(rt);
}
void update1(int p, int add, int l, int r, int rt) {
  if (l == r) {
    seg[rt].ma = seg[rt].sum = add + 0LL;
    return;
  }
  int m = (l + r) >> 1;
  if (p <= m)
    update1(p, add, l, m, rt << 1);
  else
    update1(p, add, m + 1, r, rt << 1 | 1);
  pushup(rt);
}
void update2(int L, int R, int add, int l, int r, int rt) {
  if (seg[rt].ma < add) return;
  if (l == r) {
    seg[rt].ma %= add;
    seg[rt].sum %= add;
    return;
  }
  int m = (l + r) >> 1;
  if (L <= m) update2(L, R, add, l, m, rt << 1);
  if (R > m) update2(L, R, add, m + 1, r, rt << 1 | 1);
  pushup(rt);
}
long long query(int L, int R, int l, int r, int rt) {
  if (L <= l && R >= r) return seg[rt].sum;
  int m = (l + r) >> 1;
  long long ret = 0;
  if (L <= m) ret += query(L, R, l, m, rt << 1);
  if (R > m) ret += query(L, R, m + 1, r, rt << 1 | 1);
  return ret;
}
void solve() {
  int n, q;
  scanf("%d %d", &n, &q);
  build(1, n, 1);
  for (int i = 1; i <= q; i++) {
    int op;
    scanf("%d", &op);
    if (op == 1) {
      int l, r;
      scanf("%d %d", &l, &r);
      printf("%I64d\n", query(l, r, 1, n, 1));
    } else if (op == 2) {
      int l, r, x;
      scanf("%d %d %d", &l, &r, &x);
      update2(l, r, x, 1, n, 1);
    } else {
      int k, x;
      scanf("%d %d", &k, &x);
      update1(k, x, 1, n, 1);
    }
  }
  return;
}
