結果
問題 | No.2122 黄金比で擬似乱数生成 |
ユーザー | googol_S0 |
提出日時 | 2022-08-10 15:17:42 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 697 ms / 2,000 ms |
コード長 | 957 bytes |
コンパイル時間 | 291 ms |
コンパイル使用メモリ | 87,068 KB |
実行使用メモリ | 83,772 KB |
最終ジャッジ日時 | 2023-09-26 02:48:22 |
合計ジャッジ時間 | 8,107 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge11 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 107 ms
82,224 KB |
testcase_01 | AC | 174 ms
82,720 KB |
testcase_02 | AC | 170 ms
83,172 KB |
testcase_03 | AC | 171 ms
82,496 KB |
testcase_04 | AC | 169 ms
82,980 KB |
testcase_05 | AC | 185 ms
82,580 KB |
testcase_06 | AC | 170 ms
83,192 KB |
testcase_07 | AC | 172 ms
82,480 KB |
testcase_08 | AC | 180 ms
83,004 KB |
testcase_09 | AC | 174 ms
82,784 KB |
testcase_10 | AC | 185 ms
82,544 KB |
testcase_11 | AC | 150 ms
83,048 KB |
testcase_12 | AC | 186 ms
83,388 KB |
testcase_13 | AC | 187 ms
82,876 KB |
testcase_14 | AC | 192 ms
82,740 KB |
testcase_15 | AC | 228 ms
82,968 KB |
testcase_16 | AC | 335 ms
83,408 KB |
testcase_17 | AC | 183 ms
83,156 KB |
testcase_18 | AC | 108 ms
82,132 KB |
testcase_19 | AC | 185 ms
83,508 KB |
testcase_20 | AC | 108 ms
82,144 KB |
testcase_21 | AC | 107 ms
82,132 KB |
testcase_22 | AC | 681 ms
82,948 KB |
testcase_23 | AC | 697 ms
83,240 KB |
testcase_24 | AC | 599 ms
83,772 KB |
testcase_25 | AC | 672 ms
83,484 KB |
ソースコード
S=input() P=int(input()) I=int(input()) assert(len(S)==4) for i in range(4): assert('0'<=S[i]) assert(S[i]<='9') assert(0<=P) assert(P<=(10**18)) assert(0<=I) assert(I<=(10**18)) N=int(S) mod=10**9+7 def m(a,b): r=[[0]*len(b[0]) for i in range(len(a))] for i in range(len(a)): for k in range(len(b)): for j in range(len(b[0])): r[i][j]=(r[i][j]+a[i][k]*b[k][j])%mod return r def p(a,n): r=[[0]*len(a) for i in range(len(a))] b=[] for i in range(len(a)): r[i][i]=1 b.append(a[i][:]) l=n while l>0: if l&1: r=m(b,r) b=m(b,b) l>>=1 return r def f(n): global mod if n==0: return 0 mod=10000 A=[[n,1],[1,0]] A=p(A,P) if P&1: A[0][1]-=1 return A[0][1]%mod T=[[0]*10000 for i in range(62)] T[0]=[f(i) for i in range(10000)] for i in range(61): for j in range(10000): T[i+1][j]=T[i][T[i][j]] for i in range(62): if (I>>i)&1: N=T[i][N] print(str(N).zfill(4))