結果
問題 | No.1050 Zero (Maximum) |
ユーザー |
![]() |
提出日時 | 2020-05-09 13:14:24 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,058 bytes |
コンパイル時間 | 85 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 20,736 KB |
最終ジャッジ日時 | 2024-07-05 17:48:42 |
合計ジャッジ時間 | 3,940 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | TLE * 1 -- * 14 |
ソースコード
def mat_mul(A,B):m=len(A)ans=[[0]*m for _ in range(m)]for i in range(m):for j in range(m):x=0for k in range(m):x+=A[i][k]*B[k][j]ans[i][j]=xreturn ansdef mat_mod(A,mod):m=len(A)for i in range(m):for j in range(m):A[i][j]%=modreturn Adef mat_pow(A,n,mod):m=len(A)ans=[[0]*m for _ in range(m)]for i in range(m):for j in range(m):if i==j:ans[i][j]=1while n>0:if n%2==1:ans=mat_mul(ans,A)A=mat_mul(A,A)ans=mat_mod(ans,mod)n//=2return ansfrom sys import stdindef main():#入力readline=stdin.readlineM,K=map(int,readline().split())mod=10**9+7c=[[1]*M for _ in range(M)]for i in range(M):for j in range(M):for x in range(M):if i*x%M==j:c[i][j]+=1c=mat_pow(c,K,mod)print(c[0][0])if __name__=="__main__":main()