oj.uz

Submission #970387

# Submission time Handle Problem Language Result Execution time Memory
970387 2024-04-26T13:14:54 Z steveonalex Circle selection (APIO18_circle_selection) C++17
52 / 100
 2953 ms 55128 KB 
#include <bits/stdc++.h>
 
using namespace std;
 
typedef long long ll;
typedef unsigned long long ull;
 
#define ALL(v) (v).begin(), (v).end()
#define MASK(i) (1LL << (i))
#define GETBIT(mask, i) (((mask) >> (i)) & 1)
 
// mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());
mt19937_64 rng(1);
ll rngesus(ll l, ll r){return ((ull) rng()) % (r - l + 1) + l;}
 
ll max(ll a, ll b){return (a > b) ? a : b;}
ll min(ll a, ll b){return (a < b) ? a : b;}
 
ll LASTBIT(ll mask){return mask & (-mask);}
ll pop_cnt(ll mask){return __builtin_popcountll(mask);}
ll ctz(ll mask){return __builtin_ctzll(mask);}
ll clz(ll mask){return __builtin_clzll(mask);}
ll logOf(ll mask){return 63 - clz(mask);}
 
template <class T1, class T2>
    bool minimize(T1 &a, T2 b){
        if (a > b){a = b; return true;}
        return false;
    }
template <class T1, class T2>
    bool maximize(T1 &a, T2 b){
        if (a < b){a = b; return true;}
        return false;
    }
template <class T>
    void printArr(T& a, string separator = " ", string finish = "\n", ostream& out = cout){
        for(auto i: a) out << i << separator;
        out << finish;
    }
template <class T>
    void remove_dup(vector<T> &a){
        sort(ALL(a));
        a.resize(unique(ALL(a)) - a.begin());
    }

const int N = 3e5 + 69;
const ll INF = 1e9;
int n; 
int ans[N];
vector<array<ll, 4>> a;

ll sqr(ll x){return x * x;}
ll block_sz;
map<pair<ll, ll>, vector<int>> grid;

pair<ll, ll> dumb_down(pair<ll, ll> a){
    return {a.first / block_sz, a.second / block_sz};
}

void set_up(ll arg){
    block_sz = arg;

    grid.clear();
    for(int i= 1; i<=n; ++i) if (ans[i] == 0){
        grid[dumb_down(make_pair(a[i][0], a[i][1]))].push_back(i);
    }
}

ll dis(pair<ll, ll> a, pair<ll, ll> b){return sqr(a.first - b.first) + sqr(a.second - b.second);}

bool check(pair<ll, ll> point, ll r, pair<ll, ll> des){
    if (des.first < 0 || des.second < 0) return false;
    for(int x = 0; x <= 1; ++x) for(int y = 0; y<=1; ++y){
        pair<ll, ll> ligma = {des.first * block_sz, des.second * block_sz};
        if (x) ligma.first += block_sz - 1;
        if (y) ligma.second += block_sz - 1;
        if (dis(point, ligma) <= sqr(r)) return true;
    }
    return false;
}

bool query(int i, int j){
    return sqr(a[i][0] - a[j][0]) + sqr(a[i][1] - a[j][1]) <= sqr(a[i][2] + a[j][2]);
}

int main(void){
    ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);
 
    cin >> n;
    a.resize(n+1);
    for(int i = 1; i<=n; ++i) {
        for(int j = 0; j < 3; ++j) cin >> a[i][j];
        a[i][0] += INF; a[i][1] += INF;
        a[i][3] = i;
    }

    vector<array<ll, 4>> b = a;
    sort(1 + ALL(b), [](array<ll, 4> x, array<ll, 4> y){
        if (x[2] != y[2]) return x[2] > y[2];
        return x[3] < y[3];
    });

    set_up((b[1][2] + 1) / 2);
    for(int i = 1; i<=n; ++i) if (ans[b[i][3]] == 0){
        int _i = b[i][3];
        ans[_i] = _i;

        if (b[i][2] < block_sz)
            set_up((b[i][2] + 1) / 2);

        pair<ll, ll> cur = make_pair(b[i][0], b[i][1]);
        pair<ll, ll> magnifico = dumb_down(cur);
        const int EXTEND = 10;
        for(int x = -EXTEND; x <= EXTEND; ++x) for(int y = -EXTEND; y<=EXTEND; ++y) {
            pair<ll, ll> cu = make_pair(magnifico.first + x, magnifico.second + y);
            
            if (max(abs(x), abs(y)) == EXTEND + 1 && check(cur, b[i][2] * 2, cu)) exit(1);

            if (!check(cur, b[i][2] * 2, cu) || !grid.count(cu)) continue;
            for(int j: grid[cu]){
                if (ans[j] == 0 && query(_i, j)) ans[j] = _i;
            }
        }
    }

    for(int i = 1; i<=n; ++i) 
        cout << ans[i] << " "; 
    cout << "\n";

    return 0;
}

Subtask #1
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 0 ms 348 KB Output is correct
2 Correct 0 ms 348 KB Output is correct
3 Correct 1 ms 348 KB Output is correct
4 Correct 1 ms 348 KB Output is correct
5 Correct 1 ms 348 KB Output is correct
6 Correct 0 ms 348 KB Output is correct
7 Correct 0 ms 348 KB Output is correct
8 Correct 1 ms 348 KB Output is correct
9 Correct 0 ms 348 KB Output is correct
10 Correct 0 ms 348 KB Output is correct
11 Correct 1 ms 344 KB Output is correct
12 Correct 1 ms 348 KB Output is correct
13 Correct 1 ms 348 KB Output is correct
14 Correct 1 ms 600 KB Output is correct
15 Correct 2 ms 348 KB Output is correct
16 Correct 1 ms 348 KB Output is correct
17 Correct 1 ms 344 KB Output is correct
18 Correct 1 ms 348 KB Output is correct
19 Correct 3 ms 788 KB Output is correct
20 Correct 3 ms 1048 KB Output is correct
21 Correct 2 ms 860 KB Output is correct
22 Correct 27 ms 1500 KB Output is correct
23 Correct 37 ms 1624 KB Output is correct
24 Correct 30 ms 1748 KB Output is correct
25 Correct 28 ms 1372 KB Output is correct
26 Correct 31 ms 1372 KB Output is correct

Subtask #2
0 / 12.0

# Verdict Execution time Memory Grader output
1 Correct 113 ms 23032 KB Output is correct
2 Correct 128 ms 22708 KB Output is correct
3 Correct 121 ms 23468 KB Output is correct
4 Correct 127 ms 22524 KB Output is correct
5 Incorrect 195 ms 30340 KB Output isn't correct
6 Halted 0 ms 0 KB -

Subtask #3
15.0 / 15.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 344 KB Output is correct
2 Correct 776 ms 18512 KB Output is correct
3 Correct 2953 ms 55076 KB Output is correct
4 Correct 2728 ms 55076 KB Output is correct
5 Correct 2433 ms 51700 KB Output is correct
6 Correct 796 ms 26648 KB Output is correct
7 Correct 344 ms 14392 KB Output is correct
8 Correct 54 ms 3332 KB Output is correct
9 Correct 2917 ms 54784 KB Output is correct
10 Correct 2062 ms 51364 KB Output is correct

Subtask #4
0 / 23.0

# Verdict Execution time Memory Grader output
1 Correct 2842 ms 55116 KB Output is correct
2 Correct 1670 ms 55128 KB Output is correct
3 Incorrect 583 ms 44152 KB Output isn't correct
4 Halted 0 ms 0 KB -

Subtask #5
30.0 / 30.0

# Verdict Execution time Memory Grader output
1 Correct 0 ms 348 KB Output is correct
2 Correct 0 ms 348 KB Output is correct
3 Correct 1 ms 348 KB Output is correct
4 Correct 1 ms 348 KB Output is correct
5 Correct 1 ms 348 KB Output is correct
6 Correct 0 ms 348 KB Output is correct
7 Correct 0 ms 348 KB Output is correct
8 Correct 1 ms 348 KB Output is correct
9 Correct 0 ms 348 KB Output is correct
10 Correct 0 ms 348 KB Output is correct
11 Correct 1 ms 344 KB Output is correct
12 Correct 1 ms 348 KB Output is correct
13 Correct 1 ms 348 KB Output is correct
14 Correct 1 ms 600 KB Output is correct
15 Correct 2 ms 348 KB Output is correct
16 Correct 1 ms 348 KB Output is correct
17 Correct 1 ms 344 KB Output is correct
18 Correct 1 ms 348 KB Output is correct
19 Correct 3 ms 788 KB Output is correct
20 Correct 3 ms 1048 KB Output is correct
21 Correct 2 ms 860 KB Output is correct
22 Correct 27 ms 1500 KB Output is correct
23 Correct 37 ms 1624 KB Output is correct
24 Correct 30 ms 1748 KB Output is correct
25 Correct 28 ms 1372 KB Output is correct
26 Correct 31 ms 1372 KB Output is correct
27 Correct 5 ms 1360 KB Output is correct
28 Correct 5 ms 1368 KB Output is correct
29 Correct 6 ms 1368 KB Output is correct
30 Correct 60 ms 2484 KB Output is correct
31 Correct 76 ms 2392 KB Output is correct
32 Correct 55 ms 2440 KB Output is correct
33 Correct 49 ms 11708 KB Output is correct
34 Correct 50 ms 11472 KB Output is correct
35 Correct 56 ms 10716 KB Output is correct
36 Correct 644 ms 21116 KB Output is correct
37 Correct 759 ms 21216 KB Output is correct
38 Correct 720 ms 21260 KB Output is correct
39 Correct 305 ms 18128 KB Output is correct
40 Correct 328 ms 17920 KB Output is correct
41 Correct 308 ms 17972 KB Output is correct
42 Correct 309 ms 20076 KB Output is correct
43 Correct 626 ms 21076 KB Output is correct
44 Correct 632 ms 21084 KB Output is correct
45 Correct 635 ms 21060 KB Output is correct
46 Correct 640 ms 21076 KB Output is correct
47 Correct 620 ms 21072 KB Output is correct
48 Correct 677 ms 21088 KB Output is correct
49 Correct 637 ms 21072 KB Output is correct
50 Correct 608 ms 21076 KB Output is correct

Subtask #6
0 / 13.0

# Verdict Execution time Memory Grader output
1 Correct 0 ms 348 KB Output is correct
2 Correct 0 ms 348 KB Output is correct
3 Correct 1 ms 348 KB Output is correct
4 Correct 1 ms 348 KB Output is correct
5 Correct 1 ms 348 KB Output is correct
6 Correct 0 ms 348 KB Output is correct
7 Correct 0 ms 348 KB Output is correct
8 Correct 1 ms 348 KB Output is correct
9 Correct 0 ms 348 KB Output is correct
10 Correct 0 ms 348 KB Output is correct
11 Correct 1 ms 344 KB Output is correct
12 Correct 1 ms 348 KB Output is correct
13 Correct 1 ms 348 KB Output is correct
14 Correct 1 ms 600 KB Output is correct
15 Correct 2 ms 348 KB Output is correct
16 Correct 1 ms 348 KB Output is correct
17 Correct 1 ms 344 KB Output is correct
18 Correct 1 ms 348 KB Output is correct
19 Correct 3 ms 788 KB Output is correct
20 Correct 3 ms 1048 KB Output is correct
21 Correct 2 ms 860 KB Output is correct
22 Correct 27 ms 1500 KB Output is correct
23 Correct 37 ms 1624 KB Output is correct
24 Correct 30 ms 1748 KB Output is correct
25 Correct 28 ms 1372 KB Output is correct
26 Correct 31 ms 1372 KB Output is correct
27 Correct 113 ms 23032 KB Output is correct
28 Correct 128 ms 22708 KB Output is correct
29 Correct 121 ms 23468 KB Output is correct
30 Correct 127 ms 22524 KB Output is correct
31 Incorrect 195 ms 30340 KB Output isn't correct
32 Halted 0 ms 0 KB -