最大子矩阵面积

http://blog.csdn.net/sunmenggmail/article/details/7709004

#include<iostream>
using namespace std;
int f[2001][2001];
int h[2001],l[2001],r[2001];
int n,i,j;
int main()
{
    cin>>n;
	int w;
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=n;j++)
        {
			cin>>w;
            f[i][j]=w;
        }
    }
    int ma=0;
    for(i=1;i<=n;i++)//一行一行处理，h[j]存储的是到第i行时，在第j列从i行开始往上的0的个数
    {
        for(j=1;j<=n;j++)
        {
            if(f[i][j]==1)
                h[j]++;
            else
                h[j]=0;
        }
        for(j=1;j<=n;j++)//l[]上一行的结果不影响这一行的
        {
            l[j]=j;//先假设第i行的第j个元素的左边界为j;
            while(l[j]>1&&h[j]<=h[l[j]-1])//将左边界向左扩展，但是要保持要扩展的边界l[j]-1的高度要大于等于j的高度
                l[j]=l[l[j]-1];
        }
        for(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(j=1;j<=n;j++)
        {
            if(ma<h[j]*(r[j]-l[j]+1))
                ma=h[j]*(r[j]-l[j]+1);
        }
    }
    cout<<ma<<endl;
    return 0;
}