[NOIP2014 普及组] 子矩阵
题目背景
NOIP2014 普及组 T4
题目描述
给出如下定义:
- 子矩阵:从一个矩阵当中选取某些行和某些列交叉位置所组成的新矩阵(保持行与列的相对顺序)被称为原矩阵的一个子矩阵。
例如,下面左图中选取第 2 2 2、 4 4 4行和第 2 2 2、 4 4 4、 5 5 5列交叉位置的元素得到一个 2 × 3 2 \times 3 2×3的子矩阵如右图所示。
9 3 3 3 9
9 4 8 7 4
1 7 4 6 6
6 8 5 6 9
7 4 5 6 1
的其中一个 2 × 3 2 \times 3 2×3的子矩阵是
4 7 4
8 6 9
-
相邻的元素:矩阵中的某个元素与其上下左右四个元素(如果存在的话)是相邻的。
-
矩阵的分值:矩阵中每一对相邻元素之差的绝对值之和。
本题任务:给定一个 n n n行 m m m列的正整数矩阵,请你从这个矩阵中选出一个 r r r行 c c c列的子矩阵,使得这个子矩阵的分值最小,并输出这个分值。
(本题目为2014NOIP普及T4)
输入格式
第一行包含用空格隔开的四个整数 n , m , r , c n,m,r,c n,m,r,c,意义如问题描述中所述,每两个整数之间用一个空格隔开。
接下来的 n n n行,每行包含 m m m个用空格隔开的整数,用来表示问题描述中那个 n n n行 m m m列的矩阵。
输出格式
一个整数,表示满足题目描述的子矩阵的最小分值。
样例 #1
样例输入 #1
5 5 2 3
9 3 3 3 9
9 4 8 7 4
1 7 4 6 6
6 8 5 6 9
7 4 5 6 1
样例输出 #1
6
样例 #2
样例输入 #2
7 7 3 3
7 7 7 6 2 10 5
5 8 8 2 1 6 2
2 9 5 5 6 1 7
7 9 3 6 1 7 8
1 9 1 4 7 8 8
10 5 9 1 1 8 10
1 3 1 5 4 8 6
样例输出 #2
16
提示
【输入输出样例1说明】
该矩阵中分值最小的 2 2 2行 3 3 3列的子矩阵由原矩阵的第 4 4 4行、第 5 5 5行与第 1 1 1列、第 3 3 3列、第 4 4 4列交叉位置的元素组成,为
6 5 6
7 5 6
,其分值为:
|6−5| + |5−6| + |7−5| + |5−6| + |6−7| + |5−5| + |6−6| =6 。
【输入输出样例2说明】
该矩阵中分值最小的3行3列的子矩阵由原矩阵的第 4 4 4行、第 5 5 5行、第 6 6 6行与第 2 2 2列、第 6 6 6列、第 7 7 7列交叉位置的元素组成,选取的分值最小的子矩阵为
9 7 8
9 8 8
5 8 10
【数据说明】
对于 50 % 50\% 50%的数据, 1 ≤ n ≤ 12 , 1 ≤ m ≤ 12 1 ≤ n ≤ 12,1 ≤ m ≤ 12 1≤n≤12,1≤m≤12,矩阵中的每个元素 1 ≤ a i j ≤ 20 1 ≤ a_{ij} ≤ 20 1≤aij≤20;
对于 100 % 100\% 100%的数据, 1 ≤ n ≤ 16 , 1 ≤ m ≤ 16 1 ≤ n ≤ 16,1 ≤ m ≤ 16 1≤n≤16,1≤m≤16,矩阵中的每个元素 1 ≤ a i j ≤ 1 , 000 , 1 ≤ r ≤ n , 1 ≤ c ≤ m 1 ≤ a_{ij} ≤ 1,000,1 ≤ r ≤ n,1 ≤ c ≤ m 1≤aij≤1,000,1≤r≤n,1≤c≤m。
思路
暴力-60pts
加上最优性剪枝-80pts
记忆化-100pts(当然正解就是dp)
但是学校里的老旧机器只跑了70分
预处理-100pts
时间变化:
代码
#include<cmath>
#include<cstdio>
#include<cstring>
typedef long long ll;
ll read(){
ll s=0,w=1;char c=getchar();
while(c<48||c>57){if(c=='-')w=-1;c=getchar();}
while(c>47&&c<58)s=(s<<3)+(s<<1)+(c^48),c=getchar();
return s*w;
}
void getmin(ll &a,ll b){a=((a>b)?b:a);return;}
ll n,m,l,r,ans=0x3f3f3f3f,a[20][20],hang[20],lie[20],aans[20][20][20],ljdzc[20],lh[20][20];
void ycl(){
memset(aans,0x3f,sizeof(aans));memset(ljdzc,0,sizeof(ljdzc));memset(lh,0,sizeof(lh));
for(int i=1;i<=m;i++)
for(int j=2;j<=l;j++)ljdzc[i]+=abs(a[hang[j]][i]-a[hang[j-1]][i]);
for(int i=1;i<=l;i++)
for(int j=1;j<=m;j++)
for(int k=1;k<j;k++)lh[k][j]+=abs(a[hang[i]][j]-a[hang[i]][k]);
}
void dfs2(ll wz,ll sl,ll nbans,ll pre){
if(sl>r||sl+(m-wz+1)<r||(pre!=-1&&nbans>=aans[wz][sl][pre]))return;
aans[wz][sl][pre]=nbans;
if(wz==m+1&&sl==r){getmin(ans,nbans);return;}
dfs2(wz+1,sl,nbans,pre);
lie[sl+1]=wz;
nbans+=ljdzc[wz];
if(sl!=0)nbans+=lh[pre][wz];
if(nbans>=ans)return;
dfs2(wz+1,sl+1,nbans,wz);
return;
}
void dfs1(ll wz,ll sl){
if(sl>l||sl+(n-wz+1)<l)return;
if(wz==n+1&&sl==l){ ycl();dfs2(1,0,0,-1);return;}
hang[sl+1]=wz;dfs1(wz+1,sl+1);dfs1(wz+1,sl);
return;
}
int main(){
n=read();m=read();l=read();r=read();
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)a[i][j]=read();
dfs1(1,0);printf("%lld",ans);return 0;
}