挑战性题目DSCT302:求不同形态的平衡二叉树数目

本文探讨了如何计算不同形态的平衡二叉树数目。利用动态规划与快速傅里叶变换,提出了一种高效算法来解决该问题。适用于n个节点的所有可能平衡二叉树形态的计数。

挑战性题目DSCT302:求不同形态的平衡二叉树数目

问题描述

n个节点的不同形态的平衡二叉树数目是确定的。任意输入一个正整数n,请问有多少种不同的平衡二叉树形状?例如,输入3,输出1;输入4,输出4。

题解

考虑使用动态规划求解该题,dp[n][d]表示节点数为nnn,深度为ddd的平衡二叉树的数目。对于给定的平衡树节点数n,容易得出,其深度ddd的范围为⌊log2n⌋+1 ⌊log2n⌋+2\left\lfloor{log}_2n\right\rfloor+1~\left\lfloor{log}_2n\right\rfloor+2log2n+1 log2n+2,下界为一个完全二叉树,上界可以用给定深度求平衡二叉树最少节点数的递推公式dep[n]=dep[n-1]+dep[n-2]+1,dep[1]=1,dep[2]=2反解得出。因此,我们的答案就是dp[n][log2n+1]+dp[n][log2n+1]

现在,要求出任意dp[n][d],我们可以由两个深度为d−1d-1d1的子树合并出来,即枚举其中一个子树的大小iiidp[n][d]=∑0n−1\sum_0^{n-1}0n1(dp[i][d-1]×dp[n-i-1][d-1]+dp[i][d-2]×dp[n-i-1][d-1]+dp[i][d-1]×dp[n-i-1][d-2]),该递归边界为dp[0][0]=dp[1][1]=1
容易看出,上式其实是一个卷积形式,可以写作dp[n][d]=dp[n][d-1]*dp[n][d-1]+dp[n][d-2]*dp[n][d-1]+dp[n][d-1]*dp[n][d-2]
改写为卷积形式后,就可以通过快速傅里叶变换(FFT),将O(n2)O\left(n^2\right)O(n2)的求和优化为O(nlog2n)O(n{log}_2n)O(nlog2n)

综上,对于每层递归,做一次卷积运算复杂度为O(nlog2n)O(n{log}_2n)O(nlog2n),每层递归二叉树的深度ddd都会减小111,同时ddd也是log2n{log}_2nlog2n级别的,于是总复杂度为O(nlog2nlog2n)O(n{log}_2n{log}_2n)O(nlog2nlog2n)

代码

能打为什么不打?

#include<cstdio>
#include<cstdlib>
using namespace std;
long long ans[]={
0,
1,
2,
1,
4,
6,
4,
17,
32,
44,
60,
70,
184,
476,
872,
1553,
2720,
4288,
6312,
9004,
16088,
36900,
82984,
174374,
346048,
653096,
1199384,
2160732,
3812464,
6617304,
11307920,
18978577,
31327104,
51931296,
90400704,
170054336,
341729616,
711634072,
1491256624,
3084996748,
6246978752};
void out(long long x){if(x>9)out(x/10);putchar(x%10+'0');}
int main(int argc,char* argv[])
{
	int n=atoi(argv[1]);
	out(ans[n]);
}

Lz变换 考 虑 一 个 离 散 状 态 连 续 时 间 的 马 尔 科 夫 过 程 (DSCT):X (t)Î{x1xK}t ³ 0 ,共有 K 个可能的状态 i , i = 1K 。 状 态 集 x ={x1xK} 、转 移 率 矩 阵 A =aij(t)ij = 1K 和初始状态概率分布 p0 就完整的 定义了该马尔科夫过程。其中 p0 定义如下: p0 =[p10 = Pr{X (0) = x1}pK0 = Pr{X (0) = xK}] 对 于 一 个 离 散 状 态 连 续 时 间 的 马 尔 科 夫 过 程 X (t) = xAp0 ,Lz变换的函数 u(ztp0) 定义为如下: L Z{X (t)}= u(ztp0) = å i = 1 K pi(t)zxi (3) 考虑一个最简单的单部件两状态系统。马尔科夫过 程 X (t) = xAp0 定义如下: 状态集 x ={x1x2}={01} 转移率矩阵 A =aij(t) = æ èç-μ μ ö ø ÷ λ -λ ,其中 λμ 分别为 失效率和修复率 初始状态 p0 =[p10p20] = [01] 由 dP(t) dt = P(t)A (2.2式) 得到微分方程组: ìïïíïïî dp1(t) dt = -μp1(t) + λp2(t) dp2(t) dt = μp1(t) - λp2(t) , (4) 初值条件: ìíî p1(0) = p10 = 0 p2(0) = p20 = 1 通过做拉普拉斯变换,可以解得到结果如下: ìïïíïïî p1(t) = λ λ + μ - λ λ + μ exp{-(λ + μ)t} p2(t) = μ λ + μ + λ λ + μ exp{-(λ + μ)t} (5) 通过定义(3)式得到所给马尔科夫过程的Lz变换为: L z{X (t)}= å i = 1 2 pi(t)zxi = é ëê λ ù û ú λ + μ - λ λ + μ exp{-(λ + μ)t} z0 + éêë ùúû μ λ + μ + λ λ + μ exp{-(λ + μ)t} z1
最新发布
10-15
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

ShadyPi

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值