Submission #545497

# Submission time Handle Problem Language Result Execution time Memory
545497 2022-04-04T15:56:03 Z rainboy Space Pirate (JOI14_space_pirate) C
47 / 100
2000 ms 84412 KB
#include <stdio.h>
#include <stdlib.h>
#include <string.h>

#define N	100000
#define L	60	/* L = ceil(log2(10^18)) */

int *ej[N], eo[N], pp[N], pp_[L + 1][N], ii[N + 1], tt[N], tt_[N], n; long long cc[N], cc_[L + 1][N], k;

void append(int i, int j) {
	int o = eo[i]++;

	if (o >= 2 && (o & o - 1) == 0)
		ej[i] = (int *) realloc(ej[i], o * 2 * sizeof *ej[i]);
	ej[i][o] = j;
}

void dfs1(int i, int t) {
	int o;

	if (tt_[i] != -1)
		return;
	tt_[i] = t;
	for (o = eo[i]; o--; ) {
		int j = ej[i][o];

		dfs1(j, t);
	}
}

void update(int i, long long k, int x) {
	int l;

	for (l = L; l >= 0; l--)
		if (k >= 1LL << l)
			cc_[l][i] += x, k -= 1LL << l, i = pp_[l][i];
}

int qu[N * 2];

void dfs2(int p, int i, int cnt, int h_, int h1) {
	int o, h;

	qu[cnt++] = i;
	if (tt[i] != -1)
		h1 = cnt - 1;
	for (h = h1; h >= h_; h--)
		cc[qu[cnt - 1 - (k + cnt - (tt[qu[h]] + h) - 1) % (cnt - h)]]++;
	for (o = eo[i]; o--; ) {
		int j = ej[i][o];

		if (j != p)
			dfs2(i, j, cnt, h_, h1);
	}
}

int main() {
	int h, i, i_, j, t, l, cnt;
	long long k_;

	scanf("%d%lld", &n, &k);
	for (i = 0; i < n; i++)
		ej[i] = (int *) malloc(2 * sizeof *ej[i]);
	for (i = 0; i < n; i++) {
		scanf("%d", &pp[i]), pp[i]--;
		append(pp[i], i);
	}
	for (i = 0; i < n; i++)
		pp_[0][i] = pp[i];
	for (l = 1; l <= L; l++)
		for (i = 0; i < n; i++)
			pp_[l][i] = pp_[l - 1][pp_[l - 1][i]];
	memset(tt, -1, n * sizeof *tt), memset(tt_, -1, n * sizeof *tt_);
	for (i = 0, t = 0; tt[i] == -1; i = pp[i], t++)
		ii[tt[i] = t] = i, dfs1(i, t);
	ii[t] = i_ = i;
	for (i = 0; i < n; i++)
		if (tt_[i] == -1)
			tt_[i] = t;
	j = i_, k_ = (k - t) % (t - tt[i_]);
	while (k_--)
		j = pp[j];
	for (i = 0; i < n; i++)
		if (tt[i] == -1)
			cc[j] += n;
	for (i = 0; i < n; i++)
		update(i, k, 1), update(i, k - tt_[i], -1);
	for (l = L; l > 0; l--)
		for (i = 0; i < n; i++)
			cc_[l - 1][i] += cc_[l][i], cc_[l - 1][pp_[l - 1][i]] += cc_[l][i];
	for (i = 0; i < n; i++)
		cc[i] += cc_[0][i];
	cnt = 0;
	for (h = t - 1; h > tt[i_]; h--)
		qu[cnt++] = ii[h];
	for (h = t; h > tt[i_]; h--)
		dfs2(ii[h - 1], ii[h], cnt, cnt - (t - tt[i_]) + 1, cnt), qu[cnt++] = ii[h];
	for (i = 0; i < n; i++)
		printf("%lld\n", cc[i]);
	return 0;
}

Compilation message

space_pirate.c: In function 'append':
space_pirate.c:13:23: warning: suggest parentheses around '-' in operand of '&' [-Wparentheses]
   13 |  if (o >= 2 && (o & o - 1) == 0)
      |                     ~~^~~
space_pirate.c: In function 'main':
space_pirate.c:61:2: warning: ignoring return value of 'scanf' declared with attribute 'warn_unused_result' [-Wunused-result]
   61 |  scanf("%d%lld", &n, &k);
      |  ^~~~~~~~~~~~~~~~~~~~~~~
space_pirate.c:65:3: warning: ignoring return value of 'scanf' declared with attribute 'warn_unused_result' [-Wunused-result]
   65 |   scanf("%d", &pp[i]), pp[i]--;
      |   ^~~~~~~~~~~~~~~~~~~
# Verdict Execution time Memory Grader output
1 Correct 1 ms 980 KB Output is correct
2 Correct 1 ms 980 KB Output is correct
3 Correct 1 ms 980 KB Output is correct
4 Correct 1 ms 980 KB Output is correct
5 Correct 1 ms 980 KB Output is correct
6 Correct 1 ms 980 KB Output is correct
7 Correct 1 ms 980 KB Output is correct
8 Correct 1 ms 980 KB Output is correct
9 Correct 1 ms 980 KB Output is correct
10 Correct 1 ms 980 KB Output is correct
11 Correct 1 ms 980 KB Output is correct
12 Correct 1 ms 980 KB Output is correct
13 Correct 1 ms 980 KB Output is correct
14 Correct 1 ms 980 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 980 KB Output is correct
2 Correct 1 ms 980 KB Output is correct
3 Correct 1 ms 980 KB Output is correct
4 Correct 1 ms 980 KB Output is correct
5 Correct 1 ms 980 KB Output is correct
6 Correct 1 ms 980 KB Output is correct
7 Correct 1 ms 980 KB Output is correct
8 Correct 1 ms 980 KB Output is correct
9 Correct 1 ms 980 KB Output is correct
10 Correct 1 ms 980 KB Output is correct
11 Correct 1 ms 980 KB Output is correct
12 Correct 1 ms 980 KB Output is correct
13 Correct 1 ms 980 KB Output is correct
14 Correct 1 ms 980 KB Output is correct
15 Correct 8 ms 3156 KB Output is correct
16 Correct 5 ms 3156 KB Output is correct
17 Correct 9 ms 3256 KB Output is correct
18 Correct 158 ms 3420 KB Output is correct
19 Correct 49 ms 3328 KB Output is correct
20 Correct 101 ms 3388 KB Output is correct
21 Correct 103 ms 3460 KB Output is correct
22 Correct 37 ms 3348 KB Output is correct
23 Correct 98 ms 3452 KB Output is correct
24 Correct 26 ms 3360 KB Output is correct
25 Correct 4 ms 3136 KB Output is correct
26 Correct 178 ms 3440 KB Output is correct
27 Correct 42 ms 3308 KB Output is correct
28 Correct 65 ms 3316 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 465 ms 77720 KB Output is correct
2 Execution timed out 2088 ms 84412 KB Time limit exceeded
3 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 1 ms 980 KB Output is correct
2 Correct 1 ms 980 KB Output is correct
3 Correct 1 ms 980 KB Output is correct
4 Correct 1 ms 980 KB Output is correct
5 Correct 1 ms 980 KB Output is correct
6 Correct 1 ms 980 KB Output is correct
7 Correct 1 ms 980 KB Output is correct
8 Correct 1 ms 980 KB Output is correct
9 Correct 1 ms 980 KB Output is correct
10 Correct 1 ms 980 KB Output is correct
11 Correct 1 ms 980 KB Output is correct
12 Correct 1 ms 980 KB Output is correct
13 Correct 1 ms 980 KB Output is correct
14 Correct 1 ms 980 KB Output is correct
15 Correct 8 ms 3156 KB Output is correct
16 Correct 5 ms 3156 KB Output is correct
17 Correct 9 ms 3256 KB Output is correct
18 Correct 158 ms 3420 KB Output is correct
19 Correct 49 ms 3328 KB Output is correct
20 Correct 101 ms 3388 KB Output is correct
21 Correct 103 ms 3460 KB Output is correct
22 Correct 37 ms 3348 KB Output is correct
23 Correct 98 ms 3452 KB Output is correct
24 Correct 26 ms 3360 KB Output is correct
25 Correct 4 ms 3136 KB Output is correct
26 Correct 178 ms 3440 KB Output is correct
27 Correct 42 ms 3308 KB Output is correct
28 Correct 65 ms 3316 KB Output is correct
29 Correct 465 ms 77720 KB Output is correct
30 Execution timed out 2088 ms 84412 KB Time limit exceeded
31 Halted 0 ms 0 KB -