Asymptotic Notation and Recurrences

⚠ 转载请注明出处:作者:ZobinHuang,更新日期:Mar.3 2021

    注:在本节中我们将花比较多的篇幅来讨论分析算法所使用的数学工具。

1. 渐近符号

(1) 近似符号 O:近似上界

    数学定义O(g(n)) 表示了一个函数集合,这些函数 f(n) 满足:c>0,n0>0,使得 0f(n)cg(n) 对于所有 nn0 成立,即 O(g(n))={f(n):c>0,n0>0,使0f(n)cg(n)nn0}。当我们写出 f(n)=O(g(n)) 的时候,= 在这里的含义应该理解为 , 即等号左边是等号右边的一个子集。在含有渐近符号的等式中,应该意识到等号两边是 不等价 的。有些教材/论文中会写成 , 后面会介绍写成 = 的一些好处。
    加深理解:现在我们知道了,当等式中出现 O(g(n)) 这样的渐近符号时,应该把它理解为一个集合,而 = 应该理解为 。这里更准确的说法是 宏替换(Macro Substitution),即等式中的集合代表了一个在这个集合中的匿名函数(anonymous function),我们在对带有集合的等式做分析时,我们需要将等式中的集合符号替换为匿名函数,从而对等式的含义做出解释。下面我们通过一些等式的例子加深对于 符号 O 的理解。
        (a) 2n2=O(n3)̲
        2n2 属于 O(n3),用宏替换来说就是在集合 O(n3) 中存在一个函数 h(n),使得 2n2=h(n) 成立,注意此处的 = 是真正的相等关系。因此,当集合符号出现在等式右边时,它应该被解释为 “在集合中存在一个函数,使得等号条件成立”。
        (b) f(n)=n3+O(n2)̲
        在集合 O(n2) 中存在一个函数 h(n),使得 f(n)=n3+h(n) 等号条件成立。这个式子的意义在于,它对 f(n) 做出了更加精确的定义,我们当然可以简单地说 f(n) 是一个以 n3 为上界的一个函数,但是这样定义是模糊的,O(n2) 的加入说明 f(n) 除了一个 n3 的主项之外,还存在一个上界为 n2 的误差项。
        (c) n2+O(n)=O(n2)̲
        这个等式在等号左边出现了集合,基于我们把 = 当作 的理解方式,运用宏替换,其解释方法为:对于任何的 h1(n)O(n),在集合 O(n2) 中一定存在一个函数 h2(n),使得 n2+h1(n)=h2(n) 等号条件成立。因此,当集合符号出现在等式左边时,它应该被解释为 “任何在集合中的函数都能使得等号条件成立”。在 (c) 例子中在等号两侧都出现了集合符号,所以综合 (a) 例子的解释,我们可以得到,当等号两边都存在集合符号时,应该被解释为 “任何在等号左边集合中的函数,都一定存在一个属于等号右边集合的函数,使得等号条件成立”。

(2) 近似符号 Ω:近似下界

    数学定义:我们有了 O 来表示上界,同样地,我们也需要能表示下界的数学定义。这里我们给出数学符号 Ω 的数学定义:Ω(g(n)) 表示了一个函数集合,这些函数 f(n) 满足:c>0,n0>0,使得 0cg(n)f(n) 对于所有 nn0 成立,即 Ω(g(n))={f(n):c>0,n0>0,使0cg(n)f(n)nn0}
    加深理解:下面我们分析一个例子。
        (a) n=Ω(lgn)̲
        运用基本的逻辑理解,可以把这个等式解释为:只要 n 足够大,n 至少会是 Ω(lgn) 的常数倍。当然我们在真正解释这个式子的数学含义的时候还是得套用上面所用的宏替换的解释,即:在集合 Ω(lgn) 中存在一个函数,使得 n=Ω(lgn) 等号条件成立。

(3) 近似符号 Θ:近似相等

    当我们有了可以类比于 O 和可以类比于 Ω 后,我们希望也有可以类比于 = 的数学定义,这里我们给出数学符号 Θ 的数学定义:Θ(g(n))=O(g(n))Ω(g(n))。简单地理解就像是我们在 Lecture 1 中提出来的那样,对于 f(n)=Θ(g(n))g(n)= “把 f(n) 的低阶项和高阶项常数忽略的部分”,即系数为1的高阶项。例如:12n2-2n=Θ(n2)

(4) 严格符号 o:严格上界

    相比于近似符号 O(g(n)) 要求集合中的函数 f(n) 满足: “存在c>0n0>0 使得 0f(n)cg(n) 对于所有 n>n0 成立,严格上界符号 o(g(n)) 要求集合中的函数 f(n) 满足: 对于 “任意c>0,一定存在 n0>0 使得 0f(n)cg(n) 对于所有 n>n0 成立。可以发现这里对系数 c存在任意 的区别,对于 f(n)=O(g(n))O(g(n)) 要求你只要存在 c 使得在n 足够大之后 cg(n)f(n) 的上界就可以了,而对于 f(n)=o(g(n))o(g(n)) 则要求不管你在 g(n) 前面加多小的系数,cg(n) 一定得在 n 足够大之后继续超过 f(n),即要求 g(n)f(n) 的严格上界。例如我们分析 2n2=o(n3)(*) 这个式子,我们可以得到 cn32n2,故 n0=2c,即对于所有的 c,都一定有:对于所有 nn0cn32n2,因此 (*) 式是成立的,证毕。

(5) 严格符号 ω:严格下界

    同样地,我们定义严格下界符号 ω。严格下界符号 ω(g(n)) 要求集合中的函数 f(n) 满足: 对于 “任意c>0,一定存在 n0>0 使得 0cg(n)f(n) 对于所有 n>n0 成立。即不管你在 g(n) 前面加多大的系数,f(n) 一定得在 n 足够大之后继续超过 cg(n),即要求 g(n)f(n) 的严格下界。例如我们分析 n=ω(lgn)(*),我们可以得到 nclgn,这个式子的答案会是 n0=1+1c (怎么求的?),即对于所有的 c,都一定有:对于所有 n1+1cnclgn,因此 (*) 式是成立的,证毕。

2. 解递归式

    目前我们没有一种通用的方法用于解出递归式,我们只能运用各种各样的技巧来解递归式,就像求解积分那样。但是检查我们是否正确地解出递归式是简单的。下面我们介绍三种解递归式的技巧。

(1) 代换法(数学归纳法)

    我们看这样一个例子:T(n)=4T(n2)+n,求解这个递归式,即求出 T(n)
    我们计划使用数学归纳法的方式把 T(n) 求出来。我们初步猜测:

T(n)=O(n3) ...... ①

    插曲:

    注意在数学归纳法的使用过程中,我们不能让我们的式子带有集合符号(e.g. O, Ω, etc)。比如:n=O(1) (即函数 f(n)=n 的渐近上界是 O(1)),这是一个伪命题(如果这是真的,那么所有算法的复杂度都将是 O(1)),但是如果我们使用数学归纳法强加证明,我们会看到:

已知基本情况: 1=O(1),这是一个真命题
如果假设 n-1=O(1)
那么我们可以得到 n=(n-1)O(1)+(1)O(1)=O(1)

    但是很明显上面的出来的结论是错误的。证明过程中的错误就在于我们保留了 O(1)O(1)实际上就是一堆常数的集合。如果我们只是对常数做有限次叠加,那问题不大,它依然是常数。但是在我们的证明过程中我们将其进行了n次叠加,这时这个常数就是依赖于n变化的了。也就是说,对于一条有限长度的 O 关系链,首尾的包含关系是可以成立的,例如 n2=O(n2)=O(n3)=O(n4) 从而推出 n2=O(n4) 是允许的。但是一条无限长度的 O 关系链就没有上面这样的链式法则。


    因此,为了防止使用 O 符号导致的常数被忽略的易犯错误,我们在数学归纳法中必须使用常数把 ① 式展开。下面我们开始进行数学归纳法的证明。
    首先我们分析基本情况:

T(1)=Θ(1)cn3=c ...... ②

    只要 c 足够大,② 式是一定成立的。然后我们假设当 k<n 时,有:

T(k)ck3 (k<n) ...... ③

    然后我们证明 k=n 时的情况,有:

T(n)=4T(n2)+n4c(n2)3+n=c2n3+n ...... ④

    我们最终的证明目标是 T(n)cn3,因此我们进行构造:

T(n)=c2n3+n=cn3-(12cn3-n)cn3,当且仅当 12cn3-n0 ...... ⑤

    易证 ⑤ 式的这个条件是很容易达成的,因为这个 c 是可以由我们控制的。我们只要控制 c 足够大,使得 ② 式和 ⑤ 式的条件都成立,这个命题就是成立的。但是一般来说,算法的基本情况(i.e. n=1时的情况) 对后续并没有太大的影响,所以我们常常会忽略对基本情况的讨论。
    至此,我们根据这个递推式,证明了 T(n) 的上界就是 n3,但注意 n3 并不就一定是严格的上界,我们只是在一开始猜测 n3T(n) 的上界,然后将我们的猜测是用数学归纳法进行证明而已。

    下面,还是上述的递推式 T(n)=4T(n2)+n,我们尝试继续将上界压缩,证明 T(n)=O(n2),即证:T(n)cn2 ...... ①。
    首先我们假设当 k<n 时,有:

T(k)ck2 (k<n) ...... ②

    然后我们得到:

T(n)=4T(n2)+n=4c(n2)2+n=cn2+n ...... ③

    我们发现通过 ③ 式我们无法证明 T(n)cn2,因此无法完成归纳。此时回顾 ① 式,我们发现我们限制得太死,我们仅尝试用二次项 n2 去约束 T(n),但是实际上一次项也是允许出现在上界条件中的。我们可以稍微将条件放宽,将证明目标变成:

T(n)c1n2-c2n ...... ①'

    然后我们假设当 k<n 时,有:

T(k)c1k2-c2k (k<n) ...... ②'

    然后我们得到:

T(n)=4T(n2)+n=4(c1(n2)2-c2(n2))+n=c1n2-2c2n+n ...... ③'

    结合我们的证明目标 ①',我们构造:

T(n)=c1n2-2c2n+n=c1n2-c2n-(-1+c2)nc1n2-c2n
当且仅当 (-1+c2)n0 ...... ④'

    综上,我们证明了对于 c1>0c2>1T(n)=O(n2)。特别注意 c1 必须要足够大才能满足基本情况 (n=1的情况),这里不再证明。

(2) 递归树法 (Recursion-tree)

    我们看一个更复杂的例子:T(n)=T(n4)+T(n2)+n2,求解这个递归式。

    首先我们将递推式子写成如上所示的奇怪的树状形式,我们把递推项放到叶子节点上,然后继续将这些递推项递归地拆成这样的树状结构,我们可以得到:

    于是我们得到了最终会为分解一颗叶子节点为 Θ(1) 的树。注意上图有可能产生误解,这是棵树应该是一颗半满二叉树,即叶子节点是不在同一行上的,明显的树的左边被递归分解得比较快,值得注意。
    现在我们把重点放在如何计算出这棵树的和上面来。我们按行来看,如上图所示,我们将每一行的所有数相加完写在右边。仔细观察我们会发现这是一个几何级数(等比数列)的形式,因此我们可以得到:

T(n)=n2(1+516+25256+......)n2((516)0+(516)1+(516)2+......) ........ ①

    对于几何级数,有这样一个性质:

1+12+14+18...=2

    实际上上式代表的就是二进制数 1.111111......,当小数点后为无穷位时,这个数就等于2了。
    运用这个性质,我们继续处理我们的 ① 式,有:

T(n)n2((516)0+(516)1+(516)2+......)<2n2 ........ ②

    至此,我们证出了 T(n)=O(n2),同时从 T(n)=T(n4)+T(n2)+n2 也可以得到 T(n)=Ω(n2) 的结论。

(3) 主方法 (Master method)

    首先我们给出主方法的使用条件:主方法适用于形如:T(n)=aT(nb)+f(n),a1,b>1,f(n)n0n>n0f(n)>0 的递归式子。观察这个递归式子,其规律就在于 T(n) 包含的是 a 个相同的子问题,且这些子问题的输入规模都是 nb,即 T(nb),外加一个非递归项 f(n)a1 的条件保证了至少会递归一次,b>1 的条件保证了子问题的规模会比当前的问题小。
    只要满足主方法的使用条件,主方法就能够帮助我们求解递归式。下面我们首先给出主方法给出的solution,然后我们使用递归数直观但不严谨地证明为什么主方法是正确的,以便于我们理解主方法。
    主方法以非递归项 f(n) 作为分类依据,给出了三种情况下的solution:
    CASE 1. 当 ε>0 使得 f(n)=O(nlogba-ε),也即 f(n) 的增长速度以多项式级别地缓慢于 nlogba 时(受 n-ε 的控制),T(n)=Θ(nlogba)
    CASE 2. 当 k0 使得 f(n)=O(nlogbalgkn),也即 f(n) 的增长速度接近于 nlogba 时,T(n)=Θ(nlogbalgk+1n)
    CASE 3. 当 ① ε>0 使得 f(n)=O(nlogba+ε),也即 f(n) 的增长速度以多项式级别地快于 nlogba 时(受 nε 的控制),同时 ② c<1 使得 af(nb)cf(n),也即保证每次递归时分解出来的 a 个子问题 f(nb) 的规模一定比 f(n) 要小时,T(n)=Θ(f(n))

    下面我们使用递归树对主方法做一个直观的观察。

    首先我们可以得到递归树的高度为 h=logbn,因此可得叶子节点的个数为 #LEAVES=ah=alogbn=nlogba怎么转化的?
    我们首先分析 CASE 3 的情况。在 CASE 3 中,我们发现的是 f(n) 的增长速度以多项式级别地快于 nlogba,也即递归树的总和受第一层 f(n) 项的主导,因此我们有 T(n)=Θ(f(n))
    对于 CASE 1 的情况,我们发现的是 f(n) 的增长速度以多项式级别地缓慢于 nlogba,因此递归树的总和受最后一层叶子节点 nlogbaT(1) 的主导,因此我们有T(n)=Θ(nlogba)
    对于 CASE 2 的情况,我们发现的是 f(n) 的增长速度接近于 nlogba,我们假设 f(n)=nlogbalgkn,因此我们有 T(n)=f(n)h=nlogbalgknlogbn=Θ(nlogbalgk+1n)

Related