斜率优化

斜率优化

---

前置技能点：

• DP初步
• 单调队列
• 求凸包

闲话：

前面写这么多都是为了写这个做准备的…本来还应该写一下凸包，但凸包算法比较多，迟点吧…

问题：

给定非负整数序列$a$长度为$n$，状态转移方程为$DP[i]={\min }_{k=1}^i(DP[k]+(sum[i]-sum[k]{)}^2+C)$，求DP[n]

思路：

最朴素的做法，有状态方程就对着做，$for(i)$套$for(k)$，复杂度$O(n^2)$

概念引入：

设$D{P}_k[i]=DP[k]+(sum[i]-sum[k]{)}^2+C$，$j<p$且$D{P}_p[i]$≤$D{P}_j[i]$

---

$D{P}_j[i]\ge D{P}_p[i]$

$DP[j]+{(sum[i]-sum[j])}^2+C\ge DP[p]+{(sum[i]-sum[p])}^2+C$ 拆开

$DP[j]-2\ast sum[i]\ast sum[j]+{sum[j]}^2\ge DP[p]-2\ast sum[i]\ast sum[p]+{sum[p]}^2$ 拆开平方项并消去常数$sum[i{]}^2$和$C$

$(DP[j]+{sum[j]}^2)-(DP[p]+{sum[p]}^2)\ge 2\ast sum[i]\ast (sum[j]-sum[p])$ 移项

---

设$f(x)=DP[x]+{sum[x]}^2$，$g(x)=2\ast sum[x]$，$h(x)=sum[x]$

则原式等价于$f(j)-f(p)\ge h(i)\ast (g(j)-g(p))$

由于g(x)单调不降

两边同除以$(g(j)-g(p))$（事实上可能为0，不能除，这里是为了看上去直观一点）得到：

$\frac{f(j)-f(p)}{g(j)-g(p)}\le h(i)$

而由于$h(i)$单调不降，$g(j)-g(p)\le 0$故当$i$增大时，不等式恒成立

也即是一旦对于$DP[i]$，不等式成立，则对于任意$i^{\prime }>i$不等式恒成立，即决策点p永远优于决策点j

定义使得$D{P}_k[x]$取得最小值的点$k$为$p[x]$(perfect)。

这意味着，对于要得到的$DP[i]$，我们只需要将$p[i-1]$到$i-1$的所有元素遍历一遍即可得到$p[i]$。

目前还属于常数剪枝，不影响复杂度…

---

现在，以$g(x)$为$x$轴，$f(x)$为$y$轴建立直角坐标系，由于$g(x)$递增，故输入的每个元素的顺序与其在x轴上的顺序相同。

将所有下标为$p[i-1]$到$i-1$的点画到图上，方便起见，就画三个点$l,m,r$。

回顾之前的式子

$\frac{f(j)-f(p)}{g(j)-g(p)}\le h(i)$

不难发现左边部分在这个直角坐标系上就是$j,p$两点间的斜率。

设$k_1=k(l,m)$，$k_2=k(m,r)$，$k_3=k(l,r)$。

一共有两种情况：

• $k_1\ge k_2$

• $k_1<k2$
  
  我们分开来讨论

• 情况 1
  $k_1\ge k_3\ge k_2$
  对于一个给定的$h(i)$，若
  $h(i)\ge k_1$，则$m$点比$l$点更优，但同时$h(i)\ge k_2$，故$r$点在三点中最优。
  $k_1>h(i)\ge k_3$，则$r$点比$l$点更优，同时$h(i)\ge k_2$，故$r$点在三点中最优。
  $k_3>h(i)\ge k_2$，则$r$点比$m$点更优，但$h(i)<k_3$，故$l$点在三点中最优。
  $k_2>h(i)$，则$l$点在三点中最优。
  也即是当出现这种情况时，$m$点根本没有判定的必要，只需要判定$k(l,r)$与$h(i)$的关系即可。

• 情况 2
  $k_2＞k_3＞k_1$
  对于一个给定的$h(i)$，若
  $h(i)\ge k_2$，则$r$点比$m$点更优，同时$h(i)\ge k_1$，故$r$点在三点中最优。
  $k_2>h(i)\ge k_3$，则$r$点比$l$点更优，但$h(i)<k_2$，故$m$点在三点中最优。
  $k_3>h(i)\ge k_1$，则$m$点比$l$点更优，同时$h(i)<k_2$，故$m$点在三点中最优。
  $k_1>h(i)$，则$l$点在三点中最优。
  也即是当出现这种情况时，我们不需要考虑$k(l,r)$，只需要考虑相邻两点之间的斜率即可。

得到了上面的结论后，我们可以用水平法求下凸包进行优化，但是点$p[i-1]$可能会因为$i-1$的进入而在做下凸包时被删除。
所以我们需要知道新加的元素导致凸包尾被删除到了哪个点，在进行$p[i]$的求值时，以此判断是从$p[i-1]$还是直接从单调栈中最后两个元素开始对斜率进行判定。

此外还有另一个办法，我们注意到即使存在点$j<p[i-1]<i-1$使得出现之前讨论过的情况1，但实际上从点$p[i-1]$到点$i-1$的斜率一样足以让我们判定点$i-1$优于点$p[i-1]$（$h(i)\ge h(i-1)\ge k_1$）

所以我们可以通过直接移动栈底，即单调队列的弹出队头操作来维护得到$p[i]$到$i-1$之间的下凸包，而非从$1$到$i-1$之间的下凸包。

两个方法的代码几乎没有不同，第二个可能快点吧，但第一种方法能顺便得到下凸包。

复杂度：

时间复杂度$O(n)$，空间复杂度$O(n)$

例题

• hdu 3507
  就是上面的问题，给定非负整数序列$a$长度为$n$，状态转移方程为$DP[i]={\min }_{k=1}^i(DP[k]+(sum[i]-sum[k]{)}^2+C)$，求DP[n]
  解题报告