結果

問題 No.1043 直列大学
ユーザー m4c_physics99
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;
}
0