#include <stdio.h>
#include <stdlib.h>
#define N 70000
unsigned int Z = 12345;
int rand_() {
return (Z *= 3) >> 1;
}
int *ej[N], *ew[N], eo[N], n, d_;
void append(int i, int j, int w) {
int o = eo[i]++;
if (o >= 2 && (o & o - 1) == 0) {
ej[i] = (int *) realloc(ej[i], o * 2 * sizeof *ej[i]);
ew[i] = (int *) realloc(ew[i], o * 2 * sizeof *ew[i]);
}
ej[i][o] = j, ew[i][o] = w;
}
void detach(int i, int j) {
int o;
for (o = eo[i]; o--; )
if (ej[i][o] == j) {
eo[i]--, ej[i][o] = ej[i][eo[i]], ew[i][o] = ew[i][eo[i]];
return;
}
}
int sz[N], n_, c;
void dfs1(int p, int i) {
int o, centroid;
sz[i] = 1, centroid = 1;
for (o = eo[i]; o--; ) {
int j = ej[i][o];
if (j != p) {
dfs1(i, j);
sz[i] += sz[j];
if (sz[j] * 2 > n_)
centroid = 0;
}
}
if ((n_ - sz[i]) * 2 > n_)
centroid = 0;
if (centroid)
c = i;
}
int ii[N], qu[N], cnt;
int pp[N], rr[N]; long long dd[N]; int kk[N];
void dfs2(int p, int i, int r, long long d) {
int o, h, lower, upper;
qu[cnt++] = i;
ii[n_++] = i;
pp[i] = p, rr[i] = r, dd[i] = d, sz[i] = 1, kk[i] = 1;
for (o = eo[i]; o--; ) {
int j = ej[i][o], w = ew[i][o];
if (j != p) {
dfs2(i, j, p == -1 ? j : r, d + w);
sz[i] += sz[j];
}
}
lower = -1, upper = cnt - 1;
while (upper - lower > 1) {
h = (lower + upper) / 2;
if (dd[i] - dd[qu[h]] > d_)
lower = h;
else
upper = h;
}
if (lower != -1)
kk[qu[upper]] += kk[i];
cnt--;
}
void sort(int *ii, int l, int r) {
while (l < r) {
int i = l, j = l, k = r, i_ = ii[l + rand_() % (r - l)], tmp;
while (j < k) {
int c = dd[ii[j]] != dd[i_] ? (dd[ii[j]] < dd[i_] ? -1 : 1) : sz[i_] - sz[ii[j]];
if (c == 0)
j++;
else if (c < 0) {
tmp = ii[i], ii[i] = ii[j], ii[j] = tmp;
i++, j++;
} else {
k--;
tmp = ii[j], ii[j] = ii[k], ii[k] = tmp;
}
}
sort(ii, l, i);
l = k;
}
}
long long ans[N]; int kk_[N], ss[N];
void dfs3(int p, int i) {
int o, h, lower, upper;
qu[cnt++] = i;
lower = -1, upper = cnt - 1;
while (upper - lower > 1) {
h = (lower + upper) / 2;
if (dd[i] - dd[qu[h]] > d_)
lower = h;
else
upper = h;
}
kk_[i] += ss[qu[upper]];
ss[i] = (p == -1 ? 0 : ss[p]) + kk_[i];
for (o = eo[i]; o--; ) {
int j = ej[i][o];
if (j != p)
kk_[j] -= ss[qu[upper]], dfs3(i, j);
}
if (p != -1)
ans[p] += (long long) kk_[i] * sz[i];
cnt--;
}
void cd(int n, int i) {
static int kk1[N];
int r, h, h1, h2, j, o, k;
n_ = n, dfs1(-1, i), r = c;
n_ = 0, dfs2(-1, r, r, 0);
for (h = 1; h < n; h++) {
i = ii[h];
ans[i] += (long long) (kk[i] - 1) * (n - sz[rr[i]]);
}
sort(ii, 0, n);
kk1[r] = 0;
for (o = eo[r]; o--; ) {
i = ej[r][o];
kk1[i] = 0;
}
for (h1 = 0, h2 = n - 1, k = 0; h1 < n; h1++) {
i = ii[h1];
while (h2 >= 0 && dd[i] + dd[j = ii[h2]] > d_) {
if (dd[j] <= d_)
k += kk[j], kk1[rr[j]] += kk[j];
h2--;
}
kk_[i] = k - kk1[rr[i]];
}
for (h = n - 1; h > 0; h--) {
i = ii[h];
kk_[i] -= kk_[pp[i]];
}
dfs3(-1, r);
for (o = eo[r]; o--; ) {
i = ej[r][o];
detach(i, r), cd(sz[i], i);
}
}
int main() {
int h, i, j, w;
scanf("%d%d", &n, &d_);
for (i = 0; i < n; i++) {
ej[i] = (int *) malloc(2 * sizeof *ej[i]);
ew[i] = (int *) malloc(2 * sizeof *ew[i]);
}
for (h = 0; h < n - 1; h++) {
scanf("%d%d%d", &i, &j, &w);
append(i, j, w), append(j, i, w);
}
cd(n, 0);
for (i = 0; i < n; i++)
printf("%lld\n", ans[i]);
return 0;
}
Compilation message
Main.c: In function 'append':
Main.c:17:23: warning: suggest parentheses around '-' in operand of '&' [-Wparentheses]
17 | if (o >= 2 && (o & o - 1) == 0) {
| ~~^~~
Main.c: In function 'main':
Main.c:174:2: warning: ignoring return value of 'scanf' declared with attribute 'warn_unused_result' [-Wunused-result]
174 | scanf("%d%d", &n, &d_);
| ^~~~~~~~~~~~~~~~~~~~~~
Main.c:180:3: warning: ignoring return value of 'scanf' declared with attribute 'warn_unused_result' [-Wunused-result]
180 | scanf("%d%d%d", &i, &j, &w);
| ^~~~~~~~~~~~~~~~~~~~~~~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
2396 KB |
Output is correct |
2 |
Correct |
0 ms |
2396 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
2396 KB |
Output is correct |
2 |
Correct |
0 ms |
2396 KB |
Output is correct |
3 |
Correct |
1 ms |
2396 KB |
Output is correct |
4 |
Correct |
1 ms |
2396 KB |
Output is correct |
5 |
Correct |
1 ms |
2396 KB |
Output is correct |
6 |
Correct |
2 ms |
2652 KB |
Output is correct |
7 |
Correct |
1 ms |
2652 KB |
Output is correct |
8 |
Correct |
0 ms |
2396 KB |
Output is correct |
9 |
Correct |
1 ms |
2396 KB |
Output is correct |
10 |
Correct |
1 ms |
2396 KB |
Output is correct |
11 |
Correct |
1 ms |
2396 KB |
Output is correct |
12 |
Correct |
1 ms |
2396 KB |
Output is correct |
13 |
Correct |
2 ms |
2396 KB |
Output is correct |
14 |
Correct |
1 ms |
2396 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
2396 KB |
Output is correct |
2 |
Correct |
268 ms |
15164 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
2396 KB |
Output is correct |
2 |
Correct |
0 ms |
2396 KB |
Output is correct |
3 |
Correct |
1 ms |
2396 KB |
Output is correct |
4 |
Correct |
268 ms |
15164 KB |
Output is correct |
5 |
Correct |
284 ms |
15444 KB |
Output is correct |
6 |
Correct |
246 ms |
14932 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
2396 KB |
Output is correct |
2 |
Correct |
0 ms |
2396 KB |
Output is correct |
3 |
Correct |
177 ms |
10924 KB |
Output is correct |
4 |
Correct |
285 ms |
14932 KB |
Output is correct |
5 |
Correct |
237 ms |
15404 KB |
Output is correct |
6 |
Correct |
1 ms |
2484 KB |
Output is correct |
7 |
Correct |
200 ms |
11100 KB |
Output is correct |
8 |
Correct |
195 ms |
11092 KB |
Output is correct |
9 |
Correct |
190 ms |
11092 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
2396 KB |
Output is correct |
2 |
Correct |
0 ms |
2396 KB |
Output is correct |
3 |
Correct |
177 ms |
10924 KB |
Output is correct |
4 |
Correct |
285 ms |
14932 KB |
Output is correct |
5 |
Correct |
237 ms |
15404 KB |
Output is correct |
6 |
Correct |
1 ms |
2484 KB |
Output is correct |
7 |
Correct |
200 ms |
11100 KB |
Output is correct |
8 |
Correct |
195 ms |
11092 KB |
Output is correct |
9 |
Correct |
190 ms |
11092 KB |
Output is correct |
10 |
Correct |
107 ms |
11604 KB |
Output is correct |
11 |
Correct |
130 ms |
11020 KB |
Output is correct |
12 |
Correct |
117 ms |
11092 KB |
Output is correct |
13 |
Correct |
121 ms |
11204 KB |
Output is correct |
14 |
Correct |
121 ms |
11164 KB |
Output is correct |
15 |
Correct |
118 ms |
11092 KB |
Output is correct |
16 |
Correct |
30 ms |
11492 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
2396 KB |
Output is correct |
2 |
Correct |
0 ms |
2396 KB |
Output is correct |
3 |
Correct |
1 ms |
2396 KB |
Output is correct |
4 |
Correct |
1 ms |
2396 KB |
Output is correct |
5 |
Correct |
1 ms |
2396 KB |
Output is correct |
6 |
Correct |
2 ms |
2652 KB |
Output is correct |
7 |
Correct |
1 ms |
2652 KB |
Output is correct |
8 |
Correct |
0 ms |
2396 KB |
Output is correct |
9 |
Correct |
1 ms |
2396 KB |
Output is correct |
10 |
Correct |
1 ms |
2396 KB |
Output is correct |
11 |
Correct |
1 ms |
2396 KB |
Output is correct |
12 |
Correct |
1 ms |
2396 KB |
Output is correct |
13 |
Correct |
2 ms |
2396 KB |
Output is correct |
14 |
Correct |
1 ms |
2396 KB |
Output is correct |
15 |
Correct |
1 ms |
2396 KB |
Output is correct |
16 |
Correct |
268 ms |
15164 KB |
Output is correct |
17 |
Correct |
284 ms |
15444 KB |
Output is correct |
18 |
Correct |
246 ms |
14932 KB |
Output is correct |
19 |
Correct |
177 ms |
10924 KB |
Output is correct |
20 |
Correct |
285 ms |
14932 KB |
Output is correct |
21 |
Correct |
237 ms |
15404 KB |
Output is correct |
22 |
Correct |
1 ms |
2484 KB |
Output is correct |
23 |
Correct |
200 ms |
11100 KB |
Output is correct |
24 |
Correct |
195 ms |
11092 KB |
Output is correct |
25 |
Correct |
190 ms |
11092 KB |
Output is correct |
26 |
Correct |
107 ms |
11604 KB |
Output is correct |
27 |
Correct |
130 ms |
11020 KB |
Output is correct |
28 |
Correct |
117 ms |
11092 KB |
Output is correct |
29 |
Correct |
121 ms |
11204 KB |
Output is correct |
30 |
Correct |
121 ms |
11164 KB |
Output is correct |
31 |
Correct |
118 ms |
11092 KB |
Output is correct |
32 |
Correct |
30 ms |
11492 KB |
Output is correct |
33 |
Correct |
158 ms |
10852 KB |
Output is correct |
34 |
Correct |
119 ms |
12356 KB |
Output is correct |
35 |
Correct |
132 ms |
11604 KB |
Output is correct |
36 |
Correct |
137 ms |
11604 KB |
Output is correct |
37 |
Correct |
116 ms |
11252 KB |
Output is correct |
38 |
Correct |
121 ms |
11256 KB |
Output is correct |
39 |
Correct |
142 ms |
11348 KB |
Output is correct |
40 |
Correct |
43 ms |
12004 KB |
Output is correct |