1. 程式人生 > >樹狀陣列(例題:HDU1166)

樹狀陣列(例題:HDU1166)

樹狀陣列(Binary Indexed Tree(B.I.T), Fenwick Tree)是一個查詢和修改複雜度都為log(n)的資料結構。主要用於查詢任意兩位之間的所有元素之和,但是每次只能修改一個元素的值。

修改和求和的時間複雜度:logn

設節點編號為x,那麼這個節點管轄的區間為2^k(其中k為x二進位制末尾0的個數)個元素。

例如:

C1(0001) = A1                                                                 k=0,管轄區域1

C2(0010) = A1 + A2                                                        k=1,管轄區域2

C3(0011) = A3                                                                 k=0,管轄區域1

C4(0100) = A1 + A2 + A3 + A4                                        k=1,管轄區域2

C5(0101) = A5                                                                 k=0,管轄區域1

C6(0110) = A5 + A6                                                         k=1,管轄區域2

C7(0111) = A7                                                                  k=0,管轄區域1

C8 (1000)= A1 + A2 + A3 + A4 + A5 + A6 + A7 + A8      k=3,管轄區域8

求和:

step1: 令sum = 0,轉第二步;

step2: 假如n <= 0,演算法結束,返回sum值,否則sum = sum + Cn,轉第三步;

step3: 令n = n – lowbit(n),轉第二步。

lowbit操作:求二進位制末尾的1的位置,與自身負數進行與操作

//求二進位制末尾的1的位置
int lowbit(int x)
{
	return x&-x;
}

update操作:

//更新第k位置+x,如果是減就傳入-x
void update(int k,int x)
{
    //修改所有管轄到k的位置,直到大於所有的個數
	while(k<=n)
	{
		tree[k]+=x;
        //尋找下一個位置
		k+=lowbit(k);
	}
}

Query操作:這裡是求和

//其實是update的逆操作
int query(int k)
{
	int res=0;
	while(k)
	{
		res+=tree[k];
		k-=lowbit(k);
	}
	return res;
}

例題:HDU1166

敵兵佈陣

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 128070    Accepted Submission(s): 53664

 

Problem Description

C國的死對頭A國這段時間正在進行軍事演習,所以C國間諜頭子Derek和他手下Tidy又開始忙乎了。A國在海岸線沿直線佈置了N個工兵營地,Derek和Tidy的任務就是要監視這些工兵營地的活動情況。由於採取了某種先進的監測手段,所以每個工兵營地的人數C國都掌握的一清二楚,每個工兵營地的人數都有可能發生變動,可能增加或減少若干人手,但這些都逃不過C國的監視。
中央情報局要研究敵人究竟演習什麼戰術,所以Tidy要隨時向Derek彙報某一段連續的工兵營地一共有多少人,例如Derek問:“Tidy,馬上彙報第3個營地到第10個營地共有多少人!”Tidy就要馬上開始計算這一段的總人數並彙報。但敵兵營地的人數經常變動,而Derek每次詢問的段都不一樣,所以Tidy不得不每次都一個一個營地的去數,很快就精疲力盡了,Derek對Tidy的計算速度越來越不滿:"你個死肥仔,算得這麼慢,我炒你魷魚!”Tidy想:“你自己來算算看,這可真是一項累人的工作!我恨不得你炒我魷魚呢!”無奈之下,Tidy只好打電話向計算機專家Windbreaker求救,Windbreaker說:“死肥仔,叫你平時做多點acm題和看多點演算法書,現在嚐到苦果了吧!”Tidy說:"我知錯了。。。"但Windbreaker已經掛掉電話了。Tidy很苦惱,這麼算他真的會崩潰的,聰明的讀者,你能寫個程式幫他完成這項工作嗎?不過如果你的程式效率不夠高的話,Tidy還是會受到Derek的責罵的.

Input

第一行一個整數T,表示有T組資料。
每組資料第一行一個正整數N(N<=50000),表示敵人有N個工兵營地,接下來有N個正整數,第i個正整數ai代表第i個工兵營地裡開始時有ai個人(1<=ai<=50)。
接下來每行有一條命令,命令有4種形式:
(1) Add i j,i和j為正整數,表示第i個營地增加j個人(j不超過30)
(2)Sub i j ,i和j為正整數,表示第i個營地減少j個人(j不超過30);
(3)Query i j ,i和j為正整數,i<=j,表示詢問第i到第j個營地的總人數;
(4)End 表示結束,這條命令在每組資料最後出現;
每組資料最多有40000條命令

Output

對第i組資料,首先輸出“Case i:”和回車,
對於每個Query詢問,輸出一個整數並回車,表示詢問的段中的總人數,這個數保持在int以內。

Sample Input

1
10
1 2 3 4 5 6 7 8 9 10
Query 1 3
Add 3 6
Query 2 7
Sub 10 2
Add 6 3
Query 3 10
End 

Sample Output

Case 1:
6
33
59

Source

題意:很顯然,直接暴力會超時,這個用線段樹也可以,這裡用樹狀陣列來實現。

程式碼:

#include <cstdio>
#include <iostream>
#include <cstring>
#include <string>
using namespace std;

int tree[50005];
int T,n,x,y;
string s;

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

//更新
void update(int k,int x)
{
	while(k<=n)
	{
		tree[k]+=x;
		k+=lowbit(k);
	}
}

//求和
int query(int k)
{
	int res=0;
	while(k)
	{
		res+=tree[k];
		k-=lowbit(k);
	}
	return res;
}

int main()
{
	scanf("%d",&T);
	for (int cnt=1;cnt<=T;cnt++)
	{
		printf("Case %d:\n",cnt);
		scanf("%d",&n);
		memset(tree,0,sizeof(tree));
		for (int i=1;i<=n;i++)
		{
			scanf("%d",&x);
			update(i,x);
		}
		while(cin >> s && s!="End")
		{
			if (s=="Query")
			{
				scanf("%d%d",&x,&y);
				printf("%d\n",query(y)-query(x-1));
			}
			else if (s=="Add")
			{
				scanf("%d%d",&x,&y);
				update(x,y);
			}
			else
			{
				scanf("%d%d",&x,&y);
				update(x,-y);
			}				
		}
	}
	return 0;
}