Codeforces.392E.Deleting Substrings(区间DP)

题目链接

(Description)

(Solution)

合法的子序列只有三种情况：递增，递减，前半部分递增然后一直递减（下去了就不会再上去了）（当然还要都满足(|a_{i+1}-a_i|=1)）。
容易想到区间DP。(f[i][j])表示把区间([i,j])全部删除的最大收益，还需要(g[i][j])表示将区间([i,j])删成连续上升的一段（(a_isim a_j)）的最大收益，(h[i][j])表示将区间([i,j])删成连续下降的一段（(a_isim a_j)）的最大收益。
那么(g[i][j])的元素个数就是(a_j-a_i+1)，(h[i][j])的元素个数为(a_i-a_j+1)，合并(g[i][k],h[k][j])后的元素个数就是(2a_k-a_i-a_j+2-1)（减掉1个(a_k)）。
那么$$f[i][j]=max{f[i][k]+f[k+1][j], g[i][k]+h[k][j]+v[2a_k-a_i-a_j+1]}$$

其实(g,h)用一个数组就可以了，因为只需要判断一下(a_i,a_j)的大小关系，就知道是上升序列还是下降序列了。
最后的答案就是(f)的最大子段和。(n^2)求出来即可。

//78ms	1100KB
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
typedef long long LL;
const int N=405,INF=0x3f3f3f3f;

int A[N],val[N],f[N][N],g[N][N],dp[N];

inline int read()
{
	int now=0,f=1;register char c=gc();
	for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
	for(;isdigit(c);now=now*10+c-'0',c=gc());
	return now*f;
}

int main()
{
	int n=read();
	for(int i=1; i<=n; ++i) val[i]=read();
	for(int i=1; i<=n; ++i) A[i]=read();

	for(int i=1; i<=n; ++i)
	{
		f[i][i]=val[1], g[i][i]=0;
		for(int j=i+1; j<=n; ++j) f[i][j]=g[i][j]=-INF;
	}

	for(int len=1; len<n; ++len)
		for(int i=1; i+len<=n; ++i)
		{
			int j=i+len;
			for(int k=i; k<j; ++k)
			{
				f[i][j]=std::max(f[i][j],f[i][k]+f[k+1][j]);
				if((A[i]<A[j] && A[j]==A[k]+1)||(A[i]>A[j] && A[j]==A[k]-1))
					g[i][j]=std::max(g[i][j],g[i][k]+f[k+1][j-1]);
			}
			for(int k=i; k<=j; ++k)
				if(A[k]>=A[i] && A[k]>=A[j] && 2*A[k]-A[i]-A[j]+1<=n)//这东西显然不会<=0啊 
					f[i][j]=std::max(f[i][j],g[i][k]+g[k][j]+val[2*A[k]-A[i]-A[j]+1]);
		}
	for(int i=1; i<=n; ++i)
	{
		dp[i]=dp[i-1];
		for(int j=1; j<=i; ++j) dp[i]=std::max(dp[i],dp[j-1]+f[j][i]);
	}
	printf("%d
",dp[n]);

	return 0;
}