玉蟾宫(最大子矩形面积)codevs2491

来源:互联网 发布:淘宝卖家怎么赚淘金币 编辑:程序博客网 时间:2024/05/12 23:18

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