第一眼看就是跑不掉的DP,O(N^2)的算法给出状态转移方程:
f[i]=max(f[i],f[j]+(sum[i]-sum[j])-i*100);
不过N比较大,两重循环绝对不行,需要把内层循环优化到O(1),因为每次都是找max(f[j]-sum[j]),而且要花费N次,i,i+1找的区间只是错开了一个,有了这个特征,就可以可以用一种比较强大的方法
。今天就学一种动规的优化方法,单调队列优化。这道题用的是递减队列。
这个题是是参考了网上题解才AC成的,在此感谢原作者
dongdongzhang_。
代码如下:
/*2014.9.19 OJ1326《超级教主》 程序猿:卢裕东
东哥出品,必属精品!!!
版权所有,侵权必究!!!*/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
int N,M,now,sum[2000010],f[2000010];
int q[2000010],head,tail;
int GET(int i)
{
while(head<tail&&(f[q[head]]<i*100)) head++;//去掉无法到达的高度。
return q[head];
}
void PUSH(int i)//维护该单调递减的队列。
{
while(head<tail&&(f[i]-sum[i]>f[q[tail-1]]-sum[q[tail-1]])) tail--;
q[tail++]=i;
}
void init()
{
scanf("%d%d",&N,&M);
for(int i=1;i<=N;i++)
{
scanf("%d",&now);
sum[i]=sum[i-1]+now;
}
f[0]=M;
q[tail++]=0;//记录初始最高高度。
}
void DP()
{
for(int i=1;i<=N;i++)
{
int k=GET(i);
if(f[k]>=i*100) f[i]=max(f[i],f[k]+(sum[i]-sum[k])-i*100);//这里f[k]>=i*100,等于也可取!注意!!!
PUSH(i);
}
printf("%d\n",f[N]);
}
int main()
{
init();
DP();
return 0;
}
/*思路:DP+单调递减队列优化;
时间复杂度:O(N)。*/