敌兵布阵 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
这是一道典型的线段树（当然树状数组也是可以的）。因为据上次做线段树的题目的时间已经很长了，具体怎么做已经忘的差不多了，只隐约记得一些细节，所以今天做这道题的目的是为了复习。不找其他相关资料， 凭着记忆，把代码拼凑起来，因为我觉得只有这样印象才会深刻，会更好的理解，每写下一段代码，问问自己为什么这样做。
好，开始
首先我记得一共是有3个函数（目前就针对这道题而言）：
1.建树（buildUp)
2.查询(query)
3.修改点信息（change）

**线段树以二叉为模型建树
所有我们可以不严谨的说一下：
设一共n个点的信息
设m=log2（n）（取m为整数）
那么 1+2^1+^2……+2^m的和小于n*4**

static int n;
static int num[n<<2];//我们要的选段树各个节点信息存储在这

假设1到n的值都保存在 a数组

static int a[];

现在写buildUp函数

static void buildUp(int s,int L,int R)
{
    if(L==R)
    {
        num[s]=a[R];
        return;
    }
    int m=(L+R)>>1;
    buildUp(s<<1,L,m);//1
    buildUp(s<<1|1,m+1,R);//2
    num[s]=num[s<<1]+num[s<<1|1];//3
}

该怎么理解？递归函数也得递归的去理解，我们现在假设buildUp能完成这个操作，假设1 和2代码完成了，那么说明num[s<<1]和num[s<<1|1]已经完成填充，那么num[s]的值便能从上面得出，然后我们在做我们力所能及的，就是递归边界，终止函数的无限递归，写递归函数的感觉就像写一个递推关系（个人感觉）。关于 s与 R和L的关系，我们可以理解为
（R，L）是逻辑位置，s是物理位置，特定的（R，L）（为什么说特定呢，因为假设n=4，那么（1，2）有对应的s值，但是（1，3）就没有）对应一个不同的s值
理论上讲，是可以写出一个 s=f（R，L）函数，实际上 （R，L）对于那个或者多个s我们是不清楚的，除了一个，那就是（1，n）对应s=1，所以我们以这个为关键，加上递推关系，作为调用的入口值。
下面是change：

static void change(int p,int c,int s,int L,int R)//s是当前位置，
{
    if(L==R&&p==R)//说明已经到了要修改的叶子节点,可以结束了
    {
        num[s]+=c;
        return;
    }
    if(L<=m&&m<=R)//只有满足这个条件我才递归调用
    {
        int m=(L+R)>>1;
        num[s]+=c;
        change(p,c,L,m);
        change(p,c,m+1,R);//这些函数会自然遇到上面的if语句，如果不符合，就自然终止
    }

如果一个增加或减少某个值，那么只要包含这个数的区间都会改变相同的值。自顶向下去调用即可。
query：

static int query(int l,int r,int s,int L,int R)//r,l是我们所要求的区间的和的对应索引
{
    if(l<=L&&R<=r)
    {
        return num[s];R//为什么这样就直接return？，因为有一部分是为了下面的递归函数做准备，return的值可能被加到了ans
    }
    int m=(R+L)>>1;
    int ans=0;
    if(m>=L)//说明有有一部分值有可能留在做区间
        ans+=query(l.r,s<<1,L,m);
    if(m<R)//m+1<=R
        ans+=query(l,r,s<<1|1,m+1,R);
    return ans;

这道题的核心代码就写完了
完整版：

import java.util.Scanner;

public class Main 
{   
    static int num[];
    static int a[];
    static Scanner sc;
    public static void main(String[] args) 
    {
        sc=new Scanner(System.in);
        int t=sc.nextInt();
        int count=1;
        while((t--)>0)
        {
            System.out.println("Case "+count+":");
            int n=sc.nextInt();
            num=new int[n*4+1];
            /*a=new int[n+1];
            for(int i=1;i<=n;i++)
                a[i]=sc.nextInt();*/
            buildUp(1,1,n);
            int x,y;
            for(;;)
            {
                String s=sc.next();
                if(s.compareTo("End")==0)
                    break;
                else
                {
                    x=sc.nextInt();
                    y=sc.nextInt();
                    if(s.compareTo("Query")==0)
                    {
                        System.out.println(query(1,x,y,1,n));
                    }
                    else
                    {
                        if(s.compareTo("Sub")==0)
                            change(x,-y,1,1,n);
                        else
                            change(x,y,1,1,n);
                    }
                }
            }
            count++;
        }
    }
    static void buildUp(int s,int l,int r)
    {
        if(l==r)
        {
            num[s]=sc.nextInt();//想到，其实a数组是没有后续用处的，为了节省时间，想想为什么可以这样直接写？提示：下面的1，2行的代码顺序不能调换
            return;
        }
        int m=(l+r)>>1;
        buildUp(s<<1,l,m);//1
        buildUp(s<<1|1,m+1,r);//2
        num[s]=num[s<<1]+num[s<<1|1];
    }
    static void change(int p,int c,int s,int l,int r)
    {
        if(l==r&&r==p)
        {
            num[s]+=c;
            return;
        }
        if(l<=p&&p<=r)
        {
            num[s]+=c;
            int m=(l+r)>>1;
            change(p,c,s<<1,l,m);
            change(p,c,s<<1|1,m+1,r);       
        }
    }
    static int query(int s,int tl,int tr,int l,int r)
    {
        if(tl<=l&&r<=tr)
            return num[s];
        else
        {
            int m=(l+r)>>1;
            int ans=0;
            if(m>=tl)ans+=query(s<<1,tl,tr,l,m);
            if(m<tr)ans+=query(s<<1|1,tl,tr,m+1,r);
            return ans;
        }
    }
}

目前而言上述算是初级线段树，以后根据题目在详细展开，包括区间值修改，线段树的非递归实现等
今天更新一下非递归版本的线段树

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#define N 50005
#define INF 0x3f3f3f3f
typedef long long ll;
using namespace std;
int n;
int sum[N<<2];
int M;
void build()
{
    memset(sum,0,sizeof(sum));
    M=1;
    while(M<n+2)M<<=1;
    //cout<<"haha";
    for(int i=1;i<=n;i++)
        scanf("%d",&sum[M+i]);
        //cout<<"haha";
    for(int i=M-1;i>=1;i--)
        sum[i]=sum[i<<1]+sum[i<<1|1];
}
void change(int p,int c)
{
    for(int i=M+p;i;i>>=1)
        sum[i]+=c;
}
int query(int L,int R)
{
    L+=M-1;
    R+=M+1;
    int ans=0;
    for(;L^R^1;)
    {
        if(!(L&1))
            ans+=sum[L+1];
        if(R&1)
            ans+=sum[R-1];
        L>>=1;
        R>>=1;
    }
    return ans;
}
int main()
{
    int t;
    char s[10];
    int x,y;
    scanf("%d",&t);
    int cal=1;
    while(t--)
    {
        printf("Case %d:\n",cal++);
        scanf("%d",&n);
        build();
        //cout<<"haha";
        for(;;)
        {
            scanf("%s",s);
            if(s[0]=='E')
                break;
            scanf("%d%d",&x,&y);
            if(s[0]=='Q')
                 printf("%d\n",query(x,y));
            if(s[0]=='A')
                change(x,y);
            if(s[0]=='S')
                change(x,-y);
        }
    }
    return 0;
}