Problem: At the children's day, the child came to Picks's house, and messed his house up. Picks was angry at him. A lot of important things were lost, in particular the favorite sequence of Picks.

Fortunately, Picks remembers how to repair the sequence. Initially he should create an integer array a[1], a[2], ..., a[n]. Then he should perform a sequence of m operations. An operation can be one of the following:

  1. Print operation l, r. Picks should write down the value of <image>. 
  2. Modulo operation l, r, x. Picks should perform assignment a[i] = a[i] mod x for each i (l ≤ i ≤ r). 
  3. Set operation k, x. Picks should set the value of a[k] to x (in other words perform an assignment a[k] = x). 



Can you help Picks to perform the whole sequence of operations?

Input

The first line of input contains two integer: n, m (1 ≤ n, m ≤ 105). The second line contains n integers, separated by space: a[1], a[2], ..., a[n] (1 ≤ a[i] ≤ 109) — initial value of array elements.

Each of the next m lines begins with a number type <image>. 

  * If type = 1, there will be two integers more in the line: l, r (1 ≤ l ≤ r ≤ n), which correspond the operation 1. 
  * If type = 2, there will be three integers more in the line: l, r, x (1 ≤ l ≤ r ≤ n; 1 ≤ x ≤ 109), which correspond the operation 2. 
  * If type = 3, there will be two integers more in the line: k, x (1 ≤ k ≤ n; 1 ≤ x ≤ 109), which correspond the operation 3. 

Output

For each operation 1, please print a line containing the answer. Notice that the answer may exceed the 32-bit integer.

Examples

Input

5 5
1 2 3 4 5
2 3 5 4
3 3 5
1 2 5
2 1 3 3
1 1 3


Output

8
5


Input

10 10
6 9 6 7 6 1 10 10 9 5
1 3 9
2 7 10 9
2 5 10 8
1 4 7
3 3 7
2 7 9 9
1 2 4
1 6 6
1 5 9
3 1 10


Output

49
15
23
1
9

Note

Consider the first testcase:

  * At first, a = {1, 2, 3, 4, 5}. 
  * After operation 1, a = {1, 2, 3, 0, 1}. 
  * After operation 2, a = {1, 2, 5, 0, 1}. 
  * At operation 3, 2 + 5 + 0 + 1 = 8. 
  * After operation 4, a = {1, 2, 2, 0, 1}. 
  * At operation 5, 1 + 2 + 2 = 5. 

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
template <typename T>
T BigMod(T b, T p, T m) {
  if (p == 0) return 1;
  if (p % 2 == 0) {
    T s = BigMod(b, p / 2, m);
    return ((s % m) * (s % m)) % m;
  }
  return ((b % m) * (BigMod(b, p - 1, m) % m)) % m;
}
template <typename T>
T ModInv(T b, T m) {
  return BigMod(b, m - 2, m);
}
template <typename T>
T in() {
  char ch;
  T n = 0;
  bool ng = false;
  while (1) {
    ch = getchar();
    if (ch == '-') {
      ng = true;
      ch = getchar();
      break;
    }
    if (ch >= '0' && ch <= '9') break;
  }
  while (1) {
    n = n * 10 + (ch - '0');
    ch = getchar();
    if (ch < '0' || ch > '9') break;
  }
  return (ng ? -n : n);
}
template <typename T>
T POW(T B, T printf) {
  if (printf == 0) return 1;
  if (printf & 1)
    return B * POW(B, printf - 1);
  else
    return (POW(B, printf / 2) * POW(B, printf / 2));
}
template <typename T>
T Bigmod(T b, T p, T m) {
  if (p == 0)
    return 1;
  else if (!(p & 1))
    return (Bigmod(b, p / 2, m) * Bigmod(b, p / 2, m)) % m;
  else
    return ((b % m) * Bigmod(b, p - 1, m)) % m;
}
template <typename T>
T Dis(T x1, T y1, T x2, T y2) {
  return sqrt((x1 - x2 * x1 - x2) + (y1 - y2 * y1 - y2));
}
template <typename T>
T Angle(T x1, T y1, T x2, T y2) {
  return atan(double(y1 - y2) / double(x1 - x2));
}
template <typename T>
T DIFF(T a, T b) {
  T d = a - b;
  if (d < 0)
    return -d;
  else
    return d;
}
template <typename T>
T ABS(T a) {
  if (a < 0)
    return -a;
  else
    return a;
}
template <typename T>
T gcd(T a, T b) {
  if (a < 0) return gcd(-a, b);
  if (b < 0) return gcd(a, -b);
  return (b == 0) ? a : gcd(b, a % b);
}
template <typename T>
T lcm(T a, T b) {
  if (a < 0) return lcm(-a, b);
  if (b < 0) return lcm(a, -b);
  return a * (b / gcd(a, b));
}
template <typename T>
T euclide(T a, T b, T &x, T &y) {
  if (a < 0) {
    T d = euclide(-a, b, x, y);
    x = -x;
    return d;
  }
  if (b < 0) {
    T d = euclide(a, -b, x, y);
    y = -y;
    return d;
  }
  if (b == 0) {
    x = 1;
    y = 0;
    return a;
  } else {
    T d = euclide(b, a % b, x, y);
    T t = x;
    x = y;
    y = t - (a / b) * y;
    return d;
  }
}
template <typename T>
void ia(T a[], int n) {
  for (int i = 0; i < n; i++) cin >> a[i];
}
template <typename T>
void pa(T a[], int n) {
  for (int i = 0; i < n - 1; i++) cout << a[i] << " ";
  cout << a[n - 1] << endl;
}
template <typename T>
long long int isLeft(T a, T b, T c) {
  return (a.x - b.x) * (b.y - c.y) - (b.x - c.x) * (a.y - b.y);
}
int Set(int N, int pos) { return N = N | (1 << pos); }
int Reset(int N, int pos) { return N = N & ~(1 << pos); }
bool Check(int N, int pos) { return (bool)(N & (1 << pos)); }
template <class T, class first>
inline T togglebit(T a, first i) {
  T t = 1;
  return (a ^ (t << i));
}
int toInt(string s) {
  int sm;
  stringstream ss(s);
  ss >> sm;
  return sm;
}
int toLlint(string s) {
  long long int sm;
  stringstream ss(s);
  ss >> sm;
  return sm;
}
int cdigittoint(char ch) { return ch - '0'; }
bool isVowel(char ch) {
  ch = toupper(ch);
  if (ch == 'A' || ch == 'U' || ch == 'I' || ch == 'O' || ch == 'E')
    return true;
  return false;
}
bool isConst(char ch) {
  if (isalpha(ch) && !isVowel(ch)) return true;
  return false;
}
double DEG(double x) { return (180.0 * x) / (2.0 * acos(0.0)); }
double RAD(double x) { return (x * (double)2.0 * acos(0.0)) / (180.0); }
long long int Tree[4 * 100007][2];
void Build(int idx, int l, int r) {
  if (l == r) {
    scanf("%d", &Tree[idx][0]);
    Tree[idx][1] = Tree[idx][0];
    return;
  }
  int mid = (l + r) >> 1, lft = idx << 1, rgt = lft | 1;
  Build(lft, l, mid);
  Build(rgt, mid + 1, r);
  Tree[idx][0] = Tree[lft][0] + Tree[rgt][0];
  Tree[idx][1] = max(Tree[lft][1], Tree[rgt][1]);
}
void Update(int idx, int l, int r, int s, int e, int v) {
  if (Tree[idx][1] < v) return;
  if (l == r) {
    Tree[idx][0] %= v;
    Tree[idx][1] %= v;
    return;
  }
  int mid = (l + r) >> 1, lft = idx << 1, rgt = lft | 1;
  if (e <= mid)
    Update(lft, l, mid, s, e, v);
  else if (s > mid)
    Update(rgt, mid + 1, r, s, e, v);
  else
    Update(lft, l, mid, s, mid, v), Update(rgt, mid + 1, r, mid + 1, e, v);
  Tree[idx][0] = Tree[lft][0] + Tree[rgt][0];
  Tree[idx][1] = max(Tree[lft][1], Tree[rgt][1]);
}
void Modify(int idx, int l, int r, int pos, int v) {
  if (l == r) {
    Tree[idx][0] = Tree[idx][1] = v;
    return;
  }
  int mid = (l + r) >> 1, lft = idx << 1, rgt = lft | 1;
  if (pos <= mid)
    Modify(lft, l, mid, pos, v);
  else
    Modify(rgt, mid + 1, r, pos, v);
  Tree[idx][0] = Tree[lft][0] + Tree[rgt][0];
  Tree[idx][1] = max(Tree[lft][1], Tree[rgt][1]);
}
long long int Query(int idx, int l, int r, int s, int e) {
  if (l == s && r == e) return Tree[idx][0];
  int mid = (l + r) >> 1, lft = idx << 1, rgt = lft | 1;
  if (e <= mid)
    return Query(lft, l, mid, s, e);
  else if (s > mid)
    return Query(rgt, mid + 1, r, s, e);
  else
    return Query(lft, l, mid, s, mid) + Query(rgt, mid + 1, r, mid + 1, e);
}
int main() {
  int tc, cas = 1;
  int n, m;
  scanf("%d%d", &n, &m);
  Build(1, 1, n);
  while (m--) {
    int(cs);
    scanf("%d", &cs);
    if (cs == 1) {
      int x, y;
      scanf("%d%d", &x, &y);
      printf("%lld\n", Query(1, 1, n, x, y));
    }
    if (cs == 2) {
      int x, y, z;
      scanf("%d%d%d", &x, &y, &z);
      Update(1, 1, n, x, y, z);
    }
    if (cs == 3) {
      int x, y;
      scanf("%d%d", &x, &y);
      Modify(1, 1, n, x, y);
    }
  }
  return 0;
}