HDU 3507 Print Article（斜率优化dp）

题目：http://acm.split.hdu.edu.cn/showproblem.php?pid=3507
题意：N words，每个word有Ci，如果连续输出k个，那么代价就是求最少花费
这题用来入门斜率优化dp

g[i,j] = getUP(i,j) / getDOWN(i,j) < sum[i] i的决策优于j的决策 斜率优化的地方：
设k < <script type="math/tex" id="MathJax-Element-19"><</script>j < <script type="math/tex" id="MathJax-Element-20"><</script>i，如果g[i,j] < <script type="math/tex" id="MathJax-Element-21"><</script>g[j,k]，我们假设g[i,j] < <script type="math/tex" id="MathJax-Element-22"><</script>sum[i]，那么就是说i点要比j点优，排除j点。
如果g[i,j]$>=$sum[i]，那么j点此时是比i点要更优，但是同时g[j,k]$>$g[i,j]$>$sum[i]。这说明还有k点会比j点更优，同样排除j点。

思路：

于是对于这题我们对于斜率优化做法可以总结如下：
1，用一个单调队列来维护解集。
2，假设队列中从头到尾已经有元素a b c。那么当d要入队的时候，我们维护队列的上凸性质，即如果g[d,c] < <script type="math/tex" id="MathJax-Element-26"><</script>g[c,b]，那么就将c点删除。直到找到g[d,x]>=g[x,y]为止，并将d点加入在该位置中。
3，求解时候，从队头开始，如果已有元素a b c，当i点要求解时，如果g[b,a] < <script type="math/tex" id="MathJax-Element-27"><</script>sum[i]，那么说明b点比a点更优，a点可以排除，于是a出队。最后dp[i]=getDp(q[head])。

代码：

#include<bits/stdc++.h>
using namespace std;

const int N = 500005;
int n,M;
int dp[N],sum[N],q[N];
int head,tail;

int getUP(int j,int k)
{
    return dp[j] + sum[j]*sum[j] - (dp[k] + sum[k]*sum[k]);
}
int getDOWN(int j,int k)
{
    return (sum[j] - sum[k])*2;
}
int getDP(int i,int j)
{
    return dp[j] + M + (sum[i] - sum[j])*(sum[i] - sum[j]);
}
int main()
{
    while(~scanf("%d%d",&n,&M))
    {
        int x;
        for(int i = 1;i <= n;i++)
            scanf("%d",&x),sum[i] = sum[i-1] + x;
        head = tail = 0;
        q[tail++] = 0;
        dp[0] = 0;
        for(int i = 1;i <= n;i++)
        {
            while(head + 1 < tail && getUP(q[head+1],q[head]) <= getDOWN(q[head+1],q[head])*sum[i])
                head++;//g[b,a]<sum[i]，那么说明b点比a点更优，出队
            dp[i] = getDP(i,q[head]);
            while(head + 1 < tail && getUP(i,q[tail-1])*getDOWN(q[tail-1],q[tail-2]) <= getUP(q[tail-1],q[tail-2])*getDOWN(i,q[tail-1]))
                tail--;
            q[tail++] = i;//找到了更优的，入队
        }
        printf("%d\n",dp[n]);
    }
    return 0;
}