HDU 4267 A Simple Problem with Integers

树状数组。2种操作，1 a b k c ：表示把a到b区间内所有下标为i且满足(i - a) % k == 0的值加c，也就是说，下标从a开始，每次加k，直到超过b。2 a：表示询问a点的值。

当k等于1时，就是成段更新，单点查询。想到这里不妨把数字分组，对k取余相同的数进一组，那么对每组来说，都是成段更新。单点查询时只需把a在10个组内的值都加起来，再加上原来的值，就是当前a点的值了。

然而，把数字分组不是很好搞啊，算着算着就乱了。就想把所有数都放进去，然后反正有的数也用不到（余数不同），不碍事。

#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<vector>
#include<queue>
#include<cmath>
///LOOP
#define REP(i, n) for(int i = 0; i < n; i++)
#define FF(i, a, b) for(int i = a; i < b; i++)
#define FFF(i, a, b) for(int i = a; i <= b; i++)
#define FD(i, a, b) for(int i = a - 1; i >= b; i--)
#define FDD(i, a, b) for(int i = a; i >= b; i--)
///INPUT
#define RI(n) scanf("%d", &n)
#define RII(n, m) scanf("%d%d", &n, &m)
#define RIII(n, m, k) scanf("%d%d%d", &n, &m, &k)
#define RIV(n, m, k, p) scanf("%d%d%d%d", &n, &m, &k, &p)
#define RV(n, m, k, p, q) scanf("%d%d%d%d%d", &n, &m, &k, &p, &q)
#define RFI(n) scanf("%lf", &n)
#define RFII(n, m) scanf("%lf%lf", &n, &m)
#define RFIII(n, m, k) scanf("%lf%lf%lf", &n, &m, &k)
#define RFIV(n, m, k, p) scanf("%lf%lf%lf%lf", &n, &m, &k, &p)
#define RS(s) scanf("%s", s)
///OUTPUT
#define PN printf("\n")
#define PI(n) printf("%d\n", n)
#define PIS(n) printf("%d ", n)
#define PS(s) printf("%s\n", s)
#define PSS(s) printf("%s ", n)
#define PC(n) printf("Case %d: ", n)
///OTHER
#define PB(x) push_back(x)
#define CLR(a, b) memset(a, b, sizeof(a))
#define CPY(a, b) memcpy(a, b, sizeof(b))
#define display(A, n, m) {REP(i, n){REP(j, m)PIS(A[i][j]);PN;}}
#define lson l, m, rt << 1
#define rson m + 1, r, rt << 1 | 1

using namespace std;
typedef long long LL;
typedef pair<int, int> P;
const int MOD = 9901;
const int INFI = 1e9 * 2;
const LL LINFI = 1e17;
const double eps = 1e-6;
const double pi = acos(-1.0);
const int N = 50005;
const int M = 11;
const int move[8][2] = {0, 1, 0, -1, 1, 0, -1, 0, 1, 1, 1, -1, -1, 1, -1, -1};

int c[M][M][N], num[N];

void add(int x, int y, int z, int p)
{
    while(z < N)
    {
        c[x][y][z] += p;
        z += z & (-z);
    }
}

int sum(int x, int y, int z)
{
    int s = 0;
    while(z)
    {
        s += c[x][y][z];
        z -= z & (-z);
    }
    return s;
}

int main()
{
    //freopen("input.txt", "r", stdin);

    int n, m, a, b, d, k, op, s;
    while(RI(n) != EOF)
    {
        FFF(i, 1, n)RI(num[i]);
        CLR(c, 0);
        RI(m);
        while(m--)
        {
            RI(op);
            if(op == 1)
            {
                RIV(a, b, k, d);
                add(k, a % k, a, d);
                add(k, a % k, b + 1, -d);
            }
            else
            {
                RI(a);
                s = num[a];
                FFF(i, 1, 10)s += sum(i, a % i, a);
                PI(s);
            }
        }
    }
    return 0;
}