最优查找树-实验报告
- 格式:doc
- 大小:37.00 KB
- 文档页数:8
计算机算法设计与分析实验报告实验名称:最优树查找班级:08211306姓名:刘宇轩学号:08211312一、一般原理二叉排序树(Binary Sort Tree)又称二叉查找树。
它或者是一棵空树;或者是具有下列性质的二叉树:(1)若左子树不空,则左子树上所有结点的值均小于它的根结点的值;(2)若右子树不空,则右子树上所有结点的值均大于它的根结点的值;(3)左、右子树也分别为二叉排序树;问题描述给定一个有序序列K={k1<k2<k3<,……,<kn}和他们被查询的概率P={p1,p2,p3,……,pn},要求构造一棵二叉查找树T,使得查询所有元素的总的代价最小。
对于一个搜索树,当搜索的元素在树内时,表示搜索成功。
当不在树内时,表示搜索失败,用一个“虚叶子节点”来标示搜索失败的情况,因此需要n+1个虚叶子节点{d0<d1<……<dn}。
其中d0表示搜索元素小于k1的失败结果,dn表示搜索元素大于kn的失败情况。
di(0<i<n)表示搜索节点在ki和k(i+1)之间时的失败情况。
对于应di的概率序列是Q={q0,q1,……,qn}。
二、算法的复杂度下面是对于查找期望代价的解释:对于键值ki, 如果其在构造的二叉查找树里的深度(离开树根的分支数)为depthT(ki),则搜索该键值的代价= depthT(ki) +1(需要加上深度为0 的树根节点)。
由于每个键值被查找的概率分别为pi,i=1,2,3…,n。
在二叉树中T内搜索一次的期望代价为:E[T] = (depth(ki)+1)*pi //对每个i=1~n,搜索成功情况+ (depth(di)+1)*qi //对每个i=0~n,搜索失败情况时间复杂度1、穷举穷举构造最优二叉查找树,其实就是这样的一个问题:给一个拥有n 个数的已排序的节点,可以将其构造成多少种不同的BST(用来找到一个最优的二叉查找树)?设可以构造成T(n)个,那么枚举每一个元素作为根节点的情况,当第一个元素作为根节点时,其余n-1 个构成右子树,无左子树,是n-1 情况时的子问题,共T(n-1)种;当第二个元素作为根节点时,左子树有1 个元素,右子树有n-2 个元素,根据乘法原理共有T(1)T(n-2)种情况……依此类推得到:T(n)= T(0)T(n-1) + T(1)T(n-2) + T(2)T(n-3) + ...... + T(n-2)T(1) + T(n-1)T(0);此外,有T(0)=T(1)=1。
下面来求解T(n):定义函数f(x) = T(0) + T(1)*x + T(2)*x2 + ......那么有:f(x)2 = (T(0)2) + (T(0)T(1) + T(1)T(0)) ·x + (T(0)T(2) + T(1)T(1) + T(2)T(0)) ·x2 + ......= T(1) + T(2) ·x + T(3) ·x2 + ......= (f(x) - T(0)) / x= (f(x) - 1) / x这样解方程得到f(x) = [1 - (1 - 4x)1/2] / 2x右边进行泰勒展开,再与定义式比较最终得到:T(n) = (2n)! / (n!(n+1)!)然后根据Stirling 公式:n! ~(2πn)1/2*(n/e)n于是有(2n)!/n!(n+1)!~(4n1/2*2n2n)/(2n1/2*nn*(2(n+1))1/2*(n+1)n)~4n*(n+1)-3/2*(n/(n+1))n ~4n*n-3/2因此最后得到穷举方法构造最优二叉查找树的时间复杂度:T(n) = Ο(4n*n-3/2)2、递归实际上左右子树是互不影响的,不需要穷举所有左右子树的组合,所以不需要用乘法原理,加法原理就可以了,这样式子变为:T(n) = T(0) + T(n-1) + T(1) + T(n-2) + T(2) + T(n-3) + ...... + T(n-2) + T(1) + T(n-1) + T(0)= 2(T(0) + T(1) + T(2) + ...... + T(n-1))= 3T(n-1)所以得到T(n) = O(3n),还是指数级的一个算法三、动态规划上面得到指数级算法的原因在于,计算了很多重复的子树情况,一些子树的查找代价被计算了很多遍;而一棵树如果是最优二叉搜索树,那么要么它是空树,要么它的左、右子树也是最优二叉搜索树,因此只需要将子树的查找代价记录下来,采用记忆化搜索或者是自底向上的动态规划的方法,虽然需要消耗一定的空间,但可以把时间复杂度从指数级降到多项式级,这些空间消耗也是可以接受的。
步骤一:寻找最优子结构。
一个最优二叉树的子树必定包含连续范围的关键字ki~kj,1<=i<=j<=n,同时也必须含有连续的虚叶子节点di-1~dj。
如果一棵最优二叉查找树T有一棵含有关键字ki~kj的子树T',那么,T'也是一棵最优查找树,这通过剪贴思想可以证明。
现在开始构造最优子结构:在ki~kj中,选定一个r,i<=r<=j,使以kr为根,ki~k(r-1)和k(r+1)~kj为左右孩子的最优二叉树。
注意r=i 或者r=j的情况,表示左子树或右子树只有虚叶子节点。
步骤二:一个递归解定义e[i,j]为一棵包含关键字ki~kj的最优二叉树的期望代价。
当j=i-1时没有真实的关键在,只有虚叶子节点d(i-1)。
于是:当j=i-1时,e[i,i-1]=q(i-1)。
当j>=i时,需要选择合适的kr作为根节点,然后其余节点ki~K(r-1)和k(r+1)~kj构造左右孩子。
这时要考虑左右孩子这些节点成为一个节点的子树后,它的搜索代价的变化:根据E[T]的计算,得知它们的期望代价增加了“子树中所有概率的总和”w。
w[i,j] = pl // 对每个l=i~j+ql //对每个l=i-1~j于是当j>=i时,e[i,j]=pr + (e[i,r-1]+w[i,r-1])+(e[r+1,j]+w[r+1,j]) = e[i,r-1] + e[r+1,j]+w[i,j];步骤三:计算最优二叉树的期望代价e[i,j] = q(i-1) //如果j=i-1min(e[i,r-1] + e[r+1,j]+w[i,j])//如果i<=j,其中i<=r<=j w[i,j] = q(i-1) // 如果j=i-1w[i,j]=w[i,j-1]+pj+qj //如果i<=j四、算法实现#include <iostream>using namespace std;#define MAX_COUNT 100#define MAX_NUM 65536//p中为有序关键字k1到k5的搜索概率,k1<k2<k3<k4<k5double p[MAX_COUNT] = {0.22 , 0.13 , 0.03 , 0.26 , 0.36};double q[MAX_COUNT] = {0.04 , 0.06 , 0.04 , 0.06 , 0.09};void optimal_bst(double e[][MAX_COUNT] ,int root[][MAX_COUNT],double w[][MAX_COUNT],int n) {int i =0,j=0;//针对左或右孩子为空树情况初始化for(i = 1;i<=n+1;i++){e[i][i-1] = q[i-1];w[i][i-1] = q[i-1];}int l = 0;/*计算顺序如下:根据计算式:e[i,j] = e[i,r-1]+e[r+1,j首先计算节点个数为1的最优二叉树的代价e[1,1],e[2,2]……接着计算节点个数为1的最优二叉树的代价e[1,2],e[2,3]…………最后计算结点个数为n的最优二叉树的代价e[1,n],利用之前保存的较少结点最优二叉树的结果。
*/for(l = 1;l<=n;l++){for(i = 1;i<=n-l+1;i++){j = i+l-1;e[i][j] = MAX_NUM;w[i][j] = w[i][j-1] + p[j]+q[j];for(int r = i;r<=j;r++){double t = 0;t = e[i][r-1]+e[r+1][j] + w[i][j];if(t<e[i][j]){e[i][j]= t;root[i][j] = r;}}}}}int main(){double e[MAX_COUNT][MAX_COUNT];int r[MAX_COUNT][MAX_COUNT];double w[MAX_COUNT][MAX_COUNT];optimal_bst(e,r,w,5);for(int i =1;i<=6;i++) {for(int j = 0;j<=5;j++){cout << r[i][j] << " ";}cout << endl;}system("pause");return 0;}五、实验总结通过这次试验,使得我对最优查找树重新进行了认识与更加深入的了解。
虽然以往在数据结构的课程中曾接触过类似内容,但是当时并没有深究。
在一开始理解中,感觉最优查找树与哈弗曼树十分相似。
但经过仔细的比对,发现二者还是有很大的差别的,最大的不同就是最优查找树要求左叶子的值必须比右叶子小,通过这一个研究对比,我对最有查找树的理解又上升了一个层次。