[hdu多校] Maximal submatrix

单调栈+悬线法+dp解决子矩阵问题，细节还是挺多的

#include<bits/stdc++.h>
using namespace std;
const int maxn=2e3+20;

int a[maxn][maxn];
int up[maxn][maxn],l[maxn],r[maxn];
int stk[maxn],top;
int T;

int main()
{
    cin>>T; 
    while(T--)
    {
    	int n,m;
    	scanf("%d%d",&n,&m);
    	memset(up,0,sizeof(int)*(n+1)*(m+1));
    	memset(l,0,sizeof(int)*m+5);
    	memset(l,0,sizeof(int)*n+5);
    	
		for(int i=1;i<=n;i++)  //dp出up数组 
    	{
        	for(int j=1;j<=m;j++)
        	{
            	scanf("%d",&a[i][j]);
            	up[i][j]=1;
            	if (i>1 && a[i][j]>=a[i-1][j]) up[i][j]=up[i-1][j]+1;
        	}
    	}
    	int ans=0;
    	for(int i=1;i<=n;i++)  //枚举矩阵的底边 
    	{
        	top=0; 
        	up[i][m+1]=0, up[i][0]=0;	//这边考虑a[i]中的数都相等，那么第一个比ai小的数就是右边界，所以把右边界设为0，左边界同理
        	for(int j=1; j<=m+1; j++)  //单调栈分别求出l,r数组 
        	{
				while(top && up[i][j]<up[i][stk[top]])
				{
					r[stk[top]]=j-1;
					top--;
				}
				stk[++top]=j;
        	}
        	top=0;
        	for(int j=m; j>=0; j--)  
        	{
				while(top && up[i][j]<up[i][stk[top]])
				{
					l[stk[top]]=j+1;
					top--;
				}
				stk[++top]=j;
        	}
        	l[1]=1, r[m]=m; 
        	for (int j=1; j<=m; j++)
        		ans=max(ans,up[i][j]*(r[j]-l[j]+1));
    	}
    	cout<<ans<<endl;
	}
    return 0;
}