結果
問題 | No.2122 黄金比で擬似乱数生成 |
ユーザー | 👑 rin204 |
提出日時 | 2022-11-04 23:29:02 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 2,339 bytes |
コンパイル時間 | 278 ms |
コンパイル使用メモリ | 86,744 KB |
実行使用メモリ | 80,160 KB |
最終ジャッジ日時 | 2023-09-26 01:54:02 |
合計ジャッジ時間 | 8,837 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge15 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 93 ms
77,024 KB |
testcase_01 | AC | 204 ms
79,532 KB |
testcase_02 | AC | 185 ms
78,200 KB |
testcase_03 | AC | 193 ms
79,768 KB |
testcase_04 | AC | 186 ms
78,688 KB |
testcase_05 | AC | 214 ms
79,136 KB |
testcase_06 | AC | 180 ms
78,180 KB |
testcase_07 | AC | 196 ms
79,788 KB |
testcase_08 | AC | 207 ms
79,492 KB |
testcase_09 | AC | 203 ms
79,248 KB |
testcase_10 | AC | 214 ms
78,736 KB |
testcase_11 | AC | 155 ms
78,168 KB |
testcase_12 | AC | 202 ms
79,312 KB |
testcase_13 | AC | 210 ms
78,484 KB |
testcase_14 | AC | 235 ms
78,924 KB |
testcase_15 | AC | 261 ms
78,588 KB |
testcase_16 | AC | 392 ms
79,424 KB |
testcase_17 | AC | 198 ms
79,148 KB |
testcase_18 | AC | 90 ms
77,152 KB |
testcase_19 | AC | 198 ms
79,180 KB |
testcase_20 | AC | 90 ms
76,980 KB |
testcase_21 | AC | 90 ms
77,004 KB |
testcase_22 | WA | - |
testcase_23 | WA | - |
testcase_24 | AC | 754 ms
79,232 KB |
testcase_25 | AC | 770 ms
79,804 KB |
ソースコード
""" S = 0001 の遷移 わかってない """ MOD = 10000 def add(x, y): if x[0] == 0: return y elif y[0] == 0: return x x = list(x) y = list(y) if x[1] > y[1]: x, y = y, x y[0] <<= y[1] - x[1] z0 = x[0] + y[0] z1 = x[1] return f(z0, z1) def times(x, y): x = list(x) y = list(y) z0 = x[0] * y[0] z1 = x[1] + y[1] return f(z0, z1) def matpow(A, B, w): l = len(A) while w: if w & 1: C = [(0, 0)] * l for i in range(l): for j in range(l): C[i] = add(C[i], times(A[i][j], B[j])) B = C C = [[(0, 0)] * l for _ in range(l)] for i in range(l): for j in range(l): for k in range(l): C[i][j] = add(C[i][j], times(A[i][k], A[k][j])) A = C w >>= 1 return B def matpow_normal(A, B, w): l = len(A) while w: if w & 1: C = [0] * l for i in range(l): for j in range(l): C[i] += A[i][j] * B[j] C[i] %= MOD B = C C = [[0] * l for _ in range(l)] for i in range(l): for j in range(l): for k in range(l): C[i][j] += A[i][k] * A[k][j] C[i][j] %= MOD A = C w >>= 1 return B S = int(input()) m = int(input()) L = int(input()) nex = [0] * 10000 def f(x, r=0): if x == 0: return 0, 0 while x % 2 == 0: r += 1 x //= 2 return x % MOD, r for n in range(1, 10000): d = n * n + 4 B = [f(2), f(0)] A = [[f(n, -1), f(d, -1)], [f(1, -1), f(n, -1)]] ret = matpow(A, B, m) ret = ret[1] nex[n] = ret[0] * pow(2, ret[1], MOD) % MOD if m % 2 == 1: nex[n] -= 1 if nex[n] == -1: nex[n] = 9999 A = [[1, 1], [1, 0]] B = [0, 1] nex[1] = matpow_normal(A, B, m)[0] if m % 2 == 1: nex[1] -= 1 if nex[1] == -1: nex[1] = 9999 n = int(S) se = set() lst = [] while n not in se: se.add(n) lst.append(n) n = nex[n] if len(lst) >= L + 1: ans = str(lst[L]).zfill(4) print(ans) exit() ii = lst.index(n) L -= ii lst = lst[ii:] L %= len(lst) ans = str(lst[L]).zfill(4) print(ans)