結果

問題 No.1043 直列大学
コンテスト
ユーザー Ricky_pon
提出日時 2019-11-13 17:22:54
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 115 ms / 2,000 ms
コード長 1,514 bytes
コンパイル時間 1,885 ms
コンパイル使用メモリ 167,644 KB
実行使用メモリ 162,304 KB
最終ジャッジ日時 2024-12-22 19:21:34
合計ジャッジ時間 4,877 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define For(i, a, b) for(int (i)=(a); (i)<(b); ++(i))
#define rFor(i, a, b) for(int (i)=(a)-1; (i)>=(b); --(i))
#define rep(i, n) For((i), 0, (n))
#define rrep(i, n) rFor((i), (n), 0)
#define fi first
#define se second
using namespace std;
typedef long long lint;
typedef unsigned long long ulint;
typedef pair<int, int> pii;
typedef pair<lint, lint> pll;
typedef complex<double> xy_t;
typedef vector<lint> poly;

const lint mod = 1e9 + 7;
const lint INF = mod * mod;
const int MAX = 100010;

lint dpv[110][MAX], dpr[110][MAX];

int main(){
    int n, m;
    scanf("%d%d", &n, &m);
    int v[n], r[m];
    rep(i, n) scanf("%d", &v[i]);
    rep(i, m) scanf("%d", &r[i]);
    lint a, b;
    scanf("%lld%lld", &a, &b);

    rep(i, n+1)rep(j, MAX) dpv[i][j] = 0;
    rep(i, m+1)rep(j, MAX) dpr[i][j] = 0;
    dpv[0][0] = 1;
    rep(i, n)rep(j, MAX)if(dpv[i][j]){
        (dpv[i+1][j] += dpv[i][j]) %= mod;
        if(j+v[i] < MAX) (dpv[i+1][j+v[i]] += dpv[i][j]) %= mod;
    }
    dpr[0][0] = 1;
    rep(i, m)rep(j, MAX)if(dpr[i][j]){
        (dpr[i+1][j] += dpr[i][j]) %= mod;
        if(j+r[i] < MAX) (dpr[i+1][j+r[i]] += dpr[i][j]) %= mod;
    }

    lint vsum[MAX];
    vsum[0] = 0;
    rep(i, MAX-1) vsum[i+1] = (vsum[i] + dpv[n][i]) % mod;

    lint ans = 0;
    For(i, 1, MAX)if(dpr[m][i]){
        lint L = a*(lint)i, R = b*(lint)i + 1;
        if(L > MAX) break;
        (ans += (vsum[min(R, (lint)MAX-1)] - vsum[L] + mod) % mod * dpr[m][i]) %= mod;
    }
    printf("%lld\n", ans);
}
0