POJ3709 K-Anonymous Sequence(k匿名序列)（斜率优化dp）

poj3709题目传送
看到这名字第一反应是那个面具
斜率优化的动态规划，第一次做斜率优化的dp题。。。而且自己也不会，是看着书慢慢弄懂的。书是《挑战》p341 ，这一面有很多错误，刚开始自己写完后照书上一对，一模一样后发现答案很不对，后来才发现书上的状态转移方程存在错误（在下面的公式里你自己对一遍就知道了） 总之就是，书上的不一定对，自己领悟并理解才最重要。

首先不难得出，$a_0$是最小的值，为了使操作的次数较小，我们不可能把最小值拿着减。所以我们发现，对于$a_0$这样的较小的项，我们可以选择从小到大来选择需要减少到$a_0$的项。设$dp_i$为只考虑前i项的前提下最少的操作次数，则有$dp_i$=min{$dp_j$+($a[j+1]$-$a[j]$)+…+($a[i-1]$-$a[j]$)},其中$0\le j \le i-k$。

直接计算自然没问题，但是$O^3$的复杂度不太好的样子，于是考虑优化：

例如我们记一个前缀和，

$$S_i=a[0]+...+a[i-1]$$
有 $dp_i$=min{ $dp_j$+ $S[i]$- $S[j+1]$- $a_j$*(i-j-1)}
观察方程我们可以发现,若我们设一个函数（用x替换i） $f_j(x)=-a_j*x+dp[j]-S[j+1]+a_j*(j+1)$，则有 $dp[i]=S[i]+min{f_j(i)}$,于是计算dp[n]就变为了在i-k+1条直线中寻找x=i的最小值就好了。
于是我们考虑用双端队列去维护所有可能成为最小值的直线的集合，并以从小到大的顺序排列。那么我们如何判断要去掉某条直线呢？

设：

$f_1(x)=a_1*x+b_1$
$f_2(x)=a_2*x+b_2$
$f_3(x)=a_3*x+b_3$

我们有$(a_2-a_1)*(b_3-b_2)\ge(b_2-b_1)*(a_3-a_2)$
时，$f_2$为要排除的直线，为什么呢？
简证：当$f_1$与$f_3$的交点在$f_2$以下时，联立$f_1$与$f_3$
得$x_0=(\frac{b_1-b_3}{a_3-a_1})$
然后对于$f_1$,$f_3$,有$f(x_0)=a_1(\frac{b_1-b_3}{a_3-a_1})+b_1=(\frac{b_1*a_3-b_3*a_1}{a_3-a_1})$
而$f_2(x_0)=a_2*(\frac{b_1-b_3}{a_3-a_1})+b_2$
当$f_2(x_0)\le f_1(x_0)和 f_3(x_0)$时，有$a_2*(\frac{b_1-b_3}{a_3-a_1})+b_2\le (\frac{b_1-b_3}{a_3-a_1})$
拿纸来算吧 化简就有$(a_2-a_1)*(b_3-b_2)\ge(b_2-b_1)*(a_3-a_2)$了。

所以我们就可以根据上述方法来做双端队列，详见标程：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>

using namespace std;
typedef long long ll;
const ll maxn=505000;
ll cas=0;
ll n,k;
ll a[maxn];
ll deq[maxn];//虽然是双端队列 但存的是值
ll dp[maxn*10];
ll sh[maxn];//前缀和

bool check(ll f1,ll f2,ll f3)
{
    ll a1=-a[f1];
    ll b1=dp[f1]-sh[f1+1]+a[f1]*(f1+1);
    ll a2=-a[f2];
    ll b2=dp[f2]-sh[f2+1]+a[f2]*(f2+1);
    ll a3=-a[f3];
    ll b3=dp[f3]-sh[f3+1]+a[f3]*(f3+1);
    return (a2-a1)*(b3-b2)>=(b2-b1)*(a3-a2);
}

ll f(ll j,ll x)
{
    return -a[j]*x+dp[j]-sh[j+1]+a[j]*(j+1);
}

int main()
{
    //freopen("1.in","r",stdin);
    //freopen("2.out","w",stdout);
    scanf("%lld",&cas);
    for(ll z=0;z<cas;z++)
    {
        memset(sh,0,sizeof(sh));
        memset(a,0,sizeof(a));
        memset(dp,0,sizeof(dp));
        memset(deq,0,sizeof(deq));
        scanf("%lld%lld",&n,&k);
        for(ll i=0;i<n;i++)
        {
            scanf("%lld",&a[i]);
        }
        for(ll i=0;i<n;i++)
        {
            sh[i+1]=sh[i]+a[i];
        }

        ll s=0;
        ll t=1;
        deq[s]=0;
        dp[0]=0;
        for(ll i=k;i<=n;i++)
        {
            if(i-k>=k)
            {
                while(s+1<t && check(deq[t-2],deq[t-1],i-k))
                {
                    t--;
                }
                deq[t]=i-k;
                t++;
            }
            while(s+1<t && f(deq[s],i)>=f(deq[s+1],i))
            {
                s++;
            }
            dp[i]=sh[i]+f(deq[s],i);
        }
        printf("%lld\n",dp[n]);
    }
    return 0;
}

书上没看懂得就来看一下这个吧QAQ