运输小猫（斜率优化）

斜率优化的题目真的是不好做呢。

历程：知道是斜率优化，但是还是不知道该怎样定义状态，以及进行状态之间的转移，理解这道题目之后，我发现dp和数学基础真的很有关系，从这道题目我有一点启发。

启发：首先对于题目的理解一定要到位，这道题目一开始说实话我都没怎么理解这里边的含义，可能还是做题太少了吧。首先是状态的定义，我们首先对题目的要求的求得的结果进行分析，这个结果的所代表的值是什么，是不是可以分解，还有就是对结果产生影响的因素是什么，从因素入手！

题目：
小S是农场主，他养了 M 只猫，雇了 P 位饲养员。

农场中有一条笔直的路，路边有 N 座山，从 1 到 N 编号。

第 i 座山与第 i-1 座山之间的距离为 Di。

饲养员都住在 1 号山。

有一天，猫出去玩。

第 i 只猫去 Hi 号山玩，玩到时刻 Ti 停止，然后在原地等饲养员来接。

饲养员们必须回收所有的猫。

每个饲养员沿着路从 1 号山走到 N 号山，把各座山上已经在等待的猫全部接走。

饲养员在路上行走需要时间，速度为1 米/单位时间。

饲养员在每座山上接猫的时间可以忽略，可以携带的猫的数量为无穷大。

例如有两座相距为 1 的山，一只猫在 2 号山玩，玩到时刻 3 开始等待。

如果饲养员从 1 号山在时刻 2 或 3 出发，那么他可以接到猫，猫的等待时间为 0 或 1。

而如果他于时刻 1 出发，那么他将于时刻 2 经过 2 号山，不能接到当时仍在玩的猫。

你的任务是规划每个饲养员从 1 号山出发的时间，使得所有猫等待时间的总和尽量小。

饲养员出发的时间可以为负。

输入格式
第一行包含三个整数N，M，P。

第二行包含n-1个整数，D2,D3,…,DN。

接下来M行，每行包含两个整数Hi和Ti。

输出格式
输出一个整数，表示所有猫等待时间的总和的最小值。

数据范围
2≤N≤10⁵,
1≤M≤10⁵,
1≤P≤100,
1≤Di<1000,
1≤Hi≤N,
0≤Ti≤10⁹
输入样例：
4 6 2
1 3 5
1 0
2 1
4 9
1 10
2 10
3 12
输出样例：
3

这道题目要求求得的是小猫的最小等待时间，我们首先观察下影响小猫等待时间的因素，首先是饲养员的出发时间，还有就是距离，还有就是小猫的玩耍时间，因为小猫必须要等到玩耍时间过了之后才可以回家，真的是调皮的小猫~找到了小猫等待时间的因素，我们就来看那些因素是我们可以控制的，发现只有饲养员的出发时间，并且题目中说了只要要求规划饲养员的出发时间，并且饲养员的出发时间可以为负数，那么进行下一步，怎么规划饲养员的出发时间可以是小猫的等待时间最短，暴力枚举每个饲养员的出发时间？这道题目的范围可是1e9
呀，不行不行，那么我们就缩小集合，想一下哪些出发时间可能是最优解，（这里是是解题的关键）在思考中应该可以发现每位饲养员走完后接到的小猫，其中一定有一只或几只小猫等待的时间为0， 怎样证明呢，我们可以用反证法，如果饲养员接到小猫等待时间都不是为0，那么我们就可以将出发的时间向前推，使得你某些小猫的等待时间减少为0，并且携带的小猫不变，并且题目中已经明确出饲养员的出发时间可以为负数。这样就证出了每位饲养员的出发时间肯定是某些小猫的玩耍恰好结束的时间和饲养员在路上走路的时间的差值，这样就缩小了枚举的范围，但是。。饲养员的个数是100个，小猫的个数是1e5，由于是组合数的复杂度，（组合数的话就有点dp那味了哈哈）我们肯定也是过不了的。由于我们已经得到饲养员的出发时间的集合，那么我们只需要设置状态f[i][j]为用i个饲养员运送j值小猫的最小值，得到表达式，可以发现这是可以用到斜率优化的方式进行优化的。在进行计算饲养员出发时间可以接走那些小猫的时候，我们可以用预先排序进行优化，因为接走的一组小猫的等待时间最小值就是0，我们可以按照小猫的玩耍恰好结束的时间和饲养员在路上走路的时间的差值进行排序，其实做到这里就差不多应该已经知道要用到斜率优化了，这是一个斜率优化很常见的模型。（当然一直以暴力的思路想的目的并不是要用暴力来解决问题，而是用一种特别的暴力手段，dp）

式子还是写一下吧：

首先是状态转移的式子：

$f[i][j]=max(f[i-1[k]+a[j]\ast (j-k)-(s[j]-s[k])$ 这里的a数组表示的是小猫的玩耍时间与距离的差值，s数组表示的是a数组的前缀和。

转化成斜率式子：

$f[i-1][k]+s[k]=a[j]\ast k+f[i][j]-a[j]\ast j+s[j]$

代码：

#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

typedef long long LL;

const int N = 1e5 + 10, M = 110;

LL d[N], t[N], f[M][N], s[N], a[N];

int q[N], n, m, p;

LL get(int i, int k)
{
    return f[i - 1] [k] + s[k];
}

int main()
{
    scanf("%d %d %d",&n, &m, &p);
    
    for (int i = 2; i <= n; i ++) 
    {
        scanf("%lld", &d[i]);
        
        d[i] += d[i - 1];
    }
    
    for (int i = 1; i <= m; i ++) 
    {
        int h;scanf("%d%lld", &h, t + i);
        
        a[i] = t[i] - d[h];
    }
    sort(a + 1, a + m + 1);
    
    for (int i = 1; i <= m; i ++) s[i] = s[i - 1] + a[i];
    
    memset(f, 0x3f, sizeof f);
    
    for (int i = 0; i <= p; i ++) f[i][0] = 0;
    for (int i = 1; i <= p; i ++)
    {
        int hh = 0, tt = 0;
        q[0] = 0;
        for (int j = 1; j <= m; j ++)
        {
            while (hh < tt && (get(i, q[hh + 1]) - get(i, q[hh])) <= a[j] * (q[hh + 1] - q[hh])) hh ++;
            int k = q[hh];
            f[i][j] = f[i - 1][k]  + a[j] * (j - k) - s[j] + s[k];
            while (hh < tt && (get(i, q[tt]) - get(i, q[tt - 1])) * (j - q[tt - 1]) >= (get(i, j) - get (i, q[tt - 1])) * ((q[tt] - q[tt - 1])))tt --;
            
            q[ ++ tt] = j;
        }
    }
    printf("%lld", f[p][m]);
    return 0;
}