結果
問題 | No.3006 ベイカーの問題 |
ユーザー |
![]() |
提出日時 | 2025-01-20 03:34:42 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 30 ms / 2,000 ms |
コード長 | 849 bytes |
コンパイル時間 | 591 ms |
コンパイル使用メモリ | 12,032 KB |
実行使用メモリ | 10,240 KB |
最終ジャッジ日時 | 2025-01-20 03:34:45 |
合計ジャッジ時間 | 2,981 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 24 |
ソースコード
# 行列の計算(numpyを使えないとき,modを使用)def prod(A,B,k,l,m):# A:k*l,B:l*mC=[[None for i in range(m)] for j in range(k)]for i in range(k):for j in range(m):ANS=0for pl in range(l):ANS=(ANS+A[i][pl]*B[pl][j])%modC[i][j]=ANSreturn Cdef plus(A,B,k,l):# a,B:k*lC=[[None for i in range(l)] for j in range(k)]for i in range(k):for j in range(l):C[i][j]=(A[i][j]+B[i][j])%modreturn Cmod=998244353x,y,n=map(int,input().split())POWA=[[[x,y,1,0],[-5*y,x,0,1],[0,0,1,0],[0,0,0,1]]]for i in range(60):POWA.append(prod(POWA[-1],POWA[-1],4,4,4)) # ベキを求めてX=[[x,y,0,0]]while n:X=prod(X,POWA[n.bit_length()-1],1,4,4) # n乗の場合n-=1<<(n.bit_length()-1)print(X[0][2],X[0][3])