51nod-1158 . 全是1的最大子矩阵 &&CODEVS-2491 玉蟾宫（枚举+单调栈 or dp）

http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1158

http://codevs.cn/problem/2491/

//这两个题都是求从给定矩阵找出一个子矩阵全部为1，并且面积最大，这题可以转换为最大矩形那题，枚举每一行，每行上方就是一个最大矩形问题，if(map[i][j]) f[j]++,else f[j]=0;

这样只要枚举+单调栈就能解决问题，还有需要注意在solve中改变了f数组，所以应该先拷贝到另外一个数组，在调用。并且拷贝不能在函数中，memcpy(ff,f,sizeof(f)),把f数组的值拷贝到ff数组中，并且f只能是数组，而不是参数，参数只是一个指针，而不是数组。注意注意，就是在这里犯错。

下面贴一份后一题的代码。

#include<cstdio>
#include<cstring>
#include<stack>
using namespace std;
int f[1001],ff[1001],map[1001][1001];
char s[2010];
int sum=0;
void solve(int f[],int m)
{
    int i,ans,temp;
    stack<int>s;
    ff[m]=-1;
   // for(i=0;i<=m;i++)
       // printf("%d\n",ff[i]);
    for(i=0;i<=m;i++)
    {
        if(s.empty()||ff[i]>ff[s.top()])
        {
            s.push(i);
        }
        else if(ff[i]<ff[s.top()])
        {
            while(!s.empty()&&ff[i]<ff[s.top()])
            {
                temp=(i-s.top())*ff[s.top()];
                if(temp>sum) sum=temp;
                ans=s.top();
                s.pop();
            }
        s.push(ans);
        ff[ans]=ff[i];
        }
    }
}
int main()
{
    int i,j,k,n,m;
    scanf("%d%d",&n,&m);
    getchar();
    for(i=0;i<n;i++)
    {
        j=0;
        gets(s);
        for(k=0;k<strlen(s);k++)
        {
            if(s[k]=='R') map[i][j++]=0;
            else if(s[k]=='F') map[i][j++]=1;
        }
    }
   /* for(i=0;i<n;i++)
    {
        for(j=0;j<m;j++)
        {
            printf("%d",map[i][j]);
        }
        printf("\n");
    }*/
    memset(f,0,sizeof(f));
    for(i=0;i<n;i++)
    {
        for(j=0;j<m;j++)
            if(map[i][j]) f[j]++;
            else f[j]=0;
        memcpy(ff,f,sizeof(f));
        solve(ff,m);
    }
    printf("%d\n",sum*3);
    return 0;
}

同样这题可以用dp写，思想跟上面一样。下面贴一份第一题的代码。

#include<cstdio>
#include<cstring>
const int N = 501;
int f[N],ff[N],l[N],r[N],map[N][N],sum=0;
void solve(int ff[],int m)
{
    int i;
    for(i=1;i<=m;i++)
    {
        l[i]=r[i]=i;
        //printf("%d\n",ff[i]);
    }
    ff[0]=ff[m+1]=-1;
    for(i=1;i<=m;i++)
    {   //预处理左边界
        while(ff[i]<=ff[l[i]-1])
            l[i]=l[l[i]-1];
    }
    for(i=m;i>=1;i--)
    {   //预处理右边界
        while(ff[i]<=ff[r[i]+1])
            r[i]=r[r[i]+1];
    }
    for(i=1;i<=m;i++) //遍历更新最大值
        if(ff[i]*(r[i]-l[i]+1)>sum)
        sum=ff[i]*(r[i]-l[i]+1);
}
int main()
{
    int n,m,i,j;
    scanf("%d%d",&m,&n);
    for(i=1;i<=n;i++)
        for(j=1;j<=m;j++)
        scanf("%d",&map[i][j]);
    memset(f,0,sizeof(f));
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=m;j++)
            if(map[i][j]) f[j]++;
            else f[j]=0;
       // for(j=1;j<=m;j++)
            //printf("%d\n",f[j]);
        memcpy(ff,f,sizeof(f));
        //for(j=1;j<=m;j++)
            //printf("%d\n",ff[j]);
        solve(ff,m);
    }
    printf("%d\n",sum);
    return 0;
}