BZOJ 1057 ZJOI 2007 바둑판 제작 DP+ 현수법

제목의 대의: 01로 이루어진 행렬을 제시하고 이 행렬의 최대 면적의 정사각형과 사각형을 묻는다. 그 중 임의의 사각형이 서로 인접한 것은 모두 다른 칸이다.
사고방식: 사실 모든 (i+j) & 1의 위치에서 숫자가 틀리거나 틀리면 0 또는 1의 최대 정사각형과 직사각형으로 변한다.첫 번째 질문은 물DP이고, 두 번째 질문은 단조로운 창고나 현수막으로 할 수 있다.다 잘 써요.
CODE:

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define MAX 2010
using namespace std;

int m,n;
int src[MAX][MAX];
int f[MAX][MAX];

int up[MAX][MAX],_left[MAX][MAX],_right[MAX][MAX];

int main()
{
	cin >> m >> n;
	for(int i = 1; i <= m; ++i)
		for(int j = 1; j <= n; ++j) {
			scanf("%d",&src[i][j]);
			if((i + j)&1)
				src[i][j] ^= 1;
		}
	int ans = 0;
	for(int i = 1; i <= m; ++i)
		for(int j = 1; j <= n; ++j)
			if(src[i][j]) {
				f[i][j] = min(f[i - 1][j - 1],min(f[i - 1][j],f[i][j - 1])) + 1;
				ans = max(ans,f[i][j]);
			}
	memset(f,0,sizeof(f));
	for(int i = 1; i <= m; ++i)
		for(int j = 1; j <= n; ++j)
			if(!src[i][j]) {
				f[i][j] = min(f[i - 1][j - 1],min(f[i - 1][j],f[i][j - 1])) + 1;
				ans = max(ans,f[i][j]);
			}
	cout << ans * ans << endl;
	ans = 0;
	for(int i = 1; i <= m; ++i) {
		for(int j = 1; j <= n; ++j)
			_left[i][j] = src[i][j] ? _left[i][j - 1] + 1:0;
		for(int j = n; j; --j)
			_right[i][j] = src[i][j] ? _right[i][j + 1] + 1:0;
	}
	for(int i = 1; i <= m; ++i)
		for(int j = 1; j <= n; ++j)
			if(src[i][j] && src[i - 1][j]) {
				up[i][j] = up[i - 1][j] + 1;
				_left[i][j] = min(_left[i][j],_left[i - 1][j]);
				_right[i][j] = min(_right[i][j],_right[i - 1][j]);
				ans = max(ans,(_left[i][j] + _right[i][j] - 1) * (up[i][j] + 1));
			}
	memset(_left,0,sizeof(_left));
	memset(_right,0,sizeof(_right));
	memset(up,0,sizeof(up));
	for(int i = 1; i <= m; ++i) {
		for(int j = 1; j <= n; ++j)
			_left[i][j] = src[i][j] ? 0:_left[i][j - 1] + 1;
		for(int j = n; j; --j)
			_right[i][j] = src[i][j] ? 0:_right[i][j + 1] + 1;
	}
	for(int i = 1; i <= m; ++i)
		for(int j = 1; j <= n; ++j)
			if(!src[i][j] && !src[i - 1][j]) {
				up[i][j] = up[i - 1][j] + 1;
				_left[i][j] = min(_left[i][j],_left[i - 1][j]);
				_right[i][j] = min(_right[i][j],_right[i - 1][j]);
				ans = max(ans,(_left[i][j] + _right[i][j] - 1) * (up[i][j] + 1));
			}
	cout << ans << endl;
	return 0;
}