题目链接: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;
}