区间dp小结

——这篇文章主要想总结下区间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;

}