2019牛客暑期多校训练营第八场A All-one Matrices（单调栈+不被覆盖的全1矩阵）

题意：一个01矩阵，求所有不被其他矩阵（满足该条件的）覆盖的全一子矩阵的个数

题解：单调栈求该高度左右能扩大的最优位置，前缀和求该高度向上能扩大的最大位置，再枚举每行每列。这样会出现被覆盖的情况，所以还需要一个单调栈用来去除被覆盖的情况：该单调栈维护一个单调不降的h[i][j],在判断该行该高度下的L[i],R[i]的下一行的同等区间是否也都是1，是则跳过，负责cnt++.

#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn = 3e3+10;
const int INF = 0x3f3f3f3f;
const ll MOD = 1e9 + 7;
int L[maxn][maxn];int R[maxn][maxn];
int h[maxn][maxn];//每一列的浅醉 
int a[maxn][maxn];
int sum[maxn][maxn];//每一行的前缀 
stack<int> s;
int main(){
	int n,m; string t; cin>>n>>m; 
	for(int i=1;i<=n;i++){
		cin>>t;
		for(int j=1;j<=m;j++){
			a[i][j]=t[j-1]-'0';
			sum[i][j]=sum[i][j-1]+a[i][j];
			if(a[i][j]) h[i][j]=h[i-1][j]+a[i][j];
		}
	}
	int ans=0;
	for(int i=1;i<=n;i++){
		while(!s.empty()) s.pop();
		for(int j=1;j<=m;j++){
			while(s.size()&&h[i][s.top()]>=h[i][j]) s.pop();
	 	    if(s.empty())  L[i][j]=0;
	 	    else           L[i][j]=s.top();
	 	    s.push(j);
		}
		while(!s.empty()) s.pop();
		for(int j=m;j>=1;j--){
			while(s.size()&&h[i][s.top()]>=h[i][j]) s.pop();
			if(s.empty())  R[i][j]=m+1;
			else           R[i][j]=s.top();
			s.push(j);
		}
		while(!s.empty()) s.pop();
		for(int j=1;j<=m;j++){
			if(h[i][j]==0){
				while(!s.empty()) s.pop();
			}else{
				while(s.size()&&s.top()>h[i][j]) s.pop(); //维护严格不降 
				if(s.empty()||h[i][j]!=s.top()){
					if((R[i][j]-1-L[i][j])!=(sum[i+1][R[i][j]-1]-sum[i+1][L[i][j]])){
						ans++;
					} 
			    }
			}
		    s.push(h[i][j]);
		}
	}
	cout<<ans<<endl;
	return 0;
}