結果
問題 | No.1043 直列大学 |
ユーザー |
|
提出日時 | 2020-05-01 22:12:58 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 80 ms / 2,000 ms |
コード長 | 1,393 bytes |
コンパイル時間 | 1,061 ms |
コンパイル使用メモリ | 94,584 KB |
最終ジャッジ日時 | 2025-01-10 04:59:20 |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 28 |
ソースコード
#include <iostream> #include <vector> #include <algorithm> #include <queue> #include <stack> #include <deque> #include <math.h> using namespace std; typedef long long ll; constexpr ll mod=1e9+7; ll dp[2][102100]; ll dp2[2][102101]; int main(){ cin.tie(nullptr); ios::sync_with_stdio(false); int n,m; cin >> n >> m; vector<int> v(n),r(m); for(int i=0;i<n;i++){ cin >> v[i]; } for(int i=0;i<m;i++){ cin >> r[i]; } int pre=0; int nex=1; dp[0][0]=1; int a,b; cin >> a >> b; for(int i=0;i<n;i++){ for(int j=0;j<=100000;j++){ (dp[nex][j+v[i]]+=dp[pre][j])%=mod; (dp[nex][j]+=dp[pre][j])%=mod; dp[pre][j]=0; } swap(pre,nex); } int ppre=0; int nnex=1; dp2[0][0]=1; for(int i=0;i<m;i++){ for(int j=0;j<=100000;j++){ (dp2[nnex][j+r[i]]+=dp2[ppre][j])%=mod; (dp2[nnex][j]+=dp2[ppre][j])%=mod; dp2[ppre][j]=0; } swap(nnex,ppre); } for(int i=0;i<=100100;i++){ (dp[pre][i+1]+=dp[pre][i])%=mod; } ll res=0; for(int i=1;i<=100000;i++){ if(dp2[ppre][i]==0)continue; ll aa=min((ll)a*i,100001LL); ll bb=min((ll)b*i,100001LL); (res+=(dp[pre][bb]-dp[pre][aa-1])%mod*dp2[ppre][i])%=mod; } if(res<0)res+=mod; printf("%lld\n",res); }