由于每次合并可以刻画为向外延伸，那么考虑区间 \(\text{dp}\)，设 \(dp_{l,r,m,c}\) 表示考虑了 \([l,r]\) 且剩下了一个质量为 \(m \in [0,K)\) 颜色为 \(c\) 的史莱姆的答案。
状态过大且转移方程不便于优化而考虑优化状态，由于对于一个极短的需要合并成一个史莱姆的区间 \([p,q]\)，最终的颜色只能是 \(c_p\) 或 \(c_q\)，中间的颜色不可能延展出去否则会与极短不符。
所以设 \(f_{l,r}\) 表示处理掉 \([l,r]\) 的所有史莱姆的答案，\(fl_{l,r,p}\) 表示剩下一个质量为 \(p\) 且颜色为 \(c_l\) 的史莱姆的答案，\(fr_{l,r,p}\) 则是颜色为 \(c_r\) 的。
至于 \(p_i-p_{i-1}< w\) 则表示要合并的史莱姆并不会手动提升它的质量。
先考虑转移 \(fl_{l,r,p}\)，首先可以直接删除 \((l,r]\)，即 \(fl_{l,r,m_l}=f_{l+1,r}\)，否则可以先保留 \(l\)，再删除中间一段区间，再把剩下的 \(fl\) 和 \(l\) 这个史莱姆合并，即 \(fl_{l,r,p}=\min_{l< k \le r,c_l=c_k} f_{l+1,k-1}+fl_{k,r,p-m_l}\)，\(fr_{l,r,p}\) 同理可得。
至于 \(f_{l,r}\) 的转移分为是否是一个极短合并区间，如果需要合并首先需要满足 \(c_l=c_r\)，那么 \(f_{l,r}=\min_{l\le k< r,0< i< K,0< j< K}fl_{l,k,i}+fr_{k+1,r,j}+p_{i+j}\)，如果不合并而直接合并两边的 \(f\)，那么 \(f_{l,r}=\min_{l\le k< r,c_k\neq c_{r+1} \vee c_{k+1}\neq c_{l-1}}\)，\(c\) 的限制是因为不满足就会自动与 \([l,r]\) 外的史莱姆合并。
可能有写错的地方记得来打我，可能有人要问 \(fl,fr,f\) 是不是还应该有其他关于 \(c\) 的限制啊，应该是存在一种顺序满足条件吧，其实我也不大知道，加上就过不了。
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int n,s,w,c[151],m[151],p[21];
ll f[151][151],fl[151][151][21],fr[151][151][21];
void upd(ll &x,ll y){x=max(x,y);}
int main(){
memset(f,-0x3f,sizeof(f)),memset(fl,-0x3f,sizeof(fl)),memset(fr,-0x3f,sizeof(fr));
scanf("%d%d%d",&n,&s,&w);
for(int i=1;i<=n;i++) scanf("%d",&c[i]);
for(int i=1;i<=n;i++) scanf("%d",&m[i]);
for(int i=s;i<=2*s-2;i++) scanf("%d",&p[i]);
for(int i=1;i<s;i++) p[i]=p[s]-(s-i)*w;
for(int i=1;i<=n;i++) f[i][i]=p[m[i]],fl[i][i][m[i]]=fr[i][i][m[i]]=f[i][i-1]=0;
for(int t=2;t<=n;t++){
for(int l=1;l+t-1<=n;l++){
int r=l+t-1;
upd(fl[l][r][m[l]],f[l+1][r]);
for(int k=l+1;k<=r;k++) if(c[l]==c[k]) for(int i=1;i<s-m[l];i++) upd(fl[l][r][i+m[l]],f[l+1][k-1]+fl[k][r][i]);
upd(fr[l][r][m[r]],f[l][r-1]);
for(int k=l;k<r;k++) if(c[r]==c[k]) for(int i=1;i<s-m[r];i++) upd(fr[l][r][i+m[r]],f[k+1][r-1]+fr[l][k][i]);
for(int k=l;k<r;k++) if(c[k]!=c[r+1]||c[k+1]!=c[l-1]) upd(f[l][r],f[l][k]+f[k+1][r]);
if(c[l]==c[r]) for(int k=l;k<r;k++) for(int i=1;i<s;i++) for(int j=1;j<s;j++) upd(f[l][r],fl[l][k][i]+fr[k+1][r][j]+p[i+j]);
}
}
printf("%lld",f[1][n]);
return 0;
}