Submission #252756

# Submission time Handle Problem Language Result Execution time Memory
252756 2020-07-26T08:24:33 Z islingr Drzava (COCI15_drzava) C++14
64 / 160
1000 ms 5796 KB
#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(3) << sqrt(ld(r));
}
# Verdict Execution time Memory Grader output
1 Correct 2 ms 1920 KB Output is correct
2 Correct 2 ms 1920 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3 ms 1920 KB Output is correct
2 Correct 2 ms 1920 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 1920 KB Output is correct
2 Correct 379 ms 2064 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 2 ms 1920 KB Output is correct
2 Correct 305 ms 2184 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3 ms 1920 KB Output is correct
2 Correct 146 ms 2040 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 2 ms 1920 KB Output is correct
2 Correct 348 ms 2048 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 4 ms 1920 KB Output is correct
2 Correct 459 ms 2048 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 2 ms 1920 KB Output is correct
2 Correct 84 ms 2036 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 2 ms 1920 KB Output is correct
2 Execution timed out 1083 ms 4156 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 2 ms 1920 KB Output is correct
2 Execution timed out 1080 ms 5752 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 2 ms 1920 KB Output is correct
2 Execution timed out 1078 ms 5600 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 2 ms 1920 KB Output is correct
2 Execution timed out 1027 ms 5472 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 2 ms 1920 KB Output is correct
2 Execution timed out 1096 ms 5796 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 4 ms 1920 KB Output is correct
2 Execution timed out 1093 ms 5712 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 2 ms 1920 KB Output is correct
2 Execution timed out 1097 ms 5572 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 2 ms 1920 KB Output is correct
2 Execution timed out 1096 ms 5624 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 2 ms 1920 KB Output is correct
2 Execution timed out 1086 ms 5444 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 2 ms 1920 KB Output is correct
2 Execution timed out 1078 ms 5564 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 1 ms 1920 KB Output is correct
2 Execution timed out 1090 ms 5548 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 2 ms 1920 KB Output is correct
2 Execution timed out 1095 ms 5624 KB Time limit exceeded