Submission #391062

# Submission time Handle Problem Language Result Execution time Memory
391062 2021-04-17T17:50:30 Z phathnv Circle selection (APIO18_circle_selection) C++11
7 / 100
3000 ms 46264 KB
#include <bits/stdc++.h>
using namespace std;

typedef long long ll;

const int N = 3e5 + 7;

struct Circle{
    int x, y, r;
};

int n, answer[N], ord[N], ordX[N], numRow, inRow[N];
vector<int> row[N];
Circle c[N];

bool Intersect(const Circle &c1, const Circle &c2){
    int dx = c1.x - c2.x;
    int dy = c1.y - c2.y;
    int sr = c1.r + c2.r;
    return (ll) dx * dx + (ll) dy * dy <= (ll) sr * sr;
}

void Rebuild(int k){
    for(int i = 1; i <= numRow; i++)
        row[i].clear();

    numRow = 0;
    int last = -1;
    for(int i = 1; i <= n; i++){
        int ind = ordX[i];
        if (answer[ind])
            continue;
        if ((c[ind].x >> k) != last){
            last = c[ind].x >> k;
            numRow++;
        }
        row[numRow].push_back(ind);
        inRow[ind] = numRow;
    }
    for(int i = 1; i <= numRow; i++)
        sort(row[i].begin(), row[i].end(), [&](const int &a, const int &b){
                return (c[a].y >> k) < (c[b].y >> k);
             });
}

int main(){
    ios_base::sync_with_stdio(0);
    cin.tie(0);

    cin >> n;
    for(int i = 1; i <= n; i++){
        cin >> c[i].x >> c[i].y >> c[i].r;
        c[i].x += 1e9;
        c[i].y += 1e9;
        ord[i] = i;
        ordX[i] = i;
    }

    sort(ord + 1, ord + 1 + n, [&](const int &a, const int &b){
            if (c[a].r != c[b].r)
                return c[a].r > c[b].r;
            return a < b;
         });
    sort(ordX + 1, ordX + 1 + n, [&](const int &a, const int &b){
            return c[a].x < c[b].x;
         });
    int k = 30;
    Rebuild(k);
    for(int i = 1; i <= n; i++){
        int ind = ord[i];
        if (answer[ind])
            continue;
        while (c[ind].r >= (1 << (k - 1)))
            Rebuild(--k);

        for(int x = max(1, inRow[ind] - 2); x <= min(numRow, inRow[ind] + 2); x++){
            if (abs((c[row[x][0]].x >> k) - (c[ind].x >> k)) > 2)
                continue;
            int l = 0, r = row[x].size() - 1, p = row[x].size();
            while (l <= r){
                int mid = (l + r) >> 1;
                if ((c[row[x][mid]].y >> k) >= (c[ind].y >> k) - 2){
                    p = mid;
                    r = mid - 1;
                } else {
                    l = mid + 1;
                }
            }

            while (p < (int) row[x].size()){
                if (Intersect(c[ind], c[row[x][p]]) && !answer[row[x][p]])
                    answer[row[x][p]] = ind;
                if ((c[row[x][p]].y >> k) > (c[ind].y >> k) + 2)
                    break;
                p++;
            }
         }
    }

    for(int i = 1; i <= n; i++)
        cout << answer[i] << ' ';

    return 0;
}
# Verdict Execution time Memory Grader output
1 Correct 4 ms 7372 KB Output is correct
2 Correct 5 ms 7372 KB Output is correct
3 Correct 4 ms 7372 KB Output is correct
4 Correct 5 ms 7372 KB Output is correct
5 Correct 5 ms 7368 KB Output is correct
6 Correct 5 ms 7372 KB Output is correct
7 Correct 5 ms 7372 KB Output is correct
8 Correct 5 ms 7372 KB Output is correct
9 Correct 5 ms 7376 KB Output is correct
10 Correct 5 ms 7372 KB Output is correct
11 Correct 5 ms 7372 KB Output is correct
12 Correct 5 ms 7372 KB Output is correct
13 Correct 5 ms 7372 KB Output is correct
14 Correct 4 ms 7372 KB Output is correct
15 Correct 5 ms 7372 KB Output is correct
16 Correct 7 ms 7436 KB Output is correct
17 Correct 7 ms 7420 KB Output is correct
18 Correct 7 ms 7420 KB Output is correct
19 Correct 14 ms 7876 KB Output is correct
20 Correct 13 ms 7952 KB Output is correct
21 Correct 14 ms 7908 KB Output is correct
22 Correct 99 ms 7616 KB Output is correct
23 Correct 102 ms 7628 KB Output is correct
24 Correct 101 ms 7708 KB Output is correct
25 Correct 104 ms 7700 KB Output is correct
26 Correct 100 ms 7628 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 236 ms 19312 KB Output is correct
2 Correct 252 ms 18916 KB Output is correct
3 Correct 247 ms 25432 KB Output is correct
4 Correct 689 ms 46264 KB Output is correct
5 Execution timed out 3056 ms 21456 KB Time limit exceeded
6 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 5 ms 7372 KB Output is correct
2 Execution timed out 3086 ms 10724 KB Time limit exceeded
3 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 3093 ms 16108 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 4 ms 7372 KB Output is correct
2 Correct 5 ms 7372 KB Output is correct
3 Correct 4 ms 7372 KB Output is correct
4 Correct 5 ms 7372 KB Output is correct
5 Correct 5 ms 7368 KB Output is correct
6 Correct 5 ms 7372 KB Output is correct
7 Correct 5 ms 7372 KB Output is correct
8 Correct 5 ms 7372 KB Output is correct
9 Correct 5 ms 7376 KB Output is correct
10 Correct 5 ms 7372 KB Output is correct
11 Correct 5 ms 7372 KB Output is correct
12 Correct 5 ms 7372 KB Output is correct
13 Correct 5 ms 7372 KB Output is correct
14 Correct 4 ms 7372 KB Output is correct
15 Correct 5 ms 7372 KB Output is correct
16 Correct 7 ms 7436 KB Output is correct
17 Correct 7 ms 7420 KB Output is correct
18 Correct 7 ms 7420 KB Output is correct
19 Correct 14 ms 7876 KB Output is correct
20 Correct 13 ms 7952 KB Output is correct
21 Correct 14 ms 7908 KB Output is correct
22 Correct 99 ms 7616 KB Output is correct
23 Correct 102 ms 7628 KB Output is correct
24 Correct 101 ms 7708 KB Output is correct
25 Correct 104 ms 7700 KB Output is correct
26 Correct 100 ms 7628 KB Output is correct
27 Correct 12 ms 8016 KB Output is correct
28 Correct 11 ms 8004 KB Output is correct
29 Correct 12 ms 8008 KB Output is correct
30 Correct 429 ms 8016 KB Output is correct
31 Correct 382 ms 8016 KB Output is correct
32 Correct 387 ms 8012 KB Output is correct
33 Correct 223 ms 21064 KB Output is correct
34 Correct 215 ms 21300 KB Output is correct
35 Correct 147 ms 14196 KB Output is correct
36 Execution timed out 3088 ms 13152 KB Time limit exceeded
37 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 4 ms 7372 KB Output is correct
2 Correct 5 ms 7372 KB Output is correct
3 Correct 4 ms 7372 KB Output is correct
4 Correct 5 ms 7372 KB Output is correct
5 Correct 5 ms 7368 KB Output is correct
6 Correct 5 ms 7372 KB Output is correct
7 Correct 5 ms 7372 KB Output is correct
8 Correct 5 ms 7372 KB Output is correct
9 Correct 5 ms 7376 KB Output is correct
10 Correct 5 ms 7372 KB Output is correct
11 Correct 5 ms 7372 KB Output is correct
12 Correct 5 ms 7372 KB Output is correct
13 Correct 5 ms 7372 KB Output is correct
14 Correct 4 ms 7372 KB Output is correct
15 Correct 5 ms 7372 KB Output is correct
16 Correct 7 ms 7436 KB Output is correct
17 Correct 7 ms 7420 KB Output is correct
18 Correct 7 ms 7420 KB Output is correct
19 Correct 14 ms 7876 KB Output is correct
20 Correct 13 ms 7952 KB Output is correct
21 Correct 14 ms 7908 KB Output is correct
22 Correct 99 ms 7616 KB Output is correct
23 Correct 102 ms 7628 KB Output is correct
24 Correct 101 ms 7708 KB Output is correct
25 Correct 104 ms 7700 KB Output is correct
26 Correct 100 ms 7628 KB Output is correct
27 Correct 236 ms 19312 KB Output is correct
28 Correct 252 ms 18916 KB Output is correct
29 Correct 247 ms 25432 KB Output is correct
30 Correct 689 ms 46264 KB Output is correct
31 Execution timed out 3056 ms 21456 KB Time limit exceeded
32 Halted 0 ms 0 KB -