HDU1166 敌兵布阵（线段树、树状数组）

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. 题目对数组有三种操作：①查询[x,y]的总和。②某个数组元素tag加上addVal。③某个数组元素tag减上addVal。
2. 暴力求解：修改原数组元素，时间复杂度为o(1)。每次查询[x,y]都从x加到y，时间复杂度为o(n)。复杂度过高会超时。
3. 题目目的查询和修改数组的子数组的元素总和。为了节省查询时间，可以使用线段树或者树状数组来保存[x,y]的结果。

---

解题代码：

树状数组：

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

int num[50050], sum[50050], t, n;
int lowbit(int x) {
    return x&(-x);
}
int summ(int x) {
    int res = 0;
    while(x > 0) {
        res += sum[x];
        x -= lowbit(x);
    }
    return res;
}
int getSum(int L, int R) {
    return summ(R)-summ(L-1);
}
void update(int tag, int addVal) {
    while(tag <= n) {
        sum[tag] += addVal;
        tag += lowbit(tag);
    }
}
int main() {

    int a, b;
    char op[20];
    scanf("%d", &t);
    for(int Case = 1; Case <= t; Case++) {
        printf("Case %d:\n", Case);
        scanf("%d", &n);
        memset(num, 0, sizeof(num));
        memset(sum, 0, sizeof(sum));
        for(int i = 1; i <= n; i++) {
            scanf("%d", &num[i]);
            update(i, num[i]);
        }
        while(scanf("%s", op)) {
            if(!strcmp(op, "End"))
                break;
            scanf("%d %d", &a, &b);
            if(!strcmp(op, "Query")) {
                printf("%d\n", getSum(a, b));
            } else if(!strcmp(op, "Add")) {
                num[a] += b;
                update(a, b);
            } else if(!strcmp(op, "Sub")) {
                num[a] -= b;
                update(a, -b);
            }
        }
    }
    return 0;
}

线段树： 

#include <cstdio>
#include <cstring>
using namespace std;
typedef struct Node {
    int val, addMark, L, R;
    Node(){}
    Node(int L, int R, int val, int addMark) {
        this->L = L;
        this->R = R;
        this->val = val;
        this->addMark = addMark;
    }
}Node;

Node tree[200000];
int num[50050], t, n;

void build(int root, int L, int R) {
    if(L == R) {
        tree[root] = Node(L, R, num[L], 0);
        return;
    }
    int mid = (L+R)/2;
    build(root*2, L, mid);
    build(root*2+1, mid+1, R);
    tree[root] = Node(L, R, tree[root*2].val + tree[root*2+1].val, 0);
}
void updateOne(int root,int tag, int addVal) {
    tree[root].val += addVal;
    if(tree[root].L == tree[root].R) {
        return;
    }
    if(tag <= tree[root*2].R)
        updateOne(root*2, tag, addVal);
    else
        updateOne(root*2+1, tag, addVal);
}
int query(int root, int tL, int tR) {
    if(tL > tree[root].R || tR < tree[root].L)
        return 0;
    if(tL <= tree[root].L && tree[root].R <= tR)
        return tree[root].val;
    return query(root*2, tL, tR) + query(root*2+1, tL, tR);
}
int main() {

    int a, b;
    char op[20];
    scanf("%d", &t);
    for(int Case = 1; Case <= t; Case++) {
        printf("Case %d:\n", Case);
        scanf("%d", &n);
        for(int i = 1; i <= n; i++)
            scanf("%d", &num[i]);
        build(1, 1, n);
        while(scanf("%s", op)) {
            if(!strcmp(op, "End"))
                break;
            scanf("%d %d", &a, &b);
            if(!strcmp(op, "Query")) {
                printf("%d\n", query(1, a, b));
            } else if(!strcmp(op, "Add")) {
                updateOne(1, a, b);
            } else if(!strcmp(op, "Sub")) {
                updateOne(1, a, -b);
            }
        }
    }
    return 0;
}