子矩阵[DP]

传送门

暴力搜索选哪些行,再在选出来的行DP

最坏赋值度 O(C(n,n/2) * n^ 3)

---

#include<bits/stdc++.h>
#define N 20
#define inf 0x3fffffff
using namespace std;
int Map[N][N],tmp[N][N];
int n,m,r,c,ans=inf;
int cost[N] , val[N][N] , f[N][N];
int Dp(){
	memset(f,127,sizeof(f)); 
	memset(cost,0,sizeof(cost));
	memset(val,0,sizeof(val));
	for(int i=1;i<=m;i++)
		for(int j=i+1;j<=m;j++)
			for(int k=1;k<=r;k++)
				val[i][j] += abs(tmp[k][i]-tmp[k][j]);
	for(int i=1;i<=m;i++)
		for(int j=2;j<=r;j++)
			cost[i]+=abs(tmp[j][i]-tmp[j-1][i]);
	for(int i=1;i<=m;i++) f[i][1]=cost[i];
		
	for(int i=1;i<=m;i++)//第几列 
		for(int j=2;j<=min(i,c);j++)//个数 
			for(int k=1;k<i;k++)
				f[i][j] = min(f[k][j-1]+val[k][i]+cost[i],f[i][j]);
	for(int i=1;i<=m;i++) ans=min(ans,f[i][c]);
}
void dfs(int u,int x,int d,int s){
	if(u>n && x==d){
		int cnt=0;
		for(int i=1;s;s>>=1,i++)
			if(s&1) {cnt++; for(int j=1;j<=m;j++) tmp[cnt][j]=Map[i][j];}
		Dp(); return;
	}
	if(u>n) return;
	dfs(u+1,x+1,d,s|(1<<(u-1)));
	dfs(u+1,x,d,s);
}
int main(){
	scanf("%d%d%d%d",&n,&m,&r,&c);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			scanf("%d",&Map[i][j]);
	dfs(1,0,r,0); printf("%d",ans); return 0;
}