单调栈的应用

最新推荐文章于 2022-07-31 11:25:10 发布

月光下的魔术师 

最新推荐文章于 2022-07-31 11:25:10 发布

阅读量246

点赞数

分类专栏：其它文章标签： ACM 单调栈最大子矩阵

本文链接：https://blog.csdn.net/weixin_42264485/article/details/106212639

版权

其它专栏收录该内容

25 篇文章 0 订阅

订阅专栏

直方图求最大矩阵

POJ2559

题解

记录每个条柱最左和最右的位置， $h*(right-left+1)$ ，即为它能覆盖大小。
单调栈递增维护长度。遍历数组，若目前的长度 $h[i]$ 大于等于栈顶长度，则把下标入栈。
当栈顶的长度大于现在的长度，弹出栈顶cur，显然 $right=i-1$ 。那怎么求left？如果left=栈顶的下标，可否？显然不可以。给出如下图，如果此时 $i=6$ ，那么栈顶下标为5，显然 $right=5$ ，但是left却为3，因为中间有比2和5都高的部分。所以真正的操作为，弹出栈顶元素，left为弹出后栈顶元素的下标+1。为了计算方便，数组开头和末尾都加一个0。

#include <iostream>
#include <cstdio>
#include <stack>
using namespace std;
typedef long long ll;
int const N = 1e5 + 10;
int h[N],n;
stack<int>st;
int main(){
    while(~scanf("%d",&n) && n){
        while(!st.empty())  st.pop();
        for(int i=1;i<=n;i++)
            scanf("%d",&h[i]);
        h[0] = 0;   //为了处理2,1,1的情况
        h[++n] = 0;
        ll ans = 0;
        for(int i=0;i<=n;i++){
            while(!st.empty() && h[st.top()] > h[i]){
                int cur = st.top();
                st.pop();
                int left = st.top() + 1;
                int right = i - 1;
                ans = max(ans,1ll*(right - left + 1) * h[cur]);   
            }
            st.push(i);
        }
        printf("%lld\n",ans);
    }
    return 0;
}

最大全1子矩阵

POJ3494

题解

同上，先进行预处理。将这个矩阵分为以每一行为x轴的n个柱状图，对于每一个等于1的点，它的高度都等于上一行同一列的点的高度加一。注意C++自带的stack会超时。

#include <iostream>
#include <cstdio>
#include <stack>
#include <cstring>
using namespace std;
int const N = 2000 + 10;
int n,m,a[N][N];
int st[N],top;
int main(){
    while(~scanf("%d%d",&m,&n)){ 
        memset(a,0,sizeof(a));
        for(int i=1;i<=m;i++){
            for(int j=1;j<=n;j++){
                scanf("%d",&a[i][j]);
                if(a[i][j]) a[i][j] = a[i-1][j] + 1;
            }
        }
        int ans = 0;
        for(int k=1;k<=m;k++){
            top = 0;
            a[k][0] = a[k][n+1] = 0;
            int res = 0;
            for(int i=0;i<=n+1;i++){
                while(top > 0 && a[k][st[top]] > a[k][i]){
                    int cur = st[top--];
                    int right = i - 1;
                    int left = st[top] + 1;
                    res = max(res,(right - left + 1) * a[k][cur]);
                }
                st[++top] = i;
            }
            ans = max(ans,res);
        }
        printf("%d\n",ans);
    }
    return 0;
}