結果
問題 | No.2122 黄金比で擬似乱数生成 |
ユーザー | 👑 rin204 |
提出日時 | 2022-11-04 23:20:13 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,734 bytes |
コンパイル時間 | 293 ms |
コンパイル使用メモリ | 87,200 KB |
実行使用メモリ | 80,192 KB |
最終ジャッジ日時 | 2023-09-26 01:48:01 |
合計ジャッジ時間 | 9,452 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge11 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 98 ms
76,984 KB |
testcase_01 | AC | 197 ms
78,520 KB |
testcase_02 | AC | 189 ms
78,028 KB |
testcase_03 | AC | 193 ms
78,580 KB |
testcase_04 | AC | 205 ms
78,680 KB |
testcase_05 | AC | 226 ms
79,488 KB |
testcase_06 | AC | 185 ms
78,300 KB |
testcase_07 | AC | 199 ms
78,444 KB |
testcase_08 | AC | 216 ms
78,700 KB |
testcase_09 | AC | 204 ms
78,568 KB |
testcase_10 | AC | 218 ms
78,128 KB |
testcase_11 | AC | 163 ms
78,324 KB |
testcase_12 | AC | 217 ms
79,116 KB |
testcase_13 | AC | 225 ms
78,492 KB |
testcase_14 | AC | 250 ms
80,192 KB |
testcase_15 | AC | 269 ms
78,816 KB |
testcase_16 | AC | 398 ms
78,748 KB |
testcase_17 | AC | 216 ms
78,972 KB |
testcase_18 | AC | 99 ms
76,996 KB |
testcase_19 | AC | 221 ms
78,868 KB |
testcase_20 | AC | 99 ms
77,004 KB |
testcase_21 | AC | 101 ms
76,996 KB |
testcase_22 | WA | - |
testcase_23 | WA | - |
testcase_24 | AC | 774 ms
79,372 KB |
testcase_25 | AC | 806 ms
79,116 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 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(2, 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 if m == 2: nex[1] = 1 elif m == 4: nex[1] = 3 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)