求极大全1矩形个数

题意：求极大全1矩形个数

题解：用up[i][j]表示其向上连续1的个数。i枚举矩形的下边界，从前往后枚举每一列，用一个单调递增的栈维护up值，对于栈中的每个up值，都记录以它为最小值，最远的扩展到的左右边界。左边界是遍历到它时，如果有栈元素出栈，就是栈元素的左边界，否则就是它自己。右边界就是当遍历到某个元素j,栈顶元素要出栈时，右边界就是j-1。这样上下左都是不可扩展的，再看看它能否往下扩展，即它下一行的[左边界,右边界]是否全为1，如果不是，这个就是极大矩形。判断能否往下扩展用前缀和。ps：某个元素能否入栈的条件是栈空或者它比栈顶元素大。（如果和栈顶元素相等，就有可能会重复计算了）

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <stack>
using namespace std;
const int maxn = 3005;
char a[maxn][maxn];
int up[maxn][maxn],sum[maxn][maxn];
int l[maxn];
int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1; i<=n; i++)
    {
        scanf("%s",a[i]+1);
    }
    for(int i=1; i<=n; i++)
    {
        for(int j=1; j<=m+1; j++)
        {
            if(a[i][j] == '1')
            {
                up[i][j] = up[i-1][j]+1;
                sum[i][j] = sum[i][j-1]+1;
            }
            else
                sum[i][j] = sum[i][j-1];
        }
    }
    int ans = 0;
    for(int i=1; i<=n; i++)
    {
        stack<int>s;
        while(!s.empty())
            s.pop();
        for(int j=1; j<=m+1; j++)
        {
            l[j] = j;
            while(!s.empty() && up[i][s.top()]>up[i][j])
            {
                l[j] = l[s.top()];
                if(i==n || (sum[i+1][j-1]-sum[i+1][l[s.top()]-1]<j-l[s.top()]))
                    ans++;
                s.pop();
            }
            if(up[i][j] && (s.empty()||up[i][s.top()]<up[i][j]))
                s.push(j);
        }
    }
    printf("%d",ans);
}