oj.uz

Submission #100453

# Submission time Handle Problem Language Result Execution time Memory
100453 2019-03-11T10:40:45 Z JPN20 Construction of Highway (JOI18_construction) C++17
100 / 100
 1162 ms 75720 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), CC = 0;
		
		for (int i = 0; i < 20; i++) {
			M = (1 << i);
			int G = prevs(pos, M);
			int P = query(G);
			if (E != P) { R = (1 << i); CC = i + 2; maxn = min(maxn, (1 << i)); break; }
		}

		for (int i = 0; i < CC; 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)); rems++;
		pos = prevs(pos, maxn);
	}
	reverse(U.begin(), U.end());
	return U;
}

int main() {
	//FILE *in = freopen("in1.txt", "r", stdin);
	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:143:7: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
  scanf("%d", &N);
  ~~~~~^~~~~~~~~~
construction.cpp:144: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:146: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 7 ms 6656 KB Output is correct
4 Correct 8 ms 6784 KB Output is correct
5 Correct 9 ms 6912 KB Output is correct
6 Correct 8 ms 6912 KB Output is correct
7 Correct 9 ms 6912 KB Output is correct
8 Correct 9 ms 6912 KB Output is correct
9 Correct 8 ms 6912 KB Output is correct
10 Correct 8 ms 6912 KB Output is correct
11 Correct 10 ms 6912 KB Output is correct
12 Correct 8 ms 6912 KB Output is correct
13 Correct 8 ms 6848 KB Output is correct
14 Correct 9 ms 6912 KB Output is correct
15 Correct 9 ms 6912 KB Output is correct
16 Correct 9 ms 6912 KB Output is correct
17 Correct 10 ms 6912 KB Output is correct
18 Correct 9 ms 6912 KB Output is correct
19 Correct 8 ms 6912 KB Output is correct
20 Correct 9 ms 6912 KB Output is correct
21 Correct 11 ms 6912 KB Output is correct
22 Correct 11 ms 6912 KB Output is correct
23 Correct 12 ms 6912 KB Output is correct
24 Correct 9 ms 6912 KB Output is correct
25 Correct 9 ms 6912 KB Output is correct
26 Correct 9 ms 6912 KB Output is correct
27 Correct 10 ms 6912 KB Output is correct
28 Correct 9 ms 6912 KB Output is correct
29 Correct 10 ms 6868 KB Output is correct
30 Correct 9 ms 6912 KB Output is correct
31 Correct 9 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 11 ms 6912 KB Output is correct
35 Correct 9 ms 6812 KB Output is correct
36 Correct 10 ms 6912 KB Output is correct
37 Correct 12 ms 6912 KB Output is correct
38 Correct 11 ms 6912 KB Output is correct
39 Correct 10 ms 6912 KB Output is correct
40 Correct 9 ms 6912 KB Output is correct
41 Correct 9 ms 6912 KB Output is correct
42 Correct 9 ms 6912 KB Output is correct
43 Correct 9 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 7 ms 6656 KB Output is correct
4 Correct 8 ms 6784 KB Output is correct
5 Correct 9 ms 6912 KB Output is correct
6 Correct 8 ms 6912 KB Output is correct
7 Correct 9 ms 6912 KB Output is correct
8 Correct 9 ms 6912 KB Output is correct
9 Correct 8 ms 6912 KB Output is correct
10 Correct 8 ms 6912 KB Output is correct
11 Correct 10 ms 6912 KB Output is correct
12 Correct 8 ms 6912 KB Output is correct
13 Correct 8 ms 6848 KB Output is correct
14 Correct 9 ms 6912 KB Output is correct
15 Correct 9 ms 6912 KB Output is correct
16 Correct 9 ms 6912 KB Output is correct
17 Correct 10 ms 6912 KB Output is correct
18 Correct 9 ms 6912 KB Output is correct
19 Correct 8 ms 6912 KB Output is correct
20 Correct 9 ms 6912 KB Output is correct
21 Correct 11 ms 6912 KB Output is correct
22 Correct 11 ms 6912 KB Output is correct
23 Correct 12 ms 6912 KB Output is correct
24 Correct 9 ms 6912 KB Output is correct
25 Correct 9 ms 6912 KB Output is correct
26 Correct 9 ms 6912 KB Output is correct
27 Correct 10 ms 6912 KB Output is correct
28 Correct 9 ms 6912 KB Output is correct
29 Correct 10 ms 6868 KB Output is correct
30 Correct 9 ms 6912 KB Output is correct
31 Correct 9 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 11 ms 6912 KB Output is correct
35 Correct 9 ms 6812 KB Output is correct
36 Correct 10 ms 6912 KB Output is correct
37 Correct 12 ms 6912 KB Output is correct
38 Correct 11 ms 6912 KB Output is correct
39 Correct 10 ms 6912 KB Output is correct
40 Correct 9 ms 6912 KB Output is correct
41 Correct 9 ms 6912 KB Output is correct
42 Correct 9 ms 6912 KB Output is correct
43 Correct 9 ms 6912 KB Output is correct
44 Correct 10 ms 6912 KB Output is correct
45 Correct 11 ms 7168 KB Output is correct
46 Correct 24 ms 9088 KB Output is correct
47 Correct 23 ms 9060 KB Output is correct
48 Correct 22 ms 9088 KB Output is correct
49 Correct 23 ms 9260 KB Output is correct
50 Correct 18 ms 9216 KB Output is correct
51 Correct 18 ms 9216 KB Output is correct
52 Correct 20 ms 9208 KB Output is correct
53 Correct 20 ms 9208 KB Output is correct
54 Correct 20 ms 9216 KB Output is correct
55 Correct 21 ms 9216 KB Output is correct
56 Correct 23 ms 9088 KB Output is correct
57 Correct 39 ms 9088 KB Output is correct
58 Correct 31 ms 9088 KB Output is correct
59 Correct 27 ms 9088 KB Output is correct
60 Correct 27 ms 9088 KB Output is correct
61 Correct 31 ms 9216 KB Output is correct
62 Correct 22 ms 9216 KB Output is correct
63 Correct 23 ms 9208 KB Output is correct
64 Correct 23 ms 9088 KB Output is correct
65 Correct 23 ms 9088 KB Output is correct
66 Correct 22 ms 9088 KB Output is correct
67 Correct 22 ms 9088 KB Output is correct
68 Correct 18 ms 9216 KB Output is correct
69 Correct 19 ms 9216 KB Output is correct
70 Correct 17 ms 9216 KB Output is correct
71 Correct 20 ms 9188 KB Output is correct
72 Correct 31 ms 9088 KB Output is correct
73 Correct 27 ms 9088 KB Output is correct
74 Correct 21 ms 9216 KB Output is correct
75 Correct 19 ms 9088 KB Output is correct
76 Correct 19 ms 9088 KB Output is correct
77 Correct 27 ms 9088 KB Output is correct
78 Correct 24 ms 9088 KB Output is correct
79 Correct 19 ms 9088 KB Output is correct
80 Correct 20 ms 9088 KB Output is correct
81 Correct 23 ms 9088 KB Output is correct
82 Correct 25 ms 9088 KB Output is correct
83 Correct 31 ms 9088 KB Output is correct
84 Correct 24 ms 9136 KB Output is correct
85 Correct 28 ms 9088 KB Output is correct
86 Correct 21 ms 9088 KB Output is correct

Subtask #3
84.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 7 ms 6656 KB Output is correct
4 Correct 8 ms 6784 KB Output is correct
5 Correct 9 ms 6912 KB Output is correct
6 Correct 8 ms 6912 KB Output is correct
7 Correct 9 ms 6912 KB Output is correct
8 Correct 9 ms 6912 KB Output is correct
9 Correct 8 ms 6912 KB Output is correct
10 Correct 8 ms 6912 KB Output is correct
11 Correct 10 ms 6912 KB Output is correct
12 Correct 8 ms 6912 KB Output is correct
13 Correct 8 ms 6848 KB Output is correct
14 Correct 9 ms 6912 KB Output is correct
15 Correct 9 ms 6912 KB Output is correct
16 Correct 9 ms 6912 KB Output is correct
17 Correct 10 ms 6912 KB Output is correct
18 Correct 9 ms 6912 KB Output is correct
19 Correct 8 ms 6912 KB Output is correct
20 Correct 9 ms 6912 KB Output is correct
21 Correct 11 ms 6912 KB Output is correct
22 Correct 11 ms 6912 KB Output is correct
23 Correct 12 ms 6912 KB Output is correct
24 Correct 9 ms 6912 KB Output is correct
25 Correct 9 ms 6912 KB Output is correct
26 Correct 9 ms 6912 KB Output is correct
27 Correct 10 ms 6912 KB Output is correct
28 Correct 9 ms 6912 KB Output is correct
29 Correct 10 ms 6868 KB Output is correct
30 Correct 9 ms 6912 KB Output is correct
31 Correct 9 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 11 ms 6912 KB Output is correct
35 Correct 9 ms 6812 KB Output is correct
36 Correct 10 ms 6912 KB Output is correct
37 Correct 12 ms 6912 KB Output is correct
38 Correct 11 ms 6912 KB Output is correct
39 Correct 10 ms 6912 KB Output is correct
40 Correct 9 ms 6912 KB Output is correct
41 Correct 9 ms 6912 KB Output is correct
42 Correct 9 ms 6912 KB Output is correct
43 Correct 9 ms 6912 KB Output is correct
44 Correct 10 ms 6912 KB Output is correct
45 Correct 11 ms 7168 KB Output is correct
46 Correct 24 ms 9088 KB Output is correct
47 Correct 23 ms 9060 KB Output is correct
48 Correct 22 ms 9088 KB Output is correct
49 Correct 23 ms 9260 KB Output is correct
50 Correct 18 ms 9216 KB Output is correct
51 Correct 18 ms 9216 KB Output is correct
52 Correct 20 ms 9208 KB Output is correct
53 Correct 20 ms 9208 KB Output is correct
54 Correct 20 ms 9216 KB Output is correct
55 Correct 21 ms 9216 KB Output is correct
56 Correct 23 ms 9088 KB Output is correct
57 Correct 39 ms 9088 KB Output is correct
58 Correct 31 ms 9088 KB Output is correct
59 Correct 27 ms 9088 KB Output is correct
60 Correct 27 ms 9088 KB Output is correct
61 Correct 31 ms 9216 KB Output is correct
62 Correct 22 ms 9216 KB Output is correct
63 Correct 23 ms 9208 KB Output is correct
64 Correct 23 ms 9088 KB Output is correct
65 Correct 23 ms 9088 KB Output is correct
66 Correct 22 ms 9088 KB Output is correct
67 Correct 22 ms 9088 KB Output is correct
68 Correct 18 ms 9216 KB Output is correct
69 Correct 19 ms 9216 KB Output is correct
70 Correct 17 ms 9216 KB Output is correct
71 Correct 20 ms 9188 KB Output is correct
72 Correct 31 ms 9088 KB Output is correct
73 Correct 27 ms 9088 KB Output is correct
74 Correct 21 ms 9216 KB Output is correct
75 Correct 19 ms 9088 KB Output is correct
76 Correct 19 ms 9088 KB Output is correct
77 Correct 27 ms 9088 KB Output is correct
78 Correct 24 ms 9088 KB Output is correct
79 Correct 19 ms 9088 KB Output is correct
80 Correct 20 ms 9088 KB Output is correct
81 Correct 23 ms 9088 KB Output is correct
82 Correct 25 ms 9088 KB Output is correct
83 Correct 31 ms 9088 KB Output is correct
84 Correct 24 ms 9136 KB Output is correct
85 Correct 28 ms 9088 KB Output is correct
86 Correct 21 ms 9088 KB Output is correct
87 Correct 54 ms 12972 KB Output is correct
88 Correct 194 ms 25516 KB Output is correct
89 Correct 780 ms 69836 KB Output is correct
90 Correct 755 ms 69780 KB Output is correct
91 Correct 765 ms 69764 KB Output is correct
92 Correct 512 ms 74232 KB Output is correct
93 Correct 528 ms 74232 KB Output is correct
94 Correct 547 ms 74360 KB Output is correct
95 Correct 482 ms 73076 KB Output is correct
96 Correct 511 ms 73092 KB Output is correct
97 Correct 553 ms 73244 KB Output is correct
98 Correct 504 ms 73204 KB Output is correct
99 Correct 486 ms 72056 KB Output is correct
100 Correct 1034 ms 69880 KB Output is correct
101 Correct 1162 ms 71800 KB Output is correct
102 Correct 1157 ms 71824 KB Output is correct
103 Correct 1112 ms 71800 KB Output is correct
104 Correct 883 ms 74248 KB Output is correct
105 Correct 896 ms 74104 KB Output is correct
106 Correct 873 ms 74220 KB Output is correct
107 Correct 625 ms 71036 KB Output is correct
108 Correct 696 ms 71224 KB Output is correct
109 Correct 830 ms 71524 KB Output is correct
110 Correct 520 ms 75720 KB Output is correct
111 Correct 478 ms 74996 KB Output is correct
112 Correct 457 ms 74712 KB Output is correct
113 Correct 437 ms 73592 KB Output is correct
114 Correct 940 ms 71544 KB Output is correct
115 Correct 1071 ms 71312 KB Output is correct
116 Correct 900 ms 73720 KB Output is correct
117 Correct 428 ms 72696 KB Output is correct
118 Correct 454 ms 72184 KB Output is correct
119 Correct 403 ms 72184 KB Output is correct
120 Correct 495 ms 72184 KB Output is correct
121 Correct 406 ms 72032 KB Output is correct
122 Correct 389 ms 71808 KB Output is correct
123 Correct 887 ms 72772 KB Output is correct
124 Correct 882 ms 72440 KB Output is correct
125 Correct 860 ms 72264 KB Output is correct
126 Correct 892 ms 72440 KB Output is correct
127 Correct 814 ms 72024 KB Output is correct
128 Correct 816 ms 71848 KB Output is correct