1.题意:
如果我们可以从左到右从一个序列中截取长度为k,k-1,k-2,…,1的子段(是截取所以子段两两不能重叠),且这些子段的元素和单调递增,则称其“k递增”。现在给你个序列,求k最大值。
2.思路:
注意到(1+k)*k/2<=n,故k范围很小。
考虑dp做。为了方便,本题中,从左到右元素和单增,等价于从右到左元素和单减,所以我们维护 dp[i][j]:长度为i的后缀截取长度为1到j的子段,得到的最长子段元素和的最大值(其实这样的想法基于已下事实:当确定在下标i截取长度为j+1的子段的前提下,最优解法就是在剩下的部分中截取j个子段并使得最长子段元素和最大)
状态转移方程如下:dp[i][j]=max(dp[i+1][j],dp[i+j][j-1](if sum[i,i+j-1]<dp[i+j][j-1] or j==1)); 而sum可以用前缀和得到。
然后就是各种细节了…主要就是非法状态、边界、初始化等问题,详情见注释
3.代码
#include <bits/stdc++.h>
#define pb push_back
#define fi first
#define se second
#define makp make_pair
using namespace std;
typedef pair<int, int> PII;
typedef pair<int,pair<int,int>>PIII;
typedef long long LL;
const int N=1e5+10,M=450,Mod=1e9+7;
template <typename T> bool chkMax(T &x, T y) { return (y > x) ? x = y, 1 : 0; }
template <typename T> bool chkMin(T &x, T y) { return (y < x) ? x = y, 1 : 0; }
template <typename T> void inline read(T &x) {
int f = 1; x = 0; char s = getchar();
while (s < '0' || s > '9') { if (s == '-') f = -1; s = getchar(); }
while (s <= '9' && s >= '0') x = x * 10 + (s ^ 48), s = getchar();
x *= f;
}
int n,m;
LL sum[N],dp[N][M];
void fr(char &x)
{
x=0;while(x!='M'&&x!='A') x=getchar();
}
void acc()
{
ios_base::sync_with_stdio(false);
cin.tie(NULL);
}
int main()
{
acc();
int t;
cin>>t;
while(t--)
{
cin>>n;
m=sqrt(2*n);
for(int i=1;i<=n;i++)
{
cin>>sum[i];
sum[i]+=sum[i-1];
}
for(int i=1;i<=n+1;i++)//由于状态更新的特殊性,可能会更新到n+1
for(int j=1;j<=m;j++)
dp[i][j]=0;//0代表无意义
for(int i=n;i>=1;i--)
for(int j=1;j<=m;j++)
{
if(n-i+1<j*(j+1)/2)break;//先确保后缀长度足够才能进行状态转移
if(j==1 || sum[i+j-1]-sum[i-1]<dp[i+j][j-1])dp[i][j]=max(dp[i][j],sum[i+j-1]-sum[i-1]);
dp[i][j]=max(dp[i][j],dp[i+1][j]);
}
int ans=1;
for(int i=1;i<=m;i++)
if(!dp[1][i])break;
else ans=i;
cout<<ans<<"\n";
}
}
4.收获
最优子结构的是dp的影子