結果

問題 No.895 MESE
ユーザー noisy_noiminnoisy_noimin
提出日時 2019-09-27 23:15:46
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,835 bytes
コンパイル時間 1,896 ms
コンパイル使用メモリ 169,464 KB
実行使用メモリ 8,696 KB
最終ジャッジ日時 2024-09-25 01:44:55
合計ジャッジ時間 6,808 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
6,812 KB
testcase_01 AC 2 ms
6,940 KB
testcase_02 WA -
testcase_03 WA -
testcase_04 AC 3 ms
6,940 KB
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
権限があれば一括ダウンロードができます

ソースコード

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<s;++i) {
        if(i > a) ab_patterns[i] = ab_patterns[i-1];
        ab_patterns[i] = comb(a+b-i-1, a-i);
        ab_patterns[i] %= MOD;
        // cerr << ab_patterns[i] << endl;
    }

    ll ans = 0;
    ll patterns = 0;
    for(int i=s-3;i>=0;--i) {
        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