Submission #1160119

#TimeUsernameProblemLanguageResultExecution timeMemory
1160119xnqsCircle selection (APIO18_circle_selection)C++20
0 / 100
3094 ms32836 KiB
#pragma GCC optimize("Ofast")

#include <iostream>
#include <fstream>
#include <vector>
#include <queue>
#include <utility>
#include <algorithm>
#include <cmath>
#include <unordered_map>
#include <unordered_set>
#include <set>
#include <map>

struct Circle {
	int64_t x, y, r, idx;
	Circle():
		x(0), y(0), r(0), idx(0) {}
	Circle(int64_t x, int64_t y, int64_t r, int64_t idx):
		x(x), y(y), r(r), idx(idx) {}

	bool operator == (const Circle& other) const {
		if (this->x!=other.x) return 0;
		if (this->y!=other.y) return 0;
		if (this->r!=other.r) return 0;
		if (this->idx!=other.idx) return 0;
		return 1;
	}
};

struct PairHasher {
	std::size_t operator() (const std::pair<int64_t,int64_t>& p) const {
		return p.first * 1000000006 + p.second;
	}
};

int64_t n;
std::unordered_set<int64_t> active;
std::vector<Circle> circles;
std::vector<int64_t> eliminator;

int64_t grid_size;
std::unordered_map<std::pair<int64_t,int64_t>, std::vector<int64_t>, PairHasher> grid;

inline int64_t euclid_dist(const Circle& a, const Circle& b) {
	return 1LL*(a.x-b.x)*(a.x-b.x) + 1LL*(a.y-b.y)*(a.y-b.y);
}

/*
	d(a,b) - a.r <= b.r
	sqrt((a.x-b.x)^2 + (a.y-b.y)^2) - a.r <= b.r
	sqrt((a.x-b.x)^2 + (a.y-b.y)^2) <= b.r + a.r
	(a.x-b.x)^2 + (a.y-b.y)^2 <= (b.r + a.r)^2
*/
inline bool intersect(const Circle& a, const Circle& b) {
	return euclid_dist(a,b) <= 1LL*(a.r+b.r)*(a.r+b.r);
}

void init_grid(int64_t size) {
	grid.clear();
	grid_size = size;
	for (auto i : active) {
		if (eliminator[circles[i].idx]) {
			continue;
		}

		grid[std::pair<int64_t,int64_t>(circles[i].x/size,circles[i].y/size)].emplace_back(i);
	}
}

int main() {
	std::ios_base::sync_with_stdio(false);
	std::cin.tie(NULL);
	std::cout.tie(NULL);

	std::cin >> n;
	for (int64_t i = 0, x, y, r; i < n; i++) {
		active.emplace(i);
		std::cin >> x >> y >> r;
		circles.emplace_back(x,y,r,i);
	}

	std::sort(circles.begin(),circles.end(),[](const Circle& a, const Circle& b) {
		return a.r > b.r;
	});

	eliminator.assign(n, 0);
	init_grid(4LL*circles[0].r*circles[0].r);
	for (const auto& c : circles) {
		auto [x, y, r, idx] = c;
		if (4LL*r*r < grid_size/2) {
			init_grid(4*r*r);
		}
		if (eliminator[idx]) {
			active.erase(idx);
			continue;
		}

		eliminator[idx] = idx+1;
		active.erase(idx);

		//std::cout << idx << " " << grid_size << "\n";

		int64_t grid_x = x / grid_size;
		int64_t grid_y = y / grid_size;

		for (int64_t gx = grid_x-2; gx <= grid_x+2; gx++) {
			for (int64_t gy = grid_y-2; gy <= grid_y+2; gy++) {
				auto& tmp = grid[std::pair<int64_t,int64_t>(gx,gy)];
				for (int64_t i = 0; i < tmp.size(); i++) {
					if (eliminator[circles[tmp[i]].idx]) {
						continue;
					}

					if (intersect(c, circles[tmp[i]])) {
						eliminator[circles[tmp[i]].idx] = idx+1;
						//active.erase(circles[tmp[i]].idx);
					}
				}
			}
		}
	}

	for (int64_t i = 0; i < n; i++) {
		std::cout << eliminator[i] << " ";
	}
	std::cout << "\n";
}
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...