根节点存储的永远是总区间,孩子节点每次都把其父亲节点分成两部分,父亲节点表示的区间为[L,R], 其左儿子表示的区间为[L,(L+R)/2],右儿 子表示的区间为[(L+R)/2 + 1,R]。
用线段树统计的东西,必须符合区间加法:
数字之和——总数字之和 = 左区间数字之和 + 右区间数字之和;
最大公因数(GCD)——总GCD = gcd( 左区间GCD , 右区间GCD );
最大值or最小值——总最大值=max(左区间最大值,右区间最大值);
线段树适用于:单点修改,区间查询 区间修改,区间查询
维护区间最值
//维护区间最小值
void push_up(int i)
{
tr[i].minn = min(tr[i<<1].minn,tr[i<<1|1].minn);
}
//维护区间最大值
void push_up(int i)
{
tr[i].minn = max(tr[i<<1].minn,tr[i<<1|1].minn);
}
敌兵布阵 HDU - 1166 (单点修改,区间查询)
https://vjudge.net/problem/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<bits/stdc++.h>//万能头文件
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
using namespace std;
typedef long long ll;
const int MAX=5e4+5;
int arr[MAX];
struct tree
{
int L,R;//[L,R]此节点代表的区间
int sum;//区间和
}tr[MAX<<2];
void push_up(int i)
{
tr[i].sum=tr[i<<1].sum+tr[i<<1|1].sum;
}
void build(int l,int r,int i)//建树
{
tr[i].L=l;
tr[i].R=r;
if(l==r)
{
tr[i].sum=arr[l];
return ;
}
int mid=(l+r)>>1;
build(l,mid,i<<1);//i<<1==i*2 左孩子
build(mid+1,r,i<<1|1);//i<<1|1==i*2+1 右孩子
push_up(i);
}
void update(int pos,int val,int i)//将arr[pos] += val i为寻找到当前区间线段树的下标
{
if(tr[i].L==pos&&tr[i].R==pos)
{
tr[i].sum+=val;
return ;
}
int mid=(tr[i].L+tr[i].R)>>1;
if(pos>mid)
update(pos,val,i<<1|1);
else
update(pos,val,i<<1);
push_up(i);
}
ll query(int l,int r,int i)// 查询[l,r] 上的区间和
{
if(l<=tr[i].L&&r>=tr[i].R)
return tr[i].sum;
int mid=(tr[i].L+tr[i].R)>>1;
if(l>mid)
return query(l,r,i<<1|1);// 完全在右孩子区间
if(r<=mid)
return query(l,r,i<<1);// 完全在左孩子区间
return query(l,r,i<<1)+query(l,r,i<<1|1);
}
int main()
{
ios_base::sync_with_stdio(false), cin.tie(0);// 输入输出加速
int t,n,i,j,a,b;
char s[20];
cin>>t;
for(i=1;i<=t;i++)
{
cout<<"Case "<<i<<":"<<endl;
cin>>n;
for(j=1;j<=n;j++)
cin>>arr[j];
build(1,n,1);
while(cin>>s&&s[0]!='E')
{
cin>>a>>b;
if(s[0]=='Q')
cout<<query(a,b,1)<<endl;
if(s[0]=='A')
update(a,b,1);
if(s[0]=='S')
update(a,-b,1);
}
}
return 0;
}
A Simple Problem with IntegersPOJ - 3468 (区间修改,区间查询)
https://cn.vjudge.net/problem/POJ-3468
You have N integers, A1, A2, ... , AN. You need to deal with two kinds of operations. One type of operation is to add some given number to each number in a given interval. The other is to ask for the sum of numbers in a given interval.
Input
The first line contains two numbers N and Q. 1 ≤ N,Q ≤ 100000.
The second line contains N numbers, the initial values of A1, A2, ... , AN. -1000000000 ≤ Ai ≤ 1000000000.
Each of the next Q lines represents an operation.
"C a b c" means adding c to each of Aa, Aa+1, ... , Ab. -10000 ≤ c ≤ 10000.
"Q a b" means querying the sum of Aa, Aa+1, ... , Ab.
Output
You need to answer all Q commands in order. One answer in a line.
Sample Input
10 5
1 2 3 4 5 6 7 8 9 10
Q 4 4
Q 1 10
Q 2 4
C 3 6 3
Q 2 4
Sample Output
4
55
9
15
Hint
The sums may exceed the range of 32-bit integers.
-
区间修改
-
递归查找区间
-
发现现在区间完全属于需修改的区间,更新当前节点的信息,增加懒惰标记。递归结束
-
若不完全属于则继续递归
-
-
询问区间
-
发现当前区间完全属于被查询区间,直接返回当前节点值
-
或者发现当前区间不完全属于被查询区间且当前节点有懒惰标记,将懒惰标记下放至叶节点,并更新叶节点信息。取消自身的懒惰标记。
-
继续递归
-
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
const int MAX=1e5+5;
int n,m;
struct tree
{
tree()
{
sum=0;
lazy=0;
}
int L,R;//此节点代表的区间
ll sum,lazy;//区间和 懒惰标记
}tr[MAX<<2];
void push_up(int i)//单点求值
{
tr[i].sum=tr[i<<1].sum+tr[i<<1|1].sum;
}
void push_down(int i,int len)
{
if(tr[i].lazy)//将懒惰标记下放至叶节点,并更新叶节点信息。取消自身的懒惰标记
{
tr[i<<1|1].sum+=tr[i].lazy*(ll)(len>>1);
tr[i<<1].sum+=tr[i].lazy*(ll)(len-(len>>1));
tr[i<<1].lazy+=tr[i].lazy;
tr[i<<1|1].lazy+=tr[i].lazy;
tr[i].lazy=0;
}
}
void build(int l,int r,int i)//递归建树
{
tr[i].L=l;
tr[i].R=r;
if(l==r)
{
cin>>tr[i].sum;
return ;
}
int mid=(l+r)>>1;
build(l,mid,i<<1);
build(mid+1,r,i<<1|1);
push_up(i);
}
void update(int l,int r,int val,int i)//区间更新
{
if(l<=tr[i].L&&r>=tr[i].R)
{
tr[i].lazy+=(ll)val;
tr[i].sum+=(tr[i].R-tr[i].L+1)*(ll)val;
return ;
}
push_down(i,tr[i].R-tr[i].L+1);
int mid=(tr[i].L+tr[i].R)>>1;
if(l<=mid)
update(l,r,val,i<<1);
if(r>mid)
update(l,r,val,i<<1|1);
push_up(i);
}
ll query(int l,int r,int i)//区间求和
{
if(l<=tr[i].L&&r>=tr[i].R)
{
return tr[i].sum;
}
push_down(i,tr[i].R-tr[i].L+1);
int mid=(tr[i].R+tr[i].L)>>1;
if(l>mid)
return query(l,r,i<<1|1);
if(r<=mid)
return query(l,r,i<<1);
return query(l,r,i<<1)+query(l,r,i<<1|1);
}
int main()
{
ios_base::sync_with_stdio(false),cin.tie(0);
cin>>n>>m;
build(1,n,1);
while(m--)
{
char s;int a,b,c;
cin>>s>>a>>b;
if(s=='C')
{
cin>>c;
update(a,b,c,1);
}
else
cout<<query(a,b,1)<<endl;
}
return 0;
}