拿走一個區間的代價只與最大最小值有關，並且如果最後一次拿走包含區間右端點的子序列一定不會使答案更劣，於是設f[i][j][x][y]為使i~j區間剩余最小值為x最大值為y且若有數剩余一定包含j的最小代價。轉移時考慮j最後和哪一段一起拿走，有f[i][j][min(x,a[j])][max(y,a[j])]=min{f[i][d-1][x][y]+f[d][j-1][0][0]}，這樣就能處理拿走一段後區間的合並了。

　　區間dp好難啊。

#include<iostream> 
#include<cstdio>
#include<cmath>
#include
 
<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
    int x=0,f=1;char c=getchar();
    while (c<‘0‘||c>‘9‘) {if (c==‘-‘) f=-1;c=getchar();}
    while (c>=‘0‘&&c<=‘9‘) x=(x<<1)+(x<<3)+(c^48),c=getchar();
    return x*f;
}
 
#define N 55
int n,w,v,a[N],b[N],f[N][N][N][N];
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj4897.in","r",stdin);
    freopen("bzoj4897.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read(),w=read(),v=read();
    for (int i=1;i<=n;i++) b[i]=a[i]=read();
    sort(b
 
+1,b+n+1);int t=unique(b+1,b+n+1)-b-1;
    for (int i=1;i<=n;i++) a[i]=lower_bound(b+1,b+t+1,a[i])-b;
    memset(f,42,sizeof(f));
    for (int i=1;i<=n+1;i++)
    {
        f[i][i-1][0][0]=0;
        for (int x=1;x<=t;x++)
            for (int y=x;y<=t;y++)
            f[i][i-1][x][y]=0;
    }
    for (int k=1;k<=n;k++)
        for (int i=1;i<=n-k+1;i++)
        {
            int j=i+k-1;
            for (int x=1;x<=t;x++)
                for (int y=x;y<=t;y++)
                    for (int d=i;d<=j;d++)
                    f[i][j][min(x,a[j])][max(y,a[j])]=min(f[i][j][min(x,a[j])][max(y,a[j])],f[i][d-1][x][y]+f[d][j-1][0][0]);
            for (int x=1;x<=t;x++)
                for (int y=x;y<=t;y++)
                f[i][j][0][0]=min(f[i][j][0][0],f[i][j][x][y]+w+v*(b[y]-b[x])*(b[y]-b[x]));
        }
    cout<<f[1][n][0][0];
    return 0;
}

BZOJ4897 THUSC2016成績單（區間dp）