2019暑假牛客第8场-A-All-one Matrices-单调栈+前缀和

题面：


思路：

因为左右上边界已经确定了，所以只要判断下边界能否拓展即可（如果不可以拓展就说明这是一个极大矩形，++ans）.
比赛的时候想跑完然后再判断那些是极大矩形，然后各种判断，最后不出所料的写崩了。
边跑便判断，还有就是跑完，然后再判断。
掌握局部信息就判断（这样可以避免无关信息的干扰），掌握全局信息之后再判断，就容易受各种信息的干扰，也就是随着时间的运行，无用信息将有用信息覆盖了。
还有就是这种判断全是1的矩阵，无论是面积的最大值，还是数量好像都可以用单调栈解决。

#include<bits/stdc++.h>
#define per(i,a,b) for(int i = (a);i <= (b);++i)
#define rep(i,a,b) for(int i = (a;i >= (b);--i))
using namespace std;
typedef long long LL;
const int maxn = 3e3 + 10;
int n = 0,m = 0;
int pre[maxn][maxn];
char s[maxn][maxn];
struct node{
	int h,l,r;
}hei[maxn][maxn];
struct Tree{
	int h,pos;
};
void init(){
	per(i,0,n){
		per(j,0,m){
			pre[i][j] = 0;
		}
	}
}
int val(char ch){
	return ch - '0';
}
void solve(){
	per(i,1,n){
		per(j,1,m){
			pre[i][j] = pre[i][j-1] + val(s[i][j]);//利用前缀和判断矩形（左右，上边界都已经确定）下边界是否可以拓展
		}
	}
	int ans = 0;
	per(i,1,n){
		int head = 1,tail = 0;
		Tree st[maxn];
		per(j,1,m+1){//m+1设置一个哨兵，将最后剩余的栈的元素弹出
			if(j != m+1){
				if(s[i][j] == '1'){
					hei[i][j].h = hei[i-1][j].h+1;
				}else{
					hei[i][j].h = 0;
				}
			}else{
				hei[i][j].h = 0;
			}
			hei[i][j].l = hei[i][j].r = j;
			while(head <= tail && hei[i][j].h <= st[tail].h){
				int pos = st[tail].pos;
				if(hei[i][j].h != st[tail].h){//高度相等的时候不能计算，因为此时还没有到达极大值，极大值由弹出他的高度决定
					int l = hei[i][pos].l,r = j-1;
					if(pre[i+1][r] - pre[i+1][l-1] != (r-l+1)){//不可以向下拓展
						++ans;
					}
				}
				hei[i][j].l = hei[i][pos].l;
				--tail;
			}
			st[++tail].h = hei[i][j].h;st[tail].pos = j;
		}
	}
	printf("%d\n",ans);
}
int main(){

	while(~scanf("%d %d",&n,&m)){
		init();
		per(i,1,n){
			scanf("%s",s[i]+1);
		}
		solve();
	}
	return 0;
}