結果

問題 No.1073 無限すごろく
ユーザー brthyyjpbrthyyjp
提出日時 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
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 36 ms
53,264 KB
testcase_01 AC 35 ms
52,696 KB
testcase_02 AC 43 ms
61,076 KB
testcase_03 AC 37 ms
52,744 KB
testcase_04 AC 35 ms
53,156 KB
testcase_05 AC 35 ms
52,600 KB
testcase_06 AC 36 ms
53,552 KB
testcase_07 AC 39 ms
58,652 KB
testcase_08 AC 40 ms
58,780 KB
testcase_09 AC 40 ms
58,584 KB
testcase_10 AC 40 ms
59,328 KB
testcase_11 AC 42 ms
58,696 KB
testcase_12 AC 43 ms
59,628 KB
testcase_13 AC 44 ms
61,500 KB
testcase_14 AC 44 ms
61,356 KB
testcase_15 AC 43 ms
61,576 KB
testcase_16 AC 44 ms
59,928 KB
testcase_17 AC 43 ms
60,052 KB
testcase_18 AC 43 ms
61,952 KB
testcase_19 AC 43 ms
60,520 KB
testcase_20 AC 43 ms
60,492 KB
testcase_21 AC 46 ms
60,144 KB
testcase_22 AC 42 ms
60,972 KB
testcase_23 AC 45 ms
61,948 KB
testcase_24 AC 46 ms
61,768 KB
testcase_25 AC 44 ms
61,520 KB
testcase_26 AC 46 ms
61,996 KB
testcase_27 AC 46 ms
62,112 KB
testcase_28 AC 44 ms
61,196 KB
testcase_29 AC 45 ms
62,320 KB
testcase_30 AC 43 ms
61,520 KB
testcase_31 AC 45 ms
61,564 KB
testcase_32 AC 46 ms
62,440 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