/**
* author: wxhtzdy
* created: 27.06.2022 18:47:13
**/
#include <bits/stdc++.h>
using namespace std;
struct point {
int x, y, c;
bool operator < (const point &p) const {
return x < p.x || (x == p.x && y < p.y);
};
};
long long Dist(point a, point b) {
long long dx = abs(a.x - b.x);
long long dy = abs(a.y - b.y);
return dx * dx + dy * dy;
}
int readint(){
int x=0,f=1; char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
point a[50000];
int main() {
ios::sync_with_stdio(false);
cin.tie(0);
int n, k;
n = readint();
k = readint();
for (int i = 0; i < n; i++) {
a[i].x = readint();
a[i].y = readint();
a[i].c = readint();
a[i].c %= k;
}
sort(a, a + n);
auto Check = [&](long long x) {
int d = (double) (sqrt(x) + 1);
vector<vector<int>> g(n);
auto Add = [&](int x, int y) {
g[x].push_back(y);
g[y].push_back(x);
};
set<pair<int, int>> s;
int ptr = 0;
for (int i = 0; i < n; i++) {
while (ptr < i && a[ptr].x + d < a[i].x) {
s.erase(s.find(make_pair(a[ptr].y, ptr)));
ptr += 1;
}
auto it = s.lower_bound(make_pair(a[i].y - d, -1));
int cnt = 0;
while (it != s.end()) {
int j = it->second;
if (abs(a[i].y - a[j].y) > d) {
break;
}
cnt += 1;
if (cnt >= 8 * k) {
return true;
}
if (Dist(a[i], a[j]) <= x) {
Add(i, j);
}
it++;
}
s.emplace(a[i].y, i);
}
vector<bool> was(n);
for (int i = 0; i < n; i++) {
if (!was[i]) {
vector<int> que(1, i);
was[i] = true;
for (int b = 0; b < (int) que.size(); b++) {
if ((int) que.size() >= k) {
return true;
}
int j = que[b];
for (int p : g[j]) {
if (!was[p]) {
was[p] = true;
que.push_back(p);
}
}
}
vector<bool> dp(k + 1, false);
dp[0] = true;
for (int j : que) {
vector<bool> new_dp(k + 1);
for (int p = k; p >= 0; p--) {
new_dp[p] = (dp[p] | dp[(((p - a[j].c) % k) + k) % k]);
}
swap(new_dp, dp);
if (dp[k]) {
return true;
}
}
}
}
return false;
};
long long low = 0, high = 1e18, ans = 0;
while (low <= high) {
long long mid = low + high >> 1;
if (Check(mid)) {
ans = mid;
high = mid - 1;
} else {
low = mid + 1;
}
}
cout << fixed << setprecision(3) << sqrt(ans * 1.0) << '\n';
return 0;
}
Compilation message
drzava.cpp: In function 'int main()':
drzava.cpp:112:25: warning: suggest parentheses around '+' inside '>>' [-Wparentheses]
112 | long long mid = low + high >> 1;
| ~~~~^~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
212 KB |
Output is correct |
2 |
Correct |
1 ms |
340 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
212 KB |
Output is correct |
2 |
Correct |
1 ms |
212 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
212 KB |
Output is correct |
2 |
Correct |
22 ms |
980 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
212 KB |
Output is correct |
2 |
Correct |
11 ms |
584 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
212 KB |
Output is correct |
2 |
Correct |
33 ms |
724 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
35 ms |
1280 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
21 ms |
968 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
212 KB |
Output is correct |
2 |
Correct |
15 ms |
724 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
143 ms |
4304 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
2 ms |
340 KB |
Output is correct |
2 |
Correct |
980 ms |
51336 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
212 KB |
Output is correct |
2 |
Execution timed out |
1000 ms |
50468 KB |
Time limit exceeded |
3 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
655 ms |
26104 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
212 KB |
Output is correct |
2 |
Correct |
676 ms |
29916 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
950 ms |
54188 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
350 ms |
6744 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
212 KB |
Output is correct |
2 |
Correct |
507 ms |
18396 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
212 KB |
Output is correct |
2 |
Correct |
457 ms |
21504 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
212 KB |
Output is correct |
2 |
Correct |
510 ms |
21668 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
410 ms |
12208 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
212 KB |
Output is correct |
2 |
Correct |
352 ms |
8420 KB |
Output is correct |