Submission #155223

# Submission time Handle Problem Language Result Execution time Memory
155223 2019-09-27T07:03:58 Z KCSC Drzava (COCI15_drzava) C++14
48 / 160
1000 ms 31776 KB
#include <bits/stdc++.h>
using namespace std;

const int DIM = 50005;
const int MOD = 666013;
const int CNS = 1000000007;
const double EPS = 0.000001;

int fth[DIM], msk[DIM];

vector<int> hsh[MOD];
map<pair<long long, long long>, int> mmp;

struct Point {
    int x, y, c;
} pts[DIM];

int _(pair<long long, long long> pr) {
    return (1LL * pr.first % MOD * CNS + pr.second) % MOD;
}

inline int getRoot(int x) {
    while (fth[x] > 0)
        x = fth[x];
    return x;
}

inline int sign(double x, double y) {
    if (fabs(x - y) < EPS)
        return 0;
    return (x < y ? -1 : 1);
}

int _k;
bool unite(int x, int y, double v) {
    if (sign(hypot(abs(pts[x].x - pts[y].x),
                   abs(pts[x].y - pts[y].y)), v * 2) > 0)
        return false;
    x = getRoot(x); y = getRoot(y);
    if (x == y)
        return false;
    if (fth[x] > fth[y])
        swap(x, y);
    fth[x] += fth[y]; fth[y] = x;
    for (int i = 0; i < _k; ++i) {
        if ((msk[y] >> i) & 1) {
            msk[x] |= ((1LL * msk[x]) << i) | (((1LL * msk[x]) << i) >> _k);
            msk[x] &= (1 << _k) - 1;
        }
    }
    if (msk[x] & 1)
        return true;
    return false;
}

bool solve(int n, int k, double v) {
    mmp.clear(); v /= 2;
    for (int i = 1; i <= n; ++i) {
        pair<long long, long long> sq = make_pair(pts[i].x / v, pts[i].y / v);
        if (++mmp[sq] == k)
            return true;
    }
    for (int i = 0; i < MOD; ++i)
        hsh[i].clear();
    for (int i = 1; i <= n; ++i) {
        fth[i] = -1;
        msk[i] = (1 << pts[i].c);
    }
    for (int i = 1; i <= n; ++i) {
        pair<long long, long long> sq = make_pair(pts[i].x / v, pts[i].y / v);
        for (int d1 = -3; d1 <= 0; ++d1) {
            for (int d2 = -3; d2 <= 3; ++d2) {
                pair<long long, long long> nsq = make_pair(max(0LL, sq.first + d1),
                                                           max(0LL, sq.second + d2));
                for (int j : hsh[_(nsq)])
                    if (unite(i, j, v))
                        return true;
            }
        }
        hsh[_(sq)].push_back(i);
    }
    return false;
}

inline bool cmp(Point p1, Point p2) {
    if (p1.x != p2.x)
        return p1.x < p2.x;
    return p1.y < p2.y;
}

int main(void) {
   // freopen("drzava.in", "r", stdin);
   // freopen("drzava.out", "w", stdout);
    int n, k;
    cin >> n >> k; _k = k;
    for (int i = 1; i <= n; ++i) {
        cin >> pts[i].x >> pts[i].y >> pts[i].c;
        pts[i].c %= k;
    }
    sort(pts + 1, pts + n + 1, cmp);
    double le = 0, ri = CNS * 2;
    for (int i = 1; i <= 60; ++i) {
        double md = (le + ri) / 2;
        if (solve(n, k, md))
            ri = md;
        else
            le = md;
    }
    cout << setprecision(3) << fixed << le;
    return 0;
}
# Verdict Execution time Memory Grader output
1 Correct 99 ms 15996 KB Output is correct
2 Correct 64 ms 16120 KB Output is correct
# Verdict Execution time Memory Grader output
1 Incorrect 85 ms 16124 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 88 ms 16120 KB Output is correct
2 Correct 88 ms 16300 KB Output is correct
# Verdict Execution time Memory Grader output
1 Incorrect 113 ms 16120 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 96 ms 16120 KB Output is correct
2 Correct 68 ms 16096 KB Output is correct
# Verdict Execution time Memory Grader output
1 Incorrect 81 ms 15992 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 77 ms 16024 KB Output is correct
2 Correct 72 ms 16248 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 91 ms 15992 KB Output is correct
2 Correct 82 ms 16220 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 99 ms 16012 KB Output is correct
2 Correct 532 ms 25480 KB Output is correct
# Verdict Execution time Memory Grader output
1 Incorrect 102 ms 16120 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 91 ms 16020 KB Output is correct
2 Execution timed out 1042 ms 31664 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Incorrect 83 ms 15992 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 78 ms 15992 KB Output is correct
2 Execution timed out 1037 ms 30884 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 79 ms 15992 KB Output is correct
2 Execution timed out 1044 ms 31264 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 100 ms 16120 KB Output is correct
2 Execution timed out 1039 ms 31184 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 93 ms 16012 KB Output is correct
2 Execution timed out 1045 ms 30556 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 100 ms 15992 KB Output is correct
2 Execution timed out 1062 ms 31776 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 80 ms 15992 KB Output is correct
2 Execution timed out 1055 ms 30060 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 89 ms 16120 KB Output is correct
2 Execution timed out 1055 ms 28732 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Incorrect 81 ms 15992 KB Output isn't correct
2 Halted 0 ms 0 KB -