hdu 4267 多维树状数组

题意：有一个序列

"1 a b k c" means adding c to each of Ai which satisfies a <= i <= b and (i - a) % k == 0. (1 <= a <= b <= N, 1 <= k <= 10, -1,000 <= c <= 1,000)
"2 a" means querying the value of Aa. (1 <= a <= N)

分析转自：http://www.cnblogs.com/Griselda/archive/2012/09/15/2686133.html

感觉是树状数组了

但是 树状数组是单点更新 而题目要求是在a b 内每隔 k 个数更新..

如果用for循环对a b 内每隔 k 个位置的数更新就会很费时间..

所以有一个方法就是 按 k 值 i%k 值还有 i/k+1 值把每隔 k 个位置的值都分组..即按 k 值和 i%k 值建树状数组 

到时候维护其中一棵 树状数组就好了..

 

维护的方法是在c[ k ][a%k][(i-a)/k+1] + 1 在c[ k ][a % k ][ (a/k) + (b-a)/k + 2 ] - 1..

 

关于建成的树状数组..

更新时有55种情况
1,2,3,4,5......
1,3,5,7,9......
2,4,6,8,10....
1,4,7,10,13...
2,5,9,12,15...
3,6,10,13,16...
.
.
.
10,20,30,40,50...

所以用55个树状数组.

建出来的树并不是普通的1,2,3,4,这样的，而是有着间隔的，其实也就是建造了许多树状数组

#include<stdio.h>

#include<queue>
#include<iostream>
#include<algorithm>
#include<string.h>
const int MAXN=50020;


int c[12][12][MAXN];
int n;

int lowbit(int x)
{
    return x&(-x);
}

void update(int t1,int t2,int i,int val)
{
    while(i<=n)
    {
        c[t1][t2][i]+=val;
        i+=lowbit(i);
    }

}
int sum(int t1,int t2,int i)
{

    int s=0;
    while(i>0)
    {
        s+=c[t1][t2][i];
        i-=lowbit(i);
    }
    return s;
}
int num[MAXN];
int main()
{
    #ifndef ONLINE_JUDGE
    freopen("1.in","r",stdin);
    #endif
    int m;
    while(scanf("%d",&n)!=EOF)
    {
        for(int i=0;i<n;i++)scanf("%d",&num[i]);
        for(int i=0;i<12;i++)
          for(int j=0;j<12;j++)
            for(int k=0;k<MAXN;k++)
              c[i][j][k]=0;
        scanf("%d",&m);
        int a,b,k,q;
        int t;
        while(m--)
        {
            scanf("%d",&t);
            if(t==1)
            {
                scanf("%d%d%d%d",&a,&b,&k,&q);
                a--;
                b--;
                int num=(b-a)/k;    
                int s=a%k-1;
                update(k,s,a/k+1,q);    //s是该数组的开头数字，a/k+1是结尾的数字
                update(k,s,a/k+num+2,-q);
            }
            else
            {
                scanf("%d",&a);
                a--;
                int ss=num[a];
                for(int i=1;i<=10;i++)
                {
                    ss+=sum(i,a%i-1,a/i+1);
                }
                printf("%d
",ss);
            }
        }
    }
    return 0;
}