結果
問題 | No.1018 suffixsuffixsuffix |
ユーザー | Kiri8128 |
提出日時 | 2020-04-04 01:49:37 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 404 ms / 2,000 ms |
コード長 | 3,802 bytes |
コンパイル時間 | 274 ms |
コンパイル使用メモリ | 87,212 KB |
実行使用メモリ | 133,940 KB |
最終ジャッジ日時 | 2023-09-16 05:12:43 |
合計ジャッジ時間 | 10,650 ms |
ジャッジサーバーID (参考情報) |
judge14 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 74 ms
71,336 KB |
testcase_01 | AC | 72 ms
71,820 KB |
testcase_02 | AC | 72 ms
71,896 KB |
testcase_03 | AC | 73 ms
71,852 KB |
testcase_04 | AC | 74 ms
71,472 KB |
testcase_05 | AC | 75 ms
71,484 KB |
testcase_06 | AC | 77 ms
71,644 KB |
testcase_07 | AC | 77 ms
71,532 KB |
testcase_08 | AC | 76 ms
71,172 KB |
testcase_09 | AC | 190 ms
90,060 KB |
testcase_10 | AC | 190 ms
91,564 KB |
testcase_11 | AC | 195 ms
92,824 KB |
testcase_12 | AC | 189 ms
92,500 KB |
testcase_13 | AC | 191 ms
91,976 KB |
testcase_14 | AC | 224 ms
96,852 KB |
testcase_15 | AC | 380 ms
117,676 KB |
testcase_16 | AC | 385 ms
122,000 KB |
testcase_17 | AC | 364 ms
117,412 KB |
testcase_18 | AC | 236 ms
100,372 KB |
testcase_19 | AC | 151 ms
85,852 KB |
testcase_20 | AC | 147 ms
86,792 KB |
testcase_21 | AC | 150 ms
86,820 KB |
testcase_22 | AC | 155 ms
86,724 KB |
testcase_23 | AC | 156 ms
87,376 KB |
testcase_24 | AC | 371 ms
126,500 KB |
testcase_25 | AC | 392 ms
130,144 KB |
testcase_26 | AC | 390 ms
124,680 KB |
testcase_27 | AC | 376 ms
127,352 KB |
testcase_28 | AC | 367 ms
120,432 KB |
testcase_29 | AC | 176 ms
96,500 KB |
testcase_30 | AC | 173 ms
93,228 KB |
testcase_31 | AC | 169 ms
90,912 KB |
testcase_32 | AC | 175 ms
96,056 KB |
testcase_33 | AC | 169 ms
94,120 KB |
testcase_34 | AC | 75 ms
71,404 KB |
testcase_35 | AC | 117 ms
91,692 KB |
testcase_36 | AC | 126 ms
90,256 KB |
testcase_37 | AC | 404 ms
133,940 KB |
ソースコード
from bisect import bisect_left as bl from bisect import bisect_right as br def SA_IS(a): a += [0] k = max(a) + 1 n = len(a) def induce_l(sa, a, n, k, stype): bucket = get_buckets(a, k, 1) for i in range(n): j = sa[i] - 1 if j >= 0 and (not stype[j]): sa[bucket[a[j]]] = j bucket[a[j]] += 1 def induce_s(sa, a, n, k, stype): bucket = get_buckets(a, k, 0) for i in range(n)[::-1]: j = sa[i] - 1 if j >= 0 and stype[j]: bucket[a[j]] -= 1 sa[bucket[a[j]]] = j def get_buckets(a, k, start = 0): bucket = [0] * k for item in a: bucket[item] += 1 s = 0 for i in range(k): s += bucket[i] bucket[i] = s - (bucket[i] if start else 0) return bucket def set_lms(a, n, k, default_order): bucket = get_buckets(a, k) sa = [-1] * n for i in default_order[::-1]: bucket[a[i]] -= 1 sa[bucket[a[i]]] = i return sa def induce(a, n, k, stype, default_order): sa = set_lms(a, n, k, default_order) induce_l(sa, a, n, k, stype) induce_s(sa, a, n, k, stype) return sa def rename_LMS_substring(sa, a, n, stype, LMS, l): sa = [_s for _s in sa if LMS[_s]] tmp = [-1] * (n//2) + [0] dupl = 0 for _ in range(1, l): i, j = sa[_-1], sa[_] for ii in range(n): if a[i+ii] != a[j+ii] or stype[i+ii] != stype[j+ii]: break if ii and (LMS[i+ii] or LMS[j+ii]): dupl += 1 break tmp[j//2] = _ - dupl tmp = [t for t in tmp if t >= 0] return tmp, dupl def calc(a, n, k): stype = [1] * n for i in range(n-1)[::-1]: if a[i] > a[i+1] or (a[i] == a[i+1] and stype[i+1] == 0): stype[i] = 0 LMS = [1 if stype[i] and not stype[i-1] else 0 for i in range(n-1)] + [1] l = sum(LMS) lms = [i for i in range(n) if LMS[i]] sa = induce(a, n, k, stype, lms) renamed_LMS, dupl = rename_LMS_substring(sa, a, n, stype, LMS, l) if dupl: sub_sa = calc(renamed_LMS, l, l - dupl) else: sub_sa = [0] * l for i in range(l): sub_sa[renamed_LMS[i]] = i lms = [lms[sub_sa[i]] for i in range(l)] sa = induce(a, n, k, stype, lms) return sa sa = calc(a, n, k) return sa N, M, Q = map(int, input().split()) S = [ord(a)-96 for a in input()] K = [int(a)-1 for a in input().split()] if S.count(S[0]) == N: print(*[N * M - k for k in K]) elif M <= 2: S *= M SA = SA_IS(S)[1:] print(*[SA[k] + 1 for k in K]) else: for i in range(1, int(N**0.5)+2): if N % i == 0 and S[i:] == S[:-i]: S = S[:i] M *= N // i N = i break S *= 2 SA = SA_IS(S)[1:] X = [] i = 0 for s in SA: if s >= N: # X.append((i, (M-1) * N + (s - N))) X.append((i << 47) + (M-1) * N + (s - N) + 1) i += 1 else: # X.append((i, s)) X.append((i << 47) + (M-2) * N + s + 1) i += M-1 def chk(a): return (a >> 47, a - (a >> 47 << 47) - 1) # print("X =", [chk(x) for x in X]) X = sorted(X) # print("SA =", SA) # print("X =", X) # print("X =", [chk(x) for x in X]) def calc(k): t = bl(X, k + 1 << 47) - 1 a, b = chk(X[t]) return b - (k - a) * N + 1 print(*[calc(k) for k in K])