题面 [USACO11OPEN] Mowing the Lawn G
题目描述
在一年前赢得了小镇的最佳草坪比赛后,Farmer John 变得很懒,再也没有修剪过草坪。现在,新一轮的最佳草坪比赛又开始了,Farmer John 希望能够再次夺冠。
然而,Farmer John 的草坪非常脏乱,因此,Farmer John 只能够让他的奶牛来完成这项工作。Farmer John 有 NNN(1≤N≤1051\le N\le 10^51≤N≤105)只排成一排的奶牛,编号为 1…N1\ldots N1…N。每只奶牛的效率是不同的,奶牛 iii 的效率为 EiE_iEi(0≤Ei≤1090\le E_i\le 10^90≤Ei≤109)。
靠近的奶牛们很熟悉,因此,如果 Farmer John安排超过 KKK 只连续的奶牛,那么,这些奶牛就会罢工去开派对 😃。因此,现在 Farmer John 需要你的帮助,计算 FJ 可以得到的最大效率,并且该方案中没有连续的超过 KKK 只奶牛。
输入格式
第一行:空格隔开的两个整数 NNN 和 KKK。
第二到 N+1N+1N+1 行:第 i+1i+1i+1 行有一个整数 EiE_iEi。
输出格式
第一行:一个值,表示 Farmer John 可以得到的最大的效率值。
样例 #1
样例输入 #1
5 2
1
2
3
4
5
样例输出 #1
12
动态规划
发现这题的限制条件类似于“没有上司的舞会”,即前面的选不选对后面造成“能不能选”的影响。
于是尝试设有与没有两种状态:
dpi,1dp_{i,1}dpi,1表示仅考虑前iii头奶牛时选择第iii头奶牛的最大值
dpi,0dp_{i,0}dpi,0表示仅考虑前iii头奶牛时不选第iii头奶牛的最大值。
接下来计算递推式:
dpi,0=maxj∈[1,i){dpj,0,dpj,1}dp_{i,0}=\mathop{max}\limits_{j\in[1,i)}\{dp_{j,0},dp_{j,1}\}dpi,0=j∈[1,i)max{dpj,0,dpj,1}
dpi,1=maxj∈[i−k,i){dpj,0+∑t=j+1iEt}dp_{i,1}=\mathop{max}\limits_{j\in[i-k,i)}\{dp_{j,0}+\sum\limits_{t=j+1}^{i}E_t\}dpi,1=j∈[i−k,i)max{dpj,0+t=j+1∑iEt}
发现这样递推有O(n3)O(n^3)O(n3)的复杂度(一层iii循环,一层jjj循环,一层ttt循环),同时n≤105n\le10^5n≤105,所以要优化。
对于dpi,0dp_{i,0}dpi,0,容易发现是单调不减的,所以可以直接写作dpi,0=max{dpi−1,0,dpi−1,1}dp_{i,0}=max\{dp_{i-1,0},dp_{i-1,1}\}dpi,0=max{dpi−1,0,dpi−1,1}
接下来看dpi,1dp_{i,1}dpi,1的优化:
首先发现∑t=j+1iEt\sum\limits_{t=j+1}^{i}E_tt=j+1∑iEt部分可以用前缀和sumi−sumjsum_i-sum_jsumi−sumj替换,时间复杂度变为O(n2)O(n^2)O(n2)
其次对式子进行处理:
dpi,1=maxj∈[i−k,i){dpj,0+sumi−sumj} =maxj∈[i−k,i){dpj,0−sumj}+sumidp_{i,1} =\mathop{max}\limits_{j\in[i-k,i)}\{dp_{j,0}+sum_i-sum_j\} \\ \ =\mathop{max}\limits_{j\in[i-k,i)}\{dp_{j,0}-sum_j\}+sum_idpi,1=j∈[i−k,i)max{dpj,0+sumi−sumj} =j∈[i−k,i)max{dpj,0−sumj}+sumi
于是变为维护dpj,0−sumjdp_{j,0}-sum_jdpj,0−sumj的最大值,一眼滑动窗口问题,单调队列优化即可。
另外,第二维是可以被优化掉的,即dpi=max(dpi,1,dpi,0)dp_i=max(dp_{i,1},dp_{i,0})dpi=max(dpi,1,dpi,0)
代码如下:
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
const int N = 1e5+5;
long long sum[N],n,m,f[N],q[N];
int t=0,h=0;
int main()
{
cin>>n>>m;
long long maxx = 0;
for(int i=1;i<=n;i++)
{
cin>>sum[i];
sum[i]+=sum[i-1];
}
for(int i=1;i<=n;i++)
{
while(t>=h && f[i-1]-sum[i]>=f[q[t]-1]-sum[q[t]])
t--;
t++;
q[t] = i;
while(t>=h && q[h]<i-m)
h++;
int j = q[h];
f[i] = sum[i]+f[j-1]-sum[j];
maxx = max(maxx,f[i]);
}
cout<<maxx<<endl;
return 0;
}