結果
問題 |
No.1043 直列大学
|
ユーザー |
![]() |
提出日時 | 2020-05-02 00:39:22 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 73 ms / 2,000 ms |
コード長 | 1,293 bytes |
コンパイル時間 | 1,863 ms |
コンパイル使用メモリ | 170,296 KB |
実行使用メモリ | 5,504 KB |
最終ジャッジ日時 | 2024-12-22 19:57:43 |
合計ジャッジ時間 | 3,371 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 28 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define rep(i,n) for(int i=0; i<(n); i++) using ll = long long; using P = pair<int,int>; const int mod = 1000000007; int main() { int n,m; cin >>n >> m; vector<ll> v(n),r(m); ll V = 0; rep(i,n){ cin >> v[i]; V += v[i]; } ll R = 0; rep(i,m){ cin >> r[i]; R += r[i]; } ll a,b; cin >> a >> b; vector<ll> dpv(V+1),dpr(R+1); dpv[0] = 1; dpr[0] = 1; rep(i,n)for(ll j = V;j >= 0;--j){ if(j+ v[i] > V)continue; dpv[j+v[i]] += dpv[j]; dpv[j+v[i]] %= mod; if(dpv[j+v[i]] <0)dpv[j+v[i]] += mod; } rep(i,m)for(ll j = R;j >= 0;--j){ if(j+ r[i] > R)continue; dpr[j+r[i]] += dpr[j]; dpr[j+r[i]] %= mod; if(dpr[j+r[i]] <0)dpr[j+r[i]] += mod; } vector<ll> smv(V+2,0); rep(i,V+1){ smv[i+1] = smv[i] + dpv[i]; smv[i+1] %= mod; if (smv[i+1] < 0)smv[i+1] += mod; } ll ans = 0; rep(i,R+1){ if(!i)continue; if(!dpr[i])continue; ll vmin = a*i; ll vmax = b*i; ans += (smv[min(V+1,vmax+1)] - smv[min(V+1,vmin)])*dpr[i]; ans %= mod; if( ans < 0)ans += mod; } cout << ans << endl; return 0; }