Codeforces 1198D（DP）

题目链接：Rectangle Painting 1

题意：给定一个N*N的01矩阵（1<=n<=50），每次可以选择一个大小为h*w的矩形，并把该矩形内部1全变成0，花费为max（h，w），问将该01矩阵变成全0矩阵最小花费

题解:DP，dp[x1][y1][x2][y2]代表以(x1,y1)为左下角，（x2,y2）为右上角矩形的最小花费，因为数据n<=50,外面可以直接枚举h,w,x1,y1和每次转移的位置,而当前矩形的花费为直接选择即max（h,w）或者由两个更小的矩阵的花费拼凑起来，选花费最小的即可，因为可以横着或竖着分成两个矩阵，分别枚举下即可

不懂和细节见代码及注释：

#include<iostream>
#include<stack>
#include<list>
#include<set>
#include<vector>
#include<algorithm>
#include<math.h>
#include<numeric>
#include<map>
#include<cstring>
#include<queue>
#include<iomanip>
#include<cmath>
#include<queue>
#include <bitset>
#include<unordered_map>
	#ifndef local
	#define endl '\n'
#endif */
#define mkp make_pair
using namespace std;
using std::bitset;
typedef long long ll;
typedef long double ld;
const int inf=0x3f3f3f3f;
const ll MAXN=2e6+10;
const ll N=1e5+100;
const ll mod=1e9+7;
const ll hash_p1=1610612741;
const ll hash_p2=805306457;
const ll hash_p3=402653189;
//-----------------------------------------------------------------------------------------------------------------*/
// ll head[MAXN],net[MAXN],to[MAXN],edge[MAXN]/*流量*/,cost[MAXN]//费用;
/* 
void add(ll u,ll v,ll w,ll s){
	to[++cnt]=v;net[cnt]=head[u];edge[cnt]=w;cost[cnt]=s;head[u]=cnt;
	to[++cnt]=u;net[cnt]=head[v];edge[cnt]=0;cost[cnt]=-s;head[v]=cnt;
}
struct elemt{
	int p,v;
};
-----------------------------------
求[1,MAXN]组合式和逆元 
ll mi(ll a,ll b){
	ll res=1;
	while(b){
		if(b%2){
			res=res*a%mod;
		}	
		a=a*a%mod;
	}
	return res;
}
ll fac[MAXN],inv[MAXN]
fac[0]=1;inv[0]=1;
for(int i=1;i<=MAXN;i){
	fac[i]=(fac[i-1]*i)%mod;
	inv[i]=mi(fac[i],mod-2);
}
ll C(int m,int n){//组合式C(m,n); 
	if(!n){
		return 1;
	}
	return fac[m]*(inv[n]*inv[m*-n]%mod)%mod;
}
---------------------------------
 unordered_map<int,int>mp;
//优先队列默认小顶堆 ， greater<int> --小顶堆  less<int> --大顶堆  
priority_queue<elemt,vector<elemt>,comp>q;
struct comp{
	public:
		bool operator()(elemt v1,elemt v2){
			return v1.v<v2.v;
		}
};
	set<int>::iterator it=st.begin();
*/
// vector<vector<int>>edge; 二维虚拟储存坐标 
//-----------------------------------------------------------------------------------------------------------------*/
  //map<int,bool>mp[N]; 
int dp[60][60][60][60];//存矩形的左下角和右上角坐标，代表当前矩形所需最小花费
int main(){
/*cout<<setiosflags(ios::fixed)<<setprecision(8)<<ans<<endl;//输出ans（float）格式控制为8位小数（不含整数部分）*/
/*cout<<setprecision(8)<<ans<<endl;//输出ans（float）格式控制为8位小数（含整数部分）*/
	ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);//同步流
	int n ;
	cin>>n;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			char d;
			cin>>d;
			if(d=='#'){
				dp[i][j][i][j]=1;//当前最基础矩形需要1点花费，为后面dp做准备
			}
		}
	}
	for(int h=1;h<=n;h++){//枚举矩形高
		for(int w=1;w<=n;w++){//枚举矩形宽
			for(int i=1;i+h-1<=n;i++){//枚举左下角左边
				for(int j=1;j+w-1<=n;j++){
					if(h==1&&w==1){
						continue;
					}
					dp[i][j][i+h-1][j+w-1]=max(h,w);//不分割的花费
					for(int k=i+1;k<=i+h-1;k++){//横着分割，枚举分割点
						dp[i][j][i+h-1][j+w-1]=min(dp[i][j][i+h-1][j+w-1],dp[i][j][k-1][j+w-1]+dp[k][j][i+h-1][j+w-1]);
					}
					for(int k=j+1;k<=j+w-1;k++){//竖着分割，枚举分割点
						dp[i][j][i+h-1][j+w-1]=min(dp[i][j][i+h-1][j+w-1],dp[i][j][i+h-1][k-1]+dp[i][k][i+h-1][j+w-1]);
					}
					//cout<<"i="<<i<<" j="<<j<<" ix="<<i+h-1<<" jy="<<j+w-1<<" dp="<<dp[i][j][i+h-1][j+w-1]<<endl;
				}
			}
		}
	}
	cout<<dp[1][1][n][n]<<endl;
	return 0;
}