結果

問題 No.895 MESE
ユーザー noisy_noiminnoisy_noimin
提出日時 2019-09-27 23:20:48
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 86 ms / 2,000 ms
コード長 1,755 bytes
コンパイル時間 1,707 ms
コンパイル使用メモリ 168,536 KB
実行使用メモリ 9,136 KB
最終ジャッジ日時 2023-10-25 06:36:15
合計ジャッジ時間 4,855 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
7,680 KB
testcase_01 AC 3 ms
7,680 KB
testcase_02 AC 3 ms
7,680 KB
testcase_03 AC 3 ms
7,680 KB
testcase_04 AC 3 ms
7,680 KB
testcase_05 AC 3 ms
7,680 KB
testcase_06 AC 3 ms
7,680 KB
testcase_07 AC 3 ms
7,680 KB
testcase_08 AC 3 ms
7,700 KB
testcase_09 AC 3 ms
7,692 KB
testcase_10 AC 3 ms
7,704 KB
testcase_11 AC 3 ms
7,708 KB
testcase_12 AC 3 ms
7,696 KB
testcase_13 AC 50 ms
8,612 KB
testcase_14 AC 50 ms
9,028 KB
testcase_15 AC 55 ms
9,060 KB
testcase_16 AC 46 ms
8,780 KB
testcase_17 AC 46 ms
8,420 KB
testcase_18 AC 86 ms
9,136 KB
testcase_19 AC 85 ms
9,136 KB
testcase_20 AC 85 ms
9,136 KB
testcase_21 AC 85 ms
9,136 KB
testcase_22 AC 84 ms
9,136 KB
testcase_23 AC 85 ms
9,136 KB
testcase_24 AC 85 ms
9,136 KB
testcase_25 AC 85 ms
9,136 KB
testcase_26 AC 85 ms
9,136 KB
testcase_27 AC 85 ms
9,136 KB
testcase_28 AC 84 ms
9,136 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

using ll =  long long;
using Pll = pair<ll, ll>;
using Pii = pair<int, int>;

constexpr ll MOD = 1000000007;
constexpr long double EPS = 1e-10;
constexpr int dyx[4][2] = {
    { 0, 1}, {-1, 0}, {0,-1}, {1, 0}
};
constexpr int S_MAX = 300000;

ll fact[S_MAX+1], rfact[S_MAX+1];

ll perm(ll n, ll r){
    return (fact[n] * rfact[r]) % MOD;
}

ll comb(ll n, ll r){
    return (perm(n, r) * rfact[n-r]) % MOD;
}

void init(ll n){
    fact[0] = fact[1] = 1;
    rfact[0] = rfact[1] = 1;
    for(int i=2;i<=n;++i) {
        fact[i] = (fact[i-1] * (ll)i) % MOD;
        rfact[i] = 1;
        ll k = MOD-2;
        ll a = fact[i];
        while(k > 0){
            if(k & 1){
                rfact[i] *= a;
                rfact[i] %= MOD;
            }
            a *= a;
            a %= MOD;
            k  >>= 1;
        }
    }
}

ll modpow(ll a, ll t) {
    ll ret = 1LL;
    while(t){
        if(t & 1LL){
            ret *= a;
            ret %= MOD;
        }
        a *= a;
        a %= MOD;
        t >>= 1;
    }
    return ret;
}

int main() {
    ios::sync_with_stdio(false); cin.tie(nullptr);
    ll a, b, c;
    cin >> a >> b >> c;
    ll s = a + b + c;
    init(s);

    ll ab_patterns[s]; // 先頭 i 文字以下が a
    ab_patterns[0] = 0;
    for(int i=1;i<=a;++i) {
        ab_patterns[i] = comb(a+b-i-1, a-i);
        ab_patterns[i] %= MOD;
    }

    ll ans = 0;
    ll patterns = 0;
    for(int i=s-3;i>=0;--i) {
        if(s-i-2 <= a) patterns += (ab_patterns[s-i-2] * comb(s-min(a, s-i-2)-1-1, c-1)) % MOD;
        patterns %= MOD;
        ans += (patterns * modpow(2LL, i)) % MOD;
        // cerr << i << ": " << patterns << endl;
        ans %= MOD;
    }
    cout << ans << endl;
}
0