线段树

线段树可以解决的问题：

1、区间查询

--询问某段区间的某些性质（极值、求和、etc）

2、区间更新

--某些操作影响了某段区间（统一加一个值...）

3、三个问题

--更新点，区间查询

--更新区间，查询点

--更新区间，查询区间

问题：一个长度为N的一维数组（a[1]~a[N]），我们每次对该数组有一些操作：

·1、修改数组中某个元素的值

-----【1，5，4，1，6】-----（a[2]=3）----->【1，3，4，1，6】

·2、询问数组中某段区间的最大值

-----【1，5，4，1，6】-----（max(1,4)=？）----->5

·3、询问数组中某段区间的和

-----【1，5，4，1，6】-----（sum(3,5)=?）----->11

如果只有一次询问？

-----枚举相应区间内的元素，输出答案O(N)

-----更多询问？若是q次询问，则时间复杂度为O(qN)，too slow

这里的主要问题是询问是针对区间的，而我们维护的信息是针对每个元素的；线段树可以做到的是在O(logn)的时间内完成每次操作，则时间复杂度是O(qlogn)。

开始线段树讲解：

线段树的本质是一棵二叉树，不同于其他二叉树，线段树的每一个节点记录的是一段区间的信息。

e.g. 对于长度为5的数组a[1]~a[5]

对于任一非叶子节点，若该区间为[L,R]，则左儿子为[L,(L+R)/2]，右儿子为[(L+R)/2+1,R]

如何建树？如何记录信息？

每一个节点记录的信息？

const int N = 100005;
struct Segtree{
    int l, r;//这个节点对应区间的端点
    int mx, sum;//维护的信息，视具体题目而定
    int lazy;//区间操作
}tr[N * 4];

如何记录左儿子和右儿子

#define lid (id << 1)
#define rid (id << 1 | 1)

我们用一个数组a记录节点，且根节点的下标为1，对于任一节点a[k]，它的左儿子为a[2*k]，它的右儿子为a[2*k+1]。

建立一棵线段树，并记录原数组信息

int a[N];
void push_up(int id){
    tr[id].sum = tr[rid].sum + tr[lid].sum;
    tr[id].mx = max(tr[lid], tr[rid]);
}

void build(int id, int l, int r){
    tr[id].l = l;
    tr[id].r = r;
    tr[id].lazy = 0;//区间查询时候的lazy操作
    if(l == r){
        tr[id].sum = tr[id].mx = a[l];
        return ;
    }
    int mid = (l + r) >> 1;
    build(lid, l, mid);
    build(rid, mid + 1, r);
    push_up(id);
}

如果原数组从a[1]~a[N]，调用build(1,1,n)即可。

更新某个区间，并维护相关点信息，需要引进lazy的思想。

void push_down(int id){
    if(tr[id].lazy){
        tr[lid].lazy += tr[id].lazy;
        tr[rid].lazy += tr[id].lazy;
        tr[lid].sum += tr[id].lazy * (tr[lid].r - tr[lid].l + 1);
        tr[rid].sum += tr[id].lazy * (tr[rid].r - tr[rid].l + 1);
        tr[id].lazy = 0;
    }
}

void update(int id, int l, int r, int z){
    if(tr[id].l == l && tr[id].r == r){
        tr[id].lazy += z;
        tr[id].sum += z * (tr[id].r - tr[id].l + 1);
        return ;
    }
    push_down(id);
    int mid = (tr[id].l + tr[id].r) >> 1;
    if(mid >= r)
        update(lid, l, r, z);
    else if(mid < l)
        update(rid, l, r, z);
    else{
        update(lid, l, mid, z);
        update(rid, mid + 1, r, z);
    }
    push_up(id);
}

更新某个点的数值，并维护相关点的信息

void update(int id, int x, int v){
    if(tr[id].l == tr[id].r){
        tr[id].sum = tr[id].mx = v;
        return ;
    }
    int mid = (tr[id].l + tr[id].r) >> 1;
    update(x <= mid ? lid : rid, x, v);
    push_up(id);
}

若更新a[k]的值为val，调用update(1, k, val)即可。

查询某区间内元素的和或最大值（以总和为例）

int query(int id, int l, int r){
    if(tr[id].l == l && tr[id].r == r)
        return tr[id].sum;
    push_down(id);//区间查询的时候需要更新子节点
    int mid = (tr[id].l + tr[id].r) >> 1;
    if(r <= mid)
        return query(lid, l, r);
    if(l > mid)
        return query(rid, l, r);
    return query(lid, l, mid) + query(rid, mid + 1, r);
}

调用query(1,l,r)即可查询[l, r]区间内元素的总和。

 

时间复杂度：

更新操作：由于总是一条路径从跟到某个叶子，而树的深度为logn，因而为O(logn)

查询操作：每层被访问的节点不超过4个，因而同样为O(logn)

小结：

1、线段树可以做很多很多与区间有关的事情

2、空间复杂度为O(N*4)，每次更新和查询操作的复杂度都是O(logn)

3、在更新和查询区间[l, r]的时候，为了保证复杂度是严格的O(logn)必须在达到被[l, r]覆盖的区间的节点时就立即返回。而为了保证这样做的正确的，需要在这两个过程中做一些相关的“懒”操作

题目：

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

 

裸的线段树：单点更新，区间求和

代码：

#include<bits/stdc++.h>
using namespace std;

#define lid (id << 1)
#define rid (id << 1 | 1)
const int N = 50005;

int a[N];

struct Segtree{
    int l, r;
    int sum, mx;
}tr[N * 4];

void push_up(int id){
    tr[id].sum = tr[rid].sum + tr[lid].sum;
    //tr[id].mx = max(tr[rid].mx, tr[lid].mx);
}
void build(int id, int l, int r){
    tr[id].l = l;
    tr[id].r = r;
    if(l == r){
        tr[id].sum  = a[l];
        return;
    }
    int mid = (l + r) >> 1;
    build(lid, l, mid);
    build(rid, mid + 1, r);
    push_up(id);
}
void update(int id, int x, int v){
    if(tr[id].l == tr[id].r){
        tr[id].sum += v;
        return ;
    }
    int mid = (tr[id].l + tr[id].r) >> 1;
    update(x <= mid ? lid : rid, x, v);
    push_up(id);
}

int query(int id, int l, int r){
    if(tr[id].l == l && tr[id].r == r)
        return tr[id].sum;
    int mid = (tr[id].l + tr[id].r) >> 1;
    if(r <= mid)
        return query(lid, l, r);
    if(l > mid)
        return query(rid, l, r);
    return query(lid, l, mid) + query(rid, mid + 1, r);
}
int main(){
    int T;
    scanf("%d", &T);
    int nCase = 1;
    while(T--){
        int n;
        scanf("%d", &n);
        for(int i = 1; i <= n; i++)
            scanf("%d", &a[i]);
        build(1, 1, n);
        printf("Case %d:\n", nCase++);
        while(1){
            char s[10];
            scanf("%s", s);
            if(s[0] == 'Q'){
                int x, y;
                scanf("%d%d", &x, &y);
                printf("%d\n", query(1, x, y));
            }else if(s[0] == 'A'){
                int pos, v;
                scanf("%d%d", &pos, &v);
                update(1, pos, v);
            }else if(s[0] == 'S'){
                int poos, vv;
                scanf("%d%d", &poos, &vv);
                update(1, poos, -vv);
            }else
                break;
        }
    }
    return 0;
}