【codevs 2491】玉蟾宫（悬线法）

2491 玉蟾宫
时间限制: 1 s 空间限制: 64000 KB 题目等级 : 大师 Master
题目描述 Description
　　有一天，小猫rainbow和freda来到了湘西张家界的天门山玉蟾宫，玉蟾宫宫主蓝兔盛情地款待了它们，并赐予它们一片土地。

　　这片土地被分成N*M个格子，每个格子里写着’R’或者’F’，R代表这块土地被赐予了rainbow，F代表这块土地被赐予了freda。
　　现在freda要在这里卖萌。。。它要找一块矩形土地，要求这片土地都标着’F’并且面积最大。
　　但是rainbow和freda的OI水平都弱爆了，找不出这块土地，而蓝兔也想看freda卖萌（她显然是不会编程的……），所以它们决定，如果你找到的土地面积为S，它们每人给你S两银子。

输入描述 Input Description
　　第一行两个整数N,M，表示矩形土地有N行M列。
　　接下来N行，每行M个用空格隔开的字符’F’或’R’，描述了矩形土地。

输出描述 Output Description
　　输出一个整数，表示你能得到多少银子，即(3*最大’F’矩形土地面积)的值。

样例输入 Sample Input
5 6
R F F F F F
F F F F F F
R R R F F F
F F F F F F
F F F F F F

样例输出 Sample Output
45

数据范围及提示 Data Size & Hint
　　对于50%的数据，1<=N,M<=200
　　对于100%的数据，1<=N,M<=1000

【题解】最大子矩阵第二种算法模板题

普通版代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
int a[1010][1010],n,m,maxr,maxl,ans;
int h[1010][1010],l[1010][1010],r[1010][1010];
int main()
{
    int i,j;
    scanf("%d%d",&n,&m);
    for(i=1;i<=n;++i)
     for(j=1;j<=m;++j)
      {
        char c;
        cin>>c;
        if(c=='R') a[i][j]=1;
         else a[i][j]=0;
      }
    for(i=1;i<=m;++i) r[0][i]=m;
    for(i=1;i<=n;++i)
     {
        maxr=m; maxl=1;
        for(j=1;j<=m;++j)
         if(a[i][j])
          {
            maxl=j+1;
            h[i][j]=l[i][j]=0;
           } 
          else
           {
            h[i][j]=h[i-1][j]+1;
            l[i][j]=max(maxl,l[i-1][j]);
           }
        for(j=m;j>0;--j)
         if(a[i][j])
          {
            maxr=j-1;
            r[i][j]=m;
          }
          else
           {
            r[i][j]=min(maxr,r[i-1][j]);
            ans=max(ans,(r[i][j]-l[i][j]+1)*h[i][j]);
           }
     }
    printf("%d\n",3*ans);
    return 0;
}

滚动数组版代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
int a[1010][1010],n,m,maxr,maxl,ans;
int h[1010],r[1010],l[1010];
int main()
{
  int i,j;
  scanf("%d%d",&n,&m);
  for(i=1;i<=n;++i)
   for(j=1;j<=m;++j)
    {
        char c;
        cin>>c;
        if(c=='F') a[i][j]=0;
         else a[i][j]=1;
        }   
  for(i=1;i<=m;++i) r[i]=m;
  for(i=1;i<=n;++i)
   {
    maxl=1; maxr=m;
    for(j=1;j<=m;++j)
     if(a[i][j])
       {
        maxl=j+1;
        h[j]=l[j]=0;
         }
      else
       {
        h[j]++;
        l[j]=max(maxl,l[j]);
       }
    for(j=m;j>0;j--)
     if(a[i][j])
       {
        maxr=j-1;
        r[j]=m;
        } 
      else
       {
        r[j]=min(maxr,r[j]);
        ans=max(ans,(r[j]-l[j]+1)*h[j]);
       }
   }
   printf("%d\n",ans*3);
   return 0;
}