听说n<=70?
我一直以为n<=400000......
区间DP
#include<cstdio> #include<algorithm> using namespace std; int n,K,F[105][105][105],b[105],Sum[105]; struct node{ int sum,id,val; }a[105]; bool cmp(node a,node b){ return a.sum<b.sum; } int main(){ scanf("%d%d",&n,&K); for (int i=1; i<=n; i++) scanf("%d",&a[i].sum); for (int i=1; i<=n; i++) scanf("%d",&a[i].id); for (int i=1; i<=n; i++) b[i]=a[i].id; for (int i=1; i<=n; i++) scanf("%d",&a[i].val); sort(a+1,a+n+1,cmp); sort(b+1,b+n+1); int N=unique(b+1,b+n+1)-b-1; for (int i=1; i<=n; i++) a[i].id=lower_bound(b+1,b+N,a[i].id)-b; for (int i=1; i<=n; i++) for (int j=i; j<=n; j++) for (int k=0; k<=N; k++) F[i][j][k]=1e9; for (int i=1; i<=n; i++) Sum[i]=Sum[i-1]+a[i].val; for (int i=1; i<=n; i++) for (int k=0; k<=N; k++){ F[i][i][k]=a[i].val; if (k>a[i].id) F[i][i][k]+=K; } for (int len=1; len<n; len++) for (int i=1; i<n; i++){ int j=i+len; if (j>n) break; for (int lim=0; lim<=N; lim++) for (int k=i; k<=j; k++){ F[i][j][lim]=min(F[i][j][lim],F[i][k-1][lim]+F[k+1][j][lim]+K+Sum[j]-Sum[i-1]); if (lim<=a[k].id) F[i][j][lim]=min(F[i][j][lim],F[i][k-1][a[k].id]+F[k+1][j][a[k].id]+Sum[j]-Sum[i-1]); } } int ans=1e9; for (int i=0; i<=N; i++) ans=min(ans,F[1][n][i]); printf("%d ",ans); return 0; }