http://poj.org/problem?id=2932
挑战程序设计竞赛P258
#include <vector>
#include <cstdio>
#include <utility> //pair
#include <iostream>
#include <algorithm>
#include <set>
#define MyPair pair<double, int>
using namespace std;
const int MAXN = 40000 + 9;
struct O
{
double x, y, r;
void read() {
scanf("%lf%lf%lf", &r, &x, &y);
}
};
O o[MAXN];
bool inside(O o1, O o2) {
double dx = o1.x - o2.x, dy = o1.y - o2.y;
return dx * dx + dy * dy <= o2.r * o2.r;
}
int main() {
//freopen("in.txt", "r", stdin);
int n;
scanf("%d", &n);
for (int i = 0; i < n; i++) {
o[i].read();
}
vector<MyPair> events;
for (int i = 0; i < n; i++) {
events.push_back(make_pair(o[i].x - o[i].r, i));
events.push_back(make_pair(o[i].x + o[i].r, i + n));
}
sort(events.begin(), events.end());
set<MyPair> outers;
vector<int> res;
for (int i = 0; i < events.size(); i++) {
int id = events[i].second % n;
if (events[i].second < n) {
set<MyPair>::iterator it = outers.lower_bound(make_pair(o[id].y, id));
if (it != outers.end() && inside(o[id], o[it->second])) continue;
if (it != outers.begin() && inside(o[id], o[(--it)->second])) continue;
res.push_back(id);
outers.insert(make_pair(o[id].y, id));
} else {
outers.erase(make_pair(o[id].y, id));
}
}
sort(res.begin(), res.end());
printf("%d\n", res.size());
printf("%d", res[0] + 1);
for (int i = 1; i < res.size(); i++) {
printf(" %d", res[i] + 1);
}
printf("\n");
}