面積最大の全1マトリクス

タイトルの説明:
1つのM*Nの行列の中で、すべての要素は0と1だけで、この行列の中から1つの面積の最大の全1サブ行列を探し出して、いわゆる最大は要素1の個数が最も多いことを指します.
入力:
入力には、複数のテストサンプルが含まれる場合があります.各テストケースについて、入力される最初の行は2つの整数m、n(1<=m、n<=1000):入力するマトリクスのサイズを表します.行列にはm行があり、各行にはn個の整数があり、それぞれ0または1であり、隣接する2つの数の間に厳密に1つのスペースで区切られている.
出力:
各テストケースに対応して、出力マトリクスの中で面積が最も大きい全1サブマトリクスの要素数.
サンプル入力:

2 2
0 0
0 0
4 4
0 0 0 0
0 1 1 0
0 1 1 0
0 0 0 0

サンプル出力:

0
4

問題解決の考え方:
最初の行から1行ずつ下に走査し、各列の現在の1の個数(0に当たるとクリアされ、高さと理解され、h[i])を記録し、各列の高さに合致する矩形の幅(j列目では幅r[j]-l[j]+1)を計算し、最後にすべての行を遍歴した最大面積が答えである.
ACコード:

#include<cstdio>
#include<cstring>

using namespace std;

int matrix[1005][1005];
int h[1005];
int r[1005];
int l[1005];

int main(int argc,char *argv[])
{
	int n,m;
	while(scanf("%d%d",&n,&m)!=EOF)
	{
		int i,j,k;
		for(i=1;i<=n;i++)
		{
			for(j=1;j<=m;j++)
			{
				scanf("%d",&matrix[i][j]);
			}
		}
		memset(h,0,sizeof(h));
		int ans=0;
		for(i=1;i<=n;i++)
		{
			for(j=1;j<=m;j++)
			{
				if(matrix[i][j]==1)
					h[j]++;
				else
					h[j]=0;
			}
			h[0]=h[m+1]=-1;
			for(j=1;j<=m;j++)//l[i] h[i] h h[i] 
			{             // l[i] i h >=h[i],r[i] 
				k=j;
				while(h[j]<=h[k-1])
					k=l[k-1];
				l[j]=k;
			}
			for(j=m;j>=1;j--)//r[i] h[i] h h[i] 
			{
				k=j;
				while(h[j]<=h[k+1])
					k=r[k+1];
				r[j]=k;
			}
			for(j=1;j<=m;j++)
			{
				if(ans<h[j]*(r[j]-l[j]+1))
					ans=h[j]*(r[j]-l[j]+1);
			}
		}
		printf("%d
",ans);
	}

	return 0;
}