HDU-1166 敌兵布阵（线段树）

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

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

分析：由于本题，数据量大且查询次数多，直接查肯定超时，这个时候就需要利用线段树进行快速查找了。利用线段树可以很方便快速的在存储数据和查找。

对于该题我们分为以下几个步骤：

步骤一：

       建立线段树，如图

        

       上图建立的是一个区间为1~10的线段树，我们输入n个阵地时，利用二分的思想建立一个区间范围（1~n）的线段树，每个根节点的区间范围（m~m），在建立过程中，不断递归最先建立最左边的根节点，即第一个阵地，输入第一个阵地敌人数量后，紧接着建立第二个根节点并输入数据，直到第n个节点为止，所有阵地数据输入完毕。

void build(int l,int r,int rt)
{
	if(l==r)//如果递归到根节点，输入数据
	{
		scanf("%d",&st[rt]);
		return;
	}
	int m=(l+r)/2;//将大区间分为两个小区间，直到区间分为l==r为止
	build(l,m,rt*2);//对左子树进行建树
	build(m+1,r,rt*2+1);//对右子树进行建树
	st[rt]=st[rt*2]+st[rt*2+1];//对父亲节点进行求和
}

步骤二：

       查询线段树

       我们在已知的查询区间内进行递归，最先归到查询区间最左边根节点，查询其值，返回到上一级。再进行右子树查询，返回上一级，最终查询到查询区间最右边的根节点，结束返回该区间所有节点内的求和。

int query(int L,int R,int l,int r,int rt)
{
	if(L<=l&&r<=R)//如果查询区间在查询的范围内 
	{
		return st[rt];//返回该部分区间节点的和 
	}
	int m=(l+r)/2;//将大查询区间分为两个小查询区间 
	int res=0;
	if(L<=m)//如果左查询区间有范围需要查询 
	{
		res+=query(L,R,l,m,rt*2);//递归左查询区间 
	}
	if(m<R)//如果右查询区间有范围需要查询 
	{
		res+=query(L,R,m+1,r,rt*2+1);//递归右查询区间 
	}
	return res;//返回最终值 
}

步骤三：

       修改线段树指定节点

       步骤三与步骤二基本一样，只要最终递归到指定节点，修改其值即可。

void update(int pos,int add,int l,int r,int rt)//pos指需要修改的节点，add指需要修改的值 
{
	if(l==r)//找到指定节点，修改其值 
	{
		st[rt]+=add;
		return;
	}
	int m=(l+r)/2;//将大查询区间分为两个小查询区间 
	if(pos<=m)//如果要修改的数在左区间，递归左区间 
	{
		update(pos,add,l,m,rt*2);
	}
	else//否则递归右区间 
	{
		update(pos,add,m+1,r,rt*2+1);
	}
	st[rt]=st[rt*2]+st[rt*2+1];//由于指定根节点的数已经变动，所以上一级节点的数也要变动 
}

最终代码：

#include<stdio.h>
const int M=1e6+5;
int st[M];
void build(int l,int r,int rt)
{
	if(l==r)
	{
		scanf("%d",&st[rt]);
		return;
	}
	int m=(l+r)/2;
	build(l,m,rt*2);
	build(m+1,r,rt*2+1);
	st[rt]=st[rt*2]+st[rt*2+1];
}
int query(int L,int R,int l,int r,int rt)
{
	if(L<=l&&r<=R)
	{
		return st[rt];
	}
	int m=(l+r)/2;
	int res=0;
	if(L<=m)
	{
		res+=query(L,R,l,m,rt*2);
	}
	if(m<R)
	{
		res+=query(L,R,m+1,r,rt*2+1);
	}
	return res;
}
void update(int pos,int add,int l,int r,int rt)
{
	if(l==r)
	{
		st[rt]+=add;
		return;
	}
	int m=(l+r)/2;
	if(pos<=m)
	{
		update(pos,add,l,m,rt*2);
	}
	else
	{
		update(pos,add,m+1,r,rt*2+1);
	}
	st[rt]=st[rt*2]+st[rt*2+1];
}
int main()
{
	int t,n,l,r,sum=0;
	char type[10];
	scanf("%d",&t);
	for(int cnt=1;cnt<=t;cnt++)
	{
		scanf("%d",&n);
		build(1,n,1);
		printf("Case %d:\n",cnt);
		while(~scanf("%s",type))
		{
			if(type[0]=='E')
			{
				break;
			}
			scanf("%d%d",&l,&r);
			if(type[0]=='Q')
			{
				printf("%d\n",query(l,r,1,n,1));
			}
			else if(type[0]=='A')
			{
				update(l,r,1,n,1);
			}
			else if(type[0]=='S')
			{
				update(l,-r,1,n,1);
			}
		}
	}
}