結果

問題 No.895 MESE
ユーザー 0w1
提出日時 2019-09-28 03:50:30
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 31 ms / 2,000 ms
コード長 1,028 bytes
コンパイル時間 2,181 ms
コンパイル使用メモリ 196,752 KB
最終ジャッジ日時 2025-01-07 19:43:26
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

const int MOD = 1e9 + 7;

int main() {
  int A, B, C;
  cin >> A >> B >> C;

  int n = A + B + C;
  vector<int> pow2(n + 1);
  for (int i = pow2[0] = 1; i <= n; ++i) {
    pow2[i] = pow2[i - 1] * 2 % MOD;
  }

  auto ipow = [](int v, int n, int m) {
    int r = 1;
    for (int i = n; i; i >>= 1) {
      if (i & 1) r = 1LL * r * v % m;
      v = 1LL * v * v % m;
    }
    return r;
  };

  vector<int> fac(n + 1);
  for (int i = fac[0] = 1; i <= n; ++i) {
    fac[i] = 1LL * fac[i - 1] * i % MOD;
  }

  vector<int> ifac(n + 1);
  ifac[n] = ipow(fac[n], MOD - 2, MOD);
  for (int i = n; i; --i) {
    ifac[i - 1] = 1LL * ifac[i] * i % MOD;
  }

  int ans = 0;
  for (int i = 1; i <= A; ++i) {
    int a = A - i;
    int b = B - 1;
    int c = C;
    int p = 1LL * fac[a + b + c] * ifac[a] % MOD * ifac[b] % MOD * ifac[c] % MOD;
    (ans += 1LL * (pow2[n - 1 - i] - 1) * ipow(n - i - 1, MOD - 2, MOD) % MOD * c % MOD * p % MOD) %= MOD;
  }

  cout << ans << endl;

  return 0;
}
0