답안 #397216

# 제출 시각 아이디 문제 언어 결과 실행 시간 메모리
397216 2021-05-01T17:42:15 Z parsabahrami 원 고르기 (APIO18_circle_selection) C++17
7 / 100
3000 ms 327040 KB
/* There's someone in my head but it's not me */ 
#include <bits/stdc++.h>
 
using namespace std;

typedef long long int ll;
typedef pair<int, int> pii;
 
#define SZ(x)                       (int) x.size()
#define F                           first
#define S                           second

const int N = 3e5 + 10, MOD = 1e9 + 7;
int n, X[N], Y[N], R[N], ord[N], xord[N], pos[N], ret[N];
struct node : set<pii> {
    node *lc, *rc; 

    node() : lc(nullptr), rc(nullptr) {}

    void build(int l = 1, int r = n + 1) {
        if (r - l < 2) return;
        int md = (l + r) >> 1;
        lc = new node(), rc = new node();
        lc->build(l, md), rc->build(md, r);
    }

    void remove(int p, int id, int l = 1, int r = n + 1) {
        this->erase({Y[id], id});
        if (r - l < 2) return;
        int md = (l + r) >> 1;
        if (p < md) lc->remove(p, id, l, md);
        else rc->remove(p, id, md, r);
    }

    int get(ll xl, ll xr, ll yl, ll yr, int l = 1, int r = n + 1) {
        if (xr <= X[xord[l]] || X[xord[r - 1]] + 1 <= xl || this->empty()) return -1;
        if (xl <= X[xord[l]] && X[xord[r - 1]] + 1 <= xr) {
            auto lb = this->lower_bound({yl, -1e9});
            if (lb != this->end() && lb->F <= yr) return lb->S;
            else return -1;
        }
        int md = (l + r) >> 1, a;
        a = lc->get(xl, xr, yl, yr, l, md);
        if (~a) return a;
        return rc->get(xl, xr, yl, yr, md, r);
    }
    
    void ins(int p, int y, int id, int l = 1, int r = n + 1) {
        this->insert({y, id});
        if (r - l < 2) return;
        int md = (l + r) >> 1;
        if (p < md) lc->ins(p, y, id, l, md);
        else rc->ins(p, y, id, md, r);
    }

};
node* seg;

int intersect(int i, int j) {
    ll x1 = X[i], x2 = X[j], y1 = Y[i], y2 = Y[j], r1 = R[i], r2 = R[j];
    return x1 * x1 + x2 * x2 + y1 * y1 + y2 * y2 - r1 * r1 - r2 * r2 <= 2 * (r1 * r2 + x1 * x2 + y1 * y2);
}

int main() {
    scanf("%d", &n);
    for (int i = 1; i <= n; i++) 
        scanf("%d%d%d", &X[i], &Y[i], &R[i]);
    iota(ord + 1, ord + n + 1, 1);
    iota(xord + 1, xord + n + 1, 1);
    sort(ord + 1, ord + n + 1, [&] (int i, int j) { return R[i] == R[j] ? i < j : R[i] > R[j]; });
    sort(xord + 1, xord + n + 1, [&] (int i, int j) { return X[i] < X[j]; });
    seg = new node(); seg->build();
    for (int i = 1; i <= n; i++) {
        seg->ins(i, Y[xord[i]], xord[i]), pos[xord[i]] = i;
    }
    for (int _ = 1; _ <= n; _++) {
        int i = ord[_]; vector<int> ids;
        if (ret[i]) continue;
        for (;;) {
            int id = seg->get(X[i] - 2ll * R[i], X[i] + 2ll * R[i] + 1, Y[i] - 2ll * R[i], Y[i] + 2ll * R[i]);
            if (~id) ids.push_back(id), seg->remove(pos[id], id);
            else break;
        }
        for (int &id : ids) {
            if (intersect(i, id)) ret[id] = i;
            else seg->ins(pos[id], Y[id], id);
        }
    }
    for (int i = 1; i <= n; i++) printf("%d ", ret[i]);
    printf("\n");
    return 0;
}

Compilation message

circle_selection.cpp: In function 'int main()':
circle_selection.cpp:65:10: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
   65 |     scanf("%d", &n);
      |     ~~~~~^~~~~~~~~~
circle_selection.cpp:67:14: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
   67 |         scanf("%d%d%d", &X[i], &Y[i], &R[i]);
      |         ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# 결과 실행 시간 메모리 Grader output
1 Correct 1 ms 304 KB Output is correct
2 Correct 1 ms 332 KB Output is correct
3 Correct 1 ms 332 KB Output is correct
4 Correct 1 ms 332 KB Output is correct
5 Correct 1 ms 332 KB Output is correct
6 Correct 1 ms 332 KB Output is correct
7 Correct 1 ms 332 KB Output is correct
8 Correct 1 ms 332 KB Output is correct
9 Correct 1 ms 332 KB Output is correct
10 Correct 1 ms 332 KB Output is correct
11 Correct 1 ms 332 KB Output is correct
12 Correct 1 ms 332 KB Output is correct
13 Correct 1 ms 332 KB Output is correct
14 Correct 1 ms 332 KB Output is correct
15 Correct 1 ms 332 KB Output is correct
16 Correct 4 ms 972 KB Output is correct
17 Correct 4 ms 972 KB Output is correct
18 Correct 4 ms 972 KB Output is correct
19 Correct 22 ms 4556 KB Output is correct
20 Correct 22 ms 4524 KB Output is correct
21 Correct 25 ms 4556 KB Output is correct
22 Correct 28 ms 4500 KB Output is correct
23 Correct 29 ms 4520 KB Output is correct
24 Correct 27 ms 4516 KB Output is correct
25 Correct 28 ms 4556 KB Output is correct
26 Correct 27 ms 4428 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Execution timed out 3037 ms 327040 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# 결과 실행 시간 메모리 Grader output
1 Correct 1 ms 332 KB Output is correct
2 Correct 1666 ms 102560 KB Output is correct
3 Execution timed out 3084 ms 326876 KB Time limit exceeded
4 Halted 0 ms 0 KB -
# 결과 실행 시간 메모리 Grader output
1 Execution timed out 3093 ms 325984 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# 결과 실행 시간 메모리 Grader output
1 Correct 1 ms 304 KB Output is correct
2 Correct 1 ms 332 KB Output is correct
3 Correct 1 ms 332 KB Output is correct
4 Correct 1 ms 332 KB Output is correct
5 Correct 1 ms 332 KB Output is correct
6 Correct 1 ms 332 KB Output is correct
7 Correct 1 ms 332 KB Output is correct
8 Correct 1 ms 332 KB Output is correct
9 Correct 1 ms 332 KB Output is correct
10 Correct 1 ms 332 KB Output is correct
11 Correct 1 ms 332 KB Output is correct
12 Correct 1 ms 332 KB Output is correct
13 Correct 1 ms 332 KB Output is correct
14 Correct 1 ms 332 KB Output is correct
15 Correct 1 ms 332 KB Output is correct
16 Correct 4 ms 972 KB Output is correct
17 Correct 4 ms 972 KB Output is correct
18 Correct 4 ms 972 KB Output is correct
19 Correct 22 ms 4556 KB Output is correct
20 Correct 22 ms 4524 KB Output is correct
21 Correct 25 ms 4556 KB Output is correct
22 Correct 28 ms 4500 KB Output is correct
23 Correct 29 ms 4520 KB Output is correct
24 Correct 27 ms 4516 KB Output is correct
25 Correct 28 ms 4556 KB Output is correct
26 Correct 27 ms 4428 KB Output is correct
27 Correct 54 ms 9284 KB Output is correct
28 Correct 51 ms 9372 KB Output is correct
29 Correct 52 ms 9224 KB Output is correct
30 Correct 74 ms 9180 KB Output is correct
31 Correct 69 ms 9168 KB Output is correct
32 Correct 67 ms 9172 KB Output is correct
33 Correct 1008 ms 105588 KB Output is correct
34 Correct 955 ms 105688 KB Output is correct
35 Correct 1056 ms 105576 KB Output is correct
36 Correct 1553 ms 104836 KB Output is correct
37 Correct 1596 ms 105020 KB Output is correct
38 Correct 1551 ms 104976 KB Output is correct
39 Execution timed out 3079 ms 103068 KB Time limit exceeded
40 Halted 0 ms 0 KB -
# 결과 실행 시간 메모리 Grader output
1 Correct 1 ms 304 KB Output is correct
2 Correct 1 ms 332 KB Output is correct
3 Correct 1 ms 332 KB Output is correct
4 Correct 1 ms 332 KB Output is correct
5 Correct 1 ms 332 KB Output is correct
6 Correct 1 ms 332 KB Output is correct
7 Correct 1 ms 332 KB Output is correct
8 Correct 1 ms 332 KB Output is correct
9 Correct 1 ms 332 KB Output is correct
10 Correct 1 ms 332 KB Output is correct
11 Correct 1 ms 332 KB Output is correct
12 Correct 1 ms 332 KB Output is correct
13 Correct 1 ms 332 KB Output is correct
14 Correct 1 ms 332 KB Output is correct
15 Correct 1 ms 332 KB Output is correct
16 Correct 4 ms 972 KB Output is correct
17 Correct 4 ms 972 KB Output is correct
18 Correct 4 ms 972 KB Output is correct
19 Correct 22 ms 4556 KB Output is correct
20 Correct 22 ms 4524 KB Output is correct
21 Correct 25 ms 4556 KB Output is correct
22 Correct 28 ms 4500 KB Output is correct
23 Correct 29 ms 4520 KB Output is correct
24 Correct 27 ms 4516 KB Output is correct
25 Correct 28 ms 4556 KB Output is correct
26 Correct 27 ms 4428 KB Output is correct
27 Execution timed out 3037 ms 327040 KB Time limit exceeded
28 Halted 0 ms 0 KB -