题解:
不知道他们怎么看出是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;
}