Codevs 2491 玉蟾宫

题目：

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

大意：
给定R，F填充的矩阵，求全F最大子矩阵面积*3；

样例输入 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

算法：递推+单调栈
题解：
递推求数组dp[i][j]表示该点向上最多有几个连续的F；
以每一层为地平线，看成高为dp，宽为1的m个连续矩形，O(m)求最大子矩形面积;
一共n次，O(nm)

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<stack>
using namespace std;
const int N=1000+50;
struct node{
    int x,y;
};
stack<node> S;
long long ans,now;
int k,n,m;
int dp[N][N],a[N][N];
char s;
void init(){
    while(!S.empty()) S.pop();
    k=0;
    now=0;

}
void done(int x){
    init();
    S.push((node){dp[x][1],1});
    for(int i=2;i<=m;i++){
        k=0;
        while(1){
            if(S.empty()){
                S.push((node){dp[x][i],k+1});
                break;
            }
            node u=S.top();
            if(u.x<dp[x][i]){
                S.push((node){dp[x][i],k+1});
                break;
            }
            if(u.x>=dp[x][i]){
                now=u.x*(k+u.y);
                ans=max(ans,now);
                k+=u.y;
                S.pop();
            }
        }
    }
    k=0;
    while(!S.empty()){
        node u=S.top();
        S.pop();
        now=u.x*(k+u.y);
        ans=max(ans,now);
        k+=u.y;
    }
    return ;
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            cin>>s;
            if(s=='F') a[i][j]=1;
            //cin>>s;
        }
    }
    /*for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++)
        cout<<a[i][j]<<" ";
        cout<<endl;
    }
    cout<<endl<<endl;*/
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            dp[i][j]=a[i][j]*(dp[i-1][j]+1);
            //神奇的递推小技巧
        }
    }
    /*for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++)
        cout<<dp[i][j]<<" ";
        cout<<endl;
    }*/
    for(int i=1;i<=n;i++){
        done(i);
    }
    ans=ans*3;
    printf("%lld",ans);
}

redo
cnt循环掉了，最后再加一个高度为-1的，简化最后的处理。

#include<iostream>
#include<cstdio>
#include<stack>

using namespace std;
const int N=1000+50;
int n,m,ans,k;
char mp[N][N];
short int cnt[N];
struct node{
    int x,y;
};
stack<node>s;
void done(){
    while(!s.empty()) s.pop();
    for(int i=1;i<=m+1;i++){
        if(s.empty()||s.top().x<cnt[i]) {
            s.push((node){cnt[i],1});
        }
        else{
            k=0;
            while(s.size()&&s.top().x>=cnt[i]){
                node u=s.top();
                s.pop();
                ans=max(ans,u.x*(u.y+k));
                k+=u.y;
            }
            s.push((node){cnt[i],k+1});
        }
    }
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            cin>>mp[i][j];

        }
    }

    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            if(mp[i][j]=='R') {
                cnt[j]=0;
                continue;
            }
            if(mp[i-1][j]=='F') cnt[j]=cnt[j]+1;
            else cnt[j]=1;
        }
        cnt[m+1]=-1;
        done();
    }
    printf("%d",ans*3);
    return 0;
}