結果

問題 No.1043 直列大学
ユーザー Ricky_ponRicky_pon
提出日時 2019-11-14 00:47:32
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
RE  
実行時間 -
コード長 1,783 bytes
コンパイル時間 1,613 ms
コンパイル使用メモリ 169,012 KB
実行使用メモリ 175,400 KB
最終ジャッジ日時 2024-11-23 05:18:21
合計ジャッジ時間 11,711 ms
ジャッジサーバーID
(参考情報)
judge3 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 22 RE * 5 TLE * 1
権限があれば一括ダウンロードができます

ソースコード

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];
        if(dpv[i+1][j] >= mod) dpv[i+1][j] -= mod;
        if(j+v[i] < MAX){
            dpv[i+1][j+v[i]] += dpv[i][j];
            if(dpv[i+1][j+v[i]] >= mod) dpv[i+1][j+v[i]] -= mod;
        }
    }
    dpr[0][0] = 1;
    rep(i, m)rep(j, MAX)if(dpr[i][j]){
        dpr[i+1][j] += dpr[i][j];
        if(dpr[i+1][j] >= mod) dpr[i+1][j] -= mod;
        if(j+r[i] < MAX){
            dpr[i+1][j+r[i]] += dpr[i][j];
            if(dpr[i+1][j+r[i]] >= mod) dpr[i+1][j+r[i]] -= 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;
        For(j, L, R){
            ans += dpv[n][j] * dpr[m][i];
            if(ans >= mod) ans %= mod;
        }
    }
    printf("%lld\n", ans);
}
0