結果
問題 |
No.1043 直列大学
|
ユーザー |
|
提出日時 | 2020-05-08 15:48:48 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 576 ms / 2,000 ms |
コード長 | 983 bytes |
コンパイル時間 | 237 ms |
コンパイル使用メモリ | 82,288 KB |
実行使用メモリ | 238,976 KB |
最終ジャッジ日時 | 2024-12-22 20:05:54 |
合計ジャッジ時間 | 11,498 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 28 |
ソースコード
import sys input = sys.stdin.readline N, M = map(int, input().split()) V = list(map(int, input().split())) R = list(map(int, input().split())) A, B = map(int, input().split()) MOD = 10**9+7 dp1 = [[0]*(10**5+1) for _ in range(N+1)] dp1[0][0] = 1 for i in range(N): for j in range(10**5+1): dp1[i+1][j] += dp1[i][j] dp1[i+1][j] %= MOD if j+V[i]<=10**5: dp1[i+1][j+V[i]] += dp1[i][j] dp1[i+1][j+V[i]] %= MOD acc = [0]*(10**5+2) for i in range(10**5+1): acc[i+1] = acc[i]+dp1[N][i] dp2 = [[0]*(10**5+1) for _ in range(M+1)] dp2[0][0] = 1 for i in range(M): for j in range(10**5+1): dp2[i+1][j] += dp2[i][j] dp2[i+1][j] %= MOD if j+R[i]<=10**5: dp2[i+1][j+R[i]] += dp2[i][j] dp2[i+1][j+R[i]] %= MOD ans = 0 for i in range(1, 10**5+1): if i*A>10**5: break ans += dp2[M][i]*(acc[min(10**5, i*B)+1]-acc[i*A]) ans %= MOD print(ans)