考虑两种暴力解法：

1. 预处理 模 \(i\) 为 \(j\) 的下标，其中的元素之和；时间复杂度：\(O(n^2) + O(m)\)；
2. 暴力求 每次询问都遍历 \(\mathrm ki + j \space (\mathrm k \in \mathbb Z^+)\)；时间复杂度：\(O(mn)\)。

考虑优化，即将两种算法合并：

1. 模数 \(< \sqrt n\)：使用方法 \((1)\)，时间复杂度：\(O(n \sqrt n) + O(m)\)；
2. 模数 \(> \sqrt n\)：使用方法 \((2)\)，时间复杂度：\(O(m \sqrt n)\)；

因此，优化后的总时间复杂度为 \(O((n + m) \sqrt n)\)。代码如下：

#include <bits/stdc++.h>

#define ur uread()

using namespace std;

unsigned uread() {
    unsigned num = 0; char ch = getchar();
    while(!isdigit(ch)) ch = getchar();
    while(isdigit(ch)) num = num * 10 + ch - '0', ch = getchar();
    return num;
}

const int N = 1.5e5 + 10;
const int M = 390;

int arr[N];
int f[M][M];

signed main() {
    int n = ur, m = ur; int b = sqrt(n);
    for (int i = 1; i <= n; ++i) {
        arr[i] = ur;
        for (int j = 1; j < b; ++j) f[j][i % j] += arr[i];
    } while (m--) {
        char op[2]; scanf("%s", op);
        int x = ur, y = ur; if (op[0] == 'C') {
            for (int i = 1; i < b; ++i) f[i][x % i] += y - arr[x];
            arr[x] = y;
        } else if (x < b) {
            printf("%d\n", f[x][y]);
        } else {
            int sum = 0; for (int i = y; i <= n; i += x) {
                sum += arr[i];
            } printf("%d\n", sum);
        }
    }
    return 0;
}