結果
問題 | No.515 典型LCP |
ユーザー | rpy3cpp |
提出日時 | 2017-05-06 00:19:31 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
RE
|
実行時間 | - |
コード長 | 2,576 bytes |
コンパイル時間 | 493 ms |
コンパイル使用メモリ | 12,928 KB |
実行使用メモリ | 55,776 KB |
最終ジャッジ日時 | 2024-09-14 09:38:18 |
合計ジャッジ時間 | 4,458 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge6 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | RE | - |
testcase_01 | TLE | - |
testcase_02 | -- | - |
testcase_03 | -- | - |
testcase_04 | -- | - |
testcase_05 | -- | - |
testcase_06 | -- | - |
testcase_07 | -- | - |
testcase_08 | -- | - |
testcase_09 | -- | - |
testcase_10 | -- | - |
testcase_11 | -- | - |
testcase_12 | -- | - |
testcase_13 | -- | - |
testcase_14 | -- | - |
testcase_15 | -- | - |
testcase_16 | -- | - |
ソースコード
import math class RMQ(object): def __init__(self, A): self._A = A self._preprocess() def _preprocess(self): ''' RMQ の前処理。 ''' n = len(self._A) max_j = int(math.log2(n)) self._M = [list(range(n))] for j in range(0, max_j): shift = 1 << j Mj = self._M[j] Mjnext = [] for k1, k2 in zip(Mj, Mj[shift:]): k = k1 if self._A[k1] < self._A[k2] else k2 Mjnext.append(k) self._M.append(Mjnext) def query(self, i, j): if i == j: return i if i > j: i, j = j, i el = int(math.log2(j - i)) k1 = self._M[el][i] k2 = self._M[el][j - (1 << el) + 1] rmq = k1 if self._A[k1] < self._A[k2] else k2 return rmq class LCP(object): '''Longest Common Prefix を計算するクラス lcp = LCP(lst_of_strings) lcp.get(i, j) ''' def __init__(self, strings): strings_with_index = [(s, i) for i, s in enumerate(strings)] strings_with_index.sort() sorted_strings = [] self._index_converter = [0] * len(strings) j = 0 for s, i in strings_with_index: sorted_strings.append(s) self._index_converter[i] = j j += 1 self._init_lcp_neighbours(sorted_strings) self._rmq = RMQ(self._lcp_neighbours) self._lens = list(map(len, sorted_strings)) def _init_lcp_neighbours(self, ss): n = len(ss) self._lcp_neighbours = [self._calc_lcp(ss[i], ss[i + 1]) for i in range(n - 1)] def _calc_lcp(self, str0, str1): mylcp = 0 for c0, c1 in zip(str0, str1): if c0 != c1: return mylcp else: mylcp += 1 return mylcp def get(self, i, j): ''' strings[i] と strings[j] の LCP の長さを返す ''' ii = self._index_converter[i] jj = self._index_converter[j] if ii == jj: return self._lens[ii] kk = self._rmq.query(ii, jj - 1) return self._lcp_neighbours[kk] if __name__ == '__main__': N = int(input()) strings = [input().rstrip() for _ in range(N)] M, x, d = map(int, input().split()) lcp = LCP(strings) Nm1 = N - 1 NN = N * (N - 1) cum = 0 for k in range(M): i, j = divmod(x, Nm1) if (i > j): i, j = j, i else: j += 1 cum += lcp.get(i, j) x = (x + d) % NN print(cum)