jzoj 2866. 【集训队互测 2012】Bomb

Description

给你\(n\)个点，坐标分别为\((xi,yi)\)。从中取出三个点，使得其两两间曼哈顿距离和最大和最小，求最大值和最小值。
对于 100% 的数据， N<=100000 ， 0<=Xi,Yi<=10^8

Solution

看完题目后感觉要分类讨论，思考1h后果断暴力O(n^3)。
但是判断了一下n<=500才跑暴力，得了30分。（™10^9过4s？不判50分）
其实可以O(n^2)暴力的。（听说O(n^2)+优化 = 100）
对于第一个最大值，我们可以\(O(n)\)求出。
我们先扫一遍\(n\)个点，可以得到最大和最小的\(x\)以及最大和最小的\(y\)。
然后我们再枚举每个点，而其所组成的三个点的最大值便是\(2 * max(maxx - a[i].x, a[i].x - minx) + 2 * max(maxy - a[i].y, a[i].y - miny)\)
因为它所组成的矩形是最长的长和最长的宽，所以值便是最大的。
这是运用了贪心的思想，可以感性理解一下。
然后我们再求最小值。
这里使用\(O(n^2)\)的方法+优化的。
我们先将点按照\(x\)从小到大排序。
对于三个点：\((x1, y1),(x2, y2),(x3, y3)\)。
\(x1 < x2 < x3\)，所以关于\(x\)的答案就应该是：
\((x2 - x1) + (x3 - x2) + (x3 - x1)\)
\(=2 * (x3 - x1)\)
我们发现与\(x2\)无关。
所以我们只需枚举两个点，即可得到有关\(x\)的答案。
而\(y\)呢，我们便只需在线维护一个\(i\) ~ \(j\)的与\(a[i].y,a[j].y\)的差的最小值即可。
\(if (a[j].x - a[i].x >= ans2) break;\)
这个优化正确性显然。

Code

#include <cstdio>
#include <algorithm>
#define N 100010
#define fo(x, a, b) for (int x = a; x <= b; x++)
#define fd(x, a, b) for (int x = a; x >= b; x--)
using namespace std;
struct node{int x, y;}a[N];
int n, m, maxy, miny, up, down, s, s1;
int ans1, ans2 = 2e9;
bool pd;

inline int read()
{
    int x = 0; char c = getchar();
    while (c < '0' || c > '9') c = getchar();
    while (c >= '0' && c <= '9') x = (x << 1) + (x << 3) + (c ^ 48), c = getchar();
    return x;
}

inline int cmp(node x, node y) {return x.x < y.x;}

inline int max(int x, int y) {return x > y ? x : y;}

inline int min(int x, int y) {return x < y ? x : y;}

int main()
{
    freopen("bomb.in", "r", stdin);
    freopen("bomb.out", "w", stdout); 
    n = read();
    a[1].x = read();
    maxy = miny = a[1].y = read();
    fo(i, 2, n)
    {
        a[i].x = read(), a[i].y = read();
        if (a[i].y > maxy) maxy = a[i].y;
        else if (a[i].y < miny) miny = a[i].y;
    }
    sort(a + 1, a + n + 1, cmp);
    ans1 = a[n].x - a[1].x + max(max(maxy - a[1].y, a[1].y - miny), max(maxy - a[n].y, a[n].y - miny));
    fo(i, 2, n - 1) ans1 = max(ans1, max(a[n].x - a[i].x, a[i].x - a[1].x) + max(maxy - a[i].y, a[i].y - miny));
    fo(i, 1, n)
    {
        if (n >= 10000 && pd && a[i].x == a[i - 1].x) continue;
        if (a[i + 1].y < a[i].y) up = 2e9, down = a[i + 1].y;
        else up = a[i + 1].y, down = -2e9;
        pd = 0;
        fo(j, i + 2, n)
        {
            if (a[j].x - a[i].x >= ans2) {pd = 1; break;}
            s = max(up - a[i].y, abs(up - a[j].y));
            s1 = max(a[i].y - down, abs(down - a[j].y));
            s = max(min(s, s1), abs(a[i].y - a[j].y));
            ans2 = min(ans2, a[j].x - a[i].x + s);
            if (a[j].y < a[i].y) down = max(down, a[j].y);
            else up = min(up, a[j].y);
        }
    }
    printf("%d\n%d\n", ans1 << 1, ans2 << 1);
    return 0;
}