結果
問題 | No.1050 Zero (Maximum) |
ユーザー |
![]() |
提出日時 | 2020-05-08 21:56:18 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 1,405 ms / 2,000 ms |
コード長 | 825 bytes |
コンパイル時間 | 95 ms |
コンパイル使用メモリ | 12,800 KB |
実行使用メモリ | 16,896 KB |
最終ジャッジ日時 | 2024-07-04 00:36:59 |
合計ジャッジ時間 | 10,843 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 15 |
ソースコード
import sys input = sys.stdin.readline M,K=map(int,input().split()) mod=10**9+7 MAT=[[1]*M for i in range(M)] for i in range(M): for j in range(M): MAT[i*j%M][i]+=1 def prod(A,B,k,l,m): C=[[None for i in range(m)] for j in range(k)] for i in range(k): for j in range(m): ANS=0 for pl in range(l): ANS=(ANS+A[i][pl]*B[pl][j])%mod C[i][j]=ANS return C def plus(A,B,k,l): C=[[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])%mod return C X=[[0]*M] X[0][0]=1 POWA=[MAT] for i in range(60): POWA.append(prod(POWA[-1],POWA[-1],M,M,M)) while K: X=prod(X,POWA[K.bit_length()-1],1,M,M) K-=1<<(K.bit_length()-1) print(X[0][0])