UVALive 3029 City Game

题目大意就是要你找到一个矩阵的最大的全0子矩阵的面积,复杂度是O(n^2)的。

我在去年12月份在CODEVS上做过一道叫做"玉蟾宫"的题,一模一样……

当时照着题解打了一发只有90分,还以为是字符串读入的问题。结果刚刚翻了一下之前的代码。

竟然能有90分.jpg

既然是n^2,就一定有什么奇淫巧技。

首先有几个结论(现在做题怎么都先找结论):

1.答案矩阵不可能往上下左右拓展(显然)。

2.我们只枚举某一行作为最后一行是可行的(显然)。

3.答案一定在一个以(i,j)点到当前第j列能扩展到的最高的点(i',j)的高为(i-i'+1)的最宽的矩阵上。

即:确定一个点(i,j),能往上走就往上走,走到不能走之后往两边扩展,扩展不了就计算,答案是正确的。

证明的话,既然有结论1,则必定在最后一行有一个(i,j)只能扩展到这么多。然后往两边扩展,答案就出来了。

其实就有了一种很好写也可以过的O(n^2logn)的做法。

枚举点(i,j),二分找到上界,二分找到左边,二分找到右边。

用一个二维前缀和维护一下1的个数就可以实现O(1)的check了。这几个二分是并的,复杂度单独计算,所以是O(n^2logn),貌似也有人写。

但是UVALive好像有多组数据,也不知道能不能过去(只怕是不能的)。

然后来考虑怎么在O(n^2)的时间内解决它。

找到上界这个东西很简单,直接SB地DP(递推)就可以了。重点是维护左右扩展。

其实这个也很简单,从左往右枚举j,维护一下当前节点左边第一个1的位置,然后和上一行的取Max。右边类似,反过来即可。

这样就实现了O(n^2)的计算,可以通过此题了。

注意在不为0时的左边右边的赋值。因为它对下面不会有约束,要把左界设成0,右界设成m+1(就是不会影响的意思)。

思路还是很喵喵喵的。

#include    <iostream>
#include    <cstdio>
#include    <cstdlib>
#include    <algorithm>
#include    <vector>
#include    <cstring>
#include    <queue>
#include    <complex>
#include    <stack>
#define LL long long int
#define dob double
#define FILE "3029"
using namespace std;

const int N = 1010;
int n,m,map[N][N],up[N][N],le[N][N],ri[N][N],Ans;
int gi(){
  int x=0;char ch=getchar();
  while(ch>'9' || ch<'0')ch=getchar();
  while(ch>='0'&&ch<='9')x=x*10+ch-48,ch=getchar();
  return x;
}

int gec(){
  char ch=getchar();
  while(ch!='F' && ch!='R')ch=getchar();
  return ch=='F' ? 0:1;
}

inline void solve(){
  Ans=0;n=gi(),m=gi();
  for(int i=1;i<=n;++i)
    for(int j=1;j<=m;++j)
      map[i][j]=gec();
  for(int i=1;i<=n;++i){
    int ls=0,rs=m+1;
    for(int j=1;j<=m;++j)
      if(map[i][j]==1){le[i][j]=0;ls=j;up[i][j]=0;}
      else if(i==1){up[i][j]=1,le[i][j]=ls+1;}
      else up[i][j]=up[i-1][j]+1,le[i][j]=max(le[i-1][j],ls+1);
    for(int j=m;j>=1;--j)
      if(map[i][j]==1){ri[i][j]=m;rs=j;}
      else{
        if(i==1){ri[i][j]=rs-1;}
        else ri[i][j]=min(ri[i-1][j],rs-1);
        Ans=max(Ans,up[i][j]*(ri[i][j]-le[i][j]+1));
      }
  }
  printf("%d
",Ans*3);
}

int main(){
  freopen(FILE".in","r",stdin);
  freopen(FILE".out","w",stdout);
  int Case=gi();
  for(int i=1;i<=Case;++i)
    solve();
}
City Game
原文地址:https://www.cnblogs.com/fenghaoran/p/7650825.html