【DP计划】11.2——[CF]CF311B（斜率优化） EASY+

小 S 是农场主，他养了 $M$ 只猫，雇了 $P$ 位饲养员。农场中有一条笔直的路，路边有 $N$ 座山，从 $1$ 到 $N$ 编号。第 $i$ 座山与第 $i-1$ 座山之间的距离是 $D_i$​​ 。饲养员都住在 $1$ 号山上。

有一天，猫出去玩。第 iii 只猫去 $H_i$号山玩，玩到时刻 $T_i$停止，然后在原地等饲养员来接。饲养员们必须回收所有的猫。每个饲养员沿着路从 $1$ 号山走到 $N$ 号山，把各座山上已经在等待的猫全部接走。饲养员在路上行走需要时间，速度为 $1$ 米每单位时间。饲养员在每座山上接猫的时间可以忽略，可以携带的猫的数量为无穷大。

例如有两座相距为 $1$ 的山，一只猫在 $2$ 号山玩，玩到时刻 $3$ 开始等待。如果饲养员从 $1$ 号山在时刻 $2$ 或 $3$ 出发，那么他可以接到猫，猫的等待时间为 $0$ 或 $1$。而如果他于时刻 $1$ 出发，那么他将于时刻 $2$ 经过 $2$ 号山，不能接到当时仍在玩的猫。

你的任务是规划每个饲养员从 $1$ 号山出发的时间，使得所有猫等待时间的总和尽量小。饲养员出发的时间可以为负。
输入格式
第一行三个整数 $N,M,P$；
第二行 $N-1$ 个正整数 $D_i$​​ ，表示第 $i$ 座山与第 $i-1$ 座山之间的距离是 $D_i$；
接下去 $M$ 行每行两个整数 $H_i,T_i$。
输出格式
输出一个整数表示答案。
样例
样例输入

4 6 2
1 3 5
1 0
2 1
4 9
1 10
2 10
3 12
样例输出

3
数据范围与提示
对于全部数据，$2\le N\le 10^5,1\le M\le 10^5,1\le p\le 100,1\le D_i<10^4,1\le H_i\le N,0\le T_i\le 10^9$

---

我们可以先做一个处理（dis[i]表示从第1座山到第i座山的距离）

	for(int i=1;i<=m;i++){
		F[i].p=read();F[i].num=read()-dis[F[i].p];
	}

因为在第$i$座山玩到时刻$T$和在第一座山玩到时刻$T-dis[i]$是一样的（$dis[i]$就是管理员要走到的时间）

这样就变成所有猫都在第一座山，在不同的时间结束。总共P个管理员，当管理员$i$在时间点$T_i$，上一个管理员在时间点$T_{last}$时，那么$[T_{last},T_i]$之间的猫都由当前管理员$i$带走。所以我们就可以轻松的写出方程：

f[i][j]=min(f[i][j],f[p][j-1]+F[i].num*(i-p)-(sum[i]-sum[p]))
//f[p][j-1]右边这一串就表示当前区间内猫的等待时间，因为每只猫等的时间为F[i].num-F[x].num
//F[i].num代表管理员出发的时间(显然与当前区间的最后一直猫同时最优)，将这个式子求个Σ，用前缀和记录，就是上面的DP

那么一顿斜率操作之后（具体的操作可以看我之前的博客）就变成了

f[i][j]+F[i].num*p=f[p][j-1]+F[i].num*i-sum[i]+sum[p]
b    +      k     *x=                                 y

接下来用单调队列维护一个下凸包即可。

#include<bits/stdc++.h>
#define MAXN 100005
#define db double
#define ll long long
using namespace std;
ll read(){
	char c;ll x;while(c=getchar(),c<'0'||c>'9');x=c-'0';
	while(c=getchar(),c>='0'&&c<='9') x=x*10+c-'0';return x;
}
ll n,m,p,w,r,d[MAXN],dis[MAXN],sum[MAXN],f[MAXN][105],q[MAXN];
struct node{
	ll p,num;
	bool operator<(const node x)const{
		return num<x.num;
	}
}F[MAXN];
db x(ll i){return i;}
db y(ll i,ll p){return f[i][p]+sum[i];}
db k(ll i,ll j,ll p){
	if(x(i)==x(j)) return 2e9;
	return (y(i,p-1)-y(j,p-1))/(x(i)-x(j));
}
int main()
{
	n=read();m=read();p=read();
	for(int i=2;i<=n;i++) d[i]=read(),dis[i]=dis[i-1]+d[i];
	for(int i=1;i<=m;i++){
		F[i].p=read();F[i].num=read()-dis[F[i].p];
	}
	sort(F+1,F+1+m);
	for(int i=1;i<=m;i++)sum[i]=sum[i-1]+F[i].num;
	memset(f,127,sizeof(f));f[0][0]=0;
	for(int j=1;j<=p;j++){
		w=r=0;
		for(int i=1;i<=m;i++){
			while(r<w&&F[i].num>k(q[r],q[r+1],j)) r++;
			f[i][j]=f[q[r]][j-1]+F[i].num*i-F[i].num*q[r]-sum[i]+sum[q[r]];
			while(r<w&&k(q[w-1],q[w],j)>k(q[w],i,j)) w--;
			q[++w]=i;
		}
	}
	printf("%lld",f[m][p]);
	return 0;
}