結果

問題 No.1073 無限すごろく
ユーザー brthyyjpbrthyyjp
提出日時 2020-06-05 23:23:13
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 78 ms / 2,000 ms
コード長 1,037 bytes
コンパイル時間 991 ms
コンパイル使用メモリ 86,920 KB
実行使用メモリ 76,372 KB
最終ジャッジ日時 2023-08-22 17:55:04
合計ジャッジ時間 5,045 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 71 ms
71,432 KB
testcase_01 AC 70 ms
71,316 KB
testcase_02 AC 74 ms
76,184 KB
testcase_03 AC 69 ms
71,276 KB
testcase_04 AC 69 ms
71,212 KB
testcase_05 AC 67 ms
71,128 KB
testcase_06 AC 69 ms
71,472 KB
testcase_07 AC 73 ms
75,748 KB
testcase_08 AC 74 ms
75,472 KB
testcase_09 AC 72 ms
75,816 KB
testcase_10 AC 74 ms
75,660 KB
testcase_11 AC 74 ms
75,452 KB
testcase_12 AC 72 ms
75,412 KB
testcase_13 AC 75 ms
76,120 KB
testcase_14 AC 75 ms
76,372 KB
testcase_15 AC 75 ms
76,276 KB
testcase_16 AC 75 ms
76,128 KB
testcase_17 AC 75 ms
76,256 KB
testcase_18 AC 76 ms
76,280 KB
testcase_19 AC 76 ms
76,256 KB
testcase_20 AC 77 ms
75,876 KB
testcase_21 AC 76 ms
76,292 KB
testcase_22 AC 75 ms
76,228 KB
testcase_23 AC 76 ms
76,236 KB
testcase_24 AC 77 ms
76,280 KB
testcase_25 AC 77 ms
76,268 KB
testcase_26 AC 77 ms
76,296 KB
testcase_27 AC 78 ms
76,124 KB
testcase_28 AC 77 ms
76,028 KB
testcase_29 AC 76 ms
76,292 KB
testcase_30 AC 78 ms
76,060 KB
testcase_31 AC 77 ms
75,892 KB
testcase_32 AC 76 ms
76,372 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