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,这样的，而是有着间隔的，其实也就是建造了许多树状数组

 

 1 #include<stdio.h>

 2 

 3 #include<queue>

 4 #include<iostream>

 5 #include<algorithm>

 6 #include<string.h>

 7 const int MAXN=50020;

 8 

 9 

10 int c[12][12][MAXN];

11 int n;

12 

13 int lowbit(int x)

14 {

15     return x&(-x);

16 }

17 

18 void update(int t1,int t2,int i,int val)

19 {

20     while(i<=n)

21     {

22         c[t1][t2][i]+=val;

23         i+=lowbit(i);

24     }

25 

26 }

27 int sum(int t1,int t2,int i)

28 {

29 

30     int s=0;

31     while(i>0)

32     {

33         s+=c[t1][t2][i];

34         i-=lowbit(i);

35     }

36     return s;

37 }

38 int num[MAXN];

39 int main()

40 {

41     #ifndef ONLINE_JUDGE

42     freopen("1.in","r",stdin);

43     #endif

44     int m;

45     while(scanf("%d",&n)!=EOF)

46     {

47         for(int i=0;i<n;i++)scanf("%d",&num[i]);

48         for(int i=0;i<12;i++)

49           for(int j=0;j<12;j++)

50             for(int k=0;k<MAXN;k++)

51               c[i][j][k]=0;

52         scanf("%d",&m);

53         int a,b,k,q;

54         int t;

55         while(m--)

56         {

57             scanf("%d",&t);

58             if(t==1)

59             {

60                 scanf("%d%d%d%d",&a,&b,&k,&q);

61                 a--;

62                 b--;

63                 int num=(b-a)/k;    

64                 int s=a%k-1;

65                 update(k,s,a/k+1,q);    //s是该数组的开头数字，a/k+1是结尾的数字

66                 update(k,s,a/k+num+2,-q);

67             }

68             else

69             {

70                 scanf("%d",&a);

71                 a--;

72                 int ss=num[a];

73                 for(int i=1;i<=10;i++)

74                 {

75                     ss+=sum(i,a%i-1,a/i+1);

76                 }

77                 printf("%d\n",ss);

78             }

79         }

80     }

81     return 0;

82 }