結果
問題 | No.997 Jumping Kangaroo |
ユーザー | maspy |
提出日時 | 2020-03-20 02:00:08 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
AC
|
実行時間 | 17 ms / 2,000 ms |
コード長 | 944 bytes |
コンパイル時間 | 156 ms |
コンパイル使用メモリ | 10,888 KB |
実行使用メモリ | 7,952 KB |
最終ジャッジ日時 | 2023-08-20 21:20:55 |
合計ジャッジ時間 | 1,710 ms |
ジャッジサーバーID (参考情報) |
judge14 / judge11 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 17 ms
7,816 KB |
testcase_01 | AC | 17 ms
7,812 KB |
testcase_02 | AC | 17 ms
7,808 KB |
testcase_03 | AC | 16 ms
7,868 KB |
testcase_04 | AC | 16 ms
7,924 KB |
testcase_05 | AC | 16 ms
7,948 KB |
testcase_06 | AC | 16 ms
7,812 KB |
testcase_07 | AC | 16 ms
7,812 KB |
testcase_08 | AC | 16 ms
7,804 KB |
testcase_09 | AC | 16 ms
7,920 KB |
testcase_10 | AC | 16 ms
7,856 KB |
testcase_11 | AC | 16 ms
7,816 KB |
testcase_12 | AC | 16 ms
7,816 KB |
testcase_13 | AC | 16 ms
7,800 KB |
testcase_14 | AC | 16 ms
7,812 KB |
testcase_15 | AC | 16 ms
7,948 KB |
testcase_16 | AC | 16 ms
7,848 KB |
testcase_17 | AC | 16 ms
7,944 KB |
testcase_18 | AC | 16 ms
7,872 KB |
testcase_19 | AC | 16 ms
7,768 KB |
testcase_20 | AC | 17 ms
7,796 KB |
testcase_21 | AC | 15 ms
7,804 KB |
testcase_22 | AC | 17 ms
7,812 KB |
testcase_23 | AC | 17 ms
7,952 KB |
testcase_24 | AC | 16 ms
7,860 KB |
testcase_25 | AC | 15 ms
7,816 KB |
testcase_26 | AC | 16 ms
7,948 KB |
testcase_27 | AC | 16 ms
7,812 KB |
ソースコード
#!/usr/bin/ python3.8 import sys read = sys.stdin.buffer.read readline = sys.stdin.buffer.readline readlines = sys.stdin.buffer.readlines MOD = 10**9 + 7 N, W, K, *A = map(int, read().split()) A.sort() dp = [0] * (W + 1) dp[0] = 1 for i in range(W): for a in A: j = i + a if j > W: break dp[j] += dp[i] P = dp[-1] % MOD dp = [0] * (W + W + 1) dp[0] = 1 for i in range(W + W): if i == W: dp[i] = 0 continue for a in A: j = i + a if j > W + W: break dp[j] += dp[i] Q = dp[-1] % MOD # a_{n+2} = Pa_{n+1} + Qa_n # a_{-1} = 0, a_0 = 1 # a_K = ? def power(n): """x^n mod (x^2 - Px - Q)""" if n == 0: return 1, 0 a, b = power(n // 2) a, b = a * a + Q * b * b, 2 * a * b + P * b * b if n % 2 == 0: return a % MOD, b % MOD a, b = Q * b, a + P * b return a % MOD, b % MOD a, b = power(K + 1) print(b)