結果
| 問題 |
No.563 超高速一人かるた large
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2017-08-26 02:18:08 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 2,338 bytes |
| コンパイル時間 | 464 ms |
| コンパイル使用メモリ | 82,656 KB |
| 実行使用メモリ | 156,124 KB |
| 最終ジャッジ日時 | 2024-10-15 16:25:10 |
| 合計ジャッジ時間 | 11,154 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 5 TLE * 2 -- * 11 |
ソースコード
import collections
MAXN = 2001
mod = 10**9 + 7
fac = [1]
for k in range(1, MAXN):
fac.append(fac[-1] * k % mod)
nCr_memo = [[0] * MAXN for _ in range(MAXN)]
nCr_memo[0][0] = 1
for n in range(1, MAXN):
nCr_memo[n][0] = 1
nCr_memo[n][n] = 1
for r in range(1, n):
nCr_memo[n][r] = (nCr_memo[n - 1][r - 1] + nCr_memo[n - 1][r]) % mod
def nCr(n, r):
if r < 0:
return 0
return nCr_memo[n][r]
def read_data():
N = int(input())
Ss = [input() for i in range(N)]
return N, Ss
def common_prefix(S1, S2):
for i in range(min(len(S1), len(S2))):
if S1[i] != S2[i]:
return i
return min(len(S1), len(S2))
def lcp(N, Ss):
row = [common_prefix(Ss[i], Ss[i + 1]) for i in range(N - 1)]
return row
def preprocess(N, Ss):
Ss.sort()
row = lcp(N, Ss)
mat = [[1] * N for _ in range(N)]
for i in range(N):
mat[i][i] = 0
for i in range(N - 1):
lcp_ij = row[i]
for j in range(i + 1, N - 1):
mat[i][j] = lcp_ij + 1
mat[j][i] = lcp_ij + 1
lcp_ij = min(lcp_ij, row[j])
if lcp_ij == 0:
break
mat[i][N - 1] = lcp_ij + 1
mat[N - 1][i] = lcp_ij + 1
return mat
def preprocess_old(N, Ss):
Ss.sort()
mat = [[0] * N for _ in range(N)]
for i in range(N):
for j in range(i + 1, N):
mat[i][j] = common_prefix(Ss[i], Ss[j]) + 1
mat[j][i] = mat[i][j]
return mat
def preprocess2(N, mat):
fs = []
for i in range(N):
counter = collections.Counter(mat[i])
lst = list(counter.items())
lst.sort()
fL = 0
fH = N - 1
data = []
for j, fj in lst[1:]:
fH -= fj
data.append((j, fj, fL, fH))
fL += fj
fs.append(data)
return fs
def solve(N, Ss):
mat = preprocess(N, Ss)
fs = preprocess2(N, mat)
ans = [0]
for k in range(N - 1):
cum = (ans[-1] * (N - k)) % mod
for f in fs:
for j, fj, fL, fH in f:
cum = (cum + j * (nCr(fL + fj, k - fH) - nCr(fL, k - fH - fj)) * fac[k]) % mod
ans.append(cum)
cum = (ans[-1] + fac[N]) % mod
ans.append(cum)
return ans
N, Ss = read_data()
ans = solve(N, Ss)
for a in ans[1:]:
print(a)