Submission #77886

# Submission time Handle Problem Language Result Execution time Memory
77886 2018-10-01T04:09:42 Z gs14004 Construction of Highway (JOI18_construction) C++17
100 / 100
1429 ms 94672 KB
#include <bits/stdc++.h>
using namespace std;
const int MAXN = 100005;
typedef long long lint;
typedef pair<int, int> pi;
 
vector<int> gph[MAXN];
int n, s[MAXN], e[MAXN], a[MAXN];
int par[17][MAXN], dep[MAXN], din[MAXN], dout[MAXN], piv;
 
int lca(int s, int e){
	auto in = [&](int s, int e){
		return din[s] <= din[e] && dout[e] <= dout[s];
	};
	if(in(s, e)) return s;
	for(int i=16; i>=0; i--){
		if(par[i][s] && !in(par[i][s], e)) s = par[i][s];
	}
	return par[0][s];
}
 
void dfs(int x, int p){
	din[x] = ++piv;
	for(auto &i : gph[x]){
		if(i != p){
			dep[i] = dep[x] + 1;
			par[0][i] = x;
			dfs(i, x);
		}
	}
	dout[x] = piv;
}
 
struct segtree{
	int tree[270000], lim;
	void init(){
		memset(tree, 0, sizeof(tree));
		for(lim = 1; lim <= n; lim <<= 1);
	}
	void upd(int x, int v){
		x += lim;
		tree[x] = v;
		while(x > 1){
			x >>= 1;
			tree[x] = max(tree[2*x], tree[2*x+1]);
		}
	}
	int query(int s, int e){
		s += lim;
		e += lim;
		int ret = 0;
		while(s < e){
			if(s%2 == 1) ret = max(ret, tree[s++]);
			if(e%2 == 0) ret = max(ret, tree[e--]);
			s >>= 1;
			e >>= 1;
		}
		if(s == e) ret = max(ret, tree[s]);
		return ret;
	}
	int left_max(int s, int e, int thres){
		s += lim;
		e += lim;
		while(s < e){
			if(s % 2 == 1) s++;
			if(e % 2 == 0){
				if(tree[e] > thres){
					while(e < lim){
						if(tree[2*e+1] > thres) e = 2 * e + 1;
						else e = 2 * e;
					}
					return tree[e];
				}
				e--;
			}
			s >>= 1;
			e >>= 1;
		}
		if(s == e && tree[e] > thres){
			while(e < lim){
				if(tree[2*e+1] > thres) e = 2 * e + 1;
				else e = 2 * e;
			}
			return tree[e];
		}
		return 0;
	}
	int right_max(int s, int e, int thres){
		s += lim;
		e += lim;
		while(s < e){
			if(s % 2 == 1){
				if(tree[s] > thres){
					while(s < lim){
						if(tree[2*s] > thres) s = 2 * s;
						else s = 2 * s + 1;
					}
					return tree[s];
				}
				s++;
			}
			if(e % 2 == 0) e--;
			s >>= 1;
			e >>= 1;
		}
		if(s == e && tree[s] > thres){
			while(s < lim){
				if(tree[2*s] > thres) s = 2 * s;
				else s = 2 * s + 1;
			}
			return tree[s];
		}
		return 0;
	}
};
 
struct bit{
	int tree[MAXN];
	void add(int x, int v){
		while(x <= n){
			tree[x] += v;
			x += x & -x;
		}
	}
	int query(int x){
		int ret = 0;
		while(x){
			ret += tree[x];
			x -= x & -x;
		}
		return ret;
	}
}bit;
 
int main(){
	vector<int> v;
	scanf("%d",&n);
	for(int i=1; i<=n; i++){
		scanf("%d",&a[i]);
		v.push_back(a[i]);
	}
	sort(v.begin(), v.end());
	v.resize(unique(v.begin(), v.end()) - v.begin());
	for(int i=1; i<=n; i++) a[i] = lower_bound(v.begin(), v.end(), a[i]) - v.begin() + 1;
	for(int i=2; i<=n; i++){
		scanf("%d %d",&s[i],&e[i]);
		gph[s[i]].push_back(e[i]);
		gph[e[i]].push_back(s[i]);
	}
	segtree seg;
	seg.init();
	dfs(1, -1);
	for(int i=1; i<17; i++){
		for(int j=1; j<=n; j++){
			par[i][j] = par[i-1][par[i-1][j]];
		}
	}
	e[1] = 1;
	seg.upd(din[1], 1);
  int cnt = 0;
	for(int i=2; i<=n; i++){
		int pos = s[i];
		vector<pi> mem;
		int x;
		while((x = seg.query(din[pos], dout[pos])) != seg.tree[1]){
          cnt++;
			int pl = seg.left_max(0, din[pos] - 1, x);
			int pr = seg.right_max(dout[pos] + 1, seg.lim-1, x);
			int nxtpos = 1;
			if(pl > x){
				int aux = lca(pos, e[pl]);
				if(dep[aux] > dep[nxtpos]) nxtpos = aux;
			}
			if(pr > x){
				int aux = lca(pos, e[pr]);
				if(dep[aux] > dep[nxtpos]) nxtpos = aux;
			}
			mem.push_back(pi(dep[pos] - dep[nxtpos], a[e[x]]));
			pos = nxtpos;
		}
		lint ans = 0;
		for(auto &i : mem){
			ans += 1ll * bit.query(i.second - 1) * i.first;
			bit.add(i.second, i.first);
		}
		ans += 1ll * bit.query(a[e[seg.tree[1]]] - 1) * (dep[pos] + 1);
		printf("%lld\n", ans);
		for(auto &i : mem) bit.add(i.second, -i.first);
		seg.upd(din[e[i]], i);
	}
  assert(cnt <= 1500000);
}

Compilation message

construction.cpp: In function 'int main()':
construction.cpp:137:7: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
  scanf("%d",&n);
  ~~~~~^~~~~~~~~
construction.cpp:139:8: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
   scanf("%d",&a[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",&s[i],&e[i]);
   ~~~~~^~~~~~~~~~~~~~~~~~~~~
# Verdict Execution time Memory Grader output
1 Correct 5 ms 3832 KB Output is correct
2 Correct 5 ms 4096 KB Output is correct
3 Correct 5 ms 4164 KB Output is correct
4 Correct 6 ms 4312 KB Output is correct
5 Correct 6 ms 4440 KB Output is correct
6 Correct 6 ms 4592 KB Output is correct
7 Correct 6 ms 4592 KB Output is correct
8 Correct 5 ms 4592 KB Output is correct
9 Correct 5 ms 4592 KB Output is correct
10 Correct 5 ms 4592 KB Output is correct
11 Correct 5 ms 4592 KB Output is correct
12 Correct 5 ms 4592 KB Output is correct
13 Correct 5 ms 4592 KB Output is correct
14 Correct 5 ms 4592 KB Output is correct
15 Correct 6 ms 4592 KB Output is correct
16 Correct 8 ms 4628 KB Output is correct
17 Correct 6 ms 4628 KB Output is correct
18 Correct 7 ms 4628 KB Output is correct
19 Correct 6 ms 4628 KB Output is correct
20 Correct 5 ms 4628 KB Output is correct
21 Correct 6 ms 4628 KB Output is correct
22 Correct 5 ms 4628 KB Output is correct
23 Correct 8 ms 4628 KB Output is correct
24 Correct 6 ms 4628 KB Output is correct
25 Correct 6 ms 4760 KB Output is correct
26 Correct 5 ms 4812 KB Output is correct
27 Correct 6 ms 4812 KB Output is correct
28 Correct 5 ms 4812 KB Output is correct
29 Correct 5 ms 4896 KB Output is correct
30 Correct 8 ms 4896 KB Output is correct
31 Correct 6 ms 4896 KB Output is correct
32 Correct 5 ms 4896 KB Output is correct
33 Correct 5 ms 4896 KB Output is correct
34 Correct 5 ms 4896 KB Output is correct
35 Correct 5 ms 4896 KB Output is correct
36 Correct 5 ms 4896 KB Output is correct
37 Correct 5 ms 4896 KB Output is correct
38 Correct 6 ms 4896 KB Output is correct
39 Correct 7 ms 4896 KB Output is correct
40 Correct 6 ms 4896 KB Output is correct
41 Correct 7 ms 4896 KB Output is correct
42 Correct 5 ms 4896 KB Output is correct
43 Correct 5 ms 4900 KB Output is correct
44 Correct 6 ms 4924 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 5 ms 3832 KB Output is correct
2 Correct 5 ms 4096 KB Output is correct
3 Correct 5 ms 4164 KB Output is correct
4 Correct 6 ms 4312 KB Output is correct
5 Correct 6 ms 4440 KB Output is correct
6 Correct 6 ms 4592 KB Output is correct
7 Correct 6 ms 4592 KB Output is correct
8 Correct 5 ms 4592 KB Output is correct
9 Correct 5 ms 4592 KB Output is correct
10 Correct 5 ms 4592 KB Output is correct
11 Correct 5 ms 4592 KB Output is correct
12 Correct 5 ms 4592 KB Output is correct
13 Correct 5 ms 4592 KB Output is correct
14 Correct 5 ms 4592 KB Output is correct
15 Correct 6 ms 4592 KB Output is correct
16 Correct 8 ms 4628 KB Output is correct
17 Correct 6 ms 4628 KB Output is correct
18 Correct 7 ms 4628 KB Output is correct
19 Correct 6 ms 4628 KB Output is correct
20 Correct 5 ms 4628 KB Output is correct
21 Correct 6 ms 4628 KB Output is correct
22 Correct 5 ms 4628 KB Output is correct
23 Correct 8 ms 4628 KB Output is correct
24 Correct 6 ms 4628 KB Output is correct
25 Correct 6 ms 4760 KB Output is correct
26 Correct 5 ms 4812 KB Output is correct
27 Correct 6 ms 4812 KB Output is correct
28 Correct 5 ms 4812 KB Output is correct
29 Correct 5 ms 4896 KB Output is correct
30 Correct 8 ms 4896 KB Output is correct
31 Correct 6 ms 4896 KB Output is correct
32 Correct 5 ms 4896 KB Output is correct
33 Correct 5 ms 4896 KB Output is correct
34 Correct 5 ms 4896 KB Output is correct
35 Correct 5 ms 4896 KB Output is correct
36 Correct 5 ms 4896 KB Output is correct
37 Correct 5 ms 4896 KB Output is correct
38 Correct 6 ms 4896 KB Output is correct
39 Correct 7 ms 4896 KB Output is correct
40 Correct 6 ms 4896 KB Output is correct
41 Correct 7 ms 4896 KB Output is correct
42 Correct 5 ms 4896 KB Output is correct
43 Correct 5 ms 4900 KB Output is correct
44 Correct 6 ms 4924 KB Output is correct
45 Correct 7 ms 5044 KB Output is correct
46 Correct 15 ms 5576 KB Output is correct
47 Correct 21 ms 5652 KB Output is correct
48 Correct 15 ms 5652 KB Output is correct
49 Correct 10 ms 5932 KB Output is correct
50 Correct 11 ms 6008 KB Output is correct
51 Correct 9 ms 6084 KB Output is correct
52 Correct 10 ms 6084 KB Output is correct
53 Correct 10 ms 6088 KB Output is correct
54 Correct 11 ms 6168 KB Output is correct
55 Correct 10 ms 6380 KB Output is correct
56 Correct 10 ms 6380 KB Output is correct
57 Correct 24 ms 6380 KB Output is correct
58 Correct 23 ms 6380 KB Output is correct
59 Correct 24 ms 6428 KB Output is correct
60 Correct 25 ms 6632 KB Output is correct
61 Correct 12 ms 6840 KB Output is correct
62 Correct 11 ms 6840 KB Output is correct
63 Correct 15 ms 6860 KB Output is correct
64 Correct 14 ms 6860 KB Output is correct
65 Correct 14 ms 6860 KB Output is correct
66 Correct 15 ms 6904 KB Output is correct
67 Correct 16 ms 6956 KB Output is correct
68 Correct 9 ms 7268 KB Output is correct
69 Correct 10 ms 7316 KB Output is correct
70 Correct 10 ms 7316 KB Output is correct
71 Correct 10 ms 7316 KB Output is correct
72 Correct 24 ms 7316 KB Output is correct
73 Correct 25 ms 7316 KB Output is correct
74 Correct 11 ms 7448 KB Output is correct
75 Correct 10 ms 7496 KB Output is correct
76 Correct 10 ms 7496 KB Output is correct
77 Correct 10 ms 7624 KB Output is correct
78 Correct 10 ms 7688 KB Output is correct
79 Correct 10 ms 7688 KB Output is correct
80 Correct 14 ms 7688 KB Output is correct
81 Correct 11 ms 7844 KB Output is correct
82 Correct 11 ms 7896 KB Output is correct
83 Correct 14 ms 8088 KB Output is correct
84 Correct 11 ms 8088 KB Output is correct
85 Correct 11 ms 8088 KB Output is correct
86 Correct 12 ms 8088 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 5 ms 3832 KB Output is correct
2 Correct 5 ms 4096 KB Output is correct
3 Correct 5 ms 4164 KB Output is correct
4 Correct 6 ms 4312 KB Output is correct
5 Correct 6 ms 4440 KB Output is correct
6 Correct 6 ms 4592 KB Output is correct
7 Correct 6 ms 4592 KB Output is correct
8 Correct 5 ms 4592 KB Output is correct
9 Correct 5 ms 4592 KB Output is correct
10 Correct 5 ms 4592 KB Output is correct
11 Correct 5 ms 4592 KB Output is correct
12 Correct 5 ms 4592 KB Output is correct
13 Correct 5 ms 4592 KB Output is correct
14 Correct 5 ms 4592 KB Output is correct
15 Correct 6 ms 4592 KB Output is correct
16 Correct 8 ms 4628 KB Output is correct
17 Correct 6 ms 4628 KB Output is correct
18 Correct 7 ms 4628 KB Output is correct
19 Correct 6 ms 4628 KB Output is correct
20 Correct 5 ms 4628 KB Output is correct
21 Correct 6 ms 4628 KB Output is correct
22 Correct 5 ms 4628 KB Output is correct
23 Correct 8 ms 4628 KB Output is correct
24 Correct 6 ms 4628 KB Output is correct
25 Correct 6 ms 4760 KB Output is correct
26 Correct 5 ms 4812 KB Output is correct
27 Correct 6 ms 4812 KB Output is correct
28 Correct 5 ms 4812 KB Output is correct
29 Correct 5 ms 4896 KB Output is correct
30 Correct 8 ms 4896 KB Output is correct
31 Correct 6 ms 4896 KB Output is correct
32 Correct 5 ms 4896 KB Output is correct
33 Correct 5 ms 4896 KB Output is correct
34 Correct 5 ms 4896 KB Output is correct
35 Correct 5 ms 4896 KB Output is correct
36 Correct 5 ms 4896 KB Output is correct
37 Correct 5 ms 4896 KB Output is correct
38 Correct 6 ms 4896 KB Output is correct
39 Correct 7 ms 4896 KB Output is correct
40 Correct 6 ms 4896 KB Output is correct
41 Correct 7 ms 4896 KB Output is correct
42 Correct 5 ms 4896 KB Output is correct
43 Correct 5 ms 4900 KB Output is correct
44 Correct 6 ms 4924 KB Output is correct
45 Correct 7 ms 5044 KB Output is correct
46 Correct 15 ms 5576 KB Output is correct
47 Correct 21 ms 5652 KB Output is correct
48 Correct 15 ms 5652 KB Output is correct
49 Correct 10 ms 5932 KB Output is correct
50 Correct 11 ms 6008 KB Output is correct
51 Correct 9 ms 6084 KB Output is correct
52 Correct 10 ms 6084 KB Output is correct
53 Correct 10 ms 6088 KB Output is correct
54 Correct 11 ms 6168 KB Output is correct
55 Correct 10 ms 6380 KB Output is correct
56 Correct 10 ms 6380 KB Output is correct
57 Correct 24 ms 6380 KB Output is correct
58 Correct 23 ms 6380 KB Output is correct
59 Correct 24 ms 6428 KB Output is correct
60 Correct 25 ms 6632 KB Output is correct
61 Correct 12 ms 6840 KB Output is correct
62 Correct 11 ms 6840 KB Output is correct
63 Correct 15 ms 6860 KB Output is correct
64 Correct 14 ms 6860 KB Output is correct
65 Correct 14 ms 6860 KB Output is correct
66 Correct 15 ms 6904 KB Output is correct
67 Correct 16 ms 6956 KB Output is correct
68 Correct 9 ms 7268 KB Output is correct
69 Correct 10 ms 7316 KB Output is correct
70 Correct 10 ms 7316 KB Output is correct
71 Correct 10 ms 7316 KB Output is correct
72 Correct 24 ms 7316 KB Output is correct
73 Correct 25 ms 7316 KB Output is correct
74 Correct 11 ms 7448 KB Output is correct
75 Correct 10 ms 7496 KB Output is correct
76 Correct 10 ms 7496 KB Output is correct
77 Correct 10 ms 7624 KB Output is correct
78 Correct 10 ms 7688 KB Output is correct
79 Correct 10 ms 7688 KB Output is correct
80 Correct 14 ms 7688 KB Output is correct
81 Correct 11 ms 7844 KB Output is correct
82 Correct 11 ms 7896 KB Output is correct
83 Correct 14 ms 8088 KB Output is correct
84 Correct 11 ms 8088 KB Output is correct
85 Correct 11 ms 8088 KB Output is correct
86 Correct 12 ms 8088 KB Output is correct
87 Correct 37 ms 9064 KB Output is correct
88 Correct 119 ms 12460 KB Output is correct
89 Correct 545 ms 23800 KB Output is correct
90 Correct 570 ms 25820 KB Output is correct
91 Correct 536 ms 28184 KB Output is correct
92 Correct 159 ms 34080 KB Output is correct
93 Correct 151 ms 36268 KB Output is correct
94 Correct 147 ms 38396 KB Output is correct
95 Correct 224 ms 38716 KB Output is correct
96 Correct 267 ms 40848 KB Output is correct
97 Correct 233 ms 43100 KB Output is correct
98 Correct 232 ms 45164 KB Output is correct
99 Correct 238 ms 46476 KB Output is correct
100 Correct 1133 ms 46476 KB Output is correct
101 Correct 1190 ms 48432 KB Output is correct
102 Correct 1225 ms 50588 KB Output is correct
103 Correct 1429 ms 52652 KB Output is correct
104 Correct 309 ms 57056 KB Output is correct
105 Correct 318 ms 59040 KB Output is correct
106 Correct 283 ms 61452 KB Output is correct
107 Correct 509 ms 61452 KB Output is correct
108 Correct 577 ms 61564 KB Output is correct
109 Correct 620 ms 63280 KB Output is correct
110 Correct 152 ms 69060 KB Output is correct
111 Correct 238 ms 69192 KB Output is correct
112 Correct 238 ms 70812 KB Output is correct
113 Correct 245 ms 71776 KB Output is correct
114 Correct 1116 ms 71776 KB Output is correct
115 Correct 1157 ms 72732 KB Output is correct
116 Correct 270 ms 76520 KB Output is correct
117 Correct 244 ms 77164 KB Output is correct
118 Correct 252 ms 78620 KB Output is correct
119 Correct 246 ms 80060 KB Output is correct
120 Correct 219 ms 82132 KB Output is correct
121 Correct 238 ms 83460 KB Output is correct
122 Correct 226 ms 84656 KB Output is correct
123 Correct 321 ms 87444 KB Output is correct
124 Correct 352 ms 88932 KB Output is correct
125 Correct 351 ms 90600 KB Output is correct
126 Correct 318 ms 92316 KB Output is correct
127 Correct 280 ms 93412 KB Output is correct
128 Correct 291 ms 94672 KB Output is correct