#include <bits/stdc++.h>
using namespace std;
int n, q;
int a[100048];
long long c[12][12][100048];
void update(int x, int y, int pos, int delta) {
  while (pos <= n) {
    c[x][y][pos] += delta;
    pos += pos & (-pos);
  }
}
long long query(int x, int y, int pos) {
  long long res = 0;
  while (pos > 0) {
    res += c[x][y][pos];
    pos -= pos & (-pos);
  }
  return res;
}
int main() {
  int i, j, m, t, x, y, z, x1, y1, d, p;
  long long ans, coef;
  scanf("%d", &n);
  for (i = 1; i <= n; i++) scanf("%d", &a[i]);
  for (m = 2; m <= 10; m += 2)
    for (i = 1; i <= n; i++) update(m, i % m, i, a[i]);
  scanf("%d", &q);
  for (i = 1; i <= q; i++) {
    scanf("%d%d%d", &t, &x, &y);
    if (t == 1) {
      d = y - a[x];
      a[x] = y;
      for (m = 2; m <= 10; m += 2) update(m, x % m, x, d);
    } else {
      scanf("%d", &z);
      ans = 0;
      for (j = 0; j <= 2 * (z - 1) - 1; j++) {
        coef = (j == 0) ? 2 : min(j, 2 * z - j);
        ans += coef * (query(2 * (z - 1), (x + j - 1) % (2 * (z - 1)), y) -
                       query(2 * (z - 1), (x + j - 1) % (2 * (z - 1)), x - 1));
      }
      printf("%I64d\n", ans);
    }
  }
  return 0;
}
