UESTC 1170 红与蓝 计算几何、贪心、红蓝点对

D - EN TARO Artanis

Time Limit:1000MS     Memory Limit:65535KB     64bit IO Format:%lld
& %llu

Submit Status Practice UESTC
1170

Description

平面上有N个红点和N个蓝点，求红点到蓝点的最近距离

Input

第一行为一个整数N 

接下来第N行每行两个整数xi,yi，表示第i个红点的坐标

接下来第N行每行两个整数xi,yi，表示第i个蓝点的坐标(1
                                  )

Output

红点到蓝点的最短距离

Sample Input

4

0 0

0 1

1 0

1 1

2 2

2 3

3 2

3 3

Sample Output

1.414

Hint

Source

UESTC 2016 Summer Training #16 Div.2

UESTC 1170

Source

贪心的做法

预处理所有Ai到O的距离， 然后根据距离排序，

之后依次对每个Bi也求出BiO， 然后二分查找， 然后处理最近的+-100个点， 答案必定在这100或者1000个点之内

比赛的时候是+-100 共200个点Accepted的

然后比赛结束后试了一下， 我这样的方法至少要+-40 共80个点才能通过那个题目的数据测试

好像还有专业计算几何的算法 红蓝点对 ??

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
typedef long long LL;
const int maxn = 1e5 + 8;

struct p{
    int x, y;
    double dist;
} val[maxn];

inline bool cmp(const p& a, const p& b)
{
    return a.dist < b.dist;
}

int main()
{
    #ifdef LOCAL
    freopen("a.txt", "r", stdin);
    //freopen("b.txt", "w", stdout);
    int T = 2;
    while(T--){
    #endif // LOCAL

    int n, x, y, a, b;
    double dist, ans = maxn;
    scanf("%d", &n);
    for(int i = 0; i < n; i++){
        scanf("%d%d", &val[i].x, &val[i].y);
        val[i].dist = sqrt(val[i].x*val[i].x + val[i].y*val[i].y);
    }
    sort(val, val + n, cmp);
    for(int i = 0; i < n; i++){
        scanf("%d%d", &x, &y);
        dist = sqrt(x*x + y*y);
        int m;
        a = 0, b = n - 1;
        while(a < b){
            m = a + (b - a)/2;
            if(val[m].dist >= dist) b = m;
            else a  = m + 1;
        }
        for(int j = max(0, a - 100); j < min(n, a + 100); j++){
            ans = min(ans, sqrt((x - val[j].x)*(x - val[j].x) + (y - val[j].y)*(y - val[j].y)));
        }
    }
    printf("%.3f", ans);

    #ifdef LOCAL
    printf("\n");
    }
    #endif // LOCAL
    return 0;
}

Thank you!

------from ProLights