左式堆合并操作时间界的证明

阅读《数据结构与算法分析:java语言描述》一书时看到了对左式堆合并操作时间界的说明，但书中没有给出详细证明，这里自己思考并补全了证明细节，如下文(全文假定讨论的左式堆具有最小堆序，最大堆序的讨论是类似的)

证法一：

我们用 $T[i,j]$ 表示右路径(从左式堆的根节点出发抵达根节点或以其右子女为根的子树中子女数不为2的节点的最短路径)长度为i的左式堆L1和右路径长度为j的左式堆R1合并的平均时间代价。现设i>=0 且 j>=0 , L1和R1合并时存在两种可能，如果R1的根节点关键码值大于等于L1根节点关键码值，则应将L1右子树代表的左式堆LR1和左式堆R1合并，此时由于L1的右路径长度为i,故LR1根节点的零路径长度为i-1,而LR1根节点的零路径长度是LR1左路径(定义类似右路径)长度和右路径长度中的较小值，且由左式堆的性质，LR1左路径长度应当大于等于右路径长度，于是LR1的右路径长度即为i-1,这样将左式堆LR1和左式堆R1合并的平均时间代价可以表示为 $T[i-1,j]$ 。如果R1的根节点关键码值小于L1根节点关键码值，则应将左式堆L1和R1右子树代表的左式堆RR1合并，通过和上文类似的分析可知合并的平均时间代价可以表示为 $T[i,j-1]$

另外一个空左式堆和任意一个左式堆合并的时间代价显然为 $O(1)$ ,因此 $T[-1,0]=T[-1,1]=---=T[-1,j]=1$

$T[0,-1]=T[1,-1]=---=T[i,-1]=1$

现在考虑由平面上整点 $(m,k)$ $-1<=m<=i-1<=k<=j$ 以及这些整点中相邻点之间的垂直水平连线构成的矩形网格 ，由上文分析知，从整点 $(i,j)$ 出发，首先右路径长度为i的左式堆和右路径长度为j的左式堆的合并操作对应于整点 $(i,j)$ ,且在整点 $(i,j)$ 对应的合并操作中要么右路径长度为i的左式堆根节点的右子树和右路径长度为j的左式堆合并，此时该合并操作对应于整点 $(i-1,j)$ ，要么右路径长度为i的左式堆和右路径长度为j的左式堆根节点的右子树合并，此时该合并操作对应于整点 $(i,j-1)$ ,完成整点 $(i,j)$ 对应的合并操作相当于首先完成整点 $(i-1,j)$ 或 $(i,j-1)$ 的合并操作，然后再将这一合并操作的合并结果(一颗左式堆)链接至i或j的右子树并在必要时调换i,或j的左右子树，因此整点 $(i,j)$ 对应的合并操作也可视为连接 $(i-1,j)$ 或 $(i,j-1)$ 和 $(i,j)$ 的一条垂直(水平边)，更确切地说整点 $(i,j)$ 对应的合并操作完全可以看做先完成 $(i-1,j)$ 或 $(i,j-1)$ 对应的合并操作再沿这一垂直(水平边)回溯至整点 $(i,j)$ . 类似地，对代表完成整点 $(i,j)$ 对应的合并操作之前需要完成的合并操作的整点 $(i-1,j)$ 或整点 $(i,j-1)$ ,可以用和上文所述完全相同的方法继续分析，不断地抵达新的整点并选择连接这些整点中相邻整点的垂直水平边，这一过程中抵达一个新的整点后会向左或向下移动至下一个相邻整点，最终就可以得到向矩形网格左下方延伸的递降路径M，从 $(i,j)$ 出发沿M向左下方前进抵达的每个整点(除 $(i,j)$ 外)的纵坐标均小于等于其M上前驱整点的纵坐标,横坐标均小于等于其M上前驱整点的横坐标。如果令 $(i,j)$ 为M的出发点，则M的终点最终必然会在 $(-1,0)(-1,1)---(-1,j)(0,-1)(1,-1)---(i,-1)$ 之一,设其为 $(a,b)$ ,这样对右路径长度为i的左式堆和右路径长度为j的左式堆的合并操作可以用M描述为首先完成 $(a,b)$ 对应的合并操作，操作必然有一个被合并左式堆为空堆，时间代价为1，然后对 $(a,b)$在M上的前驱节点 $(a2,b2)$ 将合并得到的左式堆链接至右路径长为$a2$或$b2$(这取决于构造M时从 $(a2,b2)$ 抵达 $(a,b)$ 的方向(垂直或水平),即完成 $(a2,b2)$ 对应的合并操作前最后需要完成的合并操作是什么)的左式堆的根节点右指针域并在必要时调换根节点左右子树的位置，这一操作时间复杂度显然为 $O(1)$,该操作完成后 $(a2,b2)$ 对应的合并操作即完成，此时的累计时间复杂度为2，随后从 $(a2,b2)$出发,对 $(a2,b2)$ 在路径M上的前驱节点 $(a3,b3)$ 而言将 $(a2,b2)$ 对应的合并操作的合并结果链接至右路径长为$a3$或$b3$的左式堆的根节点右指针域并在必要是调换根节点左右子树位置这一操作时间复杂度显然为 $O(1)$,该操作完成后 $(a3,b3)$ 对应的合并操作即完成，此时的累计时间复杂度为3.随后从 $(a3,b3)$ 按照这一过程继续沿M向 $(i,j)$ 回溯，当最终回溯至 $(i,j)$ 后, $(i,j)$ 对应的合并操作即完成，此时的累计时间复杂度恰为M上整点个数，即为M的路径长度加一

若 $(a,b)=(-1,m)$ $0<=m<=j$ 则M路径长为 $i+1+j-m=i+j-m+1$

M上整点个数即M的累计时间复杂度为 $i+j-m+2$ 由于M为向左下角延伸递降路径，出发点为 $(i,j)$ ,终点为 $(-1,m)$ 故所有可能的M的数量为组合数 $C(i+j-m+1,i+1)$,理想情况下可以认为每一种路径可能出现的概率是相等的，均为 $\frac{1}{C(i+j-m+1,i+1)}$ ,于是每条路径M的期望时间复杂度为 $\frac{i+j-m+2}{C(i+j-m+1,i+1)}$ ， 对所有可能的M的期望时间复杂度求和得到出发点在 $(i,j)$ 终点在 $(a,b)=(-1,m)$ $0<=m<=j$ 的所有路径的期望时间复杂度 $\frac{i+j-m+2}{C(i+j-m+1,i+1)}C(i+j-m+1,i+1)=i+j-m+2$

可以认为当M的出发点为 $(i,j)$ 时，终点落在 $(-1,0)(-1,1)---(-1,j)(0,-1)(1,-1)---(i,-1)$ 上每一点的可能性是相等的，概率均为 $\frac{1}{i+j+2}$

于是出发点在 $(i,j)$ 终点在每一个整点 $(-1,m)$ $0<=m<=j$ 的所有路径的期望时间复杂度为 ${\sum }_{m=0}^j\frac{1}{i+j+2}i+j-m+2=j+1-\frac{j(j+1)}{2(i+j+2)}$

同理出发点在 $(i,j)$ 终点在每一个整点 $(m,-1)$ $0<=m<=i$ 的所有路径的期望时间复杂度为 ${\sum }_{m=0}^i\frac{1}{i+j+2}i+j-m+2=i+1-\frac{i(i+1)}{2(i+j+2)}$

累加以上两式得到整点 $(i,j)$ 对应的左式堆合并操作对应的期望时间复杂度 $T[i,j]==i+1-\frac{i(i+1)}{2(i+j+2)}+j+1-\frac{j(j+1)}{2(i+j+2)}=i+j+2-\frac{i^2+j^2}{2(i+j+2)}-\frac{i+j}{2(i+j+2)}$

$<=i+j+2-\frac{ij}{i+j+2}-\frac{1}{2(1+\frac{2}{i+j})}$

$<=i+j+2-\frac{1}{2(1+\frac{2}{i+j})}$

当i,j充分大时，由最后一个不等式可见 $T[i,j]<=i+j+\frac{3}{2}$

故 $T[i,j]=O(i+j)$ 设右路径长度为i的左式堆节点数为N,右路径长度为j的左式堆节点数为H，有 $i<=[lo{g}_{2}N+1]${向下取整} $j<=[lo{g}_{2}H+1]${向下取整} 故 $T[i,j]=O(lo{g}_{2}N+lo{g}_{2}H)=O(max(lo{g}_{2}N,lo{g}_{2}H))$

证毕

证法二：

$T[i,j]$的含义如证法一开头所述，这里i>=0,j>=0,T的边界条件和证法一所述相同，根据证法一开头的分析，我们可以认为L1右子树代表的左式堆LR1和左式堆R1合并以及左式堆L1和R1右子树代表的左式堆RR1合并的可能性相同，概率均为 $\frac{1}{2}$.从而可以认为 $T[i,j]=\frac{1}{2}(T[i-1,j]+T[i,j-1])+1$ $(1)$

这里1表示完成对$[i-1,j]$或$[i,j-1]$对应的合并操作后在$[i,j]$对应的合并操作中链接和调换子树位置所需的常量时间

于是 $T[i,j]=\frac{1}{2}T[i-1,j]+\frac{1}{2}T[i,j-1]+1=\frac{1}{2}T[i-1,j]+\frac{1}{2}(\frac{1}{2}T[i-1,j-1]+\frac{1}{2}T[i,j-2]+1)+1$

$=\frac{1}{2}T[i-1,j]+\frac{1}{2^2}T[i-1,j-1]+\frac{1}{2^2}T[i,j-2]+\frac{1}{2}+1$

$=\frac{1}{2}T[i-1,j]+\frac{1}{2^2}T[i-1,j-1]+\frac{1}{2^2}(\frac{1}{2}T[i-1,j-2]+\frac{1}{2}T[i,j-3]+1)+\frac{1}{2}+1$

$=\frac{1}{2}T[i-1,j]+\frac{1}{2^2}T[i-1,j-1]+\frac{1}{2^3}T[i-1,j-2]+\frac{1}{2^3}T[i,j-3]+\frac{1}{2^2}+\frac{1}{2}+1$

=-----

$=\frac{1}{2}T[i-1,j]+\frac{1}{2^2}T[i-1,j-1]+\frac{1}{2^3}T[i-1,j-2]+---+\frac{1}{2^{j+1}}T[i-1,0]+\frac{1}{2^{j+1}}T[i,-1]+\frac{1}{2^j}+---+\frac{1}{2}+1^{(2)}$

由上式知，倘若有 $T[i-1,0]---T[i-1,j-2]T[i-1,j-1]T[i-1,j]$ 均 $<=2i+1$

则 $T[i,j]=\frac{1}{2}T[i-1,j]+\frac{1}{2^2}T[i-1,j-1]+\frac{1}{2^3}T[i-1,j-2]+---+\frac{1}{2^{j+1}}T[i-1,0]$

$+\frac{1}{2^{j+1}}T[i,-1]+\frac{1}{2^j}+---+\frac{1}{2}+1$

$={\sum }_{k=0}^j\frac{1}{2^{k+1}}T[i-1,j-k]+\frac{1}{2^{j+1}}T[i,-1]+2-\frac{1}{2^j}$

$<=({\sum }_{k=0}^j\frac{1}{2^{k+1}}2i+1)+\frac{1}{2^{j+1}}+2-\frac{1}{2^j}$

$=(2i+1)(1-\frac{1}{2^{j+1}})+2-\frac{1}{2^{j+1}}$

$=2i+3-\frac{2i+2}{2^{j+1}}$

$=2i+3-\frac{i+1}{2^j}<=2i+3=2(i+1)+1$

而由边界条件知 $T[-1,0]T[-1,1]---T[-1,j]$ 均 $<=1=2\ast (-1+1)+1$

故由以上结论知 $T[0,0]T[0,1]---T[0,j]$ 均 $<=2(0+1)+1$

于是进一步可推得 $T[1,0]T[1,1]---T[1,j]$ 均 $<=2(1+1)+1$

这样层层向上递推，最终得到 $T[i,0]T[i,1]---T[i,j]均<=2(i+1)+1$

即 $T[i,j]<=2(i+1)+1$ $(3)$ $i>=0j>=0$

现在取$max(i,j)=n$,首先我们有 $T[0,-1]=T[1,-1]=---T[n,-1]=T[-1,0]=T[-1,1]=---=T[-1,n]=1$

$T[0,0]=T[0,0]T[1,1]=T[1,1]----T[n,n]=T[n,n]$

然后我们有 $T[0,1]=\frac{1}{2}T[0,0]+\frac{1}{2}T[-1,1]+1=\frac{1}{2}T[0,0]+\frac{1}{2}T[1,-1]+1=T[1,0]$

接着 $T[0,2]=\frac{1}{2}T[0,1]+\frac{1}{2}T[-1,2]+1=\frac{1}{2}T[1,0]+\frac{1}{2}T[2,-1]+1=T[2,0]$

按照这一步骤向前递推最终得到

$T[0,1]=T[1,0]T[0,2]=T[2,0]---T[0,n]=T[n,0]$$(4)$

利用(4)并采用和以上和相同的步骤又可得到 $T[1,2]=T[2,1]T[1,3]=T[3,1]----T[1,n]=T[n,1]$

继续递推，最终得到 $T[h,h+1]=T[h+1,h]T[h,h+2]=T[h+2,h]----T[h,n]=T[n,h]$

$0<=h<n$

综上,对任意 $0<=g<=n0<=p<=n$ 必有 $T[g,p]=T[p,g]$

而 $i>=0$ $j>=0$ $i<=max(i,j)=n$ $j<=max(i,j)=n$ 故 $T[i,j]=T[j,i]$

这样对 $T[j,i]$ 用(3)我们有 $T[j,i]<=2(j+1)+1$

从而 $T[i,j]=T[j,i]<=2(j+1)+1$

于是我们最终得到 $T[i,j]<=2(i+1)+1$ $T[i,j]<=2(j+1)+1$

从而得到时间界 $T[i,j]=O(i)$ $T[i,j]=O(j)$ $T[i,j]=O(max(i,j))$

即(由证法一结尾的说明) $T[i,j]=O(lo{g}_{2}N)$ $T[i,j]=O(lo{g}_{2}H)$ $T[i,j]=O(max(lo{g}_{2}H,lo{g}_{2}N))$

证完