004【算法篇】分治法的时间复杂度分析

最简单的应用分治策略的算法是归并排序,下面我们先给出归并排序的伪代码:

MERGE(A,p,q,r)
    n1=q-p+1
    n2=r-q
    //let L[1..n1+1] and R[1..n2+1] be new arrays
    for i=1 to n1
        L[i]=A[p+i-1]
    for j=1 to n2
        R[j]=A[q+j]
    L[n1+1]=∞
    R[n2+1]=∞
    i=1
    j=1
    for k=p to r
        if L[i]<=R[j]
            A[k]=L[i]
            i=i+1
        else
            A[k]=R[j]
            j=j+1
            
MERGE_SORT(A,p,r)
    if p<r
        q=(p+r)/2
        MERGE_SORT(A,p,q)
        MERGE_SORT(A,q+1,r)
        MERGE(A,p,q,r)

从上述伪代码我们可知,整段程序无非做了这么几件事情:

  • 如果已经分解到最小颗粒n=1,那么直接返回即可,此时的时间复杂度𝛩(1)
  • 如果还不是最小颗粒,那么作折半分解,也即分解为L[1..n1+1]和R[1..n2+1]后进一步递归,此时的时间复杂度为2T(n/2);其中这里忽略了分解本身的耗时,仅统计了处理分解后排序耗时
  • 针对分解后的结果作一次归并,归并过程能够在线性时间内返回,所以此时的时间复杂度为𝛩(n)

综上我们可知,归并排序的时间复杂度为:
T(n)= \begin{cases} 𝛩(1),\quad若n =1 \\ 2T(n/2)+𝛩(n),\quad若n>1 \end{cases}

教科书上告诉我们,归并排序的时间复杂度T(n)=𝛩(nlgn),那么它是怎么样计算出来的呢?下面给出几种解递归式的方法。

代换法

当我们有一定编程经验之后,其实对于算法的时间复杂度分析会有一定的直觉,而且随着经验的增长这种直觉的准确性也会越来越高。所以我们不妨先假设直觉是对的,然后对其进行验证即可,这就是假设验证法的基本思想:先猜答案,再给出验证。

先来个示例:
T(n)=4T(\frac{n}{2})+n\tag{1}

针对这样的一个递归式,如何求其时间复杂度?

先一步,我们先猜测,猜测(1)式的时间复杂度为𝛰(n^2)。

第二步,我们验证假设。如果我们的猜测是对的,那么根据上一讲的内容,针对(1)式我们就马上有如下的结论:
T(k) \leq ck^2, \quad k \lt n\tag{2}
于是我们对(1)式进行如下论证:
\begin{align} T(n) & =4T(\frac{n}{2})+n\\ & \leq 4c*(n/2)^2+n\\ & = cn^2+n\\ & = cn^2-(-n)\tag{3} \end{align}
好了,现在尴尬的事情出来了,要证明上式的最终结果小于cn^2,就必须要求(-n)大于0,很显然这是不可能的。

难道我们的假设出了问题,直觉不准了?不应该呀……

问题出在(2)式上面,我们只考虑了高阶项,却忽略了低阶项对结果的影响,所以我们对(2)式进行改造,应是如下的结论:
T(k) \leq c_1k^2-c_2k,\quad k \lt n\tag{4}
我们应用(4)式再对(1)式进行改写,得到如下结果:
\begin{align} T(n) & =4T(\frac{n}{2})+n\\ & \leq 4[c_1*(n/2)^2-c_2*(n/2)]+n\\ & = c_1n^2-2c_2n+n\\ & = c_1n^2-c_2n-[(c_2-1)n]\tag{5}\\ & \leq c_1n^2-c_2n\\ \end{align}

很明显,为使(5)式中的(c_2-1)n非负,只需要c_2\geq1即可,显然这是可以的。

因此我们通过上述公式,严格论证了(1)式的时间复杂度T(n)=𝛰(n^2)针对所有的c_1\geq c_2c_2\geq1有效。

递归树法

代换法的验证过程较为严谨,这里介绍一种我喜欢的方式,作树形图求解。一般的套路是先作图求解,而后用代换法进行验证,如果对自己比较自信的话,可以忽略验证

我们再来举个稍微复杂点的递归例子:
T(n)=T(\frac{n}{4})+T(\frac{n}{2})+n^2\tag{6}
如何求(6)式的时间复杂度?

话不多说,直接上图说话:

image-20200129190348513.png

是不是突然感觉原来如此简单?

另外最底层为什么是𝛩(1)?因为归并算法最终分解之后都是单个元素,其时间复杂度自然是𝛩(1)

我们继续求解,如下图所示:

image-20200129190440990.png

于是我们得到这样一个等比数列:
\begin{align} T(n)&=T(\frac{n}{4})+T(\frac{n}{2})+n^2\\ &=n^2+\frac{5}{16}n^2+\frac{25}{256}n^2+...+\frac{5^k}{16^k}n^2+...\\ &=\frac{1-(\frac{5}{16})^k}{1-\frac{5}{16}}*n^2\\ &\leq 2n^2\\ &=𝛰(n^2)=𝛺(n^2) \end{align}

于是问题得以求解,为求稳妥可应用代换法进行验证,此处略过。

主定理法

最后介绍一种主定理法,其根本是递归树法的应用,它有诸多限制,只能应用于形如下式的递归式上:
T(n)=aT(\frac{n}{b})+f(n),\quad a\geq1,b\gt1,f(n)渐近趋正 \tag{7}

应用主定理,我们有三种情景:

  • 当f(n)<n^{log_ba}时,T(n)=𝛩(n^{log_ba})
  • 当f(n)=n^{log_ba}时,T(n)=𝛩(n^{log_ba}*lg^{k+1}n)
  • 当f(n)>n^{log_ba},同时af(\frac{n}{b})\leq(1-𝜀)+f(n)时,T(n)=𝛩(f(n)),原因是代价逐级递减之后,f(n)开始占主导地位

这里就不详细证明了,其实可以根据上述递归树法进行论证。其中,n^{log_ba}表示递归树中的叶节点个数,其中递归树的高度为log_bn

应用此主定理,我们立马可以得知(1)式符合情况1,其时间复杂度T(n)=𝛩(n^2)

所以,现在大家知道为什么归并排序的时间复杂度T(n)=𝛩(nlgn)了吗?

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 218,036评论 6 506
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 93,046评论 3 395
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 164,411评论 0 354
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,622评论 1 293
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,661评论 6 392
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,521评论 1 304
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,288评论 3 418
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,200评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,644评论 1 314
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,837评论 3 336
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,953评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,673评论 5 346
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,281评论 3 329
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,889评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,011评论 1 269
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,119评论 3 370
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,901评论 2 355

推荐阅读更多精彩内容