oj.uz

Submission #252492

# Submission time Handle Problem Language Result Execution time Memory
252492 2020-07-25T17:11:26 Z shayan_p Circle selection (APIO18_circle_selection) C++14
87 / 100
 3000 ms 222804 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));
}

const int FX = 20; ////////

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{
    set<pair<pii, int> > st;
    void add(int f, int s, int x){
	st.insert({{f, s}, x});
    }
    void ask(int pos){
	auto it = st.lower_bound( (pair<pii, int>){ {pos + 1, -1}, -1} );
	int bad = 2;
	while(sz(stock) < FX && it != st.begin()){
	    --it;
	    if(pos <= (it->F.S))
		stock.PB(it->S);
	    else{
		if(--bad == 0)
		    break;
	    }
	}
    }
};
struct node2{
    node2 *L = 0, *R = 0;
    node* seg = 0;
    void add(int f, int s, int f2, int s2, int x, int l = 0, int r = sz(vx)-1){
	if(f <= vx[l] && vx[r] <= s){
	    if(!seg)
		seg = new node();
	    seg->add(f2, s2, 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, f2, s2, x, l, mid);
	}
	if(max(f, vx[mid]) < min(s, vx[r])){
	    if(!R)
		R = new node2();
	    R->add(f, s, f2, s2, 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();
	//	cout << "IN " << I+1 << endl;
	for(int msk = 0; msk < 4; msk++)
	    root->ask(p[I].F.F + (bit(msk, 0) ? 1 : -1) * p[I].S, p[I].F.S + (bit(msk, 1) ? 1 : -1) * p[I].S);
	for(int id : stock){
	    //	    cout << "STOCK " << id+1 << endl;
	    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 - p[I].S, p[I].F.S + p[I].S + 1, 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 0 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 0 ms 384 KB Output is correct
7 Correct 1 ms 384 KB Output is correct
8 Correct 0 ms 384 KB Output is correct
9 Correct 0 ms 384 KB Output is correct
10 Correct 0 ms 384 KB Output is correct
11 Correct 0 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 0 ms 384 KB Output is correct
16 Correct 1 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 6 ms 640 KB Output is correct
20 Correct 6 ms 640 KB Output is correct
21 Correct 6 ms 640 KB Output is correct
22 Correct 11 ms 2304 KB Output is correct
23 Correct 17 ms 3456 KB Output is correct
24 Correct 10 ms 2048 KB Output is correct
25 Correct 9 ms 1664 KB Output is correct
26 Correct 10 ms 2048 KB Output is correct

Subtask #2
12.0 / 12.0

# Verdict Execution time Memory Grader output
1 Correct 347 ms 13248 KB Output is correct
2 Correct 353 ms 13260 KB Output is correct
3 Correct 349 ms 12992 KB Output is correct
4 Correct 333 ms 13260 KB Output is correct
5 Correct 507 ms 17864 KB Output is correct
6 Correct 742 ms 55832 KB Output is correct
7 Correct 620 ms 24264 KB Output is correct
8 Correct 653 ms 33868 KB Output is correct

Subtask #3
15.0 / 15.0

# Verdict Execution time Memory Grader output
1 Correct 0 ms 384 KB Output is correct
2 Correct 456 ms 60668 KB Output is correct
3 Correct 1507 ms 178368 KB Output is correct
4 Correct 1496 ms 177740 KB Output is correct
5 Correct 1865 ms 162116 KB Output is correct
6 Correct 754 ms 65716 KB Output is correct
7 Correct 292 ms 31096 KB Output is correct
8 Correct 38 ms 4744 KB Output is correct
9 Correct 1679 ms 180292 KB Output is correct
10 Correct 2378 ms 144376 KB Output is correct

Subtask #4
23.0 / 23.0

# Verdict Execution time Memory Grader output
1 Correct 1179 ms 164296 KB Output is correct
2 Correct 709 ms 78912 KB Output is correct
3 Correct 1710 ms 33356 KB Output is correct
4 Correct 759 ms 90572 KB Output is correct
5 Correct 749 ms 84284 KB Output is correct
6 Correct 1231 ms 17728 KB Output is correct

Subtask #5
30.0 / 30.0

# Verdict Execution time Memory Grader output
1 Correct 0 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 0 ms 384 KB Output is correct
7 Correct 1 ms 384 KB Output is correct
8 Correct 0 ms 384 KB Output is correct
9 Correct 0 ms 384 KB Output is correct
10 Correct 0 ms 384 KB Output is correct
11 Correct 0 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 0 ms 384 KB Output is correct
16 Correct 1 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 6 ms 640 KB Output is correct
20 Correct 6 ms 640 KB Output is correct
21 Correct 6 ms 640 KB Output is correct
22 Correct 11 ms 2304 KB Output is correct
23 Correct 17 ms 3456 KB Output is correct
24 Correct 10 ms 2048 KB Output is correct
25 Correct 9 ms 1664 KB Output is correct
26 Correct 10 ms 2048 KB Output is correct
27 Correct 11 ms 896 KB Output is correct
28 Correct 11 ms 896 KB Output is correct
29 Correct 11 ms 896 KB Output is correct
30 Correct 29 ms 5752 KB Output is correct
31 Correct 22 ms 4088 KB Output is correct
32 Correct 23 ms 3968 KB Output is correct
33 Correct 115 ms 4676 KB Output is correct
34 Correct 132 ms 4684 KB Output is correct
35 Correct 122 ms 4868 KB Output is correct
36 Correct 317 ms 43056 KB Output is correct
37 Correct 326 ms 44132 KB Output is correct
38 Correct 356 ms 48956 KB Output is correct
39 Correct 319 ms 8676 KB Output is correct
40 Correct 324 ms 8672 KB Output is correct
41 Correct 314 ms 8548 KB Output is correct
42 Correct 280 ms 8544 KB Output is correct
43 Correct 233 ms 26208 KB Output is correct
44 Correct 252 ms 26208 KB Output is correct
45 Correct 218 ms 26084 KB Output is correct
46 Correct 235 ms 26080 KB Output is correct
47 Correct 238 ms 26208 KB Output is correct
48 Correct 233 ms 26212 KB Output is correct
49 Correct 220 ms 26208 KB Output is correct
50 Correct 212 ms 26212 KB Output is correct

Subtask #6
0 / 13.0

# Verdict Execution time Memory Grader output
1 Correct 0 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 0 ms 384 KB Output is correct
7 Correct 1 ms 384 KB Output is correct
8 Correct 0 ms 384 KB Output is correct
9 Correct 0 ms 384 KB Output is correct
10 Correct 0 ms 384 KB Output is correct
11 Correct 0 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 0 ms 384 KB Output is correct
16 Correct 1 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 6 ms 640 KB Output is correct
20 Correct 6 ms 640 KB Output is correct
21 Correct 6 ms 640 KB Output is correct
22 Correct 11 ms 2304 KB Output is correct
23 Correct 17 ms 3456 KB Output is correct
24 Correct 10 ms 2048 KB Output is correct
25 Correct 9 ms 1664 KB Output is correct
26 Correct 10 ms 2048 KB Output is correct
27 Correct 347 ms 13248 KB Output is correct
28 Correct 353 ms 13260 KB Output is correct
29 Correct 349 ms 12992 KB Output is correct
30 Correct 333 ms 13260 KB Output is correct
31 Correct 507 ms 17864 KB Output is correct
32 Correct 742 ms 55832 KB Output is correct
33 Correct 620 ms 24264 KB Output is correct
34 Correct 653 ms 33868 KB Output is correct
35 Correct 0 ms 384 KB Output is correct
36 Correct 456 ms 60668 KB Output is correct
37 Correct 1507 ms 178368 KB Output is correct
38 Correct 1496 ms 177740 KB Output is correct
39 Correct 1865 ms 162116 KB Output is correct
40 Correct 754 ms 65716 KB Output is correct
41 Correct 292 ms 31096 KB Output is correct
42 Correct 38 ms 4744 KB Output is correct
43 Correct 1679 ms 180292 KB Output is correct
44 Correct 2378 ms 144376 KB Output is correct
45 Correct 1179 ms 164296 KB Output is correct
46 Correct 709 ms 78912 KB Output is correct
47 Correct 1710 ms 33356 KB Output is correct
48 Correct 759 ms 90572 KB Output is correct
49 Correct 749 ms 84284 KB Output is correct
50 Correct 1231 ms 17728 KB Output is correct
51 Correct 11 ms 896 KB Output is correct
52 Correct 11 ms 896 KB Output is correct
53 Correct 11 ms 896 KB Output is correct
54 Correct 29 ms 5752 KB Output is correct
55 Correct 22 ms 4088 KB Output is correct
56 Correct 23 ms 3968 KB Output is correct
57 Correct 115 ms 4676 KB Output is correct
58 Correct 132 ms 4684 KB Output is correct
59 Correct 122 ms 4868 KB Output is correct
60 Correct 317 ms 43056 KB Output is correct
61 Correct 326 ms 44132 KB Output is correct
62 Correct 356 ms 48956 KB Output is correct
63 Correct 319 ms 8676 KB Output is correct
64 Correct 324 ms 8672 KB Output is correct
65 Correct 314 ms 8548 KB Output is correct
66 Correct 280 ms 8544 KB Output is correct
67 Correct 233 ms 26208 KB Output is correct
68 Correct 252 ms 26208 KB Output is correct
69 Correct 218 ms 26084 KB Output is correct
70 Correct 235 ms 26080 KB Output is correct
71 Correct 238 ms 26208 KB Output is correct
72 Correct 233 ms 26212 KB Output is correct
73 Correct 220 ms 26208 KB Output is correct
74 Correct 212 ms 26212 KB Output is correct
75 Correct 397 ms 13040 KB Output is correct
76 Correct 380 ms 13188 KB Output is correct
77 Correct 370 ms 13256 KB Output is correct
78 Correct 370 ms 13004 KB Output is correct
79 Correct 590 ms 15040 KB Output is correct
80 Correct 363 ms 13004 KB Output is correct
81 Correct 1418 ms 170620 KB Output is correct
82 Correct 1314 ms 155340 KB Output is correct
83 Correct 1246 ms 155852 KB Output is correct
84 Correct 1869 ms 193200 KB Output is correct
85 Correct 1394 ms 165576 KB Output is correct
86 Correct 1254 ms 151956 KB Output is correct
87 Execution timed out 3053 ms 222804 KB Time limit exceeded
88 Halted 0 ms 0 KB -