HDU-1166 敌兵布阵

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

#include
#include
#include
int a[50010];
int b[50010];
char op[120];
//取lowbit 
int lowbit(int x)
{
	return x&(-x);
}
//加 
void Add(int x,int c,int N)
{
	for(int i=x;i<=N;i+=lowbit(i))
	{
		b[i]+=c;
	}
}
//减 
void Sub(int x,int c,int N)
{
	for(int i=x;i<=N;i+=lowbit(i))
	{
		b[i]-=c;
	}
}
//求区间的和
int Q(int x)
{
	int ans=0; 
	for(int i=x;i>0;i-=lowbit(i))
	{
		ans+=b[i];	
	}	
	return ans;
} 
int main()
{
	int T;
	scanf("%d",&T);
	for(int i=1;i<=T;i++)
	{
		memset(b,0,sizeof(b));
		int N;
		scanf("%d",&N);
		for(int j=1;j<=N;j++)scanf("%d",&a[j]);
		for(int k=1;k<=N;k++)Add(k,a[k],N);
		printf("Case %d:\n",i);
		while(1)
		{
			int p,c;
			scanf("%s",op);
			if(op[0]=='E')break;
			else if(op[0]=='A')
			{
				scanf("%d %d",&p,&c);
				Add(p,c,N);
			}
			else if(op[0]=='S')
			{
				scanf("%d %d",&p,&c);
				Sub(p,c,N);
			}
			else if(op[0]=='Q')
			{
				scanf("%d %d",&p,&c);
				printf("%d\n",Q(c)-Q(p-1));
			}
		}
	}
	return 0;
} 

下面是线段树

#include
#include
#include
#define maxn 50010
int A[maxn];
int sum[maxn<<2];
char op[120];
//父节点等于子节点相加 
void Sumup(int rt){sum[rt]=sum[rt<<1]+sum[rt<<1|1];}
//初始化线段树 
void Build(int l,int r,int rt)
{
	if(l==r)
	{
		sum[rt]=A[l];
		return;
	}
	int mid=(l+r)>>1;
	Build(l,mid,rt<<1);
	Build(mid+1,r,rt<<1|1);
	Sumup(rt);
}
//点修改 
void Add(int L,int C,int l,int r,int rt)
{
	if(l==r)
	{
		sum[rt]+=C;
		return;
	}
	int mid=(l+r)>>1;
	if(L<=mid)Add(L,C,l,mid,rt<<1);
	else Add(L,C,mid+1,r,rt<<1|1);
	Sumup(rt);
}
//求区间
int Q(int L,int R,int l,int r,int rt)
{
	if(L<=l&&R>=r)
	{
		return sum[rt];	
	}	
	int ans=0;
	int mid=(l+r)>>1;
	if(mid>=L)ans+=Q(L,R,l,mid,rt<<1);
	if(mid<R)ans+=Q(L,R,mid+1,r,rt<<1|1);
	return ans;
} 
int main()
{
	int T,i;
	scanf("%d",&T);
	for(i=1;i<=T;i++)
	{
		memset(A,0,sizeof(A));
		memset(sum,0,sizeof(sum));
		int N;
		scanf("%d",&N);
		for(int j=1;j<=N;j++)
		{
			scanf("%d",&A[j]);
		}
		Build(1,N,1);
		printf("Case %d:\n",i);
		while(1)
		{
			int p,c;
			scanf("%s",op);
			if(op[0]=='E')break;
			else if(op[0]=='A')
			{
				scanf("%d %d",&p,&c);
				Add(p,c,1,N,1);
			}
			else if(op[0]=='S')
			{
				scanf("%d %d",&p,&c);
				Add(p,-c,1,N,1);
			}
			else if(op[0]=='Q')
			{
				scanf("%d %d",&p,&c);
				printf("%d\n",Q(p,c,1,N,1));
			}
		}
	}
	return 0;
}