poj3494 连续区域最大矩形面积 我要做的足够 强大

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<vector>
#include<string.h>
using namespace std;
const int maxn=2010;
int N,M;
int height[maxn];
int dis[maxn];
int ans;
int max(int a,int b)
{
    return a>b?a:b;
}
void get_sum(int n)
{
    vector<int>L(2*n,0);
    vector<int>R(2*n,0);
    height[0]=height[n+1]=-1;
    for(int i=0; i<=n; i++)
    {
        L[i]=i;
        R[i]=i;
    }
    for(int i=1; i<=n; i++)
        while(height[i]<=height[L[i]-1])
            L[i]=L[L[i]-1];
    for(int i=n; i>=1; i--)
        while(height[i]<=height[R[i]+1])
            R[i]=R[R[i]+1];

    for(int i=1; i<=n; i++)
        ans=max(ans,((dis[R[i]]-dis[L[i]]+1)*height[i]));
}
int main()
{
    int tem;
    while(scanf("%d%d",&N,&M)!=EOF)
    {
        ans=0;
        memset(dis,0,sizeof(dis));
        memset(height,0,sizeof(height));
        for(int i=1; i<=N; i++)
            dis[i]=1+dis[i-1];
        for(int i=1; i<=N; i++)
        {
            for(int j=1; j<=M; j++)
            {
                scanf("%d",&tem);
                if(tem)
                    height[j]++;
                else
                    height[j]=0;
            }
            get_sum(M);
        }
        printf("%d\n",ans);
    }
    return 0;
}