单调栈的应用与总结

背景问题

给定长度为n的数组\rm{A},其每个元素为非负整数,计算其所有连续子序列的最小值之和

问题分析

首先可以很直观的想到,只要遍历\rm{A}中的每个元素\rm{A[i]},计算以该元素为最后一个元素的全部连续子序列中最小值的和,记作\rm{sum[i]},最后将\rm{sum[0]}\rm{sum[n - 1]}相加就能得到所求的结果。因此解决问题的关键就在于,如何计算\rm{sum}数组?

算法设计

考虑对\rm{A[j]}计算\rm{sum[j]},有两种情况:1. 对任意的0 \leq i < j,都有\rm{A[i] \geq A[j]};2. 存在一个最大的i,满足0 \leq i < j,使得\rm{A[i] < A[j]}

  • 对任意的0 \leq i < j,都有\rm{A[i] \geq A[j]}显然\rm{sum[j] = A[j] \times (j + 1)}
  • 存在一个最大的i,满足0 \leq i < j,使得\rm{A[i] < A[j]}此时可以将以连续子序列分成两部分考虑,即包括\rm{A[i]}的连续子序列不包括\rm{A[i]}的连续子序列。两种情况下的具体分析如下,根据分析此时有\rm{sum[j] = sum[i] + A[j] \times (j - i)}
  • 包括\rm{A[i]}的连续子序列:这些连续子序列可以用\rm{A[k,j]}描述,其中0 \leq k \leq i,而\rm{A[k,j]}又可以分成\rm{A[k,i]}\rm{A[i + 1,j]}两个部分。根据i的选取规则,显然\rm{A[k,j]}的最小值只存在于\rm{A[k,i]}中,那么这一部分连续子序列的最小值之和其实就是\rm{sum[i]}
  • 不包括\rm{A[i]}的连续子序列:这些连续子序列可以用\rm{A[k,j]},其中i < k \leq j,根据i的选取规则,显然这一部分连续子序列的最小值之和就是\rm{A[j] \times (j - i)}

算法实现

根据算法分析,可知\rm{sum}的计算需要从前往后计算,且关键在于如何对任意的\rm{A[j]},找到符合条件的\rm{A[i]}。最简单的方法就是对于每个j,从j - 1开始向0逐个判断,这样做的时间复杂度为O(n^2)。实际上,只要利用一种数据结构,就能把这一操作的时间复杂度降至O(n),那就是单调栈。单调栈是一种特殊的栈结构,其栈内元素是单调递增或者单调递减的,例如单调递增的单调栈,其构建形式大致如下所示:

stack<int> stack_i;
array<int, 10> arr;
for(i = 0; i < 10; i++)
{
    // 不断弹出栈顶元素,直到栈顶元素指向的原始数组内元素小于待入栈元素指向的原始数组内元素
    while(!stack.empty() && arr[stack_i.top()] >= arr[i])
    {
        stack_i.pop();
    }

    // 待入栈元素入栈
    stack_i.push(i);
}

这种栈有一个特点,即对于每个\rm{arr[j]},在弹出栈顶元素之后,栈顶元素i若存在,则其指向的元素一定是原始数组中满足\rm{arr[i] < arr[j]}\rm{arr[i]}里最靠近\rm{arr[j]}的那个,若栈为空了,则说明\rm{arr[j]}是前j + 1个元素里最小的。也就是说,利用这一特点,可以在遍历一遍数组,且每个元素至多入栈出栈一次的情况下,对任意的\rm{A[j]},找到符合条件的\rm{A[i]}。并且由于每个元素至多入栈出栈一次,时间复杂度仅为O(n)。据此,可以写出时间复杂度为O(n)的完整算法:

int sumSubarrayMins(vector<int>& A) {
    auto size = A.size();
    vector<int> sums(size);
    stack<int> seqs;
    int i;
    int res;
    
    for(i = 0; i < size; i++)
    {
        while(!seqs.empty() && A[seqs.top()] >= A[i])
            seqs.pop();
        if(seqs.empty()) sums[i] = A[i] * (i + 1);
        else sums[i] = sums[seqs.top()] + A[i] * (i - seqs.top());
        seqs.push(i);
    }

    res = 0;
    for(i = 0; i < size; i++) res += sums[i];
    return (int)res;
}

单调栈的应用

从背景问题的算法设计和实现中可以看到,单调栈可以使得寻找数组中每个元素左边第一个小于它的元素位置的时间复杂度降低至O(n),类似地,以下几种问题都可以用单调栈解决:

  • 寻找数组中每个元素左边第一个小于它的元素位置:该算法在上面已经给出
  • 寻找数组中每个元素左边第一个大于它的元素位置:构建单调递增栈
stack<int> stack_i;
array<int, 10> arr;
array<int, 10> res;
for(i = 0; i < 10; i++)
{
    // 不断弹出栈顶元素,直到栈顶元素指向的原始数组内元素大于待入栈元素指向的原始数组内元素
    while(!stack.empty() && arr[stack_i.top()] <= arr[i]) stack_i.pop();
    
    // 保存arr[i]对应的索引res[i]
    if(stack_i.empty()) res[i] = -1;
    else res[i] = stack_i.top();

    // 待入栈元素入栈
    stack_i.push(i);
}
  • 寻找数组中每个元素右边第一个小于它的元素位置:构建单调递减栈
stack<int> stack_i;
array<int, 10> arr;
array<int, 10> res;
for(i = 10; i >= 0; i--)
{
    // 不断弹出栈顶元素,直到栈顶元素指向的原始数组内元素小于待入栈元素指向的原始数组内元素
    while(!stack.empty() && arr[stack_i.top()] >= arr[i]) stack_i.pop();
    
    // 保存arr[i]对应的索引res[i]
    if(stack_i.empty()) res[i] = 10;
    else res[i] = stack_i.top();

    // 待入栈元素入栈
    stack_i.push(i);
}
  • 寻找数组中每个元素右边第一个大于它的元素位置:构建单调递减栈
stack<int> stack_i;
array<int, 10> arr;
array<int, 10> res;
for(i = 10; i >= 0; i--)
{
    // 不断弹出栈顶元素,直到栈顶元素指向的原始数组内元素大于待入栈元素指向的原始数组内元素
    while(!stack.empty() && arr[stack_i.top()] <= arr[i]) stack_i.pop();
    
    // 保存arr[i]对应的索引res[i]
    if(stack_i.empty()) res[i] = 10;
    else res[i] = stack_i.top();

    // 待入栈元素入栈
    stack_i.push(i);
}

单调索引栈与单调元素栈

上述的四种情况里,单调栈里存储的都是原始数组中某个元素对应的索引,这种单调索引栈是单调递增还是单调递减只和索引的变化方向有关。单调索引栈用于需要索引信息的问题中。上述四类问题如果不需要找到位置,只需要找到每个元素对应的那个元素的值的话,则可以使用直接使用单调元素栈,省去从索引转换到元素值的步骤。对于某一类问题,使用索引栈和元素栈时,栈的增减情况可能会有不同:

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

推荐阅读更多精彩内容