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種形式: www.2cto.com
(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
這題是zkw線段樹中單點修改與區間查詢的一般方法:
M的計算 顯然最後一行是:
+M: 0 1 2 3....M-1
但是我們求和值是不能求(0,M-1),(zkw線段樹不支持首尾查詢)
我們查詢的應該是開區間(a,b) (這樣的目的是保證不出現當前無法判定是否取-
原理:在左邊時左邊不取,而若右邊不為t(s^t≠1),則必取.
在右邊時右邊不取,而若左邊不為s(s^t≠1),則必取.
建樹時要保證n<=M-2
[cpp]
M=1;while (M-2<n) M<<=1;
從C++的角度:
由於char s[]為數組,無法直接比較
需要用 strcmp(string a,string a)==0 判定
它是字符串的cmp 返回字典序比較情況
[cpp]
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<cctype>
#include<iostream>
using namespace std;
#define MAXN (50000+10)
#define MAXAi (50+10)
int n,m,k,M,tt,t[MAXN*5];
char s[50];
int main()
{
// freopen("Hdu1166.in","r",stdin);
scanf("%d",&tt);
for (int k=1;k<=tt;k++)
{
printf("Case %d:\n",k);
scanf("%d",&n);
memset(t,0,sizeof(t));
M=1;while (M-2<n) M<<=1;
for (int i=M+1;i<=M+n;i++) scanf("%d",&t[i]);
for (int i=M;i>=1;i--) t[i]=t[i<<1]+t[(i<<1)^1];
while (scanf("%s",s)!=EOF)
{
if (!strcmp(s,"End")) break;
int p1,p2;
scanf("%d%d",&p1,&p2);
if (s[0]=='Q')
{
int ans=0;
p1--;p2++;p1+=M;p2+=M;
while (p1^p2^1>0)
{
if (~p1&1) ans+=t[p1+1];
if (p2&1) ans+=t[p2-1];
p1>>=1;p2>>=1;
}
cout<<ans<<endl;
continue;
}
else
{
p1+=M;
if (s[0]=='A') t[p1]+=p2; else t[p1]-=p2;
for (p1/=2;p1;p1/=2)
{
t[p1]=t[p1<<1]+t[(p1<<1)^1];
}
}
}
}
return 0;
}