結果

問題 No.1073 無限すごろく
ユーザー brthyyjpbrthyyjp
提出日時 2020-06-05 23:23:13
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 49 ms / 2,000 ms
コード長 1,037 bytes
コンパイル時間 154 ms
コンパイル使用メモリ 82,280 KB
実行使用メモリ 61,312 KB
最終ジャッジ日時 2024-05-09 23:43:15
合計ジャッジ時間 2,684 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 39 ms
52,352 KB
testcase_01 AC 38 ms
52,352 KB
testcase_02 AC 46 ms
59,520 KB
testcase_03 AC 38 ms
52,096 KB
testcase_04 AC 38 ms
52,736 KB
testcase_05 AC 38 ms
52,096 KB
testcase_06 AC 38 ms
52,608 KB
testcase_07 AC 43 ms
58,240 KB
testcase_08 AC 43 ms
58,368 KB
testcase_09 AC 43 ms
58,496 KB
testcase_10 AC 42 ms
58,752 KB
testcase_11 AC 43 ms
57,984 KB
testcase_12 AC 42 ms
58,624 KB
testcase_13 AC 44 ms
59,520 KB
testcase_14 AC 46 ms
59,392 KB
testcase_15 AC 46 ms
59,748 KB
testcase_16 AC 45 ms
60,288 KB
testcase_17 AC 45 ms
59,904 KB
testcase_18 AC 45 ms
59,748 KB
testcase_19 AC 46 ms
59,776 KB
testcase_20 AC 46 ms
60,160 KB
testcase_21 AC 46 ms
59,904 KB
testcase_22 AC 45 ms
59,904 KB
testcase_23 AC 47 ms
60,928 KB
testcase_24 AC 49 ms
60,544 KB
testcase_25 AC 47 ms
60,928 KB
testcase_26 AC 47 ms
60,672 KB
testcase_27 AC 47 ms
61,184 KB
testcase_28 AC 48 ms
61,312 KB
testcase_29 AC 48 ms
60,928 KB
testcase_30 AC 46 ms
60,800 KB
testcase_31 AC 48 ms
61,056 KB
testcase_32 AC 47 ms
61,056 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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