結果

問題 No.1073 無限すごろく
ユーザー qumazakiqumazaki
提出日時 2020-06-06 03:22:19
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 63 ms / 2,000 ms
コード長 795 bytes
コンパイル時間 247 ms
コンパイル使用メモリ 82,516 KB
実行使用メモリ 70,728 KB
最終ジャッジ日時 2024-12-21 09:50:13
合計ジャッジ時間 3,376 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

import copy

n = int(input())
mod = 10**9 + 7

def mat_power(a,b):
    # aもbもm*m列の行列とする。
    m = 6
    res = [[0] *m for _ in range(m)]
    for i in range(m):
        for j in range(m):
            for k in range(m):
                res[i][j] += a[i][k] * b[k][j]
            res[i][j] %= mod

    return res

# x0 = [1,0,0,0,0,0]
a = [[0]*6 for _ in range(6)]
for i in range(6):
    a[0][i] = pow(6,mod-2,mod)

for i in range(5):
    a[i+1][i] = 1

a_exp = []
a_exp.append(a)
i = 1
while(2**i < 10**18):
    ai = mat_power(a_exp[-1],a_exp[-1])
    a_exp.append(copy.deepcopy(ai))
    i += 1

a_n = [[0]*6 for _ in range(6)]
for i in range(6):
    a_n[i][i] = 1

i = 0
while(n>0):
    if(n%2==1):
        a_n = mat_power(a_n,a_exp[i])
    i += 1
    n //= 2

print(a_n[0][0])
0