题意：
给出一张大小为(ncdot n,nleq 50)的网格，现有一些点为黑点，有一些点为白点。
现在可以执行若干次操作：选定一个大小为(wcdot h)(自选)的矩形，将矩形中的所有点都变为白点，消耗代价为(max(w,h))。
问最少需要多少代价，使得最终网格中所有点为白点。
思路：
- 考虑(dp)，对于一个矩形范围来说，其代价为(max(w,h))，若要使答案更优，那么就得减少(max)。
- 假设现在考虑(dp[x_1][y_1][x_2][y_2])，若(x_2-x_1>y_2-y_1)，那么就考虑将这个矩形横着划分为两块，通过上下两个矩形合并来得到答案。
- 此时我们枚举行，若在矩形范围内存在一行为空行，此时合并则可能会使得答案更优。
这题和常规的网格中的(dp)稍微有点不同，有点类似于二维区间(dp)的感觉。想到将一个矩形横竖进行划分得到最优解基本这个题就出来了。
/*
* Author: heyuhhh
* Created Time: 2020/3/13 11:56:12
*/
#include <iostream>
#include <algorithm>
#include <cstring>
#include <vector>
#include <cmath>
#include <set>
#include <map>
#include <queue>
#include <iomanip>
#include <assert.h>
#define MP make_pair
#define fi first
#define se second
#define pb push_back
#define sz(x) (int)(x).size()
#define all(x) (x).begin(), (x).end()
#define INF 0x3f3f3f3f
#define Local
#ifdef Local
#define dbg(args...) do { cout << #args << " -> "; err(args); } while (0)
void err() { std::cout << '
'; }
template<typename T, typename...Args>
void err(T a, Args...args) { std::cout << a << ' '; err(args...); }
template <template<typename...> class T, typename t, typename... A>
void err(const T <t> &arg, const A&... args) {
for (auto &v : arg) std::cout << v << ' '; err(args...); }
#else
#define dbg(...)
#endif
using namespace std;
typedef long long ll;
typedef pair<int, int> pii;
//head
const int N = 50 + 5;
int n;
char s[N][N];
int sum[N][N], dp[N][N][N][N];
void run() {
cin >> n;
for(int i = 1; i <= n; i++) {
cin >> (s[i] + 1);
for(int j = 1; j <= n; j++) {
sum[i][j] = (s[i][j] == '#');
}
}
for(int i = 1; i <= n; i++) {
for(int j = 1; j <= n; j++) {
sum[i][j] += sum[i][j - 1];
}
}
for(int i = 1; i <= n; i++) {
for(int j = 1; j <= n; j++) {
sum[i][j] += sum[i - 1][j];
}
}
auto query = [&](int r1, int c1, int r2, int c2) {
return sum[r2][c2] + sum[r1 - 1][c1 - 1] - sum[r2][c1 - 1] - sum[r1 - 1][c2];
};
for(int x1 = 1; x1 <= n; x1++) {
for(int y1 = 1; y1 <= n; y1++) {
for(int x2 = x1; x2 <= n; x2++) {
for(int y2 = y1; y2 <= n; y2++) {
dp[x1][y1][x2][y2] = max(x2 - x1 + 1, y2 - y1 + 1);
if(x1 == x2 && y1 == y2) dp[x1][y1][x2][y2] = (s[x1][y1] == '#');
}
}
}
}
for(int k1 = 1; k1 <= n; k1++) {
for(int k2 = 1; k2 <= n; k2++) {
for(int x1 = 1; x1 <= n; x1++) {
for(int y1 = 1; y1 <= n; y1++) {
int x2 = x1 + k1 - 1, y2 = y1 + k2 - 1;
if(x2 > n || y2 > n) break;
int& now = dp[x1][y1][x2][y2];
if(k1 >= k2) {
for(int t = x1; t <= x2; t++) {
if(query(t, y1, t, y2) == 0) {
now = min(now, dp[x1][y1][t - 1][y2] + dp[t + 1][y1][x2][y2]);
}
}
} else {
for(int t = y1; t <= y2; t++) {
if(query(x1, t, x2, t) == 0) {
now = min(now, dp[x1][y1][x2][t - 1] + dp[x1][t + 1][x2][y2]);
}
}
}
}
}
}
}
cout << dp[1][1][n][n] << '
';
}
int main() {
ios::sync_with_stdio(false);
cin.tie(0); cout.tie(0);
cout << fixed << setprecision(20);
run();
return 0;
}