#1033846 (C++17(gcc12)) No.940 ワープ ε=ε=ε=ε=ε=│;p>д＜│

• 問題: No.940 ワープ ε=ε=ε=ε=ε=│;p>д＜│
• ユーザー: Caiiiiiiii
• 提出日時: 2025-01-01 01:04:03
• 言語: C++17(gcc12) (gcc 12.3.0 + boost 1.87.0)
• 結果: AC
• 実行時間: 93 ms / 5,000 ms
• コード長: 1,152 bytes
• コンパイル時間: 6,383 ms
• コンパイル使用メモリ: 202,100 KB
• 実行使用メモリ: 42,336 KB
• 最終ジャッジ日時: 2025-01-01 01:04:13
• 合計ジャッジ時間: 5,598 ms
• ジャッジサーバーID （参考情報）: judge4 / judge3

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	55 ms; 35,528 KB
testcase_01	AC	55 ms; 36,056 KB
testcase_02	AC	55 ms; 34,968 KB
testcase_03	AC	54 ms; 34,988 KB
testcase_04	AC	55 ms; 35,864 KB
testcase_05	AC	54 ms; 35,420 KB
testcase_06	AC	54 ms; 35,076 KB
testcase_07	AC	55 ms; 36,188 KB
testcase_08	AC	54 ms; 35,388 KB
testcase_09	AC	54 ms; 35,532 KB
testcase_10	AC	54 ms; 35,496 KB
testcase_11	AC	54 ms; 35,516 KB
testcase_12	AC	54 ms; 36,740 KB
testcase_13	AC	54 ms; 35,064 KB
testcase_14	AC	54 ms; 35,128 KB
testcase_15	AC	57 ms; 35,576 KB
testcase_16	AC	63 ms; 36,132 KB
testcase_17	AC	72 ms; 38,244 KB
testcase_18	AC	78 ms; 40,292 KB
testcase_19	AC	72 ms; 38,364 KB
testcase_20	AC	84 ms; 40,288 KB
testcase_21	AC	67 ms; 38,240 KB
testcase_22	AC	82 ms; 40,288 KB
testcase_23	AC	64 ms; 38,364 KB
testcase_24	AC	83 ms; 40,292 KB
testcase_25	AC	89 ms; 40,288 KB
testcase_26	AC	93 ms; 42,336 KB

ソースコード

#include <bits/stdc++.h>

using LL = long long;

const int N = 4e6 + 7;
const int MOD = 1e9 + 7;

int p, q, r, fac[N], ifc[N];
int a[N];

int pow(int x, int y) {
  int ret = 1;
  while(y) {
    if(y & 1)
      ret = 1LL * ret * x % MOD;
    x = 1LL * x * x % MOD;
    y >>= 1;
  }
  return ret;
}

int com(int x, int y) {
  return 1LL * fac[x] * ifc[y] % MOD * ifc[x - y] % MOD;
}

int main() {

  fac[0] = 1;
  for(int i = 1; i < N; ++i)
    fac[i] = 1LL * fac[i - 1] * i % MOD;
  ifc[N - 1] = pow(fac[N - 1], MOD - 2);
  for(int i = N - 2; i >= 0; --i)
    ifc[i] = 1LL * ifc[i + 1] * (i + 1) % MOD;

  scanf("%d%d%d", &p, &q, &r);

  if(!p && !q && !r) {
    puts("1");
    return 0;
  }
  
  int m = p + q + r;
  a[0] = 1;
  for(int i = 0; i <= m + 1; ++i)
    a[i] = a[i] - (m + 1 - i & 1 ? -1 : 1) * com(m + 1, i);
  for(int i = 0; i <= m; ++i) {
    a[i] = 1LL * a[i] * ifc[2] % MOD;
    a[i + 1] = (a[i + 1] + a[i]) % MOD;
  }
  
  int ans = 0;
  for(int i = 1; i <= m; ++i)
    ans = (ans + 1LL * a[i] * com(i + p - 1, p) % MOD * com(i + q - 1, q) % MOD *
	   com(i + r - 1, r)) % MOD;
  
  printf("%d\n", (ans + MOD) % MOD);
  
  return 0;
}