結果

問題 No.895 MESE
コンテスト
ユーザー vjudge1
提出日時 2026-10-16 20:04:57
言語 C++14
(gcc 15.3.0 + boost 1.92.0 + ACL)
コンパイル:
g++-15 -O2 -lm -std=c++14 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 6 ms / 2,000 ms
+ 481µs
コード長 1,115 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,036 ms
コンパイル使用メモリ 185,004 KB
実行使用メモリ 17,240 KB
最終ジャッジ日時 2026-10-16 20:05:06
合計ジャッジ時間 3,325 ms
ジャッジサーバーID
（参考情報） 		judge5_0 / judge3_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;
const int MD = 1e9 + 7;
const int MX = 3000005;
long long f[MX], g[MX], p[MX];
long long pw(long long x, long long y) {
    long long r = 1;
    while (y) {
        if (y & 1) r = r * x % MD;
        x = x * x % MD;
        y >>= 1;
    }
    return r;
}
long long C(int n, int m) {
    if (m < 0 || m > n) return 0;
    return f[n] * g[m] % MD * g[n - m] % MD;
}
int main() {
    ios::sync_with_stdio(0);
    cin.tie(0);
    int a, b, c;
    cin >> a >> b >> c;
    int n = a + b + c;
    f[0] = 1;
    p[0] = 1;
    for (int i = 1; i <= n; ++i) {
        f[i] = f[i - 1] * i % MD;
        p[i] = p[i - 1] * 2 % MD;
    }
    g[n] = pw(f[n], MD - 2);
    for (int i = n - 1; i >= 0; --i) {
        g[i] = g[i + 1] * (i + 1) % MD;
    }
    long long ans = 0;
    long long cf = C(b + c - 2, b - 1);
    for (int k = 0; k < a; ++k) {
        int i = b + c - 1 + k;
        long long tm = C(b + c - 2 + k, b + c - 2);
        long long vl = (p[i] - 1 + MD) % MD;
        ans = (ans + tm * vl) % MD;
    }
    ans = ans * cf % MD;
    cout << ans << endl;
    return 0;
}
0