九度1497(最大矩形面积+DP)

题目描述：

在一个M * N的矩阵中，所有的元素只有0和1，从这个矩阵中找出一个面积最大的全1子矩阵，所谓最大是指元素1的个数最多。

输入：

输入可能包含多个测试样例。
对于每个测试案例，输入的第一行是两个整数m、n(1<=m、n<=1000)：代表将要输入的矩阵的大小。
矩阵共有m行，每行有n个整数，分别是0或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

 

//动态规划求最大矩形面积，分行处理，以每行每列为高，动态找在该行找左右的最大区间段，然后算出面积，记录最大值

#include<iostream>
#include<string>
#include<cstdio>
using namespace std;
#define MAX 1000+10

int height[MAX][MAX];
int right_dis[MAX];
int left_dis[MAX];


int main()
{
 int n,m,i,j,temp,max;
 while(cin>>n>>m)
 {
 // memset(height,0,sizeof(height));
  for(i=1;i<=m;i++)
   height[0][i]=0;

  for(i=1;i<=n;i++)
  {
   for(j=1;j<=m;j++)
   {
    scanf("%d",&temp);

    //求出每个位置的高度
    if(temp==1)
     height[i][j]=height[i-1][j]+1;
    else 
     height[i][j]=0;
   }
  }

  max=0;
  for(i=1;i<=n;i++)
  {
   for(j=1;j<=m;j++)
   {
    temp=j;
    while(temp>1&&height[i][j]<=height[i][temp-1])
     temp=left_dis[temp-1];
    left_dis[j]=temp;
   }

   for(j=m;j>=1;j--)
   {
    temp=j;
    while(temp<m&&height[i][j]<=height[i][temp+1])
     temp=right_dis[temp+1];
    right_dis[j]=temp;
   }

    for(j=1;j<=m;j++)
    {
     temp=(right_dis[j]-left_dis[j]+1)*height[i][j];
     if(temp>max)
      max=temp;
    }
  }

  cout<<max<<endl;
 }
 return 0;
}