제출 #252755

#	제출 시각	아이디	문제	언어	결과	실행 시간	메모리
252755		islingr	Drzava (COCI15_drzava)	C++14	0 / 160	3 ms	1920 KiB

drzava

#include <bits/stdc++.h>

using namespace std;
using ll = long long;
using ld = long double;

#define rep(i, a, b) for (auto i = (a); i < (b); ++i)
#define trav(e, x) for (auto &e : x)
#define eb(x...) emplace_back(x)
#define all(x) begin(x), end(x)
#define sz(x) int((x).size())
#define lb(x...) lower_bound(x)
#define f first
#define s second

template<class T> bool ckmin(T& a, const T& b) { return a > b ? a = b, 1 : 0; }
template<class T> bool ckmax(T& a, const T& b) { return a < b ? a = b, 1 : 0; }

const int N = 1 << 16, K = 1 << 6;
const ll inf = 2e8;
using point = complex<ll>;

point p[N];
short v[N], k;

int n, nxt[N];
vector<int> cmp[N];
bool dp[K][K];

int head(int u) { return nxt[u] != u ? nxt[u] = head(nxt[u]) : u; }
void unite(int u, int v) { nxt[head(u)] = head(v); }

bool f(ll d) {
	ld r = sqrt(ld(d));
	set<pair<ll, int>> S;

	int l = 0;
	rep(u, 0, n) nxt[u] = u;
	rep(i, 0, n) {
		if (real(p[i]) - real(p[l]) > r)
			S.erase({imag(p[l]), l}), ++l;
		for (auto it = S.lb({ceil(imag(p[i]) - r), -inf});
			it != end(S) && it->f - r <= imag(p[i]); ++it)
			if (norm(p[it->s] - p[i]) <= d)
				unite(it->s, i);
		S.insert({imag(p[i]), i});
	}
	rep(u, 0, n) cmp[u].clear();
	rep(u, 0, n) cmp[head(u)].eb(u);
	rep(z, 0, n) {
		auto &a = cmp[z];
		if (a.empty()) continue;
		if (sz(a) > k) return true;
		int n = sz(a);
		rep(i, 0, n) rep(w, 0, k) dp[i][w] = 0;
		dp[0][v[a[0]]] = 1;
		rep(i, 1, n) {
			dp[i][v[a[i]]] = 1;
			rep(w, 0, k)
				if (dp[i - 1][w])
					dp[i][(w + v[a[i]]) % k] = dp[i][w] = true;
		}
		if (dp[n - 1][0]) return true;
	}
	return false;
}

pair<point, short> tmp[N];

signed main() {
	ios::sync_with_stdio(false);
	cin.tie(nullptr);

	cin >> n >> k;
	rep(i, 0, n) {
		int x, y, c; cin >> x >> y >> c;
		tmp[i] = {{x, y}, c % k};
	}
	sort(tmp, tmp + n, [&](auto a, auto b) {
		if (real(a.f) == real(b.f)) return imag(a.f) < imag(b.f);
		return real(a.f) < real(b.f);
	});
	rep(i, 0, n) tie(p[i], v[i]) = tmp[i];
	
	ll l = 0, r = inf * inf;
	while (r - l > 1) {
		ll m = l + (r - l + 1) / 2;
		if (f(m)) r = m;
		else l = m; 
	}
	cout << fixed << setprecision(4) << sqrt(ld(r));
}

• Subtask 1: 0 / 8
• Subtask 2: 0 / 8
• Subtask 3: 0 / 8
• Subtask 4: 0 / 8
• Subtask 5: 0 / 8
• Subtask 6: 0 / 8
• Subtask 7: 0 / 8
• Subtask 8: 0 / 8
• Subtask 9: 0 / 8
• Subtask 10: 0 / 8
• Subtask 11: 0 / 8
• Subtask 12: 0 / 8
• Subtask 13: 0 / 8
• Subtask 14: 0 / 8
• Subtask 15: 0 / 8
• Subtask 16: 0 / 8
• Subtask 17: 0 / 8
• Subtask 18: 0 / 8
• Subtask 19: 0 / 8
• Subtask 20: 0 / 8