結果

問題 No.1073 無限すごろく
コンテスト
ユーザー brthyyjp
提出日時 2020-06-05 23:23:13
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 46 ms / 2,000 ms
コード長 1,037 bytes
コンパイル時間 433 ms
コンパイル使用メモリ 82,200 KB
実行使用メモリ 62,440 KB
最終ジャッジ日時 2024-12-17 18:24:46
合計ジャッジ時間 2,546 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

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

p = 1*pow(6,(mod-2),mod)%mod
#print(a)

p1 = p
p2 = (p + p1*p)%mod
p3 = (p + (p2+p1)*p)%mod
p4 = (p + (p3+p2+p1)*p)%mod
p5 = (p + (p4+p3+p2+p1)*p)%mod
p6 = (p + (p5+p4+p3+p2+p1)*p)%mod
P = [p1, p2, p3, p4, p5, p6]

if n <= 6:
    print(P[n-1])
    exit()

P.reverse()

# A*B
def mul(a, b):
    c = [[0]*len(b[0]) for i in range(len(a))]
    for i in range(len(a)):
        for k in range(len(b)):
            for j in range(len(b[0])):
                c[i][j] = (c[i][j] + a[i][k]*b[k][j])%mod

    return c

#A**n
def pow(a, n):
    b = [[0]*len(a) for i in range(len(a))]
    for i in range(len(a)):
        b[i][i] = 1
    while n > 0:
        if n & 1 == 1:
            b = mul(a, b)
        a = mul(a, a)
        n = n>>1
    return b

A = [[0]*6 for i in range(6)]
for i in range(6):
    A[0][i] = p
for i in range(1, 6):
    A[i][i-1] = 1

#print(A)

B =  pow(A, n-1)
#print(B)
ans = 0
for i in range(6):
    ans += B[5][i]*P[i]
    ans %= mod
#ans *= power(a, n-1, mod)
#ans %= mod
print(ans%mod)
0