BZOJ[3675][Apio2014]序列分割斜率优化dp

最新推荐文章于 2019-11-30 17:51:23 发布

Duan2baka

最新推荐文章于 2019-11-30 17:51:23 发布

阅读量291

点赞数

分类专栏： BZOJ 斜率优化动态规划

本文链接：https://blog.csdn.net/waduan2/article/details/79721842

版权

BZOJ 同时被 3 个专栏收录

118 篇文章 0 订阅

订阅专栏

动态规划

12 篇文章 0 订阅

订阅专栏

斜率优化

5 篇文章 0 订阅

订阅专栏

传送门ber~

感谢Orz大佬的斜率优化讲解
因为得分跟顺序无关，所以可以dp
令 $f_{i,j}$ 表示第 $i$ 次在 $j$ 处分割的最大得分， $s_i$ 为 $\sum \limits _{j=1}^{j\le i} a_j$ ，则有

f i, j = m a x {f i - 1, k + (s j - s k) * (s n - s j)}

$f_{i,j}=max\{f_{i-1,k}+(s_j-s_k)*(s_n-s_j)\}$
状态

t t $t$ 比

k

$k$ 优，当且仅当

f i - 1, t + (s j - s t) * (s n - s j) > f i - 1, k + (s j - s k) * (s n - s j)

$f_{i-1,t}+(s_j-s_t)*(s_n-s_j)>f_{i-1,k}+(s_j-s_k)*(s_n-s_j)$
令

gj=fi−1,j g j = f i − 1 , j $g_j=f_{i-1,j}$ ，拆开括号，有

g t + s j * s n - s 2 j - s n * s t + s t * s j > g k + s j * s n - s 2 j - s n * s k + s k * s j

$g_t+s_j*s_n-s_j^2-s_n*s_t+s_t*s_j>g_k+s_j*s_n-s_j^2-s_n*s_k+s_k*s_j$

g t - g k - s n * s t + s n * s k > s j * (s k - s t)

$g_t-g_k-s_n*s_t+s_n*s_k>s_j*(s_k-s_t)$

g t - g k - s n * s t + s n * s k s k - s t > s j

$\frac {g_t-g_k-s_n*s_t+s_n*s_k}{s_k-s_t} > s_j$
再令

pj=gj−sj∗sn p j = g j − s j ∗ s n $p_j=g_j-s_j*s_n$ ，那么

p t - p k s t - s k < s j

$\frac {p_t-p_k}{s_t-s_k} < s_j$
这个东西就是斜率啊！
将每个状态表示为

(pi,si) ( p i , s i ) $(p_i,s_i)$ 的点对，维护一个下凸壳，因为

s s <script type="math/tex" id="MathJax-Element-33">s</script>是单调的，所以用单调队列搞搞就可以了

代码如下：

#include<algorithm>
#include<ctype.h>
#include<cstdio>
#define INF 2147483647000000ll
#define N 100050
using namespace std;
inline int read()
{
    int x=0,f=1;char c;
    do c=getchar(),f=c=='-'?-1:f; while(!isdigit(c));
    do x=(x<<3)+(x<<1)+c-'0',c=getchar(); while(isdigit(c));
    return x*f;
}
typedef long long LL;
int n,k,l,r;
LL s[N],f[N],g[N],ans;
struct Point
{
    LL x,y;
    Point(){}
    Point(LL _,LL __):x(_),y(__){}
    double operator * (const Point &b)///算斜率
    {
        if(x==b.x) return y>b.y ?-INF:INF;
        return 1.0*(b.y-y)/(b.x-x);
    }
}q[N];
inline LL calc(int k)
{
    return g[k]-s[k]*s[n];
}
inline LL GetAns(int k)
{
    while(l<r && q[l].x+q[l].y*s[k]<q[l+1].x+q[l+1].y*s[k]) l++;
    return q[l].x-s[k]*s[k]+s[k]*s[n]+q[l].y*s[k];
}
inline void Add(Point k)
{
    while(l<r && q[r-1]*q[r]>q[r]*k) r--;
    q[++r]=k;
}
main()
{
    n=read();k=read();
    for(int i=1;i<=n;i++)
        s[i]=read()+s[i-1];
    for(int i=1;i<=k;i++)
    {
        l=r=1;
        for(int j=1;j<=n;j++)
        {
            f[j]=GetAns(j);
            Add(Point(calc(j),s[j]));
        }
        for(int j=1;j<=n;j++)
            g[j]=f[j];
    }
    for(int i=1;i<n;i++)
        ans=max(ans,g[i]);
    printf("%lld",ans);
    return 0;
}