結果

問題 No.1043 直列大学
ユーザー Shibuyap
提出日時 2020-07-31 04:54:12
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 72 ms / 2,000 ms
コード長 1,255 bytes
コンパイル時間 2,764 ms
コンパイル使用メモリ 192,536 KB
最終ジャッジ日時 2025-01-12 08:32:45
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define rep(i,n) for(int i = 0; i < (n); ++i)
#define srep(i,s,t) for (int i = s; i < t; ++i)
#define drep(i,n) for(int i = (n)-1; i >= 0; --i)
using namespace std;
typedef long long int ll;
typedef pair<int,int> P;
#define yn {puts("YES");}else{puts("NO");}
#define MAX_N 200005
ll dp2[MAX_N] = {};
ll sum[MAX_N] = {};
ll dp[MAX_N] = {};
int main() {
    int n, m;
    cin >> n >> m;
    int a[n], b[m];
    rep(i,n) cin >> a[i];
    rep(i,m) cin >> b[i];
    ll A, B;
    cin >> A >> B;
    ll MOD = 1000000007;
    
    dp[0] = 1;
    rep(i,n){
        drep(j,MAX_N){
            if(dp[j] == 0) continue;
            dp[j+a[i]] += dp[j];
            dp[j+a[i]] %= MOD;
        }
    }
    
    dp2[0] = 1; 
    rep(i,m){
        drep(j,MAX_N){
            if(dp2[j] == 0) continue;
            dp2[j+b[i]] += dp2[j];
            dp2[j+b[i]] %= MOD;
        }
    }
    
    srep(i,1,MAX_N) sum[i] = sum[i-1] + dp[i];
    ll ans = 0;
    srep(i,1,MAX_N){
        if(dp2[i] == 0) continue;
        ll l = A * i;
        ll r = B * i;
        if(MAX_N <= l) continue;
        r = min(r, (ll)MAX_N-1);
        ans += (sum[r] - sum[l-1]) % MOD * dp2[i] % MOD;
        ans %= MOD;
    }

    cout << ans << endl;
    return 0;
}
0