oj.uz

Submission #100449

# Submission time Handle Problem Language Result Execution time Memory
100449 2019-03-11T10:28:44 Z JPN20 Construction of Highway (JOI18_construction) C++17
16 / 100
 2000 ms 74276 KB 
#include <bits/stdc++.h>
using namespace std;

int size_ = 1, rems = 1;
int va[1 << 18], vb[1 << 18], dat[1 << 18];
int V[1 << 17][118], VS[1 << 18];

void inits(int sz) {
	while (size_ <= sz) size_ *= 2;

	int cnt = 1;
	for (int i = size_; i >= 1; i /= 2) {
		for (int j = 0; j < size_; j += i) {
			va[cnt] = j;
			vb[cnt] = j + i;
			cnt++;
		}
	}
}
void update(int pos, int x) {
	pos += size_;
	dat[pos] = x;

	while (pos >= 2) {
		pos >>= 1;
		dat[pos] = max(dat[pos * 2], dat[pos * 2 + 1]);
	}
}

int query(int pos) {
	int ret = -(1 << 30);
	for (int i = 0; i < VS[pos]; i++) ret = max(ret, dat[V[pos][i]]);
	return ret;
}

void query_init_(int pos, int l, int r, int a, int b, int u) {
	if (l <= a && b <= r) { V[pos][VS[pos]] = u; VS[pos]++; return; }
	if (b <= l || r <= a) return;
	
	query_init_(pos, l, r, a, (a + b) >> 1, u * 2);
	query_init_(pos, l, r, (a + b) >> 1, b, u * 2 + 1);
}

void query_init(int pos, int l, int r) {
	query_init_(pos, l, r, 0, size_, 1);
}

class BIT {
public:
	vector<int>bit; int size1_;

	void init(int sz) {
		size1_ = sz + 2;
		bit.resize(size1_ + 1, 0);
	}
	void add(int pos, int x) {
		pos++;
		while (pos <= size1_) {
			bit[pos] += x; pos += (pos & -pos);
		}
	}
	int sum(int pos) {
		int s = 0; pos++;
		while (pos >= 1) {
			s += bit[pos]; pos -= (pos & -pos);
		}
		return s;
	}
};

long long check_inversions(vector<pair<int, int>> vec) {
	vector<int> A;
	for (int i = 0; i < (int)vec.size(); i++) A.push_back(vec[i].first);

	sort(A.begin(), A.end());
	A.erase(unique(A.begin(), A.end()), A.end());
	for (int i = 0; i < (int)vec.size(); i++) vec[i].first = lower_bound(A.begin(), A.end(), vec[i].first) - A.begin();

	BIT Z; Z.init(A.size() + 1);

	long long ans = 0;
	for (int i = 0; i < (int)vec.size(); i++) {
		ans += 1LL * vec[i].second * (Z.sum(A.size()) - Z.sum(vec[i].first));
		Z.add(vec[i].first, vec[i].second);
	}
	return ans;
}

int N, C[1 << 18], A[1 << 18], B[1 << 18], dist[1 << 18], dp[1 << 18][20];
int cl[1 << 18], cr[1 << 18], cnts;
vector<int> X[1 << 18];

void dfs(int pos, int dep) {
	cnts++; cl[pos] = cnts; dist[pos] = dep;
	for (int i = 0; i < (int)X[pos].size(); i++) {
		if (dist[X[pos][i]] != -1) continue;
		dp[X[pos][i]][0] = pos;
		dfs(X[pos][i], dep + 1);
	}
	cr[pos] = cnts;
}

int prevs(int pos, int x) {
	for (int i = 19; i >= 0; i--) {
		if (x >= (1 << i)) { x -= (1 << i); pos = dp[pos][i]; }
	}
	if (pos == 0) pos = 1;
	return pos;
}

vector<pair<int, int>> get_inverse(int pos) {
	vector<pair<int, int>> U;
	while (true) {
		// Binary Search.
		int L = 0, R = dist[pos] + 1, M, maxn = (1 << 30);
		int E = query(pos);

		for (int i = 0; i < 20; i++) {
			M = (L + R) / 2;
			int G = prevs(pos, M);
			int P = query(G);
			if (E != P) { maxn = min(maxn, M); R = M; }
			else { L = M; }
		}

		if (maxn == (1 << 30)) { U.push_back(make_pair(C[B[E]], dist[pos] + 1)); break; }
		U.push_back(make_pair(C[B[E]], maxn));
		pos = prevs(pos, maxn);
	}
	reverse(U.begin(), U.end());
	return U;
}

int main() {
	scanf("%d", &N);
	for (int i = 1; i <= N; i++) scanf("%d", &C[i]);
	for (int i = 1; i <= N - 1; i++) {
		scanf("%d%d", &A[i], &B[i]);
		X[A[i]].push_back(B[i]);
		X[B[i]].push_back(A[i]);
	}

	for (int i = 1; i <= N; i++) dist[i] = -1;
	dfs(1, 0);

	for (int i = 0; i < 19; i++) {
		for (int j = 1; j <= N; j++) dp[j][i + 1] = dp[dp[j][i]][i];
	}

	inits(N + 2);
	for (int i = 1; i <= N; i++) query_init(i, cl[i], cr[i] + 1);
	
	update(cl[1], 0);
	for (int i = 1; i <= N - 1; i++) {
		vector<pair<int, int>> T = get_inverse(A[i]);

		printf("%lld\n", check_inversions(T));
		update(cl[B[i]], i);
	}
	return 0;
}

Compilation message

construction.cpp: In function 'int main()':
construction.cpp:135:7: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
  scanf("%d", &N);
  ~~~~~^~~~~~~~~~
construction.cpp:136:36: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
  for (int i = 1; i <= N; i++) scanf("%d", &C[i]);
                               ~~~~~^~~~~~~~~~~~~
construction.cpp:138:8: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
   scanf("%d%d", &A[i], &B[i]);
   ~~~~~^~~~~~~~~~~~~~~~~~~~~~

Subtask #1
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 9 ms 6528 KB Output is correct
2 Correct 7 ms 6528 KB Output is correct
3 Correct 8 ms 6572 KB Output is correct
4 Correct 10 ms 6784 KB Output is correct
5 Correct 10 ms 6912 KB Output is correct
6 Correct 12 ms 6912 KB Output is correct
7 Correct 9 ms 6912 KB Output is correct
8 Correct 12 ms 6912 KB Output is correct
9 Correct 10 ms 6912 KB Output is correct
10 Correct 9 ms 6912 KB Output is correct
11 Correct 9 ms 6912 KB Output is correct
12 Correct 9 ms 6912 KB Output is correct
13 Correct 10 ms 6912 KB Output is correct
14 Correct 9 ms 6912 KB Output is correct
15 Correct 11 ms 6884 KB Output is correct
16 Correct 11 ms 6912 KB Output is correct
17 Correct 11 ms 6912 KB Output is correct
18 Correct 11 ms 6912 KB Output is correct
19 Correct 9 ms 6912 KB Output is correct
20 Correct 11 ms 6972 KB Output is correct
21 Correct 9 ms 6912 KB Output is correct
22 Correct 10 ms 6912 KB Output is correct
23 Correct 10 ms 6912 KB Output is correct
24 Correct 10 ms 6912 KB Output is correct
25 Correct 11 ms 6912 KB Output is correct
26 Correct 10 ms 6948 KB Output is correct
27 Correct 10 ms 6936 KB Output is correct
28 Correct 6 ms 6912 KB Output is correct
29 Correct 9 ms 6912 KB Output is correct
30 Correct 13 ms 6912 KB Output is correct
31 Correct 12 ms 6912 KB Output is correct
32 Correct 9 ms 6912 KB Output is correct
33 Correct 9 ms 6912 KB Output is correct
34 Correct 9 ms 6912 KB Output is correct
35 Correct 12 ms 6912 KB Output is correct
36 Correct 10 ms 6908 KB Output is correct
37 Correct 9 ms 6936 KB Output is correct
38 Correct 11 ms 6912 KB Output is correct
39 Correct 10 ms 6912 KB Output is correct
40 Correct 10 ms 6912 KB Output is correct
41 Correct 9 ms 6912 KB Output is correct
42 Correct 10 ms 6912 KB Output is correct
43 Correct 10 ms 6912 KB Output is correct
44 Correct 10 ms 6912 KB Output is correct

Subtask #2
9.0 / 9.0

# Verdict Execution time Memory Grader output
1 Correct 9 ms 6528 KB Output is correct
2 Correct 7 ms 6528 KB Output is correct
3 Correct 8 ms 6572 KB Output is correct
4 Correct 10 ms 6784 KB Output is correct
5 Correct 10 ms 6912 KB Output is correct
6 Correct 12 ms 6912 KB Output is correct
7 Correct 9 ms 6912 KB Output is correct
8 Correct 12 ms 6912 KB Output is correct
9 Correct 10 ms 6912 KB Output is correct
10 Correct 9 ms 6912 KB Output is correct
11 Correct 9 ms 6912 KB Output is correct
12 Correct 9 ms 6912 KB Output is correct
13 Correct 10 ms 6912 KB Output is correct
14 Correct 9 ms 6912 KB Output is correct
15 Correct 11 ms 6884 KB Output is correct
16 Correct 11 ms 6912 KB Output is correct
17 Correct 11 ms 6912 KB Output is correct
18 Correct 11 ms 6912 KB Output is correct
19 Correct 9 ms 6912 KB Output is correct
20 Correct 11 ms 6972 KB Output is correct
21 Correct 9 ms 6912 KB Output is correct
22 Correct 10 ms 6912 KB Output is correct
23 Correct 10 ms 6912 KB Output is correct
24 Correct 10 ms 6912 KB Output is correct
25 Correct 11 ms 6912 KB Output is correct
26 Correct 10 ms 6948 KB Output is correct
27 Correct 10 ms 6936 KB Output is correct
28 Correct 6 ms 6912 KB Output is correct
29 Correct 9 ms 6912 KB Output is correct
30 Correct 13 ms 6912 KB Output is correct
31 Correct 12 ms 6912 KB Output is correct
32 Correct 9 ms 6912 KB Output is correct
33 Correct 9 ms 6912 KB Output is correct
34 Correct 9 ms 6912 KB Output is correct
35 Correct 12 ms 6912 KB Output is correct
36 Correct 10 ms 6908 KB Output is correct
37 Correct 9 ms 6936 KB Output is correct
38 Correct 11 ms 6912 KB Output is correct
39 Correct 10 ms 6912 KB Output is correct
40 Correct 10 ms 6912 KB Output is correct
41 Correct 9 ms 6912 KB Output is correct
42 Correct 10 ms 6912 KB Output is correct
43 Correct 10 ms 6912 KB Output is correct
44 Correct 10 ms 6912 KB Output is correct
45 Correct 14 ms 7220 KB Output is correct
46 Correct 32 ms 9088 KB Output is correct
47 Correct 34 ms 9224 KB Output is correct
48 Correct 29 ms 9088 KB Output is correct
49 Correct 20 ms 9264 KB Output is correct
50 Correct 23 ms 9216 KB Output is correct
51 Correct 26 ms 9216 KB Output is correct
52 Correct 32 ms 9244 KB Output is correct
53 Correct 25 ms 9224 KB Output is correct
54 Correct 25 ms 9228 KB Output is correct
55 Correct 26 ms 9208 KB Output is correct
56 Correct 33 ms 9088 KB Output is correct
57 Correct 66 ms 9088 KB Output is correct
58 Correct 59 ms 9088 KB Output is correct
59 Correct 82 ms 9208 KB Output is correct
60 Correct 65 ms 9208 KB Output is correct
61 Correct 23 ms 9088 KB Output is correct
62 Correct 27 ms 9216 KB Output is correct
63 Correct 32 ms 9088 KB Output is correct
64 Correct 35 ms 9208 KB Output is correct
65 Correct 41 ms 9208 KB Output is correct
66 Correct 33 ms 9088 KB Output is correct
67 Correct 31 ms 9088 KB Output is correct
68 Correct 28 ms 9216 KB Output is correct
69 Correct 25 ms 9216 KB Output is correct
70 Correct 26 ms 9208 KB Output is correct
71 Correct 24 ms 9216 KB Output is correct
72 Correct 51 ms 9088 KB Output is correct
73 Correct 78 ms 9084 KB Output is correct
74 Correct 27 ms 9216 KB Output is correct
75 Correct 20 ms 9088 KB Output is correct
76 Correct 21 ms 9108 KB Output is correct
77 Correct 22 ms 9088 KB Output is correct
78 Correct 24 ms 9132 KB Output is correct
79 Correct 20 ms 9088 KB Output is correct
80 Correct 26 ms 9088 KB Output is correct
81 Correct 23 ms 9088 KB Output is correct
82 Correct 24 ms 9088 KB Output is correct
83 Correct 35 ms 9088 KB Output is correct
84 Correct 25 ms 9088 KB Output is correct
85 Correct 24 ms 9088 KB Output is correct
86 Correct 25 ms 9088 KB Output is correct

Subtask #3
0 / 84.0

# Verdict Execution time Memory Grader output
1 Correct 9 ms 6528 KB Output is correct
2 Correct 7 ms 6528 KB Output is correct
3 Correct 8 ms 6572 KB Output is correct
4 Correct 10 ms 6784 KB Output is correct
5 Correct 10 ms 6912 KB Output is correct
6 Correct 12 ms 6912 KB Output is correct
7 Correct 9 ms 6912 KB Output is correct
8 Correct 12 ms 6912 KB Output is correct
9 Correct 10 ms 6912 KB Output is correct
10 Correct 9 ms 6912 KB Output is correct
11 Correct 9 ms 6912 KB Output is correct
12 Correct 9 ms 6912 KB Output is correct
13 Correct 10 ms 6912 KB Output is correct
14 Correct 9 ms 6912 KB Output is correct
15 Correct 11 ms 6884 KB Output is correct
16 Correct 11 ms 6912 KB Output is correct
17 Correct 11 ms 6912 KB Output is correct
18 Correct 11 ms 6912 KB Output is correct
19 Correct 9 ms 6912 KB Output is correct
20 Correct 11 ms 6972 KB Output is correct
21 Correct 9 ms 6912 KB Output is correct
22 Correct 10 ms 6912 KB Output is correct
23 Correct 10 ms 6912 KB Output is correct
24 Correct 10 ms 6912 KB Output is correct
25 Correct 11 ms 6912 KB Output is correct
26 Correct 10 ms 6948 KB Output is correct
27 Correct 10 ms 6936 KB Output is correct
28 Correct 6 ms 6912 KB Output is correct
29 Correct 9 ms 6912 KB Output is correct
30 Correct 13 ms 6912 KB Output is correct
31 Correct 12 ms 6912 KB Output is correct
32 Correct 9 ms 6912 KB Output is correct
33 Correct 9 ms 6912 KB Output is correct
34 Correct 9 ms 6912 KB Output is correct
35 Correct 12 ms 6912 KB Output is correct
36 Correct 10 ms 6908 KB Output is correct
37 Correct 9 ms 6936 KB Output is correct
38 Correct 11 ms 6912 KB Output is correct
39 Correct 10 ms 6912 KB Output is correct
40 Correct 10 ms 6912 KB Output is correct
41 Correct 9 ms 6912 KB Output is correct
42 Correct 10 ms 6912 KB Output is correct
43 Correct 10 ms 6912 KB Output is correct
44 Correct 10 ms 6912 KB Output is correct
45 Correct 14 ms 7220 KB Output is correct
46 Correct 32 ms 9088 KB Output is correct
47 Correct 34 ms 9224 KB Output is correct
48 Correct 29 ms 9088 KB Output is correct
49 Correct 20 ms 9264 KB Output is correct
50 Correct 23 ms 9216 KB Output is correct
51 Correct 26 ms 9216 KB Output is correct
52 Correct 32 ms 9244 KB Output is correct
53 Correct 25 ms 9224 KB Output is correct
54 Correct 25 ms 9228 KB Output is correct
55 Correct 26 ms 9208 KB Output is correct
56 Correct 33 ms 9088 KB Output is correct
57 Correct 66 ms 9088 KB Output is correct
58 Correct 59 ms 9088 KB Output is correct
59 Correct 82 ms 9208 KB Output is correct
60 Correct 65 ms 9208 KB Output is correct
61 Correct 23 ms 9088 KB Output is correct
62 Correct 27 ms 9216 KB Output is correct
63 Correct 32 ms 9088 KB Output is correct
64 Correct 35 ms 9208 KB Output is correct
65 Correct 41 ms 9208 KB Output is correct
66 Correct 33 ms 9088 KB Output is correct
67 Correct 31 ms 9088 KB Output is correct
68 Correct 28 ms 9216 KB Output is correct
69 Correct 25 ms 9216 KB Output is correct
70 Correct 26 ms 9208 KB Output is correct
71 Correct 24 ms 9216 KB Output is correct
72 Correct 51 ms 9088 KB Output is correct
73 Correct 78 ms 9084 KB Output is correct
74 Correct 27 ms 9216 KB Output is correct
75 Correct 20 ms 9088 KB Output is correct
76 Correct 21 ms 9108 KB Output is correct
77 Correct 22 ms 9088 KB Output is correct
78 Correct 24 ms 9132 KB Output is correct
79 Correct 20 ms 9088 KB Output is correct
80 Correct 26 ms 9088 KB Output is correct
81 Correct 23 ms 9088 KB Output is correct
82 Correct 24 ms 9088 KB Output is correct
83 Correct 35 ms 9088 KB Output is correct
84 Correct 25 ms 9088 KB Output is correct
85 Correct 24 ms 9088 KB Output is correct
86 Correct 25 ms 9088 KB Output is correct
87 Correct 98 ms 13072 KB Output is correct
88 Correct 236 ms 25524 KB Output is correct
89 Correct 1104 ms 69820 KB Output is correct
90 Correct 1215 ms 69740 KB Output is correct
91 Correct 1162 ms 69752 KB Output is correct
92 Correct 637 ms 74276 KB Output is correct
93 Correct 600 ms 74264 KB Output is correct
94 Correct 602 ms 74232 KB Output is correct
95 Correct 679 ms 73196 KB Output is correct
96 Correct 690 ms 73196 KB Output is correct
97 Correct 709 ms 73076 KB Output is correct
98 Correct 753 ms 73176 KB Output is correct
99 Correct 841 ms 72184 KB Output is correct
100 Execution timed out 2053 ms 69936 KB Time limit exceeded
101 Halted 0 ms 0 KB -