数学（斜率优化dp）

数学

$\ $
$\ $
$\ $

$\ $$\ $

$\ $
$\ $
$\ $

解：

还是数学太弱了-_-~
虽然是一眼斜率优化，然而并不会推式子，啊啊啊啊！！！
首先dp方程是很好推的，$f_{i,k}=min\{f_{j,k-1}+\frac {a_{j+1}}{a_{j+1}}+\frac {a_{j+1}+a_{j+2}} {a_{j+2}} +…\frac {\sum^i_{n=j+1} a_n}{a_i}\}$
其实这一道题的dp方程是很简单的，照样子推推式子就行，只不过我们要完成的就是O1求：
$\frac {b_1}{b_1}+\frac {b_1+b_2} {b_2} +…\frac {\sum^m_{i=1} b_i}{b_m}$
首先我们是可以$O(n)$求从1开始的这一串式子，但是我们还要求从每一个起点开始的这一串。
所以设这一串式子到$i$点的值为$ans(i)$，假设我们要求$q(j+1,i)$
$ans(i)-ans(j)=\frac {\sum^{j+1}_{n=1} b_n}{b_j+1} +…\frac {\sum^i_{n=1} b_i}{b_i}$
那么我们多求的就是$\frac {\sum^j_{n=1} b_n}{b_j+1} +…\frac {\sum^j_{n=1} b_n}{b_i}$
于是我们把分子提出来：
$\sum^j_{n=1} b_n*(\frac 1 {b_{j+1}}+\frac 1 {b_{j+2}}+...+\frac 1 {b_{i}})$
于是我们想到预处理倒数前缀和($num(i)$)和前缀和($p(i)$)。
所以$f_{i,k}=min\{f_{j,k-1}+ans(i)-ans(j)-num(j)*(p(i)-p(j)\}$
然后就可以推斜率优化了。
附上代码（50个单调队列，容易写错，思路一定要清晰）：

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;

double data[200005];
double num[200005];
double p[200005];
double f[55][200005];
double ans[200005];
int n,k;
int tail[55],head[55];
int d[55][200005];
int now;

double X(int id)
{
    return f[now][id]-ans[id]+num[id]*p[id];
}
double slop(int x,int y)
{
    return (X(x)-X(y))/(num[x]-num[y]);
}

int main()
{
    scanf("%d%d",&n,&k);
    for(int i=0;i<=k;i++)
      head[i]=1;
    for(int i=1;i<=n;i++)
      scanf("%lf",&data[i]);
    for(int i=1;i<=n;i++)
      num[i]=num[i-1]+data[i],p[i]=p[i-1]+1/data[i];
    for(int i=1;i<=n;i++)
      ans[i]=ans[i-1]+num[i]/data[i];
    tail[0]++;
    for(int i=1;i<=n;i++)
      for(int j=k;j>=1;j--)
      {
        now=j-1;
        if(tail[now]-head[now]<0) continue;
        while(tail[now]-head[now]>0&&slop(d[now][head[now]],d[now][head[now]+1])<p[i])
          head[now]++;
        f[j][i]=f[now][d[now][head[now]]]+ans[i]-ans[d[now][head[now]]]-num[d[now][head[now]]]*(p[i]-p[d[now][head[now]]]);
        now=j;
        while(tail[j]-head[j]>0&&slop(d[j][tail[j]-1],d[j][tail[j]])>slop(d[j][tail[j]],i))
          tail[j]--;
        d[j][++tail[j]]=i;
      }
    printf("%lf",f[k][n]);
}

下次斜率优化别爆炸了！！！（其实还是数学太菜了）