树状数组(数据结构)

本文详细介绍了树状数组的概念及其实现方法,包括如何通过lowbit函数确定数组元素的影响范围、如何更新树状数组中的元素值以及如何求取区间和。

   对于序列a,我们设一个数组C,这个C可以称之为a[]的树状数组:
    ◦C[i] = a[i – 2^k + 1] + … + a[i]
    ◦k为i在二进制下末尾0的个数
   ◦2k就是i 保留最右边的1,其余位全变0
   ◦i从1开始算!

简单地说就是: 

树状数组
a[]  a[]数组下标  a[]数组的下标用二进制表示c[]
a[1]101c[1]=a[1]
a[2]210c[2]=a[1]+a[2]
a[3]311c[3]=a[3]
a[4]4100c[4]=a[1]+a[2]+a[3]+a[4]
a[5]5101c[5]=a[5]
a[6]6110c[6]=a[5]+a[6]
a[7]7111c[7]=a[7]
a[8]81000c[8]=a[1]+a[2]+a[3]+a[4]+a[5]+a[6]+a[7]+a[8]

 

 

 

 

 

 

 

 

 

 

所以就有:

C[i] = a[i-lowbit(i)+1] + …  + a[i],这总共有多少项,是由a[]数组的下标二进制中,从后往前看,第一个1后面有K个零,2^k决定的。

比如:a[12]  1100 第二个1后面有两个零,所以是2^2共四项。c[12]=a[9]+a[10]+a[11]+a[12]

那下面的事情就是如何找到a[]数组下标中,最后一个1所在的位置了。可以用预备函数:

1.预备函数:

int lowbit(int k)
{
	return k&(-k);
}

k=6,将返回2,k=8将返回8.

2.修改原始数据某一元素的值:

将a[k]的值增加v,需要把c[k]以及后继节点的值全部增加k,一直到最大值。

#define namx 100
int c[nmax];
void add(int k,int v)
{
	while(k<nmax)
	{
		c[k]=c[k]+v;
		k=k+lowbit(k);
	}
}

3.求原始数据前k项的和:

只需要把c[k]及其前驱节点的值累加即可。

int sum(int k)
{
	int cnt=0;
	while(k>0)
	{
		cnt=cnt+c[k];
		k=k-lowbit(k); 
	}
	return cnt;
}

 

基于可靠性评估序贯蒙特卡洛模拟法的配电网可靠性评估研究(Matlab代码实现)内容概要:本文围绕“基于可靠性评估序贯蒙特卡洛模拟法的配电网可靠性评估研究”,介绍了利用Matlab代码实现配电网可靠性的仿真分析方法。重点采用序贯蒙特卡洛模拟法对配电网进行长时间段的状态抽样与统计,通过模拟系统元件的故障与修复过程,评估配电网的关键可靠性指标,如系统停电频率、停电持续时间、负荷点可靠性等。该方法能够有效处理复杂网络结构与设备时序特性,提升评估精度,适用于含分布式电源、电动汽车等新型负荷接入的现代配电网。文中提供了完整的Matlab实现代码与案例分析,便于复现和扩展应用。; 适合人群:具备电力系统基础知识和Matlab编程能力的高校研究生、科研人员及电力行业技术人员,尤其适合从事配电网规划、运行与可靠性分析相关工作的人员; 使用场景及目标:①掌握序贯蒙特卡洛模拟法在电力系统可靠性评估中的基本原理与实现流程;②学习如何通过Matlab构建配电网仿真模型并进行状态转移模拟;③应用于含新能源接入的复杂配电网可靠性定量评估与优化设计; 阅读建议:建议结合文中提供的Matlab代码逐段调试运行,理解状态抽样、故障判断、修复逻辑及指标统计的具体实现方式,同时可扩展至不同网络结构或加入更多不确定性因素进行深化研究。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值