BZOJ4897 [Thu Summer Camp2016]成绩单

挺厉害的题

考虑区间DP，f[i][j]表示把[i,j]都拿光的最小代价

f[i][j]怎么算呢？我们考虑最后一次发成绩单

最后一次发成绩单发的一定是原序列的一个子序列，那么我们用g[i][j][k][l]表示以一个起始点>=i，最后一个元素是j，最小值为k，最大值为l的子序列为最后一次操作的代价

那么g[i][j][k][l]+f[j+1][u-1]就可以转移到g[i][u][min(k,a[u)][max(l,a[u])]

那么计算f[i][j]f的时候，枚举一下子序列的最后一个元素的位置，再加上最后那段的代价，即可算出f[i][j]

对f和g一起转移即可

复杂度n五次方

这数据范围这复杂度，一脸TC气息

#include<iostream>
#include<cstring>
#include<ctime>
#include<cmath>
#include<algorithm>
#include<iomanip>
#include<cstdlib>
#include<cstdio>
#include<map>
#include<bitset>
#include<set>
#include<stack>
#include<vector>
#include<queue>
using namespace std;
#define MAXN 52
#define MAXM 1010
#define ll long long
#define eps 1e-8
#define MOD 1000000007
#define INF 1000000000
int f[MAXN][MAXN][MAXN][MAXN];
bool visf[MAXN][MAXN][MAXN][MAXN];
int g[MAXN][MAXN];
bool visg[MAXN][MAXN];
int n;
int tls[MAXN],mx;
map<int,int>h;
int gg[MAXN];
int a[MAXN],b[MAXN];
int A,B;
int main(){
	int L,l,r,i,j,k;
	scanf("%d%d%d",&n,&A,&B);
	for(i=1;i<=n;i++){
		scanf("%d",&a[i]);
		tls[i]=a[i];
	}
	sort(tls+1,tls+n+1);
	for(i=1;i<=n;i++){
		if(tls[i]!=tls[i-1]||i==1){
			gg[h[tls[i]]=++mx]=tls[i];
		}
	}
	for(i=1;i<=n;i++){
		b[i]=h[a[i]];
	}
	memset(f,0x3f,sizeof(f));
	memset(g,0x3f,sizeof(g));
	for(i=1;i<=n;i++){
		f[i][i][b[i]][b[i]]=0;
		g[i][i]=A;
		g[i][i-1]=0;
	}
	g[n+1][n]=0;
	for(L=2;L<=n;L++){
		for(l=1;(r=l+L-1)<=n;l++){
			f[l][r][b[r]][b[r]]=g[l][r-1];
			for(i=l;i<r;i++){
				for(j=1;j<=mx;j++){
					for(k=j;k<=mx;k++){
						int tj=min(j,b[r]);
						int tk=max(k,b[r]);
						f[l][r][tj][tk]=min(f[l][r][tj][tk],f[l][i][j][k]+g[i+1][r-1]);
					}
				}
			}
			for(i=l;i<=r;i++){
				for(j=1;j<=mx;j++){
					for(k=j;k<=mx;k++){
						g[l][r]=min(g[l][r],f[l][i][j][k]+g[i+1][r]+A+B*(gg[k]-gg[j])*(gg[k]-gg[j]));
					}
				}
			}
		}
	}
	printf("%d\n",g[1][n]);
	return 0;
}

/*
10
3 1
7 10 9 10 6 7 10 7 1 2

*/