洛谷P1158 导弹拦截[NOIP2010 普及组]

题面：

输入#1（输出18）：
0 0 10 0
2
-3 3
10 0

输入#2（输出30）
0 0 6 0
5
-4 -2
-2 3
4 0
6 -2
9 1

解题：

分析：

数据范围：N≤1e5，时间复杂度大约在O（nlogn），

sort排序可通过，而双层1~N循环枚举每种可能解不太有希望通过……

预处理（复杂度O（n））：

创建D类，代表导弹；拥有变量s1、s2，分别对应与系统1、系统2距离²；

读入时，求出每个导弹的s1、s2；

排序（复杂度O（nlogn））

用sort（）函数，对s1距离降序排序。

计算（复杂度O（n））

首先，我们特判，系统1、系统2覆盖了全部点，即：将ans初始化为s1、s2的最大值；

下面打出排序好的表：

我们知道，如果系统1覆盖了距离远的点，那么距离近的点也一并被系统1覆盖，

即：如果取40，后面20、16、13的点都被系统1覆盖，系统2不需要计算这些点

则：系统2需要的最大工作半径为：当前未被覆盖的点中，s2的最大值！

有点小困惑？没关系，下面我们一起来对样例#2进行模拟

（系统1范围：k1，系统2：k2）：

• 1：k1=40，k2=10，

系统1覆盖了距离40的C点，以及比它离自己更近的、后面的（D E F）点

而系统2需要做的，就是把系统1没有覆盖到的B点覆盖，即k2=10；

• 2：k1=20，k2=10，

系统1覆盖了距离20的D点，以及比它离自己更近的、后面的（E F）点

而系统2需要做的，就是把系统1没有覆盖到的B、C点覆盖，因此，k2必须取max（4 ，10）=10；

• 3：k1=16，k2=104，

• 4：k1=13，k2=104，

如此，最优答案ans除了s1、s2的最大值（一个系统覆盖所有点，另一个系统不工作）：82、104外，

便只能在我们所枚举的k1+k2中产生，在计算k1、k2时不断更新ans即可；

❀❀❀❀❀❀❀完结撒花❀❀❀❀❀❀❀

AC代码奉上

#include<iostream>
#include<algorithm>
#include<vector>
#define MAXN 1e6+5

using namespace std;

int xa, ya, xb, yb, n;
long long ans = MAXN, s2Maxn = 0;

class D
{
public:
    int x; int y;
    long long s1, s2;//与1/2发射器的距离

    void cal() //计算与系统1、系统2的距离
    {
        this->s1 = (x - xa) * (x - xa) + (y - ya) * (y - ya);
        this->s2 = (x - xb) * (x - xb) + (y - yb) * (y - yb);
    }
};

class MyCompare  //给s1排序的仿函数
{
public:
    bool operator ()(D d1,D d2)
    {
        return d1.s1 > d2.s1;
    }
};

vector<D>v; //装"导弹"的容器（

int main()
{
    cin >> xa >> ya >> xb >> yb >> n;
    for (int i = 1; i <= n; i++)
    {
        D d;
        cin >> d.x >> d.y;
        d.cal();
        s2Maxn = max(s2Maxn, d.s2); //记一下s2的最大值，因为没有对s2排序
        v.push_back(d);
    }

    if (v.size() == 1) { cout << min(v[0].s1, v[0].s2) << endl; return 0; } //对只有一个点的情况特判

    sort(v.begin(), v.end(), MyCompare());  //排序排序
    ans = min(min(v[0].s1, s2Maxn), v[0].s2 + v[1].s1); //初始化ans，为一个系统全覆盖，另一个不工作情况
    for (int i = 1; i <= v.size() - 2; i++) 
    {
        v[i].s2 = max(v[i].s2, v[i - 1].s2);   //前面的s2更大，用大的s2把这个替换掉，因为我们只要最大的s2 
        ans = min(ans, v[i + 1].s1 + v[i].s2); //更新答案
    }
    cout << ans << endl;
    return 0;  //溜咯~
}