玉蟾宫（最大子矩形面积）codevs2491

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

Q

话说这是一道单调栈的题，求最大子矩形面积；
但暴力可做？？？
加几个优化，n^3能过；
可是我太弱，调了n年80分，最后wzhd’神牛orz一语道醒梦中人；
总之还是要勤练，代码可以不对，但思路不能出现误差；

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int n,m,dp[2001][2011],ma[2001][2011],ans=-1,d[2001],inf=1e9;
char s[2333];
void solve()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            scanf("%s",s);
            if(s[0]=='R') dp[i][j]=ma[i][j]=0;
            else dp[i][j]=ma[i][j]=1;
        }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            dp[i][j]=ma[i][j]*(dp[i-1][j]+1),d[i]=max(d[i],dp[i][j]);//这里很巧妙;
    for(int i=1;i<=n;i++) ans=max(ans,d[i]);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            if(!ma[i][j]) continue;
            if((m-j+1)*d[i]<=ans) break;
            int h=inf;
            for(int k=j;k<=m;k++)
            {
                h=min(dp[i][k],h);
                ans=max((k-j+1)*h,ans);//ans要在里面更新；
                if((m-j+1)*h<=ans) break;
            }
        }
    }
    cout<<ans*3;
}
int main()
{
    solve();
    return 0;
}

来个单调栈版本的
依矩形高度单调上升；

注意：
1.sq ：指栈顶元素高度的矩形最左边的坐标；
2.对于高度相同的矩形，只入栈一个即可，因为我们已经记录了坐标，不会丢失；
3.弹栈时，边弹边更新面积；

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<stack>
using namespace std;
const int MAXN=2001;
struct hh {int h,pos;};
int a[MAXN][MAXN],dp[MAXN][MAXN],n,m,b[MAXN],ans=-1,sq;
char ss[233];
stack<hh>s;
hh x,q;
void solve()
{
    cin>>n>>m;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            scanf("%s",ss);
            if(ss[0]=='F') a[i][j]=1;
        }
    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+1;j++)
        {
            q.h=dp[i][j],q.pos=j;
            if(s.empty()) s.push(q);
            else 
            {
                if(s.top().h<q.h) s.push(q);
                else if(s.top().h>q.h)
                {
                    sq=j;
                    while(!s.empty() && s.top().h>=q.h)
                    {
                        x=s.top(),s.pop();
                        int xx=(j-x.pos)*x.h;
                        ans=max(ans,xx),sq=x.pos;
                    }
                    q.pos=sq,s.push(q);
                }
            }
        }
    cout<<ans*3;
}
int main()
{
    solve();
    return 0;
}