分析

#include<bits/stdc++.h>
#define N 55
#define inf 0x3fffffff
using namespace std;
int x[N],w[N];
int n,a,b;
int g[N][N][N][N],f[N][N];
int read(){
  int cnt=0; char ch=0;
  while(!isdigit(ch)) ch=getchar();
  while(isdigit(ch))cnt=cnt*10+(ch-'0'),ch=getchar();
  return cnt;
}
int main(){
  n=read(),a=read(),b=read();
  for(int i=1;i<=n;i++) x[i]=w[i]=read();
  sort(x+1,x+n+1);
  int siz=unique(x+1,x+n+1)-(x+1);
  for(int i=1;i<=n;i++) w[i]=lower_bound(x+1,x+siz+1,w[i])-x;
  for(int i=1;i<=n;i++) f[i][i]=a;
  for(int len=2;len<=n;len++){
    for(int l=1,r=len;r<=n;l++,r++){
      /*----------------------------初始化----------------------------*/
      for(int i=l;i<=r;i++)
        for(int j=siz;j;j--)
          for(int k=siz;k>=j;k--)
            g[l][i][j][k]=inf;
      g[l][l][w[l]][w[l]]=0;
      /*----------------------------更新G------------------------------*/
      for(int i=l;i<r;i++)
        for(int j=siz;j;j--)
          for(int k=siz;k>=j;k--){
            if(g[l][i][j][k]==inf) continue;
            int &tmp=g[l][i+1][min(j,w[i+1])][max(k,w[i+1])];
            tmp=min(tmp,g[l][i][j][k]);
            for(int t=i+1;t<=r;t++) 
              g[l][t][j][k]=min(g[l][t][j][k],g[l][i][j][k]+f[i+1][t]); 
          }
      /*---------------------------更新F------------------------------*/
      f[l][r]=inf;
      for(int i=siz;i;i--)
        for(int j=siz;j>=i;j--)
          f[l][r]=min(f[l][r],g[l][r][i][j]+a+b*(x[j]-x[i])*(x[j]-x[i])); 
    }
  }
  printf("%d\n",f[1][n]); return 0;
}