伸展树时间复杂度分摊时间复杂度

声明：本文原创，转载请注明出处

前言

本文写的超级详细，仔细阅读一定会读懂~有疑问的同学欢迎留言！
我们都知道，splay tree可达到分摊O(logn)的时间复杂度，但这个复杂度是怎么得到的呢？下面我们通过引入势能的概念加以分析。

势能

节点v的势能定义为：以节点v为根的子树的大小取对数，即$\varphi$(v) = log(size(v))，本文里面的log默认以2为底。整棵树的势能是每个节点的势能之和$\varphi$=${\sum }_v\varphi (v)$

直观上来讲，一棵树越“倾斜”，势能越高；反之，一棵树越“平衡”，势能就越低。下面证明，最倾斜的情况下（树退化为链表），势能为O(logn)；最平衡的情况下（满树），势能为为O(n)。

• 链表情况：${\sum }_{v}log(size(v))$ = ${\sum }_{k=1}^{n}log(k)$ = $log{\prod }_{k=1}^{n}k$ = $logn!$ = $O(nlogn)$

• 满树情况：设树的高度为h，节点v的深度为d，则size(v) = $2^{d+1}-1$，和v深度相同的节点共有$2^{h-d}$个。${\sum }_{v}log(size(v))={\sum }_{d=0}^h{2}^{h-d}log(2^{d+1}-1)<{\sum }_{d=0}^h{2}^{h-d}log{2}^{d+1}={\sum }_{d=0}^h{2}^{h-d}(d+1)=2^{h+2}-h-3=O(n)$

核心公式$A_{(k)}$ = $T_{(k)}$ + $\Delta {\varphi }_{(k)}$

假设我们一共有m次操作，我们的目的是证明这m次操作花费的总时间T=O(mlogn)。我们接下来的分析围绕着一个公式：$A_{(k)}$ = $T_{(k)}$ + $\Delta {\varphi }_{(k)}$，k=1,2…m。其中k表示的是第k次操作。把这k个公式累加，$A={\sum }_{k=1}^m{A}_k$，$T={\sum }_{k=1}^m{T}_k$，得到公式：$A=T+{\varphi }_m-{\varphi }_0$.
下面我解释一下A,T,$\varphi$的含义：

• T：花费的时间；
• $\varphi$：整棵树的势能；
• A：“代价”（又称摊还代价）

所以$A_{(k)}$ = $T_{(k)}$ + $\Delta {\varphi }_{(k)}$这个公式的含义就是：每次操作的“代价”就是这次操作所花费的时间和势能增量之和。

我们已经知道，${\varphi }_m-{\varphi }_0=O(nlogn)$，代入$A=T+{\varphi }_m-{\varphi }_0$得到：$A=T+O(nlogn)$。现在我们的证明已经完成了一半。我们想要证明T=O(mlogn)，接下来只需证明A=O(mlogn)。实际上，我们可以证明$A_{(k)}=O(logn)$，然后对m个$A_{(k)}$求和即可。

证明$A_{(k)}=O(logn)$

求子步骤代价$A_{(k)}^i$

为了描述方便，设$rank(v)=log(size(v))$

假设我们要对节点v做splay操作，即把v从原位置伸展至根。整个splay操作可以被分割成多个子操作：除了最后一个子操作外，每个子操作把v提升两层，最后一个子操作有可能把v提升一层或者两层（与v的深度的奇偶有关）。每个子操作不外乎属于以下几种之一：zig,zag,zig-zig,zag-zag,zig-zag,zag-zig。下面我们分情况讨论：

• zig/zag
  我们详细说明zig，zag的情况和zig是对称的。
  
  $A_{(k)}^i=T_{(k)}^i+\Delta {\varphi }_{(k)}^i$ (1)

  在子操作i中，整棵树里只有v和root两个节点的势能可能变化。故有：

  $A_{(k)}^i=T_{(k)}^i+\Delta ran{k}_{(k)}^i(v)+\Delta ran{k}_{(k)}^i(root)$ (2)

  而root子树的size减小了，因此$\Delta {\varphi }_{(k)}^i(root)<0$，故：

  $A_{(k)}^i<T_{(k)}^i+\Delta ran{k}_{(k)}^i(v)$ (3)

  因为被提升点v只上升了一层，故$T_{(k)}^i=1$:

  $A_{(k)}^i<1+\Delta ran{k}_{(k)}^i(v)$ (4)

• zig-zag/zag-zig
  
  $A_{(k)}^i=T_{(k)}^i+\Delta {\varphi }_{(k)}^i$ (1)

  在子操作i中，整棵树里只有v,p,g三个节点的势能可能变化。故有：

  $A_{(k)}^i=T_{(k)}^i+\Delta ran{k}_{(k)}^i(v)+\Delta ran{k}_{(k)}^i(p)+\Delta ran{k}_{(k)}^i(g)$ (2)

  $A_{(k)}^i=T_{(k)}^i+ran{k}_{(k)}^i(v)-ran{k}_{(k)}^{i-1}(v)+ran{k}_{(k)}^i(p)-ran{k}_{(k)}^{i-1}(p)+ran{k}_{(k)}^i(g)-ran{k}_{(k)}^{i-1}(g)$ (3)

  容易看出$ran{k}_{(k)}^i(v)=ran{k}_{(k)}^{i-1}(g)$,上式化简为：

  $A_{(k)}^i=2-ran{k}_{(k)}^{i-1}(v)+ran{k}_{(k)}^i(p)-ran{k}_{(k)}^{i-1}(p)+ran{k}_{(k)}^i(g)$ (4)

  因为$ran{k}_{(k)}^{i-1}(p)>ran{k}_{(k)}^{i-1}(v)$:

  $A_{(k)}^i<2-2ran{k}_{(k)}^{i-1}(v)+ran{k}_{(k)}^i(p)+ran{k}_{(k)}^i(g)$ (5)

  因为：$ran{k}_{(k)}^i(p)+ran{k}_{(k)}^i(g)=log(size(p_i))+log(size(g_i))<2log((size(p_i)+size(g_i))/2)=2log((size(v_i)-1)/2)<2log(size(v_i))-2=2ran{k}_{(k)}^i(v)-2$

  $A_{(k)}^i<2ran{k}_{(k)}^i(v)-2ran{k}_{(k)}^{i-1}(v)$ (6)

• zig-zig/zag-zag
  
  $A_{(k)}^i=T_{(k)}^i+\Delta {\varphi }_{(k)}^i$ (1)

  在子操作i中，整棵树里只有v,p,g三个节点的势能可能变化。故有：

  $A_{(k)}^i=T_{(k)}^i+\Delta ran{k}_{(k)}^i(v)+\Delta ran{k}_{(k)}^i(p)+\Delta ran{k}_{(k)}^i(g)$ (2)

  $A_{(k)}^i=T_{(k)}^i+ran{k}_{(k)}^i(v)-ran{k}_{(k)}^{i-1}(v)+ran{k}_{(k)}^i(p)-ran{k}_{(k)}^{i-1}(p)+ran{k}_{(k)}^i(g)-ran{k}_{(k)}^{i-1}(g)$ (3)

  容易看出$ran{k}_{(k)}^i(v)=ran{k}_{(k)}^{i-1}(g)$,上式化简为：
  $A_{(k)}^i=2-ran{k}_{(k)}^{i-1}(v)+ran{k}_{(k)}^i(p)-ran{k}_{(k)}^{i-1}(p)+ran{k}_{(k)}^i(g)$ (4)

  因为$ran{k}_{(k)}^{i-1}(p)>ran{k}_{(k)}^{i-1}(v)$:
  $A_{(k)}^i<2-2ran{k}_{(k)}^{i-1}(v)+ran{k}_{(k)}^i(p)+ran{k}_{(k)}^i(g)$ (5)

  到（5）为止，zig-zig的证明过程均与zig-zag相同。
  $A_{(k)}^i<2-3ran{k}_{(k)}^{i-1}(v)+ran{k}_{(k)}^i(p)+ran{k}_{(k)}^i(g)+ran{k}_{(k)}^{i-1}(v)$ （6）

  因为$ran{k}_{(k)}^i(p)<ran{k}_{(k)}^i(v)$:
  $A_{(k)}^i<2-3ran{k}_{(k)}^{i-1}(v)+ran{k}_{(k)}^i(v)+ran{k}_{(k)}^i(g)+ran{k}_{(k)}^{i-1}(v)$ （7）

  因为$ran{k}_{(k)}^i(g)+ran{k}_{(k)}^{i-1}(v)=log(size(t2+t3+1))+log(size(t0+t1+1))<2log(size(t0+t1+t2+t3+2)/2)<2log(size(v_i))-2=2ran{k}_{(k)}^i(v)-2$：

  $A_{(k)}^i<3ran{k}_{(k)}^i(v)-3ran{k}_{(k)}^{i-1}$ (8)

总结一下，到目前为止，我们已经得到了三种子步骤的代价公式：

• zig/zag: $A_{(k)}^i<1+ran{k}_{(k)}^i(v)-ran{k}_{(k)}^{i-1}(v)$ （1）
• zig-zag/zag-zig:$A_{(k)}^i<2ran{k}_{(k)}^i(v)-2ran{k}_{(k)}^{i-1}(v)$ （2）
• zig-zig/zag-zag:$A_{(k)}^i<3ran{k}_{(k)}^i(v)-3ran{k}_{(k)}^{i-1}$ (3)

每个操作的代价$A_{(k)}={\sum }_i{A}_{(k)}^i$

我们已经得到了每个子步骤的“代价”，欲得到整个操作的代价，只需把这些子步骤的代价求和。

case 1

我们首先假设v在splay至根时，上升的层数是偶数层，也就是说，我们只用到了子步骤zig-zag/zag-zig/zig-zig/zag-zag（只使用了公式(2)和(3))。

对于(2),我们知道：$2ran{k}_{(k)}^i(v)-2ran{k}_{(k)}^{i-1}(v)<3ran{k}_{(k)}^i(v)-3ran{k}_{(k)}^{i-1}(v)$

因此，(2)和(3)都可以写为：$A_{(k)}^i<3ran{k}_{(k)}^i(v)-3ran{k}_{(k)}^{i-1}$

假设一共有x个子步骤，我们对这x个式子求和：$A_{(k)}={\sum }_{i=1}^x{A}_{(k)}^i<3ran{k}_{(k)}^x(v)-3ran{k}_{(k)}^0(v)<3ran{k}_{(k)}^x(v)=O(logn)$

case 2

v在splay至根时，上升的层数是奇数层。此时，我们用到了公式(1)(2)(3)。

同样地，我们可以把公式1放缩：$A_{(k)}^i<1+ran{k}_{(k)}^i(v)-ran{k}_{(k)}^{i-1}(v)<1+3ran{k}_{(k)}^i(v)-3ran{k}_{(k)}^{i-1}$

现在我们把这x个式子求和，注意(1)只使用了一次：

$A_{(k)}={\sum }_{i=1}^x{A}_{(k)}^i<1+3ran{k}_{(k)}^x(v)-3ran{k}_{(k)}^0(v)<3ran{k}_{(k)}^x(v)=O(logn)$

总结

全文围绕公式$A=T+\Delta \varphi$

我们先证明了$\varphi =O(nlogn)$,故$\Delta \varphi =O(nlogn)$;

我们证明了每个操作的代价$A_{(k)}=O(logn)$，对m个操作求和，总代价为$A=O(mlogn)$；

对于m>>n的情况,$T=A-\Delta \varphi =O(mlogn)$，每次操作的均摊时间复杂度为O(logn)