poj3714 Raid

哈，就这样水过去了？？

不科学呀。

好吧，平面最近点对(伪)。

就是一个二分的事：

1.以x排序。

2.二分，递归求解。

具体来讲，主要难在如何求 横跨两个区域的最近点对 这一环节。

那么是这样做的：首先我们会有个ans是左右的最小值对吧。

我们只考虑在划分线左右两边不超过ans的那些点(即为[_l, _r])

然后把这些点暴力求解即可AC按照y排序，这样的情况下，再只考虑与一个点上下不超过ans的那些点即可......

那么下面是我的暴力版写法：

/// poj 3714
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long LL;
const int N = 100010;
const double INF = (double)(0x7f7f7f7f);

struct Point {
    LL x, y;
    bool f;
    bool operator < (const Point &k) const {
        return x < k.x;
    }
}p[N << 1];

inline double dis(int i, int j) {
    return sqrt((double)((p[i].x - p[j].x) * (p[i].x - p[j].x) + (p[i].y - p[j].y) * (p[i].y - p[j].y)));
}

double twsolve(int l, int r) {
    if(l == r) return INF;
    int mid = (l + r) >> 1;
    double k = min(twsolve(l, mid), twsolve(mid + 1, r));
    int _l, _r;
    for(_l = mid; _l >= l; _l--) {
        if(p[_l].x <= p[mid + 1].x - k) {
            break;
        }
    }
    for(_r = mid + 1; _r <= r; _r++){
        if(p[_r].x >= p[mid].x + k) {
            break;
        }
    }
    if(_l < l) _l = l;
    if(_r > r) _r = r;
    double ans = k;
    for(int i = _l; i <= mid; i++) {
        for(int j = mid + 1; j <= _r; j++) {
            if(!(p[i].f ^ p[j].f)) continue;
            ans = min(ans, dis(i, j));
        }
    }
    return ans;
}

int main() {
    int T, n;
    scanf("%d", &T);
    while(T--) {
        scanf("%d", &n);
        for(int i = 1; i <= n; i++) {
            scanf("%I64d %I64d", &p[i].x, &p[i].y);
            p[i].f = 0;

        }
        for(int i = 1; i <= n; i++) {
            scanf("%I64d %I64d", &p[i + n].x, &p[i + n].y);
            p[i + n].f = 1;
        }
        sort(p + 1, p + (n << 1) + 1);
        printf("%.3lf
", twsolve(1, n << 1));
        //cout << twsolve(1, n << 1) << endl;
    }
    return 0;
}

附赠数据生成器：

#include<cstdio>
#include<cstdlib>
#include<ctime>
using namespace std;
int main()
{
    freopen("in.in","w",stdout);
    srand(time(0));
    printf("1
");
    int n = rand() % 100000 + 1;
    printf("%d
", n);
    n = n << 1;
    for(int i = 1; i <= n; i++) {
        printf("%d %d
", rand() % 1000010 + 1, rand() % 1000010 + 1);
    }
    return 0;
}

[update]2018.10.20

这回写的是正规的nlogn，洛谷P1429 平面最近点对

y的排序使用归并，少一个排序的log。

然后发现一点小问题：我写的错误的程序A了，但是正确的却WA了......以后写了对的正解再发上来吧。