题目链接
本题是ICPC亚洲区域赛Daejeon赛区2010年的 H 题。
分析
先来证明一个定理:
按照 di 进行排序所得到的方案,一定使得罚时最大服务项的罚时最小。
证明:
对于任意两项服务 Ja、Jb,设 da < db,当 Ja 在前时,罚时为 t1 = max(C+sa-da, C+T+sa+sb-db),C、T 均为非负常数; 当 Jb 在前时,罚时为 t2 = max(C+sb-db, C+T+sb+sa-da);
而 C+T+sb+sa-da > C+sb-db,故 t2 = C+T+sb+sa-da
又 C+T+sb+sa-da > C+T+sa+sb-db 且 C+T+sb+sa-da > C+sa-da,故 t2 > t1
既然对于任意的两项服务,di 比较小的都应该在前面,那么对于整个序列, 就应该是按 di 递增排序时,序列的罚时最大项的罚时最小。
然后看这道题
先按照 di 为第一关键字,si 为第二关键字排序,先扫描一遍整个序列,记下最大罚时的位置和次大罚时的位置,序号较小的为 p1,序号较大的为 p2,这样序 列被分成了 A,B,C 三个区间:[0,p1]、(p1,p2)、[p2,n)。
根据上述定理,A,B,C 三个区间内最大罚时已经是最小了,如果随便更改一 个区间内部的顺序,得到新序列的最大罚时不会更小,如果这个数值大于了所记 录的 p1,p2 处的罚时,比原先的情况更劣的一组解,所以说,在 A,B,C 三个区间 内部打乱顺序无助于寻找最优解。
把[0,p2),(p1,n)分别记为 D,E,如果这两个区间的顺序被扰乱,同样不可能 构造出更优的解。
因此只有一种方法可以构造出更优的解,就是改变 A∪C 中服务项的顺序, 并且一定不能只扰乱一个区间的顺序。
如果把一个 di 比较大的从 C 扔到 A 里去(p2 除外),那么对于此项来讲,它 的罚时变小了,而对于 p1、p2 来讲,罚时都变大了,这种解是更劣的,那如果把 多个 C 里的服务项扔到 A 里去,显然答案会变得更劣,因此 C 中的服务项没必要 放到A里去。
那么 A 中的项可以放到 C 里(p2 之后)吗?
如果放一个过去的话,显然它本身的实际结束时间变大了,而且它的 di 不 大于 dp2,因此它的罚时比原先的最大值还要大,而原先的最大值变小了,甚至可 能比它前面的某些项罚时还要小,因此最大值和次大值的和可能变小,可以更新 一下答案。
不用考虑继续从 A 往 C 扔更多的服务项,因为从区间 A([0,p1])中依次枚举 一项 Ji 放到 C 中作为最大罚时的服务项时,其余[0,n)-Ji 的项仍然是按照 di 为第一 关键字,si 为第二关键字排序的,根据定理可知[0,n)-Ji 的最大罚时(也就是第二大罚时)一定最小。 经过以上论述,排序后唯一能够构造出更优解的方法,就是从[0,p1]中选出一项,放到 p2 之后。将[0,p1]中的每项都尝试一下,最终一定能够得到最优解。
AC代码
#include <iostream>
#include <algorithm>
using namespace std;
#define N 520
struct node {
int s, d;
bool operator< (const node& rhs) const {
return d<rhs.d || (d==rhs.d && s<rhs.s);
}
} b[N];
int solve() {
short n; cin >> n;
for (short i=0; i<n; ++i) cin >> b[i].s >> b[i].d;
sort(b, b+n);
int s = 0, px = 0, py = 0; short x, y;
for (short i=0; i<n; ++i) {
s += b[i].s;
int v = s - b[i].d;
if (v > 0) {
if (px < 1) {
px = v; x = i;
} else if (v >= px) {
if (px > py) y = x;
py = px; px = v; x = i;
} else if (v >= py) {
py = v; y = i;
}
}
}
int ans = px + py;
if (py > 0) {
int pr = px;
if (y > x) {
int t = y; y = x; x = t; pr = py;
}
for (int i=n-1, pp=0; i>=x; s -= b[i--].s) {
py = pp;
for (int j=0, ss = 0; j<y; ++j) {
ans = min(ans, s - b[j].d + max(py, px - b[j].s));
py = max(py, (ss += b[j].s) - b[j].d);
}
ans = min(ans, s - b[y].d + max(py, pr - b[y].s));
pp = max(pp, s - b[i].d);
}
}
return ans;
}
int main() {
int t; cin >> t;
while (t--) cout << solve() << endl;
return 0;
}