結果
問題 | No.1043 直列大学 |
ユーザー |
![]() |
提出日時 | 2019-11-14 00:47:32 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
RE
|
実行時間 | - |
コード長 | 1,783 bytes |
コンパイル時間 | 1,613 ms |
コンパイル使用メモリ | 169,012 KB |
実行使用メモリ | 175,400 KB |
最終ジャッジ日時 | 2024-11-23 05:18:21 |
合計ジャッジ時間 | 11,711 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 22 RE * 5 TLE * 1 |
ソースコード
#include <bits/stdc++.h> #define For(i, a, b) for(int (i)=(a); (i)<(b); ++(i)) #define rFor(i, a, b) for(int (i)=(a)-1; (i)>=(b); --(i)) #define rep(i, n) For((i), 0, (n)) #define rrep(i, n) rFor((i), (n), 0) #define fi first #define se second using namespace std; typedef long long lint; typedef unsigned long long ulint; typedef pair<int, int> pii; typedef pair<lint, lint> pll; typedef complex<double> xy_t; typedef vector<lint> poly; const lint mod = 1e9 + 7; const lint INF = mod * mod; const int MAX = 100010; lint dpv[110][MAX], dpr[110][MAX]; int main(){ int n, m; scanf("%d%d", &n, &m); int v[n], r[m]; rep(i, n) scanf("%d", &v[i]); rep(i, m) scanf("%d", &r[i]); lint a, b; scanf("%lld%lld", &a, &b); rep(i, n+1)rep(j, MAX) dpv[i][j] = 0; rep(i, m+1)rep(j, MAX) dpr[i][j] = 0; dpv[0][0] = 1; rep(i, n)rep(j, MAX)if(dpv[i][j]){ dpv[i+1][j] += dpv[i][j]; if(dpv[i+1][j] >= mod) dpv[i+1][j] -= mod; if(j+v[i] < MAX){ dpv[i+1][j+v[i]] += dpv[i][j]; if(dpv[i+1][j+v[i]] >= mod) dpv[i+1][j+v[i]] -= mod; } } dpr[0][0] = 1; rep(i, m)rep(j, MAX)if(dpr[i][j]){ dpr[i+1][j] += dpr[i][j]; if(dpr[i+1][j] >= mod) dpr[i+1][j] -= mod; if(j+r[i] < MAX){ dpr[i+1][j+r[i]] += dpr[i][j]; if(dpr[i+1][j+r[i]] >= mod) dpr[i+1][j+r[i]] -= mod; } } lint vsum[MAX]; vsum[0] = 0; rep(i, MAX-1) vsum[i+1] = (vsum[i] + dpv[n][i]) % mod; lint ans = 0; For(i, 1, MAX)if(dpr[m][i]){ lint L = a*(lint)i, R = b*(lint)i + 1; if(L >= MAX) break; For(j, L, R){ ans += dpv[n][j] * dpr[m][i]; if(ans >= mod) ans %= mod; } } printf("%lld\n", ans); }