单调栈和单调队列

单调栈:进栈元素单调递增(减)的栈,如果碰到比栈顶元素大的元素就进栈,否则不断把栈顶元素弹出直到栈顶元素小于等于要进栈的元素或者栈为空。
应用:
已知一个数列,求每一个元素后面第一个比它大的元素到它的距离,如果没有,输出0。此类问题仅和元素下标差以及元素的值有关,而且问题具有对下标和对值的单调性,可以用这种方法把O(N^2)的算法变成O(N)的算法。

#include<iostream>
#include<stack>
using namespace std;
int n;
int a[100010];
int res[100010];
typedef struct {
    int t;
    int pos;
}T;
stack<T> s;
T now, temp;
int main() {
    while (cin >> a[n++]);
    n--;
    for (int i = 0; i < n; i++) {
        now.pos = i;
        now.t = a[i];
        if (s.empty() || s.top().t > a[i]) {
            s.push(now);
        }
        else if (s.top().t == a[i]) continue;
        else {
            while (!s.empty() && s.top().t <= a[i]) {
                temp = s.top();
                s.pop();
                res[temp.pos] = i - temp.pos;
            }
            s.push(now);
        }
    }
    while (!s.empty()) {
        res[s.top().pos] = 0;
        s.pop();
    }
    for (int i = 0; i < n; i++) cout << res[i] << ' ';
    return 0;
}

例题:
最大长方形:http://poj.org/problem?id=2559(原题的网站上不去)
先看最大正方形:https://www.luogu.org/problem/P1387
已知一个0-1矩阵a[M][N],求只包含1的最大正方形的面积。
dp[i][j]为以(i,j)为右下角的最大正方形边长,易知第一行和第一列的0方格的dp值为0,1方格的dp值为1。
对于其他的方格,如果它是0方格,dp值为0;如果它是1方格,则考虑它的左边,上边和左上边,如果有任意一个方格为0,则它的dp值为1(只有它自己),否则它的dp值为相邻三个方格dp值的最小值再加1,因为作为最小值min,说明它左上方的方格的左上方至少有一个min*min的正方形,它上方和左方的方格一定能向上和向左扩展出min,再加上它自己,边长为min+1,即:
dp[i][j]=min(dp[i-1][j],dp[i][j-1],dp[i-1][j-1])+1
动态规划即可。

#include<iostream>
#include<algorithm>
using namespace std;
int h, w;
bool a[1405][1405];
int dp[1405][1405];
int res;
int main() {
    cin >> h >> w;
    for (int i = 0; i < h; i++) {
        for (int j = 0; j < w; j++) {
            cin >> a[i][j];
            a[i][j] = !a[i][j];
        }
    }
    for (int i = 0; i < h; i++) {
        if (!a[i][0]) dp[i][0] = 1;
    }
    for (int j = 0; j < w; j++) {
        if (!a[0][j]) dp[0][j] = 1;
    }
    for (int i = 1; i < h; i++) {
        for (int j = 1; j < w; j++) {
            if (a[i][j]) dp[i][j] = 0;
            else {
                dp[i][j] = min(dp[i - 1][j], min(dp[i][j - 1], dp[i - 1][j - 1])) + 1;
                res = max(res, dp[i][j]);
            }
        }
    }
    cout << res;
    return 0;
}

时间复杂度:O(MN)
最大长方形:把上题的正方形改成长方形。
动态规划不对了,因为长方形长宽不一样。
最朴素的想法:枚举所有点,找它们能构成的最大长方形面积的最大值。联想关灯问题,可以按顺序做,枚举每一行以[i,j]为底的高最短的长方形面积,对于高可以用预处理算出来,即从上往下算最多能连续多少个1,显然有dp[i]=dp[i-1]+1,用动态规划容易算得,问题转化为求直方图上的最大长方形面积,预处理的时间复杂度为O(MN)
然后遍历每一行,每一行再遍历每一组[i,j],去最小高乘宽度再取最大值,时间复杂度为O(MN^2)
下面用单调栈把复杂度降到O(MN)
栈储存直方图上长方形的高和左端的位置pos。
对每一行,从左往右枚举元素:
若栈为空或者栈顶元素的高小于要入栈元素的高,就入栈。
否则出栈到栈顶元素的高小于要入栈元素的高为止,设置pos为最后一个出栈元素的pos,入栈,在这个过程中,对于每一个出栈的元素,都计算它和已出栈元素所能形成的最大长方形的面积:(i-pos)*h,直到所有元素都入过栈,把栈内的元素清空。
这样实际上每一行只遍历了一遍,时间复杂度为O(MN)
部分代码(求柱状图的最大矩形,能过poj2559):

#pragma warning(disable:4996)
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<stack>
using namespace std;
typedef long long ll;
const int maxn = 100010;
int n;
ll a[maxn];
ll res;
typedef struct {
    ll pos;
    ll h;
}rect;
stack<rect> s;
int main() {
    rect now, temp;
    while (1) {
        scanf("%d", &n);
        if (n == 0) break;
        for (int i = 0; i < n; i++) scanf("%I64d", &a[i]);
        res = 0;
        a[n] = -1;
        for (int i = 0; i <= n; i++) {
            now.pos = i;
            now.h = a[i];
            if (s.empty() || a[i] > s.top().h) {
                s.push(now);
            }
            else if (a[i] == s.top().h) continue;
            else {
                temp.pos = i;
                while (!s.empty() && a[i] <= s.top().h) {
                    temp = s.top();
                    s.pop();
                    res = max(res,ll((i - temp.pos)*temp.h));
                }
                now.pos = temp.pos;
                s.push(now);
            }
        }
        printf("%I64d\n", res);
    }

    return 0;
}

单调队列:https://www.luogu.org/problem/P1886
连续移动区间的最值问题。
首先想到可以用树状数组,时间复杂度为O(NlogN),似乎很好,但是此题N=10^6NlogN=2*10^7,这样还是会超时。
1s的时间限制,O(N)对应10^6及以上,O(NlogN)对应10^5O(N^2)对应10^3
用单调队列可以把时间复杂度降到O(N)
如果队列为空,从队尾进队;否则比较要进队的元素和队尾元素,如果它比队尾元素大,就进队,否则出队直到队尾元素比它小,进队。如果要进队的元素和队首元素下标差大于等于k,队首元素从队首出队(相当于滑动)。
STL实现:

#include<iostream>
#include<deque>
using namespace std;
int n, k;
const int maxn = 1000010;
struct S {
    int pos;
    int val;
}a[maxn];
deque<S> q;
int main() {
    cin >> n >> k;
    for (int i = 1; i <= n; i++) {
        cin >> a[i].val;
        a[i].pos = i;
    }
    for (int i = 1; i <= n; i++) {
        while (!q.empty() && q.back().val >= a[i].val) q.pop_back();
        q.push_back(a[i]);
        if (a[i].pos - q.front().pos >= k) q.pop_front();
        if (i >= k) cout << q.front().val << ' ';
    }
    cout << endl;
    q.clear();
    for (int i = 1; i <= n; i++) {
        while (!q.empty() && q.back().val <= a[i].val) q.pop_back();
        q.push_back(a[i]);
        if (a[i].pos - q.front().pos >= k) q.pop_front();
        if (i >= k) cout << q.front().val << ' ';
    }
    return 0;
}

数组实现:一定要记住设tail=0,head=1,判断非空的条件是head<=tail,不这样的话会出现++tail首元素没读,tail++从队首出队错误。

#include<iostream>
using namespace std;
const int maxn = 1000010;
int n, k;
struct S {
    int pos;
    int val;
}a[maxn], q[maxn];
int head, tail;
int main() {
    cin >> n >> k;
    for (int i = 1; i <= n; i++) {
        cin >> a[i].val;
        a[i].pos = i;
    }
    head = 1;
    for (int i = 1; i <= n; i++) {
        while (head <= tail && q[tail].val >= a[i].val) tail--;
        q[++tail] = a[i];//注意是++tail,不是tail++
        if (a[i].pos - q[head].pos >= k) head++;
        if (i >= k) cout << q[head].val << ' ';
    }
    cout << endl;
    head = 1;
    tail = 0;
    for (int i = 1; i <= n; i++) {
        while (head <= tail && q[tail].val <= a[i].val) tail--;
        q[++tail] = a[i];
        if (a[i].pos - q[head].pos >= k) head++;
        if (i >= k) cout << q[head].val << ' ';
    }
    return 0;
}
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 159,835评论 4 364
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 67,598评论 1 295
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 109,569评论 0 244
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 44,159评论 0 213
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 52,533评论 3 287
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 40,710评论 1 222
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 31,923评论 2 313
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 30,674评论 0 203
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 34,421评论 1 246
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 30,622评论 2 245
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 32,115评论 1 260
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 28,428评论 2 254
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 33,114评论 3 238
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 26,097评论 0 8
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 26,875评论 0 197
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 35,753评论 2 276
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 35,649评论 2 271

推荐阅读更多精彩内容