結果

問題 No.1018 suffixsuffixsuffix
ユーザー Kiri8128Kiri8128
提出日時 2020-04-04 01:32:34
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 3,626 bytes
コンパイル時間 273 ms
コンパイル使用メモリ 87,028 KB
実行使用メモリ 130,640 KB
最終ジャッジ日時 2023-09-16 05:11:43
合計ジャッジ時間 11,649 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 73 ms
71,324 KB
testcase_01 AC 74 ms
71,208 KB
testcase_02 AC 76 ms
71,564 KB
testcase_03 AC 76 ms
71,220 KB
testcase_04 AC 80 ms
71,272 KB
testcase_05 AC 77 ms
71,352 KB
testcase_06 AC 79 ms
71,356 KB
testcase_07 AC 74 ms
71,408 KB
testcase_08 AC 76 ms
71,256 KB
testcase_09 AC 187 ms
89,940 KB
testcase_10 AC 190 ms
90,704 KB
testcase_11 AC 191 ms
92,684 KB
testcase_12 AC 186 ms
92,384 KB
testcase_13 AC 188 ms
91,804 KB
testcase_14 WA -
testcase_15 AC 372 ms
117,460 KB
testcase_16 AC 377 ms
122,484 KB
testcase_17 AC 355 ms
117,556 KB
testcase_18 WA -
testcase_19 AC 145 ms
85,512 KB
testcase_20 AC 143 ms
86,120 KB
testcase_21 AC 154 ms
86,080 KB
testcase_22 AC 151 ms
86,248 KB
testcase_23 AC 152 ms
86,916 KB
testcase_24 AC 373 ms
126,544 KB
testcase_25 AC 378 ms
129,860 KB
testcase_26 AC 356 ms
120,512 KB
testcase_27 AC 373 ms
125,912 KB
testcase_28 AC 360 ms
120,940 KB
testcase_29 WA -
testcase_30 WA -
testcase_31 WA -
testcase_32 WA -
testcase_33 WA -
testcase_34 AC 73 ms
71,172 KB
testcase_35 AC 119 ms
91,248 KB
testcase_36 AC 125 ms
90,264 KB
testcase_37 AC 393 ms
130,640 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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(*[s + 1 for s in SA])
else:
    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])
0