【ybt金牌导航1-5-4】【luogu CF311B】【LOJ 10187】猫的运输 / Cats Transport（斜率优化）

猫的运输 / Cats Transport

题目链接：ybt金牌导航1-5-4 / luogu CF311B / LOJ 10187

题目大意

有一些猫，它们各在一条链的某条位置上，会在某个时刻出现，之后就会等待。
然后有一些人在链的最左边，以一个单位长度每单位时间往右走，然后他会带走他所在位置出现的所有猫。每个人的开始时间任意，可以是负数都行。
然后问你所有猫最小的等待时间。

思路

首先我们看到这个小猫它很烦，又有位置又有出现时间。
然后我们看到饲养员都在最左边，且开始行走时间任意，那我们考虑它每个小猫的要记录的值即为饲养员要在什么时刻或之后出发才能接到它。
容易得出，设此值为 $a_i$，则有 $a_i=t_i-d_i$。
（$t_i$ 为它出现时间，$d_i$ 为它到最左边的距离，这个可以通过前缀和预处理出）

那我们会发现，如果我们把猫按 $a_i$ 从小到大排，每个饲养员能拿到的猫就是从某个开始往后的全部。
那越靠后的猫如果让他来收就等的时间越长，那就要在后面继续放饲养员。
那每个饲养员拿猫的范围就是它到下一个饲养员的前面一直猫。
然后你就会觉得它可以 DP，于是我们来搞搞看：

设 $f_{i,j}$ 为 $i$ 个饲养员处理前 $j$ 值小猫的最小时间。
然后可以列出转移方程：
f i , j = min ⁡ k = 0 j − 1 { f i − 1 , k + a j ∗ ( j − k ) − ( S j − S k ) } f_{i,j}=\min_{k=0}^{j-1}\{f_{i-1,k}+a_j*(j-k)-(S_j-S_k)\} fi,j​=mink=0j−1​{fi−1,k​+aj​∗(j−k)−(Sj​−Sk​)}
（$S_i$ 是 $a_i$ 数组的前缀和）

然后你会发现这样枚举会超时，那我们考虑优化，先拆括号看看：
f i , j = min ⁡ k = 0 j − 1 { f i − 1 , k + a j ∗ j − a j ∗ k − S j + S k } f_{i,j}=\min_{k=0}^{j-1}\{f_{i-1,k}+a_j*j-a_j*k-S_j+S_k\} fi,j​=mink=0j−1​{fi−1,k​+aj​∗j−aj​∗k−Sj​+Sk​}
看到这个 $\min$ 和 $a_j\ast k$，我们可以先到斜率优化，那就试试咯。

那枚举好 $i,j$，当 $x$ 比 $y$ 优。（$a>b$）
$f_{i-1,x}+a_j\ast j-a_j\ast x-S_j+S_x<f_{i-1,y}+a_j\ast j-a_j\ast y-S_j+S_y$
$f_{i-1,x}-f_{i-1,y}+S_x-S_y<a_j\ast x-a_j\ast y$
$\frac{(f_{i-1,x}+S_x)-(f_{i-1,y}+S_y)}{x-y}<a_j$

然后你就可以开搞了。
由于你前面把 $a_j$ 数组从小到大排序了，那你就直接用单调队列搞就可以了。

---

发现之前是枚举猫再枚举人的。
这样的话要开 $q$ 个单调队列，其实不太优。
（因为今天做了一题 $q=1000$ 的，$f$ 可以滚动但是 $q$ 不能）

所以你可以枚举人，再枚举猫，这样就就可以只用一个单调队列，然后就可以滚动 $f$，然后就可以做 $q=1000$ 了。

代码

之前

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long 

using namespace std;

ll n, m, p, d[100001];
ll h, t, a[100001], s[100001];
ll l[100001], r[100001], q[101][100001], f[201][100001];

bool cmp(ll x, ll y) {
	return x < y;
}

double xl(int i, int j1, int j2) {//计算斜率
	return (1.0 * f[i][j1] - f[i][j2] + s[j1] - s[j2]) / (j1 - j2);
}

int main() {
	scanf("%lld %lld %lld", &n, &m, &p);
	for (ll i = 1; i < n; i++) scanf("%lld", &d[i]), d[i] += d[i - 1];
	
	for (ll i = 1; i <= m; i++) {
		scanf("%lld %lld", &h, &t);
		a[i] = t - d[h - 1];
	}
	
	sort(a + 1, a + m + 1, cmp);//按出发时间排序
	for (ll i = 1; i <= m; i++)
		s[i] = s[i - 1] + a[i];
	
	for (ll j = 1; j <= m; j++) {
		for (ll i = 1; i <= min(j, p); i++) {//饲养员的个数一定要套在里面
			while (l[i - 1] < r[i - 1] && xl(i - 1, q[i - 1][l[i - 1]], q[i - 1][l[i - 1] + 1]) <= a[j]) l[i - 1]++;
			f[i][j] = f[i - 1][q[i - 1][l[i - 1]]] + a[j] * j - a[j] * q[i - 1][l[i - 1]] - s[j] + s[q[i - 1][l[i - 1]]];
			while (l[i] < r[i] && xl(i, q[i][r[i]], j) < xl(i, q[i][r[i] - 1], q[i][r[i]])) r[i]--;
			q[i][++r[i]] = j;//从饲养员个数是i-1转移到个数是i
		}
	}
	
	printf("%lld", f[p][m]);
	
	return 0;
}

新的

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long 

using namespace std;

const int N = 1e5 + 1000;
ll n, m, p, d[N];
ll h, t, a[N], s[N];
ll l, r, f[2][N];
ll sta[N];

bool cmp(ll x, ll y) {
	return x < y;
}

int main() {
	scanf("%lld %lld %lld", &n, &m, &p);
	for (ll i = 1; i < n; i++) scanf("%lld", &d[i]), d[i] += d[i - 1];
	
	for (ll i = 1; i <= m; i++) {
		scanf("%lld %lld", &h, &t);
		a[i] = t - d[h - 1];
	}
	
	sort(a + 1, a + m + 1, cmp);
	for (ll i = 1; i <= m; i++)
		s[i] = s[i - 1] + a[i];
	
	memset(f, 0x7f, sizeof(f)); f[0][0] = 0;
	for (int i = 1; i <= p; i++) {
		l = r = 1; sta[1] = 0;
		for (int j = 1; j <= m; j++) {
			while (l < r && f[(i - 1) & 1][sta[l + 1]] + s[sta[l + 1]] - f[(i - 1) & 1][sta[l]] - s[sta[l]] <= a[j] * (sta[l + 1] - sta[l])) l++;
			f[i & 1][j] = f[(i - 1) & 1][sta[l]] + a[j] * (j - sta[l]) - (s[j] - s[sta[l]]);
			while (l < r && (f[(i - 1) & 1][j] + s[j] - f[(i - 1) & 1][sta[r]] - s[sta[r]]) * (sta[r] - sta[r - 1]) <= (f[(i - 1) & 1][sta[r]] + s[sta[r]] - f[(i - 1) & 1][sta[r - 1]] - s[sta[r - 1]]) * (j - sta[r])) r--;
			sta[++r] = j;
		}
	}
	
	printf("%lld", f[p & 1][m]);
	
	return 0;
}