poj 3017 Cut the Sequence（单调队列优化DP）
Cut the Sequence
$ solution: $
这道题出的真的很好，奈何数据水啊！
这道题当时看得一脸懵逼，说二分也不像二分，说贪心也不像贪心，说搜索吧这题数据范围怎么这么大？而且这题看起来也实在不好DP，当时是真的满头雾水。只能说是各个都尝试一下。最后还是选了DP来做第一步突破，因为这道题可以用最优子结构来推出最优答案，也符合常规DP套路即设 $ F[i][j] $ 表示将前面 $ i $ 个数分成 $ j $ 份，但是这一道题没有说具体的份数，而且数据范围很大，所以我们直接设 $ F[i] $ 表示前 $ i $ 个数的最优答案。这个的转移方程还是比较好想的：
$ F[i]=^{\quad\quad min}_{0<j<i,\sum^{k<i}_{k=j+1}A_k\leq M}{F[j]+max{A_k } } $
这样我们发现我们得到了一种 $ n^2 $ 的做法，然后数据范围告诉我们这样还不够。我们还得优化，但是 $ wch $ 是真的不会优化了，这个式子很有单调队列的样子，但是这个 $ max{A_k } $ ，这个式子似乎无懈可击啊！（于是他看了看书，表示：考这道题我直接爆零离场算了）
我们知道一个优秀的搜索或DP一定要尽可能去避免一些没有用的状态和决策，书上说：要保持决策集合的高度有效性和秩序性。所以我们看一看我们会不会枚举了一些没有什么用的 $ F[j] $ 。我们可以分析得出只有两种 $ F[j] $ 对我们有用：
- $ j $ 是满足 $ \sum^{k<i}_{k=j+1}A_k\leq M $ 的最小的 $ j $
- 我们的 $ A_j $ 是从 $ j $ 到 $ i $ 中间最大的权值
这个其实我们静下心来想想就可以明白，如果一个数不是从它开始到第 $ i $ 个数中最大的而且它不是可枚举的最小的第 $ j $ 个数，那么它一定不是最优解，它的左边那一个 $ F[j-1] $ 一定更优，因为它比 $ F[j ] $ 小，而且最大值一样。于是我们用维护一个包含所有的高效决策的单调队列，显然 $ j $ 单调递增 $ A_j $ 单调递减。然后我们要支持及时删去对头失效的决策，然后更新队尾加元素。
然后是这一题非常神奇的一个地方：我们要得到答案还需要 $ max{A_k } $ ，这个我们可以倍增 $ ST $ 表，但是我们分析题目性质就能发现一个跟奇妙的地方，我们的单调队列里的下一个元素就是我们要求得 $ max{A_k } $ ，这个只要想一想都知道以为我们维护的是（所有的高效决策的单调队列，并且显 $ j $ 单调递增 $ A_j $ 单调递减）。
然后还有一个我想吐槽的地方（奈何这一题数据水），我们的单调队列的对头不一定就是我们的最佳决策，这也是说我们队列里 $ A_K $ 单调递增，但是我们求值的式子里还有一个 $ F[j] $ 这个东西是单调递减的！！！！所以我们不能妄下决定，我们还需要映射一个取值的数据结构，也就是我们的 $ set $ ，我们要让这两个东西同时增加元素又同时减少相应的元素。
$ code: $
#include<iostream>
#include<cstdio>
#include<iomanip>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<ctime>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#define ll long long
#define db double
#define inf 0x7fffffff
#define rg register int
using namespace std;
int n; ll m;
int a[100005];
ll f[100005];
int q[100005];
multiset<ll >s;
inline ll qr(){
register char ch; register bool sign=0; ll res=0;
while(!isdigit(ch=getchar())) if(ch=='-')sign=1;
while(isdigit(ch)) res=res*10+(ch^48),ch=getchar();
return sign?-res:res;
}
int main(){
//freopen("in.in","r",stdin);
//freopen("cpp.out","w",stdout);
n=qr(); m=qr();
for(rg i=1;i<=n;++i) a[i]=qr();
rg l=1,r=0,k=1; ll tot=0;
for(rg i=1;i<=n;++i){
tot+=a[i];
while(tot>m)tot-=a[k++];
if(k>i){puts("-1");return 0;}
while(l<=r&&a[q[r]]<=a[i]){
if(r>l)s.erase(f[q[r-1]]+a[q[r]]);; --r;
}q[++r]=i; if(r>l)s.insert(f[q[r-1]]+a[q[r]]);
while(l<=r&&q[l]<k){
if(r>l)s.erase(f[q[l]]+a[q[l+1]]);; ++l;
}f[i]=f[k-1]+a[q[l]];
if(l<r)f[i]=min(f[i],*s.begin());
}printf("%lld\n",f[n]);
return 0;
}
