上个学期看了很多复杂的数据结构,但是由于之前对简单数据结构的了解不够深入,所以,没办法真正了解复杂数据结构的内涵,最近几天又看了一下线段树,总结一下线段树单点更新。
相信很多初学线段树的人,看了别人的模板代码都是只知其一,不知其二,只知道这道题求区间和是这么写,换一道题换一个题目背景自己就不知道要怎么修改了。
其实,可以这么理解,线段树在建树的过程中,叶子结点上储存的就是你所输入的那些数值,那么两个叶子结点的父节点储存的是什么呢?这个要取决于题意,如果要求的是区间和,那么存的就是两个子节点的和,如果要求最大值,就存两个子节点较大的一个,最小值则相反。建树的过程是从1,1,n开始的,也就是说从根节点开始建树,接下来用递归的方式,分别建左右子树,直到叶子结点为止。
那么我们不难写出建树的代码(下面的代码是求区间和的)
int build(int rt,int l,int r)
{
t[rt].r=r;
t[rt].l=l;
if(l==r)
return t[rt].sum=num[l];
int mid=(l+r)>>1;
int a=build(rt<<1,l,mid);
int b=build(rt<<1|1,mid+1,r);
return t[rt].sum=a+b;
}
接下来我们来说update函数。如果所给的点超出了我们的区间范围,那么直接返回根节点的数值即可,如果在我们所给的区间范围内,那么就开始查找,查找的过程也一样用到了递归的方法,直到我们找到所要修改的那个点,修改即可,但是我们要注意,修改的过程中,比如我们要修改第三个点,不是只有3,3这个区间需要修改,所有包含3这个节点的区间都是需要修改的,如果我们的区间是1-5,那么(1,5),(3,5),(3,4),(3,3)这四个范围都是需要修改的,所以在递归的过程中也要加上这些节点的修改。代码不难写出:
int update(int rt,int pos,int val)
{
if(pos<t[rt].l||pos>t[rt].r)
return t[rt].sum;
if(pos==t[rt].l&&pos==t[rt].r)
return t[rt].sum+=val;
int a=update(rt<<1,pos,val);
int b=update(rt<<1|1,pos,val);
t[rt].sum=a+b;
return t[rt].sum;
}
最后我们要来说一下query函数。也就是最终的查找函数,查找的过程中同样需要用到递归的方法,我们也要找到这个区间范围然后输出这个范围内所存储的值,但是如果,我们的范围是1-5,我们要查找的范围是2-5怎么办?这个范围并不是一个已经存储好值的区间,这个时候,我们就需要把(2,2)和(3,5)区间加到一起然后输出,也就是说在递归的时候有三种情况。代码同样不难写出:
int query(int rt,int l,int r)
{
if(r<t[rt].l||l>t[rt].r)
return 0;
if(l<=t[rt].l&&r>=t[rt].r)
return t[rt].sum;
int a=query(rt<<1,l,r);
int b=query(rt<<1|1,l,r);
return a+b;
}
下边贴出两道HDU线段树单点更新的代码:
HDU 1166 敌兵布阵:
敌兵布阵
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 57695 Accepted Submission(s): 24379
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的责骂的.
中央情报局要研究敌人究竟演习什么战术,所以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条命令
每组数据第一行一个正整数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以内。
对于每个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 <cstdio>
#include <iostream>
#include <cstring>
#include <string>
#include <algorithm>
#define N 50010
int num[N];
using namespace std;
struct node
{
int l,r;
int sum;
}t[N*4];
int build(int rt,int l,int r)
{
t[rt].r=r;
t[rt].l=l;
if(l==r)
return t[rt].sum=num[l];
int mid=(l+r)>>1;
int a=build(rt<<1,l,mid);
int b=build(rt<<1|1,mid+1,r);
return t[rt].sum=a+b;
}
int query(int rt,int l,int r)
{
if(r<t[rt].l||l>t[rt].r)
return 0;
if(l<=t[rt].l&&r>=t[rt].r)
return t[rt].sum;
int a=query(rt<<1,l,r);
int b=query(rt<<1|1,l,r);
return a+b;
}
int update(int rt,int pos,int val)
{
if(pos<t[rt].l||pos>t[rt].r)
return t[rt].sum;
if(pos==t[rt].l&&pos==t[rt].r)
return t[rt].sum+=val;
int a=update(rt<<1,pos,val);
int b=update(rt<<1|1,pos,val);
t[rt].sum=a+b;
return t[rt].sum;
}
int main()
{
int T,i,cas=1,n,x,y;
string b;
scanf("%d",&T);
while(T--)
{
scanf("%d",&n);
for(i=1;i<=n;i++)
scanf("%d",&num[i]);
build(1,1,n);
printf("Case %d:\n",cas++);
while(cin>>b)
{
if(b[0]=='Q')
{
scanf("%d%d",&x,&y);
cout<<query(1,x,y)<<endl;
}
if(b[0]=='A')
{
scanf("%d%d",&x,&y);
num[x]+=y;
update(1,x,y);
}
if(b[0]=='S')
{
scanf("%d%d",&x,&y);
num[x]-=y;
update(1,x,-y);
}
if(b[0]=='E')
break;
}
}
return 0;
}
HDU 1754 I hate It:
I Hate It
Time Limit: 9000/3000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 49540 Accepted Submission(s): 19451
Problem Description
很多学校流行一种比较的习惯。老师们很喜欢询问,从某某到某某当中,分数最高的是多少。
这让很多学生很反感。
不管你喜不喜欢,现在需要你做的是,就是按照老师的要求,写一个程序,模拟老师的询问。当然,老师有时候需要更新某位同学的成绩。
这让很多学生很反感。
不管你喜不喜欢,现在需要你做的是,就是按照老师的要求,写一个程序,模拟老师的询问。当然,老师有时候需要更新某位同学的成绩。
Input
本题目包含多组测试,请处理到文件结束。
在每个测试的第一行,有两个正整数 N 和 M ( 0<N<=200000,0<M<5000 ),分别代表学生的数目和操作的数目。
学生ID编号分别从1编到N。
第二行包含N个整数,代表这N个学生的初始成绩,其中第i个数代表ID为i的学生的成绩。
接下来有M行。每一行有一个字符 C (只取'Q'或'U') ,和两个正整数A,B。
当C为'Q'的时候,表示这是一条询问操作,它询问ID从A到B(包括A,B)的学生当中,成绩最高的是多少。
当C为'U'的时候,表示这是一条更新操作,要求把ID为A的学生的成绩更改为B。
在每个测试的第一行,有两个正整数 N 和 M ( 0<N<=200000,0<M<5000 ),分别代表学生的数目和操作的数目。
学生ID编号分别从1编到N。
第二行包含N个整数,代表这N个学生的初始成绩,其中第i个数代表ID为i的学生的成绩。
接下来有M行。每一行有一个字符 C (只取'Q'或'U') ,和两个正整数A,B。
当C为'Q'的时候,表示这是一条询问操作,它询问ID从A到B(包括A,B)的学生当中,成绩最高的是多少。
当C为'U'的时候,表示这是一条更新操作,要求把ID为A的学生的成绩更改为B。
Output
对于每一次询问操作,在一行里面输出最高成绩。
Sample Input
5 6 1 2 3 4 5 Q 1 5 U 3 6 Q 3 4 Q 4 5 U 2 9 Q 1 5
Sample Output
5 6 5 9HintHuge input,the C function scanf() will work better than cin
Author
linle
Source
题目链接:HDU1754 I Hate It
#include <cstdio>
#include <iostream>
#include <string>
#include <cstring>
#define N 200010
using namespace std;
struct tree
{
int l,r;
int max;
}t[N*3];
int num[N];
int build(int rt,int l,int r)
{
t[rt].l=l;
t[rt].r=r;
if(l==r)
return t[rt].max=num[l];
int mid=(l+r)>>1;
int a=build(rt<<1,l,mid);
int b=build(rt<<1|1,mid+1,r);
return t[rt].max=max(a,b);
}
int update(int rt,int pos,int val)
{
if (pos < t[rt].l|| t[rt].r< pos)
return t[rt].max ;
if (t[rt].l== pos && t[rt].r== pos)
return t[rt].max = val;
int a, b ;
a = update (2 * rt, pos, val) ;
b = update (2 * rt + 1, pos, val) ;
t[rt].max = max (a, b) ;
return t[rt].max ;
}
int query(int rt,int l,int r)
{
int mid ;
if (t[rt].l> r|| t[rt].r < l)
return 0 ;
if (l<= t[rt].l&& t[rt].r<=r)
return t[rt].max;
int a, b ;
a = query(2 * rt, l, r) ;
b = query(2 * rt + 1, l, r) ;
return max (a, b) ;
}
int main()
{
int n,m,i,x,y;
char b;
while(~scanf("%d%d",&n,&m))
{
for(i=1;i<=n;i++)
scanf("%d",&num[i]);
build(1,1,n);
for(i=0;i<m;i++)
{
scanf("%c",&b);
if(b=='Q')
{
scanf("%d%d",&x,&y);
cout<<query(1,x,y)<<endl;
}
if(b=='U')
{
scanf("%d%d",&x,&y);
num[x]=y;
update(1,x,y);
}
}
}
return 0;
}
看了上面这些代码我们会发现其实线段树单点更新只要明白了思想,写出代码是很容易的,代码的写法有很多种,希望大家也能写出属于自己的线段树~
接下来几天可能会更一些其他数据结构的文章,明天会更区间更新~