#include <bits/stdc++.h>
using namespace std;
void tempwj() {
  freopen("P3380_2.in", "r", stdin);
  freopen("hash.out", "w", stdout);
}
long long gcd(long long a, long long b) { return b == 0 ? a : gcd(b, a % b); }
long long qpow(long long a, long long b, long long mod) {
  a %= mod;
  long long ans = 1;
  while (b > 0) {
    if (b & 1) ans = ans * a % mod;
    a = a * a % mod;
    b >>= 1;
  }
  return ans % mod;
}
struct cmp {
  bool operator()(const pair<int, long long>& a,
                  const pair<int, long long>& b) {
    return a.second > b.second;
  }
};
int lb(int x) { return x & -x; }
const long long inf = 0x3f3f3f;
const double esp = 1e-9;
const long long INF = 0x3f3f3f3f3f3f3f;
const long long mod = 1e9 + 7;
const int maxn = 1e6 + 10;
struct Node {
  int l, r, num;
  long long sum;
} node[maxn << 2];
int a[maxn];
void build(int l, int r, int no) {
  node[no].l = l;
  node[no].r = r;
  if (l == r) {
    node[no].num = a[l];
    node[no].sum = a[l];
    return;
  }
  int mid = l + r >> 1;
  build(l, mid, no << 1);
  build(mid + 1, r, no << 1 | 1);
  node[no].num = max(node[no << 1].num, node[no << 1 | 1].num);
  node[no].sum = node[no << 1].sum + node[no << 1 | 1].sum;
}
void update(int l, int r, int no, int num) {
  if (node[no].l > r || node[no].r < l) return;
  if (node[no].num < num) return;
  if (node[no].l == node[no].r) {
    node[no].num %= num;
    node[no].sum %= num;
    return;
  }
  update(l, r, no << 1, num);
  update(l, r, no << 1 | 1, num);
  node[no].num = max(node[no << 1].num, node[no << 1 | 1].num);
  node[no].sum = node[no << 1].sum + node[no << 1 | 1].sum;
}
void update2(int pos, int no, int num) {
  if (node[no].l > pos || node[no].r < pos) return;
  if (node[no].l == node[no].r) {
    node[no].num = num;
    node[no].sum = num;
    return;
  }
  update2(pos, no << 1, num);
  update2(pos, no << 1 | 1, num);
  node[no].num = max(node[no << 1].num, node[no << 1 | 1].num);
  node[no].sum = node[no << 1].sum + node[no << 1 | 1].sum;
}
long long query(int l, int r, int no) {
  if (node[no].l > r || node[no].r < l) return 0;
  if (node[no].l >= l && node[no].r <= r) return node[no].sum;
  return query(l, r, no << 1) + query(l, r, no << 1 | 1);
}
int main() {
  int n, m;
  scanf("%d%d", &n, &m);
  for (int i = 1; i <= n; i++) scanf("%d", &a[i]);
  build(1, n, 1);
  while (m--) {
    int f, l, r;
    scanf("%d%d%d", &f, &l, &r);
    if (f == 1) {
      long long ans = query(l, r, 1);
      printf("%lld\n", ans);
    } else if (f == 2) {
      int x;
      scanf("%d", &x);
      update(l, r, 1, x);
    } else {
      update2(l, 1, r);
    }
  }
}
