vijos 1283 佳佳的魔杖 - 前缀和优化dp

传送门

题解：

不知道他们怎么看出是LCS的。

首先考虑转化成图，选了[l,r]就不能选择[s,t]，l<=s<=t<=r，|[s,t]|!=|[l,r]|；

所以可以用网络流做最小割，也可以在这张无环图上做dp，但是边数是O(n^4)的。

发现如果左端点相同至多只会选择一根。

考虑dp，记dp[l][r]表示选了[l,r]选了左端点最靠左的区间为[l,r]（当然有可能不合法，此时记为0）

这样dp[l][r]是从dp[s][t]转移而来，其中s>l, t>r。发现这是一个矩阵。

用类似前缀和的方法维护一个f数组。f[i][j]=dp[s][t]，s>=i, t>=j。

f[i][j]=max{f[i+1][j],f[i][j+1],dp[i][j]}。

这样dp[i][j]=f[i+1][j+1]+calc(i,j);

然后发现干脆连dp数组也省去好了，这样就和网上的程序长的一模一样了。

好像是有一点LCS的味道：

//vijos 1283
#include<iostream>
#include<cstring>
#include<cstdio>
#include<climits>
#define INF INT_MAX/2
#define N 1010
#define lint long long
using namespace std;
int minr[N],maxr[N];
lint f[N][N];
int s[N],len[N],val[N];
int main()
{
	int n,minl,maxl;scanf("%d%d%d",&n,&minl,&maxl);
	for(int i=1;i<=n;i++) scanf("%d",&len[i]);
	for(int i=1;i<=n;i++) s[i]=s[i-1]+len[i];
	for(int i=1;i<=n;i++) scanf("%d",&val[i]);
	for(int i=1;i<=n;i++) val[i]+=val[i-1];
	int mr=0;
	for(int i=1;i<=n;i++)
	{
		while(mr<=n&&s[mr]-s[i-1]<minl) mr++;
		minr[i]=mr;
	}
	s[mr=n+1]=INF;
	for(int i=n;i>=1;i--)
	{
		while(mr>=1&&s[mr]-s[i-1]>maxl) mr--;
		maxr[i]=mr;
	}
//	for(int i=1;i<=n;i++) cout<<minr[i]<<" "<<maxr[i]<<endl;
	lint ans=0;
	for(int i=n;i>=1;i--)
		for(int j=n;j>=1;j--)
		{
			if(j>=minr[i]&&j<=maxr[i]) f[i][j]=f[i+1][j+1]+val[j]-val[i-1];
			f[i][j]=max(max(f[i][j+1],f[i+1][j]),f[i][j]);
			ans=max(ans,f[i][j]);
			//cout<<i<<" "<<j<<" "<<dp[i][j]<<endl;
		}
	printf("%lld\n",ans);return 0;
}