POJ 2932 Coneology——扫描线

最新推荐文章于 2020-10-02 19:30:05 发布

hao_zong_yin

最新推荐文章于 2020-10-02 19:30:05 发布

阅读量282

点赞数 1

分类专栏：几何数据结构

本文链接：https://blog.csdn.net/hao_zong_yin/article/details/79340032

版权

数据结构同时被 2 个专栏收录

159 篇文章 0 订阅

订阅专栏

几何

15 篇文章 0 订阅

订阅专栏

乍一看还以为线段树+扫描线，仔细想想虽然不是一个东西但是相似性很高，只要维护一下与扫描线相交的最外层圆就可以了，这里用了set，第一次写写得则乱，后来看了挑战程序设计才简化了代码，骚操作学到了

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <vector>
#include <set>
using namespace std;
const int maxn = 40005;
int n;
double x[maxn], y[maxn], r[maxn];

bool inside(int i, int j) {
    double dx = x[i]-x[j], dy = y[i]-y[j];
    return dx*dx+dy*dy <= r[j]*r[j];
}

void solve() {
    vector<pair<double, int> > pos;
    for (int i = 0; i < n; i++) {
        pos.push_back(make_pair(x[i]-r[i], i));
        pos.push_back(make_pair(x[i]+r[i], i+n));
    }
    sort(pos.begin(), pos.end());
    set<pair<double, int> > s;
    vector<int> ans;
    for (int i = 0; i < pos.size(); i++) {
        int id = pos[i].second % n;
        if (pos[i].second < n) {
            set<pair<double, int> >::iterator it = s.lower_bound(make_pair(y[id], id));
            if (it != s.end() && inside(id, it->second)) continue;
            if (it != s.begin() && inside(id, (--it)->second)) continue;
            ans.push_back(id);
            s.insert(make_pair(y[id], id));
        }
        else {
            s.erase(make_pair(y[id], id));
        }
    }
    sort(ans.begin(), ans.end());
    printf("%d\n", ans.size());
    for (int i = 0; i < ans.size(); i++) {
        printf("%d%c", ans[i]+1, (i+1) == ans.size()?'\n':' ');
    }
}

int main() {
    while (~scanf("%d", &n)) {
        for (int i = 0; i < n; i++) {
            scanf("%lf %lf %lf", &r[i], &x[i], &y[i]);
        }
        solve();
    }
    return 0;
}