题目大意
给一个长度为 n ( n < = 500 ) n(n<=500) n(n<=500)的数组,每次可以把相邻的两个相同元素 a i a_i ai和 a i + 1 a_{i_+1} ai+1合并为 a i + 1 a_i+1 ai+1,问数组在若干次操作后可以取得的最小长度为多少?
分析思路
关于这道题目网上已经有了很多题解博客,但是基本上都没有把分析思路写的很详细,于是我写下这篇博客,记录一下自己的分析过程。
这个题目的数据范围为
(
n
<
=
500
)
(n<=500)
(n<=500),所以可以考虑用区间DP搞一下。设
d
p
[
i
]
[
j
]
dp[i][j]
dp[i][j]为对应区间
[
i
,
j
]
[i,j]
[i,j]能够经过合并取得的最短长度。首先考虑对任一区间
[
i
,
j
]
[i,j]
[i,j]来说,要么最终可以合并为1个元素,要么不能合并为1个元素。如果不能合并为1个元素,那么必然存在一个最优分割点
k
k
k使得
d
p
[
i
]
[
j
]
=
d
p
[
i
]
[
k
]
+
d
p
[
k
+
1
]
[
j
]
dp[i][j]=dp[i][k]+dp[k+1][j]
dp[i][j]=dp[i][k]+dp[k+1][j]如果能合并为同一个元素,那么必然存在最优分割点
k
k
k使得
d
p
[
i
]
[
j
]
=
d
p
[
i
]
[
k
]
+
d
p
[
k
+
1
]
[
j
]
+
1
dp[i][j]=dp[i][k]+dp[k+1][j]+1
dp[i][j]=dp[i][k]+dp[k+1][j]+1而取到这个最优分割点的条件便是此分割点两侧区间均能合并为同一个值,且它们相等。因此,区间
[
i
,
j
]
[i,j]
[i,j]的状态转移方程为
d
p
[
i
]
[
j
]
=
min
k
=
i
j
−
1
{
d
p
[
i
]
[
k
]
+
d
p
[
k
+
1
]
[
j
]
−
(
d
p
[
i
]
[
k
]
=
=
1
&
&
d
p
[
k
+
1
]
[
j
]
=
=
1
&
&
m
[
i
]
[
k
]
=
=
m
[
k
+
1
]
[
j
]
)
}
dp[i][j]=\min\limits_{k=i}^{j-1}\{dp[i][k]+dp[k+1][j]-(dp[i][k]==1\&\&dp[k+1][j]==1\&\&m[i][k]==m[k+1][j])\}
dp[i][j]=k=iminj−1{dp[i][k]+dp[k+1][j]−(dp[i][k]==1&&dp[k+1][j]==1&&m[i][k]==m[k+1][j])}其中,
m
[
i
]
[
j
]
m[i][j]
m[i][j]表示区间
[
i
,
j
]
[i,j]
[i,j]最终合并为同一个元素的元素值。
AC代码
#include<bits/stdc++.h>
using namespace std;
const int maxn = 507;
const int MAXINT = 6000000;
int n, dp[maxn][maxn], m[maxn][maxn], a[maxn];
void solve(){
int i, j, k, len;
//初始化
for(i=1;i<=n;++i){
for(j=1;j<=n;++j){
dp[i][j] = MAXINT;
}
dp[i][i] = 1;
m[i][i] = a[i];
}
for(len=2;len<=n;++len){
for(i=1;i<=n-len+1;++i){
j = i + len - 1;
for(k=i;k<=j-1;++k){
dp[i][j] = min(dp[i][j], dp[i][k] + dp[k+1][j]);
if(dp[i][k] == 1 && dp[k+1][j] == 1 && m[i][k] == m[k+1][j]){ //如果可以合并为一个元素
dp[i][j] = 1;
m[i][j] = m[i][k] + 1;
}
}
}
}
cout<<dp[1][n];
}
int main(){
int t, i, j;
ios::sync_with_stdio(false);
cin>>n;
for(i=1;i<=n;++i) cin>>a[i];
solve();
return 0;
}