結果

問題 No.997 Jumping Kangaroo
ユーザー nehan_der_thal
提出日時 2020-02-21 23:48:33
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 556 ms / 2,000 ms
コード長 596 bytes
コンパイル時間 301 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 44,348 KB
最終ジャッジ日時 2024-10-09 02:51:46
合計ジャッジ時間 18,303 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff # 

import numpy as np
mod = 10**9+7
N, W, K = map(int, input().split())
A = list(map(int, input().split()))
dp = [0]*(2*W+1)
dp[0] = 1
for i in range(1, 2*W+1):
    for a in A:
        if i >= a:
            dp[i] += dp[i-a]
    dp[i] %= mod
a = dp[W]
b = (dp[2*W] - dp[W]*dp[W])%mod
#print(a, b)
X = [[b, a],[b*a, a**2+b]]
#X = [[1, 0],[1, 0]]
X = np.array(X, dtype="object")
S = [1, a]
S = np.array(S, dtype="object")
#print(S)
#print("aa", np.dot(X, S))

tmp = K%2
K //= 2
for _ in range(60):
    if K%2:
        S = np.dot(X, S)%mod
    K = K//2
    X = np.dot(X, X)%mod
#print(S)
print(S[tmp])
0