結果
問題 | No.1043 直列大学 |
ユーザー |
![]() |
提出日時 | 2019-11-13 17:22:54 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 115 ms / 2,000 ms |
コード長 | 1,514 bytes |
コンパイル時間 | 1,885 ms |
コンパイル使用メモリ | 167,644 KB |
実行使用メモリ | 162,304 KB |
最終ジャッジ日時 | 2024-12-22 19:21:34 |
合計ジャッジ時間 | 4,877 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 28 |
ソースコード
#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 secondusing 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]) %= mod;if(j+v[i] < MAX) (dpv[i+1][j+v[i]] += dpv[i][j]) %= mod;}dpr[0][0] = 1;rep(i, m)rep(j, MAX)if(dpr[i][j]){(dpr[i+1][j] += dpr[i][j]) %= mod;if(j+r[i] < MAX) (dpr[i+1][j+r[i]] += dpr[i][j]) %= 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;(ans += (vsum[min(R, (lint)MAX-1)] - vsum[L] + mod) % mod * dpr[m][i]) %= mod;}printf("%lld\n", ans);}