結果

問題 No.194 フィボナッチ数列の理解(1)
ユーザー qqqqqqqq
提出日時 2019-07-08 17:17:57
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 336 ms / 5,000 ms
コード長 1,418 bytes
コンパイル時間 195 ms
コンパイル使用メモリ 82,408 KB
実行使用メモリ 78,080 KB
最終ジャッジ日時 2024-04-17 08:00:08
合計ジャッジ時間 7,049 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 36 ms
52,224 KB
testcase_01 AC 36 ms
51,968 KB
testcase_02 AC 336 ms
76,760 KB
testcase_03 AC 85 ms
76,224 KB
testcase_04 AC 164 ms
76,284 KB
testcase_05 AC 147 ms
76,288 KB
testcase_06 AC 158 ms
76,056 KB
testcase_07 AC 244 ms
76,268 KB
testcase_08 AC 78 ms
75,904 KB
testcase_09 AC 192 ms
76,172 KB
testcase_10 AC 112 ms
76,416 KB
testcase_11 AC 110 ms
76,032 KB
testcase_12 AC 151 ms
76,504 KB
testcase_13 AC 93 ms
76,288 KB
testcase_14 AC 65 ms
74,164 KB
testcase_15 AC 268 ms
76,544 KB
testcase_16 AC 241 ms
76,544 KB
testcase_17 AC 111 ms
76,384 KB
testcase_18 AC 248 ms
76,288 KB
testcase_19 AC 324 ms
76,672 KB
testcase_20 AC 62 ms
74,368 KB
testcase_21 AC 64 ms
78,080 KB
testcase_22 AC 61 ms
74,880 KB
testcase_23 AC 48 ms
62,848 KB
testcase_24 AC 56 ms
69,248 KB
testcase_25 AC 53 ms
68,352 KB
testcase_26 AC 53 ms
68,608 KB
testcase_27 AC 53 ms
67,584 KB
testcase_28 AC 49 ms
63,104 KB
testcase_29 AC 62 ms
75,392 KB
testcase_30 AC 331 ms
76,540 KB
testcase_31 AC 44 ms
59,776 KB
testcase_32 AC 135 ms
76,268 KB
testcase_33 AC 179 ms
76,160 KB
testcase_34 AC 150 ms
76,268 KB
testcase_35 AC 133 ms
76,328 KB
testcase_36 AC 265 ms
76,240 KB
testcase_37 AC 77 ms
76,032 KB
testcase_38 AC 305 ms
76,492 KB
testcase_39 AC 163 ms
76,416 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

mod = 10**9+7
k, n = map(int, input().split())
a = list(map(int, input().split()))

# 累積和
if n <= 10**6:
    aa = [0]*n
    # 1-indexed
    s = [0]*(n+1)
    for i in range(k):
        aa[i] = a[i]
        s[i+1] = s[i] + a[i]
    for i in range(k, n):
        ai = s[i] - s[i-k]
        aa[i] = ai
        aa[i] %= mod
        s[i+1] = s[i] + ai
        s[i+1] %= mod
    print(aa[-1], s[-1])
# 行列累乗
else:
    def mul(a, b):
        n = len(a)
        c = [[0]*n for i in range(n)]
        for i in range(n):
            for j in range(n):
                for k in range(n):
                    c[i][j] += (a[i][k]*b[k][j])%mod
        return c

    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:
            if n&1:
                b = mul(a, b)
            a = mul(a, a)
            n >>=1
        return b

    A = [[0]*k for i in range(k)]
    A[0] = [1]*k
    for i in range(k-1):
        A[i+1][i] = 1
    A = pow(A, n-k)
    f = 0
    for i in range(k):
        f += A[0][i]*a[k-i-1]
        f %= mod

    B = [[0]*(k+1) for i in range(k+1)]
    B[0][0] = 2
    B[0][-1] = -1
    for i in range(k):
        B[i+1][i] = 1
    B = pow(B, n-k)
    from itertools import accumulate
    s = [0] + list(accumulate(a))
    g = 0
    for i in range(k+1):
        g += B[0][i]*s[k-i]
        g %= mod
    print(f, g)
0