51Nod 1022 石子归并 V2（区间DP+四边形优化）

题目链接：http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1022

题目大意：

N堆石子摆成一个环。现要将石子有次序地合并成一堆。规定每次只能选相邻的2堆石子合并成新的一堆，并将新的一堆石子数记为该次合并的代价。计算将N堆石子合并成一堆的最小代价。

例如： 1 2 3 4，有不少合并方法

1 2 3 4 => 3 3 4(3) => 6 4(9) => 10(19)

1 2 3 4 => 1 5 4(5) => 1 9(14) => 10(24)

1 2 3 4 => 1 2 7(7) => 3 7(10) => 10(20)

括号里面为总代价可以看出，第一种方法的代价最低，现在给出n堆石子的数量，计算最小合并代价。

解题思路：经典的石子合并问题，较原来不同的是石子是环形摆放的，而且石子数目n的范围有100增加到了1000，原本O（n^3）的算法肯定是会超时的。所以需要用四边形不等式优化将复杂度降为O（n^2），并且将数组倍增把环变为链。

粗略介绍一下四边形优化的作用，具体证明看这里《动态规划加速原理之四边形不等式》。

我们原本的状态转移方程为dp[i][j]=min{dp[i][k]+dp[k+1][j]+w[i][j]}（i<j，i<=k<=j）。

上式在动态规划的状态转移方程中是很常见的，对于上式中的w(i,j)
如果符合w(i`,j) <= w(i,j`) i<i`<j<j`        那么我们称函数w满足关于区间包含的单调性。
如果符合w(i,j)+w(i`,j`) <= w(i`,j)+w(i,j`)    那么我们称函数w满足四边形不等式。

那么就可以使用两个定理（图片来源）

于是，我们可以使用s[i][j]记录使得dp[i][j]最优的分割点（k点），并且满足s[i][j-1]<=s[i][j]<=s[i+1][j+1]，那么我们的k的枚举范围就是s[i][j-1]<=s[i][j]<=s[i+1][j+1]。

复杂度证明:

代码：

 #include<iostream>
 #include<cstring>
 #include<cstdio>
 #include<algorithm>
 using namespace std;
 const int N=2e3+;
 const int INF=0x3f3f3f3f;

 int dp[N][N],s[N][N],sum[N],a[N];//s[i][j]为使dp[i][j]最优的分割点 

 int main(){
     int n;
     scanf("%d",&n);
     memset(dp,0x3f,sizeof(dp));
     for(int i=;i<=n;i++){
         scanf("%d",&a[i]);
         a[i+n]=a[i];
     }
     for(int i=;i<=*n;i++){
         dp[i][i]=;
         sum[i]=sum[i-]+a[i];
         s[i][i]=i;
     }
     for(int d=;d<n;d++){
         for(int i=;i<=*n-d;i++){
             int j=i+d;
             for(int k=s[i][j-];k<=s[i+][j];k++){
                 int tmp=dp[i][k]+dp[k+][j]+sum[j]-sum[i-];
                 if(tmp<dp[i][j]){
                     dp[i][j]=tmp;
                     s[i][j]=k;
                 }
             }
         }
     }
     int ans=INF,d=n-;
     for(int i=;i<=*n-d;i++){
         int j=i+d;
         ans=min(ans,dp[i][j]);
     }
     printf("%d\n",ans);
     return ;
 }