UVa1467/LA4850 Installations

题目链接

本题是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;
}