HDU6957 Maximal submatrix（DP）

题目链接

转化一下，标记a[i][j]>=a[i-1][j]的位置为1 ，代表可以取。然后求最大全1子矩阵面积。用单调栈/悬线法dp 都是N*M复杂度。
对于单调栈，我们用一个m+1 制造一个虚拟矩阵，为了防止RE以及清空栈。
每次出栈更新的是比当前连续1数目多的那个位置（相当于以第i行作为最底层，用h[j] 标记第j列中当前已经连续的高度）
更新范围为【栈顶，当前）以栈顶为高度的矩阵面积

ll a[2010][2010];
ll b[2010][2010];
ll h[2222];
ll que[22222];
signed main()
{


    ll t;
    read(t);
    while(t--)
    {
        ll n,m;
        read(n);
        read(m);

        for(int i=1; i<=n; i++)
        {
            for(int j=1; j<=m; j++)
            {
                read(a[i][j]);
            }

        }
        for(int i=1; i<=n; i++)
        {

            for(int j=1; j<=m; j++)
            {
                if(a[i][j]>=a[i-1][j])
                {
                    b[i][j]=1;

                }
                else
                {
                    b[i][j]=0;
                }
            }

        }

        for(int i=1; i<=m; i++)
        {
            h[i]=0;
        }
        int ans=0;
        for(int i=1; i<=n; i++)
        {
            for(int j=1; j<=m; j++)
            {
                if(b[i][j]==0)
                {
                    h[j]=1;
                }
                else
                {
                    h[j]++;
                }
            }
            // h  [j]  第j列   已经连续了  x
            int tot=0;
            h[m+1]=0;
            for(int j=1; j<=m+1; j++)
            {
                while(tot&&h[que[tot]]>h[j])
                {
                    ans=max(ans,(j-que[tot-1]-1)*h[que[tot]]);
                    tot--;
                }
                que[++tot]=j;
            }
        }
        printf("%lld\n",ans);
    }


}

感谢老铁送来的资料