——这篇文章主要想总结下区间DP,同时给自己复习巩固上一次课刚讲的这方面知识点。
区间DP
一、定义
区间DP,顾名思义是在区间上DP,它的主要思想就是先在小区间进行DP得到最优解,然后再利用小区间的最优解合并求大区间的最优解。
二、实现思路
下面给出区间DP最简单形式的伪代码(具体要根据题目修改)
//mst(dp,0) 初始化DP数组
for(int i=1;i<=n;i++)
{
dp[i][i]=初始值
}
for(int len=2;len<=n;len++) //区间长度
for(int i=1;i<=n;i++) //枚举起点
{
int j=i+len-1; //区间终点
if(j>n) break; //越界结束
for(int k=i;k<j;k++) //枚举分割点,构造状态转移方程
{
dp[i][j]=max(dp[i][j],dp[i][k]+dp[k+1][j]+w[i][j]);
}
}
区间DP,区间上做动态规划。一个大区间必然由长度不等的小区间合并而来,在合并过程中,最基本也是动态规划的必须原则:满足最优化原理和无后效性原则,所以在确定状态转移方程时,得特别留心这两点。
区间DP区别于其它线性DP题的特点是:区间DP通过枚举区间分界点实现状态转移。
词穷了,还是写题来感受吧~~,有坑自己再补上。
NOI1995 石子合并
在一个圆形操场的四周摆放N堆石子,现要将石子有次序地合并成一堆.规定每次只能选相邻的2堆合并成新的一堆,并将新的一堆的石子数,记为该次合并的得分。
试设计出1个算法,计算出将N堆石子合并成1堆的最小得分。
区间dp常用的一个状态就是dp[i][j]表示i~j这个区间的最优值是多少?
emmmm 我们可以看出题目中这个合并过程,很像是区间dp的一种合并,也就是说dp[i][j]可以由dp[i][k]和dp[k+1][j]转移过来
那么我们自然而然的想到要枚举上一次合并是在哪个位置,也就是断点k,然后进行转移
不过,由于这个题目的特殊限制,我们需要记录g[i][j]表示i到j这个区间有多少石子,来进行辅助转移
好啦!至此,我们的dp方程也就解决了
f[i][j]=min(f[i][j],f[i][k]+f[k+1][j]+g[i][k]+g[k+1][j]) ( i<=k<=j)
对于环,我们只需要把整个数组复制一遍,然后在统计答案的时候,把1-n开头的长度为n的区间的答案求一个min 就可以啦
上代码:
#include<bits/stdc++.h>
using namespace std;
const int maxn = 1010;
int g[maxn][maxn],f[maxn][maxn],i,j,k,m,a[maxn],n,g1[maxn][maxn],f1[maxn][maxn];
int ans = 0xfffffff;
int main ()
{
scanf("%d",&n);
for (i=1;i<=2*n;i++)
for (j=1;j<=2*n;j++)
f[i][j]=0xfffffff;
for (i=1;i<=n;i++)
scanf("%d",&a[i]);
for (i=n+1;i<=2*n;i++)
a[i]=a[i-n];
for (i=1;i<=2*n;i++)
{
f[i][i]=0;
g[i][i]=a[i];
}
for (i=1;i<=n;i++)
for (j=1;j<=2*n-i;j++)
{
for (k=j;k<j+i;k++)
{
int cnt=g[j][k]+g[k+1][j+i];
if (f[j][j+i]>cnt+f[j][k]+f[k+1][j+i])
{
f[j][j+i]=cnt+f[j][k]+f[k+1][j+i];
g[j][i+j]=cnt;
}
}
}
for (i=1;i<=n+1;i++)
ans=min(ans,f[i][i+n-1]);
cout<<ans<<endl;
return 0;
}