【以前的空间】bzoj [ZJOI2007]棋盘制作

具体可以去跪《浅谈用极大化思想解决最大子矩形问题》（p.s. 蒟蒻跪了还是很晕，不过想到之前usaco好像是最后一章的一道题……看了下代码顿然醒悟）

也就是如果用o（nm）的方法维护一个极大矩阵？其实很简单，按行处理，维护u[],l[],r[](向上，向左，向右）每一行先预处理处理出这行中每一个点向右能延展多长(right[j])，向左能延展多长(left[j])，然后再开始处理这一行的点，具体做法是比较这个点（i，j）与上行同一列的点（i-1，j）是否可以连在一起，可以的话，那么u[j]++，l[j]=min(left[j],l[j]),r[j]:=min(right[j],r[j]),如果不可以，那么u[j]=1,l[j]=left[j],r[j]=right[j]。记住u[],l[],r[]的作用，他们表示一个如果以这一行为底边，那么包括进j点的最大矩形的形状u[j]*(r[j]+l[j]-1)。然后好像就是这样了吧。

var
  map,f:array[..,..]of longint;
  left,right,u,l,r:array[..]of longint;
  i,j,k,n,m,ans:longint;

function min(x,y:longint):longint;
begin
  if x>y then exit(y);
  exit(x);
end;

function max(x,y:longint):longint;
begin
  if x>y then exit(x);
  exit(y);
end;

procedure work1;
var
  i,j,k:longint;
begin
  ans:=;
  for i:= to n do
    for j:= to m do begin
      if (i=)or(j=) then begin
        f[i,j]:=;
        continue;
      end;
      if (map[i-][j]=map[i][j])
        or (map[i,j-]=map[i,j])
          or (map[i-,j-]<>map[i][j]) then begin
            f[i,j]:=;
            continue;
          end;
      f[i,j]:=min(f[i-][j],min(f[i,j-],f[i-,j-]))+;
      if f[i,j]>ans then ans:=f[i,j];
    end;
  writeln(ans*ans);
end;

procedure work2;
var
  i,j,k:longint;
begin
  ans:=;
  for i:= to n do begin
    left[]:=;
    right[m]:=;
    for j:= to m do
      if map[i,j-]<>map[i,j] then left[j]:=left[j-]+ else left[j]:=;
    for j:=m- downto  do
      if map[i,j+]<>map[i,j] then right[j]:=right[j+]+ else right[j]:=;
    for j:= to m do begin
      if i= then begin
        u[j]:=;
        l[j]:=left[j];
        r[j]:=right[j];
      end
        else
          if map[i-,j]<>map[i,j] then begin
            inc(u[j]);
            l[j]:=min(l[j],left[j]);
            r[j]:=min(r[j],right[j]);
          end
            else begin
              u[j]:=;
              l[j]:=left[j];
              r[j]:=right[j];
            end;
      ans:=max(ans,(r[j]+l[j]-)*u[j]);
    end;
  end;
  writeln(ans);
end;

begin
  readln(n,m);
  for i:= to n do
    for j:= to m do
      read(map[i,j]);
  work1;
  work2;
end.

p.s. 其实看了丽洁姐的代码，发现自己真是傻掉了，两问明明是可以并在一起求得嘛，(r[j]+l[j]-1)、u[j]不就是矩形的长和高，那么两者中肯定存在一个短边构成的正方形……