結果
問題 |
No.1043 直列大学
|
ユーザー |
![]() |
提出日時 | 2020-06-02 09:10:15 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,299 bytes |
コンパイル時間 | 1,596 ms |
コンパイル使用メモリ | 170,456 KB |
実行使用メモリ | 9,344 KB |
最終ジャッジ日時 | 2024-11-23 18:21:11 |
合計ジャッジ時間 | 4,399 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 22 WA * 6 |
ソースコード
#include <bits/stdc++.h> #define rep(i,n) for(int i=(0);i<(n);i++) using namespace std; typedef long long ll; int main(){ cin.tie(0); ios::sync_with_stdio(false); ll MOD = 1e9 + 7; ll n, m; cin >> n >> m; vector<ll> v(n), r(m); rep(i, n) cin >> v[i]; rep(i, m) cin >> r[i]; ll a, b; cin >> a >> b; int max_n = 151515; ll dpv[2][max_n]; rep(i, 2) rep(j, max_n) dpv[i][j] = 0; dpv[0][0] = 1; rep(i, n) { rep(j, max_n) dpv[(i+1)%2][j] = 0; rep(j, max_n){ dpv[(i+1)%2][j] += dpv[i%2][j]; dpv[(i+1)%2][j] %= MOD; if(j + v[i] < max_n){ dpv[(i+1)%2][j+v[i]] += dpv[i%2][j]; dpv[(i+1)%2][j+v[i]] %= MOD; } } } ll dpr[2][max_n]; rep(i, 2) rep(j, max_n) dpr[i][j] = 0; dpr[0][0] = 1; rep(i, m) { rep(j, max_n) dpr[(i+1)%2][j] = 0; rep(j, max_n){ dpr[(i+1)%2][j] += dpr[i%2][j]; dpr[(i+1)%2][j] %= MOD; if(j + r[i] < max_n){ dpr[(i+1)%2][j+r[i]] += dpr[i%2][j]; dpr[(i+1)%2][j+r[i]] %= MOD; } } } vector<ll> sv(max_n+1, 0); sv[1] = dpv[n%2][1]; for(int i = 1; i < max_n; i++) { sv[i+1] = sv[i] + dpv[n%2][i+1]; sv[i+1] %= MOD; } ll ans = 0; for(int r = 1; b * r < max_n; r++){ ll num = sv[b * r] + MOD - sv[a * r - 1]; num %= MOD; ans += dpr[m%2][r] * num % MOD; ans %= MOD; } cout << ans << endl; }