bzoj3675: [Apio2014]序列分割（Dp）

题目传送门
。

解法：
emmmm一开始理解错题意了写了个错误的方程。
亮老师走过来发现我的方程错了。
赶紧指导一波嘻嘻嘻才理解题意。。
其实仔细想下他不过是分成几段然后每段之间的乘积加上每段的和。
f[i][k]表示前i个数切k次的最大获利。
那么f[i][k]=max(f[j][k]+(s[i]-s[j])*s[j]就行了。

斜率优化：
若j>t
且j相对于i来说更优。
那么有
f[j][k-1]+(s[i]-s[j])*s[j]>f[t][k-1]+(s[i]-s[t])*s[t]
展开后移项，懒得写详细了

(f[j][k-1]-f[t][k-1]+s[t]^2-s[j]^2)/(s[t]-s[j])＜s[i]
然后单调队列维护下凸包即可

代码实现：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
ll s[110000],f[2][110000];int k,head,tail,list[110000];
ll slop(int j,int t) {return (f[1-k][j]-f[1-k][t]+s[t]*s[t]-s[j]*s[j])/(s[t]-s[j]);}
int main() {
    int n,K;scanf("%d%d",&n,&K);s[0]=0ll;
    for(int i=1;i<=n;i++) {ll x;scanf("%lld",&x);s[i]=s[i-1]+x;}
    k=1;memset(f,0,sizeof(f));
    for(int t=1;t<=K;t++) {
        k=1-k;head=1;tail=0;list[1]=0;
        for(int i=t;i<=n;i++) {
            while(head<tail&&(s[list[head]]==s[list[head+1]]||slop(list[head+1],list[head])<s[i]))head++;
            int l=list[head];f[k][i]=f[1-k][l]+(s[i]-s[l])*s[l];
            while(head<tail&&(s[i]==s[list[tail]]||slop(i,list[tail])<slop(list[tail],list[tail-1])))tail--;
            list[++tail]=i;
        }
    }printf("%lld\n",f[k][n]);
    return 0;
}