oj.uz

Submission #100448

# Submission time Handle Problem Language Result Execution time Memory
100448 2019-03-11T10:27:38 Z JPN20 Construction of Highway (JOI18_construction) C++17
16 / 100
 1151 ms 138736 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);

		rems++;
		assert(rems <= 500000);
	}
	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:138:7: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
  scanf("%d", &N);
  ~~~~~^~~~~~~~~~
construction.cpp:139: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:141: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 7 ms 6528 KB Output is correct
2 Correct 8 ms 6528 KB Output is correct
3 Correct 8 ms 6656 KB Output is correct
4 Correct 9 ms 6784 KB Output is correct
5 Correct 13 ms 6912 KB Output is correct
6 Correct 10 ms 6912 KB Output is correct
7 Correct 13 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 9 ms 6912 KB Output is correct
14 Correct 8 ms 6912 KB Output is correct
15 Correct 10 ms 6912 KB Output is correct
16 Correct 12 ms 6912 KB Output is correct
17 Correct 16 ms 6912 KB Output is correct
18 Correct 12 ms 6912 KB Output is correct
19 Correct 10 ms 6912 KB Output is correct
20 Correct 10 ms 6912 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 10 ms 6912 KB Output is correct
26 Correct 11 ms 6844 KB Output is correct
27 Correct 12 ms 6912 KB Output is correct
28 Correct 10 ms 6912 KB Output is correct
29 Correct 10 ms 6912 KB Output is correct
30 Correct 12 ms 6912 KB Output is correct
31 Correct 15 ms 6912 KB Output is correct
32 Correct 12 ms 6912 KB Output is correct
33 Correct 9 ms 6912 KB Output is correct
34 Correct 10 ms 6988 KB Output is correct
35 Correct 9 ms 6912 KB Output is correct
36 Correct 10 ms 6912 KB Output is correct
37 Correct 9 ms 6912 KB Output is correct
38 Correct 9 ms 6912 KB Output is correct
39 Correct 10 ms 7040 KB Output is correct
40 Correct 9 ms 6912 KB Output is correct
41 Correct 10 ms 6912 KB Output is correct
42 Correct 9 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 7 ms 6528 KB Output is correct
2 Correct 8 ms 6528 KB Output is correct
3 Correct 8 ms 6656 KB Output is correct
4 Correct 9 ms 6784 KB Output is correct
5 Correct 13 ms 6912 KB Output is correct
6 Correct 10 ms 6912 KB Output is correct
7 Correct 13 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 9 ms 6912 KB Output is correct
14 Correct 8 ms 6912 KB Output is correct
15 Correct 10 ms 6912 KB Output is correct
16 Correct 12 ms 6912 KB Output is correct
17 Correct 16 ms 6912 KB Output is correct
18 Correct 12 ms 6912 KB Output is correct
19 Correct 10 ms 6912 KB Output is correct
20 Correct 10 ms 6912 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 10 ms 6912 KB Output is correct
26 Correct 11 ms 6844 KB Output is correct
27 Correct 12 ms 6912 KB Output is correct
28 Correct 10 ms 6912 KB Output is correct
29 Correct 10 ms 6912 KB Output is correct
30 Correct 12 ms 6912 KB Output is correct
31 Correct 15 ms 6912 KB Output is correct
32 Correct 12 ms 6912 KB Output is correct
33 Correct 9 ms 6912 KB Output is correct
34 Correct 10 ms 6988 KB Output is correct
35 Correct 9 ms 6912 KB Output is correct
36 Correct 10 ms 6912 KB Output is correct
37 Correct 9 ms 6912 KB Output is correct
38 Correct 9 ms 6912 KB Output is correct
39 Correct 10 ms 7040 KB Output is correct
40 Correct 9 ms 6912 KB Output is correct
41 Correct 10 ms 6912 KB Output is correct
42 Correct 9 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 13 ms 7168 KB Output is correct
46 Correct 31 ms 9088 KB Output is correct
47 Correct 31 ms 9088 KB Output is correct
48 Correct 42 ms 9204 KB Output is correct
49 Correct 21 ms 9216 KB Output is correct
50 Correct 20 ms 9216 KB Output is correct
51 Correct 21 ms 9216 KB Output is correct
52 Correct 22 ms 9188 KB Output is correct
53 Correct 25 ms 9216 KB Output is correct
54 Correct 32 ms 9216 KB Output is correct
55 Correct 33 ms 9180 KB Output is correct
56 Correct 24 ms 9216 KB Output is correct
57 Correct 55 ms 9088 KB Output is correct
58 Correct 47 ms 9088 KB Output is correct
59 Correct 55 ms 9088 KB Output is correct
60 Correct 72 ms 9088 KB Output is correct
61 Correct 22 ms 9088 KB Output is correct
62 Correct 25 ms 9088 KB Output is correct
63 Correct 24 ms 9088 KB Output is correct
64 Correct 30 ms 9088 KB Output is correct
65 Correct 33 ms 9088 KB Output is correct
66 Correct 31 ms 9088 KB Output is correct
67 Correct 31 ms 9088 KB Output is correct
68 Correct 23 ms 9216 KB Output is correct
69 Correct 22 ms 9216 KB Output is correct
70 Correct 24 ms 9084 KB Output is correct
71 Correct 23 ms 9216 KB Output is correct
72 Correct 53 ms 9088 KB Output is correct
73 Correct 47 ms 9088 KB Output is correct
74 Correct 23 ms 9088 KB Output is correct
75 Correct 22 ms 9088 KB Output is correct
76 Correct 21 ms 9088 KB Output is correct
77 Correct 23 ms 9088 KB Output is correct
78 Correct 20 ms 9088 KB Output is correct
79 Correct 21 ms 9060 KB Output is correct
80 Correct 29 ms 9088 KB Output is correct
81 Correct 25 ms 9088 KB Output is correct
82 Correct 24 ms 9088 KB Output is correct
83 Correct 26 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 27 ms 9088 KB Output is correct

Subtask #3
0 / 84.0

# Verdict Execution time Memory Grader output
1 Correct 7 ms 6528 KB Output is correct
2 Correct 8 ms 6528 KB Output is correct
3 Correct 8 ms 6656 KB Output is correct
4 Correct 9 ms 6784 KB Output is correct
5 Correct 13 ms 6912 KB Output is correct
6 Correct 10 ms 6912 KB Output is correct
7 Correct 13 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 9 ms 6912 KB Output is correct
14 Correct 8 ms 6912 KB Output is correct
15 Correct 10 ms 6912 KB Output is correct
16 Correct 12 ms 6912 KB Output is correct
17 Correct 16 ms 6912 KB Output is correct
18 Correct 12 ms 6912 KB Output is correct
19 Correct 10 ms 6912 KB Output is correct
20 Correct 10 ms 6912 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 10 ms 6912 KB Output is correct
26 Correct 11 ms 6844 KB Output is correct
27 Correct 12 ms 6912 KB Output is correct
28 Correct 10 ms 6912 KB Output is correct
29 Correct 10 ms 6912 KB Output is correct
30 Correct 12 ms 6912 KB Output is correct
31 Correct 15 ms 6912 KB Output is correct
32 Correct 12 ms 6912 KB Output is correct
33 Correct 9 ms 6912 KB Output is correct
34 Correct 10 ms 6988 KB Output is correct
35 Correct 9 ms 6912 KB Output is correct
36 Correct 10 ms 6912 KB Output is correct
37 Correct 9 ms 6912 KB Output is correct
38 Correct 9 ms 6912 KB Output is correct
39 Correct 10 ms 7040 KB Output is correct
40 Correct 9 ms 6912 KB Output is correct
41 Correct 10 ms 6912 KB Output is correct
42 Correct 9 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 13 ms 7168 KB Output is correct
46 Correct 31 ms 9088 KB Output is correct
47 Correct 31 ms 9088 KB Output is correct
48 Correct 42 ms 9204 KB Output is correct
49 Correct 21 ms 9216 KB Output is correct
50 Correct 20 ms 9216 KB Output is correct
51 Correct 21 ms 9216 KB Output is correct
52 Correct 22 ms 9188 KB Output is correct
53 Correct 25 ms 9216 KB Output is correct
54 Correct 32 ms 9216 KB Output is correct
55 Correct 33 ms 9180 KB Output is correct
56 Correct 24 ms 9216 KB Output is correct
57 Correct 55 ms 9088 KB Output is correct
58 Correct 47 ms 9088 KB Output is correct
59 Correct 55 ms 9088 KB Output is correct
60 Correct 72 ms 9088 KB Output is correct
61 Correct 22 ms 9088 KB Output is correct
62 Correct 25 ms 9088 KB Output is correct
63 Correct 24 ms 9088 KB Output is correct
64 Correct 30 ms 9088 KB Output is correct
65 Correct 33 ms 9088 KB Output is correct
66 Correct 31 ms 9088 KB Output is correct
67 Correct 31 ms 9088 KB Output is correct
68 Correct 23 ms 9216 KB Output is correct
69 Correct 22 ms 9216 KB Output is correct
70 Correct 24 ms 9084 KB Output is correct
71 Correct 23 ms 9216 KB Output is correct
72 Correct 53 ms 9088 KB Output is correct
73 Correct 47 ms 9088 KB Output is correct
74 Correct 23 ms 9088 KB Output is correct
75 Correct 22 ms 9088 KB Output is correct
76 Correct 21 ms 9088 KB Output is correct
77 Correct 23 ms 9088 KB Output is correct
78 Correct 20 ms 9088 KB Output is correct
79 Correct 21 ms 9060 KB Output is correct
80 Correct 29 ms 9088 KB Output is correct
81 Correct 25 ms 9088 KB Output is correct
82 Correct 24 ms 9088 KB Output is correct
83 Correct 26 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 27 ms 9088 KB Output is correct
87 Correct 80 ms 12980 KB Output is correct
88 Correct 292 ms 25472 KB Output is correct
89 Correct 1151 ms 69852 KB Output is correct
90 Correct 1128 ms 69736 KB Output is correct
91 Correct 1059 ms 69708 KB Output is correct
92 Correct 662 ms 74208 KB Output is correct
93 Correct 683 ms 74196 KB Output is correct
94 Correct 677 ms 74296 KB Output is correct
95 Correct 786 ms 73196 KB Output is correct
96 Correct 777 ms 73152 KB Output is correct
97 Correct 871 ms 73276 KB Output is correct
98 Correct 833 ms 73220 KB Output is correct
99 Correct 780 ms 72160 KB Output is correct
100 Runtime error 981 ms 138736 KB Execution killed with signal 11 (could be triggered by violating memory limits)
101 Halted 0 ms 0 KB -