详解树状数组-程序员宅基地

技术标签:   算法  深度优先  c语言  数据结构  

前言

树状数组二叉索引树(Binary Indexed Tree),又以其发明者命名为 Fenwick 树。其初衷是解决数据压缩里的累积频率的计算问题,现多用于高效计算数列的前缀和区间和。它可以以 O(logn) 的时间得到任意前缀和。并同时支持在 O(logn) 时间内支持动态单点值的修改。空间复杂度 O(n)


一、树状数组概括

树状数组是一个查询修改复杂度都为log(n)数据结构。主要用于数组的单点修改&&区间求和,另外一个拥有类似功能的是线段树

具体区别和联系如下:

1.两者在复杂度上同级, 但是树状数组常数明显优于线段树, 其编程复杂度也远小于线段树.

2.树状数组的作用被线段树完全涵盖, 凡是可以使用树状数组解决的问题, 使用线段树一定可以解决, 但是线段树能够解决的问题树状数组未必能够解决.

3. 树状数组的突出特点是其编程的极端简洁性, 使用lowbit技术可以在很短的几步操作中完成树状数组的核心操作,其代码效率远高于线段树


二、树状数组的应用

1.单点修改+区间查询

代码示例


int lowbit(int x)
{
  return x & (-x);//表示求数组下标二进制的非0最低位所表示的值
}

//查找1~x的和
int find_sum(int x)
{
  int ans = 0;
  while(x)
  {
    ans += c[x];//从右往左累加求和
    x -= lowbit(x);
  }
  return x;
}

//单点修改
void gexi(int x,int v)
{
  a[x] += v;
  while(x <= n)
  {
    c[x] += v;
    x += lowbit(x);//由叶子节点向上更新树状数组C,从左往右更新
  }
}

实现原理

模板中最常见的三个函数:

①取数组下标二进制非0最低位所表示的值;
②单点更新;
③区间查询。

树状数组,顾名思义是树状的数组,我们首先引入二叉树,叶子节点代表A[1]~A[8]。

现在变形一下:

现在定义每一列的顶端节点C数组(其实C数组就是树状数组),如图:

理解树状数组的重点

C[i]代表子树的叶子节点的权值之和,如图可以知道:

C[1]=A[1];

C[2]=A[1]+A[2];

C[3]=A[3];

C[4]=A[1]+A[2]+A[3]+A[4];

C[5]=A[5];

C[6]=A[5]+A[6];

C[7]=A[7];

C[8]=A[1]+A[2]+A[3]+A[4]+A[5]+A[6]+A[7]+A[8];

首先是区间查询(求和):

利用C[i]数组,求A数组前i项和,举两个栗子:

①i=7

前7项和:sum[7]=A[1]+A[2]+A[3]+A[4]+A[5]+A[6]+A[7];

而C[4]=A[1]+A[2]+A[3]+A[4];C[6]=A[5]+A[6];C[7]=A[7];

可以得到:sum[7]=C[4]+C[6]+C[7]。

数组下标写成二进制:sum[(111)]=C[(100)]+C[(110)]+C[(111)];

②i=5

前5项和:sum[5]=A[1]+A[2]+A[3]+A[4]+A[5];

而C[4]=A[1]+A[2]+A[3]+A[4];C[5]=A[5];

可以得到:sum[5]=C[4]+C[5];

数组下标写成二进制:sum[(101)]=C[(100)]+C[(101)];

细细观察二进制,树状数组追其根本就是二进制的应用,结合代码演示一下代码过程:


//查找1~x的和
int find_sum(int x)
{
  int ans = 0;
  while(x)
  {
    ans += c[x];
    x -= lowbit(x);
  }
  return x;
}
代码推演

对于i=7进行演示:


7(111)  ans+=C[7]

lowbit(7)=001 7-lowbit(7)=6(110) ans+=C[6]

lowbit(6)=010 6-lowbit(6)=4(100) ans+=C[4]

lowbit(4)=100 4-lowbit(4)=0(000) break;

对于i=5进行演示:


5(101)  ans+=C[5]

lowbit(5)=001 5-lowbit(5)=4(100) ans+=C[4]

lowbit(4)=100 4-lowbit(4)=0(000) break;

然后单点更新:

当我们修改A数组中某个值时,应当如何更新C数组呢?回想一下,区间查询的过程,再看一下上文中列出的过程。这里声明一下:单点更新实际上是不修改A数组的,而是修改树状数组C,向上更新区间长度为lowbit(i)所代表的节点的值。


//单点修改
void gexi(int x,int v)
{
  a[x] += v;
  while(x <= n)
  {
    c[x] += v;
    x += lowbit(x);
  }

如图:当在A[1]加上值val,即更新A[1]时,需要向上更新C[1],C[2],C[4],C[8],这个时候只需将这4个节点每个节点的值加上val即可。这里为了方便大家理解,人为添加了个A数组表示每个叶子节点的值,事实上A数组并不用修改,实际运用中也可不设置A数组,单点更新只需修改树状数组C即可。下标写成二进制:C[(001)],C[(010)],C[(100)],C[(1000)];

lowbit(1)=001 1+lowbit(1)=2(010) C[2]+=val;

lowbit(2)=010 2+lowbit(2)=4(100) C[4]+=val;

lowbit(4)=100 4+lowbit(4)=8(1000) C[8]+=val;

由于c[1] c[2] c[4] c[8] 都包含有A[1],所以在更新A[1]时实际上就是更新每一个包含A[1]的节点。

总结

树状数组的重点就是利用二进制的变化动态地更新树状数组

树状数组的每一个节点并不是代表原数组的值,而是包含了原数组多个节点的值。

所以在更新A[1]时需要将所有包含A[1]的C[i]都加上val这也就利用到了二进制的神奇之处。

如果是更新A[i]的值,则每一次对C[i] 中的 i 向上更新,即每次i+=lowbit(i),这样就能C[i] 以及C[i] 的所有父节点都加上val。

反之求区间和也是和更新节点值差不多,只不过每次 i-=lowbit(i)。

结语

如果这篇文章对您有帮助的话,请记得点赞收藏加关注吖(●'◡'●)

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/weq2011/article/details/128835460

智能推荐

c++路径转义字符_c++文件名中文转义-程序员宅基地

文章浏览阅读2.3k次。root_dir = R"(E:\data\aa\aaa\aaa\)"_c++文件名中文转义

Android Fragment简介_android fragment包括-程序员宅基地

文章浏览阅读407次。Fragment概述介绍Fragment是一种可以嵌入在活动中的UI片段,能够让程序更加合理和充分地利用大屏幕的空间,出现的初衷是为了适应大屏幕的平板电脑,可以将其看成一个小型Activity,又称作Activity片段。使用Fragment可以把屏幕划分成几块,然后进行分组,进行一个模块化管理。Fragment不能够单独使用,需要嵌套在Activity中使用,其生命周期也受到宿主Activity的生命周期的影响。2.特点·Fragment依赖于Activity,不能独立存在·一个Act_android fragment包括

jQuery 个人之前的笔记_jquery详细笔记-程序员宅基地

文章浏览阅读2.7k次,点赞3次,收藏5次。jQuery 基础篇概述jQuery是一款优秀的JavaScript库 ,从命名可以看出jQuery最主要的用途是用来做查询(iQuery=js+Query) ,正如jQuery官方Logo副标题所说(write less, domore)使用jQuery能 上我们对HTML文档遍历和操作、事件处理、动画以及Ajax变得更加简单入门和原生JS区别1、原生JS: 等DOM元素加载完毕,并且图片也加载完毕才会执行 编写多个入口函数,后面会覆盖前面2、jQuery 等Dom元素加载完毕,但_jquery详细笔记

C++实现Delaunay三角网生长算法_增量法生成delaunay三角网-程序员宅基地

文章浏览阅读8.4k次,点赞25次,收藏142次。一、概述1.1 三角网的介绍三角网是由一系列连续三角形构成的网状的平面控制图形,是三角测量中布设连续三角形的两种主要扩展形式,同时向各方向扩展而构成网状.适用于地势起伏大,通视条件比较好的场地。三角网是实现地形三维可视化,数字地面模型(Digital Terrain Model,简称DTM)是一种很有效的途径。DTM主要是由栅格和不规则三角网(Triangulated Irregular Network,简称TIN)两种数据格式来表示,相比于栅格TIN具有许多..._增量法生成delaunay三角网

TCGA数据下载和整理工具----GDCRNATools_gdcrnatools软件包进行差异基因分析-程序员宅基地

文章浏览阅读1.9w次,点赞7次,收藏62次。TCGA数据下载和整理的网站及软件发表很多了,比如Broad GDAC Firehose, Oncomine, TCGAbiolinks,TCGA-Assembler, TCGA2STAT,RTCGAToolbox等等,这些网站或软件要么使用的是TCGA更新前的数据,要么运行起来比较繁琐。当然各个工具都有其优势所在。之前在论坛里分享了自己下载和整理TCGA数据的Python代码。最近忙里偷_gdcrnatools软件包进行差异基因分析

win7更改计算机时间,win7系统自动更改日期时间是怎么回事-程序员宅基地

文章浏览阅读1.9k次。工具/原料硬件:计算机操作系统:Windows7方法/步骤1.Windows7系统不能更改日期和时间的解决方法2.在本地组策略编辑器窗口,展开Windows设置 - 安全设置 - 本地策略;3.在本地策略中找到:用户权限分配,左键点击:用户权限分配,在用户权限分配对应的右侧窗口找到:更改系统时间,并左键双击:更改系统时间;4.在打开的更改系统时间 属性窗口,我们点击:添加用户或组(U);5.在选择..._win7系统时间老是自己跳变

随便推点

Java给JPanel添加彩虹边框_给jpanel加边框-程序员宅基地

文章浏览阅读685次。使用多线程RGB变换为JPanel添加五彩斑斓的特效_给jpanel加边框

华为p8升级android8系统资源,华为新系统广受好评 目前EMUI10升级用户数已破亿-程序员宅基地

文章浏览阅读844次。华为消费者BG软件部总裁王成录昨日宣布,升级EMUI10的用户数已经突破1亿,涵盖35款机型。该数据不但充分展现了华为手机庞大的用户基数,也体现了华为在系统升级方面所付出的努力,正得到众多用户的肯定与支持。确定 时延引擎 与方舟编译器,解决安卓卡顿问题安卓手机曾经为人所诟病的一点要数“卡顿”,手机越用越慢,严重影响着安卓用户的用机体验。四年前,随着提出“天生快、一生快”标语的EMUI5.0发布,华..._华为p8最新系统

微服务架构,这一篇就够了!-程序员宅基地

文章浏览阅读1.9k次。所有的微服务都是独立的Java进程跑在独立的虚拟机上,所以服务间的通信就是IPC(inter process communication),已经有很多成熟的方案。原来的Monolithic方式开发,所有的服务都是本地的,UI可以直接调用,现在按功能拆分成独立的服务,跑在独立的一般都在独立的虚拟机上的 Java进程了。它通常不依赖其他服务。后台有N个服务,前台就需要记住管理N个服务,一个服务下线/更新/升级,前台就要重新部署,这明显不服务我们 拆分的理念,特别当前台是移动应用的时候,通常业务变化的节奏更快。_微服务架构

主键索引和非主键索引的区别-程序员宅基地

文章浏览阅读871次。总的来说,主键索引是表的唯一标识索引,具有唯一性和快速访问的特点;而非主键索引可以提供更多的灵活性和覆盖更多的查询场景,但可能性能略低于主键索引。_主键索引和非主键索引的区别

linux启动进入bios设置密码,通过bios怎么设置开机密码-程序员宅基地

文章浏览阅读1.2k次。电脑不想被他人乱动,来设置下BIOS管理员密码和开机密码,就让学习啦小编来告诉大家通过bios怎么设置开机密码的方法吧,希望对大家有所帮助。通过bios设置开机密码方法计算机开机以后,按键盘的Delete键进入BIOS的设置画面,如下图所示。因为开机可以按Delete键进入设置画面的时间很短,您可以在计算机一开机就慢慢的重复按Delete键,以免错过进入设置画面又要重新再开机。按键盘向下箭头键移到..._bios开机密码 画面

批处理获取所有文件、文件夹名字_bat获取文件夹下所有文件名和文件夹名称-程序员宅基地

文章浏览阅读1.6w次,点赞14次,收藏45次。已收藏下面这个链接的方法也不错excel批处理技巧:如何制作文件档案管理系统excel批处理技巧:如何制作文件档案管理系统http://www.360doc.com/content/18/0913/13/18781560_786337463.shtml有时候我们整理文件的时候需要列出文件夹里面所有的文件名或者文件夹名,生成一个文件目录,一个个重命名然后复制到word或者记事本的方法显示有点太繁琐了。网上有一些自动生成文件目录的程序,比如我之前一直在用的DirIndex.exe。但最近我发现_bat获取文件夹下所有文件名和文件夹名称