（洛谷P1387）最大正方形

题目描述

二维动态规划，设$dp\left(i,j\right)$为以$\left(i,j\right)$为右下顶点的正方形的最大边长。则容易想到$dp\left(i,j\right)$仅与$dp\left(i-1,j\right)$、$dp\left(i,j-1\right)$和$dp\left(i-1,j-1\right)$有关。
而我们知道，倘若存在以$\left(i,j\right)$为右下顶点的正方形，则其边长必在$\left[1,min\left(dp\left(i-1,j\right),dp\left(i,j-1\right)\right)\right]$。确定了边长范围以后，自然可以就通过遍历$dp\left(i-m,j-m\right)$得到最优解了。

状态转移方程如下：
$$dp(i,j)=\{\begin{array}{ll}0& \text{if }p[i][j]==0\\ 1& \text{if }i==0orj==0andp[i][j]==1\\ max(m)& \text{else}formin\left[1,min\left(dp\left(i-1,j\right),dp\left(i,j-1\right)\right)\right]andp\left(i-m,j-m\right)==1\end{array}$$

潜在优化：

1. 占个坑先。。。

代码如下：

#include<iostream>
#include<algorithm>
using namespace std;
int r, c;
int p[100][100] = { 0 };
int dp[100][100] = { 0 };

int findans() {
	int result = 0;
	for (int i = 1; i < r; ++i) {
		for (int j = 1; j < c; ++j) {
			if (p[i][j]) {
				int subl = min(dp[i - 1][j], dp[i][j - 1]);
				for (int m = 1; m <= subl; ++m) {
					if (!dp[i - m][j - m]) {
						subl = m - 1;
						break;
					}
				}
				dp[i][j] = subl + 1;
				result = max(result, dp[i][j]);
			}
		}
	}
	//return (result - 1) * 4; //求周长
	return result; //求边长
}

int main(){
	cin >> r >> c;
	for (int i = 0; i < r; ++i) {
		for (int j = 0; j < c; ++j) { 
			cin >> p[i][j]; 
			if ((i == 0 || j == 0) && p[i][j]) dp[i][j] = 1;
		}
	}
	cout << findans() << endl;
    return 0;
}