Bootstrap

P1725 琪露诺

https://www.luogu.com.cn/problem/P1725

题目描述

在幻想乡,琪露诺是以笨蛋闻名的冰之妖精。

某一天,琪露诺又在玩速冻青蛙,就是用冰把青蛙瞬间冻起来。但是这只青蛙比以往的要聪明许多,在琪露诺来之前就已经跑到了河的对岸。于是琪露诺决定到河岸去追青蛙。

小河可以看作一列格子依次编号为0到N,琪露诺只能从编号小的格子移动到编号大的格子。而且琪露诺按照一种特殊的方式进行移动,当她在格子i时,她只移动到区间[i+l,i+r]中的任意一格。你问为什么她这么移动,这还不简单,因为她是笨蛋啊。

每一个格子都有一个冰冻指数A[i],编号为0的格子冰冻指数为0。当琪露诺停留在那一格时就可以得到那一格的冰冻指数A[i]。琪露诺希望能够在到达对岸时,获取最大的冰冻指数,这样她才能狠狠地教训那只青蛙。

但是由于她实在是太笨了,所以她决定拜托你帮它决定怎样前进。

开始时,琪露诺在编号0的格子上,只要她下一步的位置编号大于N就算到达对岸。

输入格式

第1行:3个正整数N, L, R

第2行:N+1个整数,第i个数表示编号为i-1的格子的冰冻指数A[i-1]

输出格式

一个整数,表示最大冰冻指数。保证不超过2^31-1

输入输出样例

输入 #1复制

5 2 3
0 12 3 11 7 -2

输出 #1复制

11

说明/提示

对于60%的数据:N <= 10,000

对于100%的数据:N <= 200,000

对于所有数据 -1,000 <= A[i] <= 1,000且1 <= L <= R <= N


先来一波dp的朴素版

记dp[i]为到第i个点的最大和,那么转移方程为dp[i]=max(dp[j]) ( j的范围)+a[i]

那么这个j的范围是多少呢?我们看一个点在什么条件下能到j点

假设这个点是j,那么满足j+l<=i&&j+r>=i这个条件就能到i点,转化一下是

j<=i-l,j>=i-r;那么此时j的范围就知道了

这样就可以O(n^2)暴力dp获得60分

注:这两天容易意外RE,问题在访问越界,比如这个i-r.就完全可能<0,导致越界,所以注意max(0,i-r);

代码如下:

#include<iostream>
#include<vector>
#include<queue>
#include<cstring>
#include<cmath>
#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=200010;
typedef long long LL;
LL dp[maxn];//dp 
LL q[maxn];//下标 
LL a[maxn];//值 
LL n,l,r;

int main(void)
{
  cin.tie(0);std::ios::sync_with_stdio(false);
  cin>>n>>l>>r;
  for(LL i=0;i<=n;i++) cin>>a[i];
  memset(dp,-0x3f,sizeof(dp));
  dp[0]=0;LL ans=-0x3f3f3f3f;
  for(LL i=l;i<=n+r;i++)
  {
  		LL k=-0x3f3f3f3f;
  		for(LL j=max((LL)0,i-r);j<=i-l;j++)
		  {
		  		k=max(k,dp[j]);
		  }	
		dp[i]=k+a[i];
		if(i>n) ans=max(ans,dp[i]);
  }	
  cout<<ans<<endl;
return 0;
}


那么怎么优化?这个 max(dp[j]) ( j的范围) 是一个滑动窗口求最大值的优化),单调队列维护dp[i],这样最后就是O(n)的复杂度

讲一下注意点:dp的初始化,求最大值的时候先初始化无穷小,然后把起点的dp[0]=0;

这个单调队列维护的是dp的最值,所以最开始能到的点是i=l的点

然后维护的第一个元素是dp[l-l],然后不断O(n)维护的是dp[i-l](入队)

 

#include<iostream>
#include<vector>
#include<queue>
#include<cstring>
#include<cmath>
#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=200010;
typedef long long LL;
LL dp[maxn];//dp 
LL q[maxn];//下标 
LL a[maxn];//值 
LL n,l,r;
void solve()
{
	memset(dp,-0x3f,sizeof(dp));
	LL h=1;LL t=0;
	dp[0]=0;
	LL ans=-0x3f3f3f3f;
	for(LL i=l;i<=n;i++)
	{
	   while(h<=t&&q[h]+r<i) h++;
	   
	   while(h<=t&&dp[q[t]]<=dp[i-l]) t--;
	   q[++t]=i-l;
	   dp[i]=dp[q[h]]+a[i];	
	   
	}	
	for(LL i=n-r+1;i<=n;i++)
	ans=max(ans,dp[i]);
	cout<<ans<<endl; 	
}
int main(void)
{
  cin.tie(0);std::ios::sync_with_stdio(false);
  cin>>n>>l>>r;
  for(LL i=0;i<=n;i++) cin>>a[i];
  solve();	
return 0;
}
;