結果

問題 No.1050 Zero (Maximum)
コンテスト
ユーザー 已经死了
提出日時 2025-09-13 20:29:01
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 173 ms / 2,000 ms
コード長 892 bytes
コンパイル時間 449 ms
コンパイル使用メモリ 82,656 KB
実行使用メモリ 76,364 KB
最終ジャッジ日時 2025-09-13 20:29:04
合計ジャッジ時間 3,300 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 15
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import sys,math
input=lambda:sys.stdin.readline().rstrip()
MI=lambda:map(int,input().split())
ii=lambda:int(input())
MOD=10**9+7

def mat_mul(A,B,mod):
    n=len(A)
    C=[[0]*n for _ in range(n)]
    for i in range(n):
        Ai=A[i]
        for k in range(n):
            t=Ai[k]
            if t==0: continue
            Bk=B[k]
            for j in range(n):
                C[i][j]=(C[i][j]+t*Bk[j])%mod
    return C

def mat_pow(A,e,mod):
    n=len(A)
    R=[[0]*n for _ in range(n)]
    for i in range(n): R[i][i]=1
    while e:
        if e&1: R=mat_mul(R,A,mod)
        A=mat_mul(A,A,mod)
        e>>=1
    return R

def solve():
    M,K=MI()
    T=[[0]*M for _ in range(M)]
    for x in range(0,M):
        d=math.gcd(x,M)
        for y in range(0,M):
            T[x][y]=1+(d if y%d==0 else 0)
    TK=mat_pow(T,K,MOD)
    print(TK[0][0]%MOD)

if __name__=="__main__":
    solve()
0