oj.uz

Submission #252517

# Submission time Handle Problem Language Result Execution time Memory
252517 2020-07-25T17:42:29 Z shayan_p Circle selection (APIO18_circle_selection) C++14
87 / 100
 3000 ms 172544 KB 
// And you curse yourself for things you never done

#include<bits/stdc++.h>

#define F first
#define S second
#define PB push_back
#define sz(s) int((s).size())
#define bit(n,k) (((n)>>(k))&1)

using namespace std;

typedef long long ll;
typedef pair<int,int> pii;

const int maxn = 3e5 + 10, mod = 1e9 + 7, inf = 1e9 + 10;
const ll INF = 1e18;

pair<pii, int> p[maxn];
int arr[maxn], ans[maxn];

pii operator - (pii a, pii b){
    return {a.F - b.F, a.S - b.S};
}
pii operator + (pii a, pii b){
    return {a.F + b.F, a.S + b.S};
}
ll operator ^ (pii a, pii b){
    return 1ll * a.F * b.F + 1ll * a.S * b.S;
}
bool ok(int i, int j){
    return 1ll * (p[j].S + p[i].S) * (p[j].S + p[i].S) >= ((p[i].F - p[j].F) ^ (p[i].F - p[j].F));
}

vector<int> vx, vy;    

vector<int> val, fs;
vector<int> stock;

void push(int &lst, int x){
    fs.PB(lst);
    lst = sz(val);
    val.PB(x);
}

struct node{
    map<int, int> mp;
    void add(int pos, int x){
	mp[pos] = x;
    }
    void ask(int pos){
	auto it = mp.lower_bound(pos);
	if(it != mp.end())
	    stock.PB(it->S);
	if(it != mp.begin())
	    stock.PB(prev(it)->S);
    }
};
struct node2{
    node2 *L = 0, *R = 0;
    node* seg = 0;
    void add(int f, int s, int pos, int x, int l = 0, int r = sz(vx)-1){
	if(f <= vx[l] && vx[r] <= s){
	    if(!seg)
		seg = new node();
	    seg->add(pos, x);
	    return;
	}
	if(r-l == 1)
	    return;
	int mid = (l+r) >> 1;	
	if(max(f, vx[l]) < min(s, vx[mid])){
	    if(!L)
		L = new node2();
	    L->add(f, s, pos, x, l, mid);
	}
	if(max(f, vx[mid]) < min(s, vx[r])){
	    if(!R)
		R = new node2();
	    R->add(f, s, pos, x, mid, r);
	}
    }
    void ask(int pos, int pos2, int l = 0, int r = sz(vx)-1){
	if(seg)
	    seg->ask(pos2);
	if(r-l == 1)
	    return;
	int mid = (l+r) >> 1;
	if(pos < vx[mid] && L)
	    L->ask(pos, pos2, l, mid);
	if(vx[mid] <= pos && R)
	    R->ask(pos, pos2, mid, r);	    
    }
};


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

    node2* root = new node2();
    
    int n;
    cin >> n;

    for(int i = 0; i < n; i++){
	cin >> p[i].F.F >> p[i].F.S >> p[i].S;
	
	vx.PB(p[i].F.F - p[i].S);
	vx.PB(p[i].F.F + p[i].S + 1);
	
	vy.PB(p[i].F.S - p[i].S);
	vy.PB(p[i].F.S + p[i].S + 1);
	
	arr[i] = i;
    }
    auto cmp = [](int i, int j){
		   return (p[i].S == p[j].S ? i > j : p[i].S < p[j].S);
	       };
    sort(arr, arr + n, cmp);

    sort(vx.begin(), vx.end());
    vx.resize( unique(vx.begin(), vx.end()) - vx.begin() );

    sort(vy.begin(), vy.end());
    vy.resize( unique(vy.begin(), vy.end()) - vy.begin() );

    for(int i = n-1; i >= 0; i--){
	int I = arr[i];
	stock.clear();
	for(int w = -1; w <= 1; w++){
	    for(int w2 = -1; w2 <= 1; w2++)
		if((w != 0 || w2 != 0))
		    root->ask(p[I].F.F + w * p[I].S, p[I].F.S + w2 * p[I].S); // kafie?
	}
	for(int id : stock){
	    if(ok(I, id)){
		if(ans[I] == 0 || cmp(ans[I] - 1, id))
		    ans[I] = id + 1;
	    }
	}
	if(ans[I] == 0){
	    ans[I] = I + 1;
	    root->add(p[I].F.F - p[I].S, p[I].F.F + p[I].S + 1, p[I].F.S, I);
	}
    }
    for(int i = 0; i < n; i++){
	cout << ans[i] << " ";
    }
    return cout << endl, 0;
}

Subtask #1
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 384 KB Output is correct
2 Correct 0 ms 384 KB Output is correct
3 Correct 0 ms 384 KB Output is correct
4 Correct 0 ms 384 KB Output is correct
5 Correct 0 ms 384 KB Output is correct
6 Correct 1 ms 384 KB Output is correct
7 Correct 0 ms 384 KB Output is correct
8 Correct 0 ms 384 KB Output is correct
9 Correct 1 ms 384 KB Output is correct
10 Correct 0 ms 384 KB Output is correct
11 Correct 1 ms 384 KB Output is correct
12 Correct 1 ms 384 KB Output is correct
13 Correct 0 ms 384 KB Output is correct
14 Correct 1 ms 384 KB Output is correct
15 Correct 1 ms 384 KB Output is correct
16 Correct 2 ms 384 KB Output is correct
17 Correct 1 ms 384 KB Output is correct
18 Correct 1 ms 384 KB Output is correct
19 Correct 5 ms 640 KB Output is correct
20 Correct 5 ms 640 KB Output is correct
21 Correct 6 ms 640 KB Output is correct
22 Correct 11 ms 2048 KB Output is correct
23 Correct 20 ms 3072 KB Output is correct
24 Correct 11 ms 1920 KB Output is correct
25 Correct 10 ms 1664 KB Output is correct
26 Correct 11 ms 1920 KB Output is correct

Subtask #2
12.0 / 12.0

# Verdict Execution time Memory Grader output
1 Correct 368 ms 13328 KB Output is correct
2 Correct 413 ms 13260 KB Output is correct
3 Correct 374 ms 12992 KB Output is correct
4 Correct 360 ms 13404 KB Output is correct
5 Correct 652 ms 17608 KB Output is correct
6 Correct 824 ms 52528 KB Output is correct
7 Correct 705 ms 23552 KB Output is correct
8 Correct 824 ms 32256 KB Output is correct

Subtask #3
15.0 / 15.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 384 KB Output is correct
2 Correct 626 ms 53216 KB Output is correct
3 Correct 1835 ms 156628 KB Output is correct
4 Correct 1812 ms 156224 KB Output is correct
5 Correct 2241 ms 141120 KB Output is correct
6 Correct 913 ms 58072 KB Output is correct
7 Correct 333 ms 28000 KB Output is correct
8 Correct 38 ms 4468 KB Output is correct
9 Correct 1883 ms 157472 KB Output is correct
10 Correct 2677 ms 124864 KB Output is correct

Subtask #4
23.0 / 23.0

# Verdict Execution time Memory Grader output
1 Correct 1473 ms 145332 KB Output is correct
2 Correct 756 ms 74184 KB Output is correct
3 Correct 2201 ms 30796 KB Output is correct
4 Correct 873 ms 84168 KB Output is correct
5 Correct 827 ms 78780 KB Output is correct
6 Correct 1532 ms 17344 KB Output is correct

Subtask #5
30.0 / 30.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 384 KB Output is correct
2 Correct 0 ms 384 KB Output is correct
3 Correct 0 ms 384 KB Output is correct
4 Correct 0 ms 384 KB Output is correct
5 Correct 0 ms 384 KB Output is correct
6 Correct 1 ms 384 KB Output is correct
7 Correct 0 ms 384 KB Output is correct
8 Correct 0 ms 384 KB Output is correct
9 Correct 1 ms 384 KB Output is correct
10 Correct 0 ms 384 KB Output is correct
11 Correct 1 ms 384 KB Output is correct
12 Correct 1 ms 384 KB Output is correct
13 Correct 0 ms 384 KB Output is correct
14 Correct 1 ms 384 KB Output is correct
15 Correct 1 ms 384 KB Output is correct
16 Correct 2 ms 384 KB Output is correct
17 Correct 1 ms 384 KB Output is correct
18 Correct 1 ms 384 KB Output is correct
19 Correct 5 ms 640 KB Output is correct
20 Correct 5 ms 640 KB Output is correct
21 Correct 6 ms 640 KB Output is correct
22 Correct 11 ms 2048 KB Output is correct
23 Correct 20 ms 3072 KB Output is correct
24 Correct 11 ms 1920 KB Output is correct
25 Correct 10 ms 1664 KB Output is correct
26 Correct 11 ms 1920 KB Output is correct
27 Correct 12 ms 896 KB Output is correct
28 Correct 12 ms 896 KB Output is correct
29 Correct 11 ms 896 KB Output is correct
30 Correct 32 ms 5136 KB Output is correct
31 Correct 24 ms 3832 KB Output is correct
32 Correct 26 ms 3712 KB Output is correct
33 Correct 118 ms 4704 KB Output is correct
34 Correct 118 ms 4672 KB Output is correct
35 Correct 130 ms 4676 KB Output is correct
36 Correct 344 ms 39096 KB Output is correct
37 Correct 367 ms 40036 KB Output is correct
38 Correct 382 ms 44000 KB Output is correct
39 Correct 317 ms 7652 KB Output is correct
40 Correct 326 ms 7652 KB Output is correct
41 Correct 320 ms 7652 KB Output is correct
42 Correct 322 ms 7616 KB Output is correct
43 Correct 230 ms 24612 KB Output is correct
44 Correct 244 ms 24544 KB Output is correct
45 Correct 227 ms 24552 KB Output is correct
46 Correct 232 ms 24548 KB Output is correct
47 Correct 233 ms 24676 KB Output is correct
48 Correct 227 ms 24544 KB Output is correct
49 Correct 231 ms 24548 KB Output is correct
50 Correct 230 ms 24676 KB Output is correct

Subtask #6
0 / 13.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 384 KB Output is correct
2 Correct 0 ms 384 KB Output is correct
3 Correct 0 ms 384 KB Output is correct
4 Correct 0 ms 384 KB Output is correct
5 Correct 0 ms 384 KB Output is correct
6 Correct 1 ms 384 KB Output is correct
7 Correct 0 ms 384 KB Output is correct
8 Correct 0 ms 384 KB Output is correct
9 Correct 1 ms 384 KB Output is correct
10 Correct 0 ms 384 KB Output is correct
11 Correct 1 ms 384 KB Output is correct
12 Correct 1 ms 384 KB Output is correct
13 Correct 0 ms 384 KB Output is correct
14 Correct 1 ms 384 KB Output is correct
15 Correct 1 ms 384 KB Output is correct
16 Correct 2 ms 384 KB Output is correct
17 Correct 1 ms 384 KB Output is correct
18 Correct 1 ms 384 KB Output is correct
19 Correct 5 ms 640 KB Output is correct
20 Correct 5 ms 640 KB Output is correct
21 Correct 6 ms 640 KB Output is correct
22 Correct 11 ms 2048 KB Output is correct
23 Correct 20 ms 3072 KB Output is correct
24 Correct 11 ms 1920 KB Output is correct
25 Correct 10 ms 1664 KB Output is correct
26 Correct 11 ms 1920 KB Output is correct
27 Correct 368 ms 13328 KB Output is correct
28 Correct 413 ms 13260 KB Output is correct
29 Correct 374 ms 12992 KB Output is correct
30 Correct 360 ms 13404 KB Output is correct
31 Correct 652 ms 17608 KB Output is correct
32 Correct 824 ms 52528 KB Output is correct
33 Correct 705 ms 23552 KB Output is correct
34 Correct 824 ms 32256 KB Output is correct
35 Correct 1 ms 384 KB Output is correct
36 Correct 626 ms 53216 KB Output is correct
37 Correct 1835 ms 156628 KB Output is correct
38 Correct 1812 ms 156224 KB Output is correct
39 Correct 2241 ms 141120 KB Output is correct
40 Correct 913 ms 58072 KB Output is correct
41 Correct 333 ms 28000 KB Output is correct
42 Correct 38 ms 4468 KB Output is correct
43 Correct 1883 ms 157472 KB Output is correct
44 Correct 2677 ms 124864 KB Output is correct
45 Correct 1473 ms 145332 KB Output is correct
46 Correct 756 ms 74184 KB Output is correct
47 Correct 2201 ms 30796 KB Output is correct
48 Correct 873 ms 84168 KB Output is correct
49 Correct 827 ms 78780 KB Output is correct
50 Correct 1532 ms 17344 KB Output is correct
51 Correct 12 ms 896 KB Output is correct
52 Correct 12 ms 896 KB Output is correct
53 Correct 11 ms 896 KB Output is correct
54 Correct 32 ms 5136 KB Output is correct
55 Correct 24 ms 3832 KB Output is correct
56 Correct 26 ms 3712 KB Output is correct
57 Correct 118 ms 4704 KB Output is correct
58 Correct 118 ms 4672 KB Output is correct
59 Correct 130 ms 4676 KB Output is correct
60 Correct 344 ms 39096 KB Output is correct
61 Correct 367 ms 40036 KB Output is correct
62 Correct 382 ms 44000 KB Output is correct
63 Correct 317 ms 7652 KB Output is correct
64 Correct 326 ms 7652 KB Output is correct
65 Correct 320 ms 7652 KB Output is correct
66 Correct 322 ms 7616 KB Output is correct
67 Correct 230 ms 24612 KB Output is correct
68 Correct 244 ms 24544 KB Output is correct
69 Correct 227 ms 24552 KB Output is correct
70 Correct 232 ms 24548 KB Output is correct
71 Correct 233 ms 24676 KB Output is correct
72 Correct 227 ms 24544 KB Output is correct
73 Correct 231 ms 24548 KB Output is correct
74 Correct 230 ms 24676 KB Output is correct
75 Correct 424 ms 12948 KB Output is correct
76 Correct 392 ms 13256 KB Output is correct
77 Correct 413 ms 13248 KB Output is correct
78 Correct 383 ms 13140 KB Output is correct
79 Correct 769 ms 14912 KB Output is correct
80 Correct 419 ms 13000 KB Output is correct
81 Correct 1893 ms 150604 KB Output is correct
82 Correct 1618 ms 138436 KB Output is correct
83 Correct 1669 ms 138952 KB Output is correct
84 Correct 2530 ms 167836 KB Output is correct
85 Correct 1581 ms 146852 KB Output is correct
86 Correct 1454 ms 135924 KB Output is correct
87 Execution timed out 3072 ms 172544 KB Time limit exceeded
88 Halted 0 ms 0 KB -