【题解】垦田计划

题目来源： 第29次CCF计算机软件能力认证
题目链接： AcWing 5017. 垦田计划

一、题目

顿顿总共选中了 $n$ 块区域准备开垦田地，由于各块区域大小不一，开垦所需时间也不尽相同。
据估算，其中第 $i$ 块（$1\le i\le n$）区域的开垦耗时为 $t_i$ 天。
这 $n$ 块区域可以同时开垦，所以总耗时 $tTotal$ 取决于耗时最长的区域，即： t T o t a l = m a x { t 1 , t 2 , … , t n } tTotal=max\{t1,t2,…,tn\} tTotal=max{t1,t2,…,tn}为了加快开垦进度，顿顿准备在部分区域投入额外资源来缩短开垦时间。
具体来说：

• 在第 $i$ 块区域每投入 $c_i$ 单位资源，便可将其开垦耗时缩短 1 天；
• 耗时缩短天数以整数记，即第 $i$ 块区域投入资源数量必须是 $c_i$ 的整数倍；
• 在第 $i$ 块区域最多可投入 $c_i\times (t_i-k)$ 单位资源，将其开垦耗时缩短为 $k$ 天；
• 这里的 $k$ 表示开垦一块区域的最少天数，满足 0 < k ≤ m i n { t 1 , t 2 , … , t n } 0<k≤min\{t_1,t_2,…,t_n\} 0<k≤min{t1​,t2​,…,tn​}；
• 换言之，如果无限制地投入资源，所有区域都可以用 $k$ 天完成开垦。

现在顿顿手中共有 $m$ 单位资源可供使用，试计算开垦 $n$ 块区域最少需要多少天？

---

输入格式

输入共 $n+1$ 行。
输入的第一行包含空格分隔的三个正整数 $n$, $m$, $k$，分别表示待开垦的区域总数、顿顿手上的资源数量和每块区域的最少开垦天数。
接下来 $n$ 行，每行包含空格分隔的两个正整数 $t_i$ 和 $c_i$，分别表示第 $i$ 块区域开垦耗时和将耗时缩短 $1$ 天所需资源数量。

输出格式

输出一个整数，表示开垦 $n$ 块区域的最少耗时。

数据范围

$70\%$ 的测试数据满足：$0<n$, $t_i$, $c_i\le 100$ 且 $0<m\le 10^6$；
全部的测试数据满足：$0<n$, $t_i$, $c_i\le 10^5$ 且 $0<m\le 10^9$。

输入样例1：

4 9 2
6 1
5 1
6 2
7 1

输出样例1：

5

样例1解释
如下表所示，投入 $5$ 单位资源即可将总耗时缩短至 $5$ 天。此时顿顿手中还剩余 $4$ 单位资源，但无论如何安排，也无法使总耗时进一步缩短。

$i$	基础耗时	$t_i$ 缩减 $1$ 天所需资源	$c_i$ 投入资源数量	实际耗时
$1$	$6$	$1$	$1$	$5$
$2$	$5$	$1$	$0$	$5$
$3$	$6$	$2$	$2$	$5$
$4$	$7$	$1$	$2$	$5$

输入样例2：

4 30 2
6 1
5 1
6 2
7 1

输出样例2：

2

样例2解释
投入 $20$ 单位资源，恰好可将所有区域开垦耗时均缩短为 $k=2$ 天；受限于 $k$，剩余的 $10$ 单位资源无法使耗时进一步缩短。

二、思路和算法

由题意可知，开垦 $n$ 块地所需的最少天数 days为所有地耗时缩短后的最长时间，且满足 $k\le days\le 10^5$。

以所需的最少天数为分界，可以划分出两个不同的类别。当days小于最少天数时，其所需资源大于m；当days大于等于最少天数时，其所需资源小于等于m。所以，这个题中的数据具有二段性，可以通过二分的方法解决。

其中，每段缩短至x天所需的资源可以通过如下公式求得：
$$used=(t[i]-x)\times c[i]$$
由于题目中所给的参数 $n,t_i,c_i$ 的最大值都是 $10^5$，因此可能使用的最大资源是 $10^{15}$，超过了int的范围，需要使用long long存储所需的资源数量。

---

另一种解决方案是，遍历所有地的天数，找到最长的天数，使用已有的资源对其进行”裁剪“，类似于贪心的思想。

总体来看，可以计算出每天所需的资源数目，在缩短天数时，只需将减去要缩短的那天所需的资源即可。对于每天所需资源的计算，可以通过一个差分数组代表其每天所需的资源。对于差分数组s[i]，表示如果当前最少天数为i天，将其缩短为i-1天所需的资源数量。

题目要求输入两个正整数，分别代表开垦一块地的时间（设为 $t$）和缩短一天所需要的资源数目（设为 $c$）。那么，对于每次的输入，可以通过公式s[1]+c, s[t+1]-c初始化差分数组，然后求其前缀和（for循环中使用s[i] += s[i-1]）即可得到所需的差分数组。

最后，从可能的最大天数开始遍历差分数组，缩短逐个缩短天数即可。

三、代码及复杂度

方法一：二分 $O(nlo{g}_{2}n)$

#include <iostream>
using namespace std;

typedef long long LL;

const int N = 100010;
int n, m, k;
int t[N], c[N];

bool check(int mid)
{
    LL sum = 0;
    for (int i = 0; i < n; i++)
    	if (t[i] > mid)
            sum += (LL)(t[i] - mid) * c[i];
    return sum <= m;
}

int main()
{
    scanf("%d%d%d", &n, &m, &k);
    for (int i = 0; i < n; i++) scanf("%d%d", &t[i], &c[i]);
    int l = k, r = N;
    while (l < r)
    {
        int mid = l + r >> 1;
        if (check(mid)) r = mid;	// 如果mid可以满足，那么所有大于等于mid的都可以满足 答案->[l, mid]
        else l = mid + 1;	// 如果mid不能满足, 那么答案在区间[mid+1, r]中
    }
    printf("%d\n", l);
    return 0;
}

---

方法二：差分 $O(n)$

#include <iostream>
using namespace std;

typedef long long LL;

const int N = 100010;

int n, m, k;
LL s[N];

int main()
{
    scanf("%d%d%d", &n, &m, &k);
    for (int i = 0; i < n; i++)
    {
        int t, c;
        scanf("%d%d", &t, &c);
        s[1] += c, s[t + 1] -= c;
    }
    for (int i = 0; i < N; i++) s[i] += s[i - 1];
    LL cost = 0;
    int i = N - 1;
    while (i > k && cost + s[i] <= m) cost += s[i--];
    printf("%d\n", i);
    return 0;
}

---

复杂度分析

对于二分的方法，包括二分和枚举两部分。其中，二分答案部分的时间复杂度为 $O(lo{g}_{2}n)$。在每次二分结束后，都需要对二分的答案进行枚举计算，判断是否符合题意，该部分时间复杂度为 $O(n)$。故二分方法的时间复杂度为 $O(nlo{g}_{2}n)$。

对于差分的方法，包括初始化和遍历两部分，且这两部分是并列关系。在初始化时，需要输入数据和计算前缀和，其时间复杂度为 $O(n)$；在遍历阶段，从大到小一次遍历差分数组，其时间复杂度为 $O(n)$。综上，差分方法的时间复杂度为 $O(n)$。