cv2491 玉蟾宫（极大子矩阵）

题目描述 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

来源：Nescafe 20

分析：
算法二：悬线法求最大子矩阵

tip

ans*3

我发现了一个小规律：
如果给出了矩阵的具体情况，一般使用算法二
如果只给出了障碍点的坐标，一般使用算法一

这里写代码片
#include<cstdio>
#include<cstring>
#include<iostream>

using namespace std;

int la,ra,l[1010],r[1010],h[1010],mp[1010][1010];
int n,m,ans=0;

void doit()
{
    int i,j,k;
    for (i=1;i<=m;i++)
        h[i]=0,l[i]=1,r[i]=m;
    for (i=1;i<=n;i++)
    {
        la=0;ra=m+1;
        for (j=1;j<=m;j++)
            if (mp[i][j])
                h[j]=0,la=j,l[j]=1;
            else
                h[j]++,l[j]=max(l[j],la+1);
        for (j=m;j>=1;j--)
            if (mp[i][j])
                r[j]=m,ra=j;
            else
               r[j]=min(r[j],ra-1),
               ans=max(ans,h[j]*(r[j]-l[j]+1));
    }
    printf("%d",3*ans);
}

int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
        {
            char ch;
            cin>>ch;
            if (ch=='F') mp[i][j]=0;
            else mp[i][j]=1;
        }
    doit();
    return 0; 
}