最大全0子矩阵问题

题意：给定一个01矩阵，求：

1.最大全0正方形面积

2.最大全0子矩阵面积
-----------------------------------------我是分割线-----------------------------------------

对于第一问，我们有O(nm)的动态规划解法：

设f[i][j]表示以i,j为右下角的最大正方形边长

if(a[i][j]==0) f[i][j]=min(f[i-1][j-1],min(f[i-1][j],f[i][j-1]))+1;

else f[i][j]=0;

正确性可以自己YY一下

--------------------------------------------------------------------------------------------------

对于第二问有许多O(nm)的奇怪解法，具体可以参照2003年福州第三中学王知昆的论文《浅析用极大化思想解决最大子矩阵问题》。

我只会悬线法，还有二维单调栈解法orz

令h[i][j]表示点(i,j)向上最大延伸高度，可以想象为一根在(i,j-h[i][j]+1)到(i,j)的一根悬线，可以证明，最大子矩阵是某一根悬线的向左右尽量扩张而形成的

令l[i][j]表示i,j悬线向左最大扩充的位置，r[i][j]表示i,j悬线向右最大扩充的位置，那么max(h[i][j]*(r[i][j]-l[i][j]+1))就是答案。

现在我们要做到O(nm)处理出h,l,r数组

h数组显然可以简单递推。

关于l[i][j]，设其初值为j，考虑能否扩张到j-1及其扩张范围，显然如果h[i][j]<=h[i][j-1]则可以，否则不行，如果可以那么l[i][j]=l[i][l[i][j-1]]，否则就此停止；

现在考虑能否再次扩张，再进行比较h值，将之前的j-1换成l[i][j]-1即可

r数组同理。

为什么这个递推是O(nm)的呢？

显然行与行无关，我们只考虑一行内；

把l数组看成树结构，则l的指向即为连向父节点的边；

如果扩张一次完成，则新节点(j)的深度为(j-1)深度+1，这种事件出现次数<=m

每次扩张的比较过程如果成功，则深度会至少减少1，每次减少深度事件复杂度为O(1)

深度减少事件发生次数<=深度++次数<=m

所以每行内递推复杂度为O(m)

总复杂度O(nm)征毕。

至此，整个问题解决。

因为行与行无关，l与r可以只记录本行内，而h数组可以用滚动数组，故只需开一维即可，意义不变。

代码：

#include<bits/stdc++.h>
using namespace std;
const int maxn = 2017;
int a[maxn][maxn];
int h[maxn],l[maxn],r[maxn];
int n,ans;
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	   for(int j=1;j<=n;j++)
	      scanf("%d",&a[i][j]);
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		   if(a[i][j]) h[j]=0;
		   else h[j]++;
		for(int j=1;j<=n;j++)
		{
			l[j]=j;
			while(l[j]>1&&h[j]<=h[l[j]-1]) l[j]=l[l[j]-1];
		}
		for(int j=n;j>=1;j--)
		{
			r[j]=j;
			while(r[j]<n&&h[j]<=h[r[j]+1]) r[j]=r[r[j]+1];
		}
		for(int j=1;j<=n;j++)
		   ans=max(h[j]*(r[j]-l[j]+1),ans);
	}
	cout<<ans;
	return 0;
}