单调栈学习笔记

单调栈的定义：

设st为一个单调栈，那么st里的每一个元素都满足优先序大的较小（或较大）
用c++语言写成就是：

stack<int> st;
/*
此处进栈操作
*/
int a=st.top(),b;
st.pop();b=st.top();
满足：a<b

模板题：Acwing 830.单调栈
代码：

#include<iostream>
#include<stack>
#include<queue>
using namespace std;
int n,a;
stack<int> st;
queue<int> q;
int main(){
    cin>>n;
	st.push(-1);
    while(n--){
       cin>>a;
       while(st.top()>=a&&st.size()!=1) st.pop();
       cout<<st.top()<<" ";
       st.push(a);
    }
    return 0;
}

下面是一些应用：

1. 131. 直方图中最大的矩形（Acwing）

先贴代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int a[1000010],s[100010],w[100010],top;
ll ans=0,n;
int main(){
    scanf("%d",&n);
    while(n){
        ans=0;
        for(int i=1;i<=n;i++) scanf("%d",&a[i]);
        a[n+1]=0;
        memset(s,0,sizeof(s));
        memset(w,0,sizeof(w));
        for(int i=1;i<=n+1;i++){
            if(s[top]<=a[i]){
                s[++top]=a[i];w[top]=1;
            }else{
                ll width=0;
                while(s[top]>a[i]){
                    width+=w[top];
                    ans=max(ans,width*s[top]);
                    top--;
                }
                s[++top]=a[i];
                w[top]=width+1;
            }
        }
        printf("%lld\n",ans);
        scanf("%d",&n);
    }
}

这里用单调栈来维护目前高度，如果当前入栈小于栈顶高度，那么就出栈并处理
代码块：

               ll width=0;
                while(s[top]>a[i]){
                    width+=w[top];
                    ans=max(ans,width*s[top]);
                    top--;
                }
                s[++top]=a[i];
                w[top]=width+1;
            }

其中为了方便，用w来储存宽度，width记录当前宽度

2. 152.城市游戏（Acwing）

跟第一题一样的思路，用单调栈来维护：

#include<iostream>
using namespace std;
int n,m;
char Map[1001][1001];
int h[1001][1001],ans;
int main(){
    cin>>n>>m;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++) cin>>Map[i][j];
    for(int i=1;i<=n;i++){
        int st[1001],w[1001],top=0;
        for(int j=1;j<=m;j++){
            if(Map[i][j]=='F') h[i][j]=h[i-1][j]+1;
            if(st[top]<=h[i][j]){
                st[++top]=h[i][j];w[top]=1;
            }else{
                int width=0;
                while(st[top]>h[i][j]&&top>0){
                    width+=w[top];
                    ans=max(ans,width*st[top]);top--;
                }
                st[++top]=h[i][j];w[top]=width+1;
            }
        }
        int width=0;
        while(top>0){
            width+=w[top];
            ans=max(ans,width*st[top]);top--;
        }
    }
    cout<<ans*3;
    return 0;
}

简单来说就是枚举地基（把直方图看成房子），用h储存高度
其他方面是一样的
同理题：Largest Submatrix of All 1’s

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
int h[2001][2001],a[2001][2001],n,m,ans;
int main(){
	while(~scanf("%d%d",&n,&m)){
		memset(h,0,sizeof(h));
		ans=0;
		for(int i=1;i<=n;i++)
			for(int j=1;j<=m;j++) scanf("%d",&a[i][j]);
		for(int i=1;i<=n;i++){
			int st[10001],w[10001],top=0;
			for(int j=1;j<=m;j++){
				if(a[i][j]==1) h[i][j]=h[i-1][j]+1;
				if(h[i][j]>=st[top]){
					st[++top]=h[i][j];
					w[top]=1;
				}else{
					int width=0;
					while(st[top]>h[i][j]){
						width+=w[top];
						ans=max(ans,width*st[top]);top--;
					}
					st[++top]=h[i][j];w[top]=width+1;
				}
			
			}
			int width=0;
			while(top>0){
				width+=w[top];
				ans=max(ans,width*st[top]);top--;
			}
		}
		printf("%d\n",ans);
	}
	return 0;
}