HDU 6957 Maximal submatrix （单调栈，柱形图最大矩阵面积）

题目：Problem - 6957 (dingbacode.com)

        题目大意是给出一个矩阵，找出面积最大的子矩阵，这个子矩阵的每一列上的数字都要是非递减的。

        直接处理非常麻烦，所以先进行预处理，提前处理出每一列（第j列）上以第i个开始的非递减数列长度b[i][j],再用扫描线从1到n扫一遍，处理出以第i行为底边的最大矩阵面积，这样就把题目转变为求柱形图的最大矩阵面积了。

然后用单调栈维护答案即可。

代码：

#include<iostream>
#include<string.h>
using namespace std;
int a[2003][2003], b[2003][2003];
int d[4000006];
int len[4000006];
int main()
{
    int tt, n, m;
    int i, j;
    int h,last,cnt,x,ans;
    cin >> tt;
    while (tt--)
    {
        ans = 0;
        memset(d, 0, sizeof(d));
        memset(len, 0, sizeof(len));
        cin >> n >> m;
        for (i = 1;i <= n;i++)
            for (j = 1;j <= m;j++)scanf_s("%d", &a[i][j]);
        for (j = 1;j <= m;j++)
        {
            last = 0;h = 0;
            for (i = n;i >= 1;i--)
            {
                if (a[i][j] <= last)h++;
                else h = 1;
                b[i][j] = h;
                last = a[i][j];
            }
        }
        for (i = 1;i <= n;i++)
        {
            h = 1;
            d[1] = b[i][1];len[1] = 1;
            b[i][m + 1] = 0;
            for (j = 2;j <= m+1;j++)
            {
                x = b[i][j];
                if (x >= d[h])d[++h] = x, len[h] = 1;
                else
                {
                    cnt = 0;
                    while (d[h] > x)
                    {
                        cnt += len[h];
                        if (d[h] * cnt > ans)ans = d[h] * cnt;
                        d[h] = 0, len[h] = 0;
                        h--;
                    }
                    cnt++;
                    if (x * cnt > ans)ans = x * cnt;
                    d[++h] = x;len[h] = cnt;
                }
            }
        }
        printf("%d\n", ans);
    }
}