結果

問題 No.1474 かさまＪ
ユーザー se1ka2se1ka2
提出日時 2021-04-09 23:13:28
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 1,703 ms / 2,500 ms
コード長 4,372 bytes
コンパイル時間 1,071 ms
コンパイル使用メモリ 78,368 KB
実行使用メモリ 201,240 KB
最終ジャッジ日時 2023-09-07 12:59:30
合計ジャッジ時間 9,094 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
7,560 KB
testcase_01 AC 2 ms
7,404 KB
testcase_02 AC 8 ms
30,116 KB
testcase_03 AC 2 ms
7,584 KB
testcase_04 AC 4 ms
13,624 KB
testcase_05 AC 218 ms
104,616 KB
testcase_06 AC 127 ms
111,972 KB
testcase_07 AC 21 ms
95,636 KB
testcase_08 AC 37 ms
165,228 KB
testcase_09 AC 231 ms
168,136 KB
testcase_10 AC 14 ms
24,092 KB
testcase_11 AC 4 ms
5,840 KB
testcase_12 AC 4 ms
7,544 KB
testcase_13 AC 955 ms
156,292 KB
testcase_14 AC 155 ms
175,016 KB
testcase_15 AC 120 ms
163,340 KB
testcase_16 AC 32 ms
16,196 KB
testcase_17 AC 269 ms
71,580 KB
testcase_18 AC 1,396 ms
175,132 KB
testcase_19 AC 1,647 ms
193,400 KB
testcase_20 AC 1,703 ms
201,240 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <functional>
#include <iostream>
#include <vector>
using namespace std;
typedef long long ll;

const ll MOD = 1000000007;

long long modpow(long long x, long long n){
    long long  res = 1;
    long long mul = x;
    while(n){
        if(n % 2) res = res * mul % MOD;
        mul = mul * mul % MOD;
        n /= 2;
    }
    return res;
}

struct Combination
{
    int maxn;
    std::vector<long long> fact;
    std::vector<long long> ifact;
    
    Combination(int maxn) : maxn(maxn){
        fact.resize(maxn + 1);
        ifact.resize(maxn + 1);
        fact[0] = 1;
        for(int i = 1; i <= maxn; i++) fact[i] = fact[i - 1] * i % MOD;
        ifact[maxn] = modpow(fact[maxn], MOD - 2);
        for(int i = maxn - 1; i >= 0; i--) ifact[i] = ifact[i + 1] * (i + 1) % MOD;
    }

    long long com(int n, int k){
        if(n < 0 || k > n || k < 0) return 0;
        if(k == n || k == 0) return 1;
        return (fact[n] * ifact[k] % MOD) * ifact[n - k] % MOD;
    }
};

template <typename T>
struct SegmentTree
{
    using F = function<T(T, T)>;
    
    const F f;
    const T e;
    
    int n;
    vector<T> seg;
    
    SegmentTree(int nn, const F f, const T e) : f(f), e(e){
        n = 1;
        while(n <= nn) n <<= 1;
        seg.assign(n * 2, e);
    }
    
    void set(int i, T x){
        seg[i + n] = x;
    }
    
    void build(){
        for(int k = n - 1; k > 0; k--){
            seg[k] = f(seg[k * 2], seg[k * 2 + 1]);
        }
    }
    
    void update(int i, T x){
        int k = i + n;
        seg[k] = x;
        while(k >>= 1){
            seg[k] = f(seg[k * 2], seg[k * 2 + 1]);
        }
    }
    
    T query(int l, int r){
        l += n, r += n;
        T L = e, R = e;
        for(; l != r; l >>= 1, r >>= 1){
            if(l % 2) L = f(L, seg[l++]);
            if(r % 2) R = f(seg[--r], R);
        }
        return f(L, R);
    }
    
    T operator[](const int i)const {
        return seg[i + n];
    }
    
    template<typename C>
    int right_bound_sub(int k, const C &check, T x){
        while(k < n){
            if(check(f(x, seg[k * 2]))){
                k = k * 2;
            }
            else{
                x = f(x, seg[k * 2]);
                k = k * 2 + 1;
            }
        }
        return k - n;
    }
    
    template <typename C>
    int right_bound(int i, const C &check){
        T x = e;
        for(int l = i + n, r = n * 2; l != r; l >>= 1, r >>= 1){
            if(l % 2){
                if(check(f(x, seg[l]))){
                    return right_bound_sub(l, check, x);
                }
                x = f(x, seg[l]);
                l++;
            }
        }
        return -1;
    }
    
    template <typename C>
    int left_bound_sub(int k, const C &check, T x){
        while(k < n){
            if(check(f(seg[k * 2 + 1], x))){
                k = k * 2 + 1;
            }
            else{
                x = f(seg[k * 2 + 1], x);
                k = k * 2;
            }
        }
        return k - n;
    }
    
    template <typename C>
    int left_bound(int i, const C &check){
        T x = e;
        for(int l = n, r = i + n; l != r; l >>= 1, r >>= 1){
            if(r % 2){
                if(check(f(seg[--r], x))){
                    return left_bound_sub(r, check, x);
                }
                x = f(x, seg[r]);
            }
        }
        return -1;
    }
};

ll dp[42][42][20002];

int main()
{
    int n, mp, mq, l;
    cin >> n >> mp >> mq >> l;
    int s[42];
    for(int i = 0; i < n; i++) cin >> s[i];
    dp[0][0][0] = 1;
    for(int i = 1; i <= n; i++){
        for(int j = 0; j < i; j++){
            for(int k = 0; k <= mq; k++) dp[i][j][k] = dp[i - 1][j][k];
        }
        for(int j = 1; j <= i; j++){
            SegmentTree<ll> seg(mq + 1, [](ll a, ll b){return (a + b) % MOD;}, 0);
            for(int k = 0; k <= mq; k++) seg.set(k, dp[i - 1][j - 1][k]);
            seg.build();
            for(int k = 1; k <= mq; k++){
                dp[i][j][k] = (dp[i][j][k] + seg.query(max(0, k - s[i - 1]), k)) % MOD;
            }
        }
    }
    Combination com(mq + mp + n);
    ll ans = 0;
    for(int j = 0; j <= n; j++){
        for(int k = 0; k <= mq; k++){
            ll x = mp + k - j * l;
            ans = (ans + dp[n][j][k] * com.com(x + n - 1, n - 1)) % MOD;
        }
    }
    cout << ans << endl;
}
0