結果
問題 | No.1662 (ox) Alternative |
ユーザー | 👑 ygussany |
提出日時 | 2021-08-09 10:10:00 |
言語 | C (gcc 12.3.0) |
結果 |
AC
|
実行時間 | 96 ms / 2,000 ms |
コード長 | 1,018 bytes |
コンパイル時間 | 199 ms |
コンパイル使用メモリ | 29,952 KB |
実行使用メモリ | 8,116 KB |
最終ジャッジ日時 | 2023-08-13 07:26:06 |
合計ジャッジ時間 | 948 ms |
ジャッジサーバーID (参考情報) |
judge13 / judge12 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 7 ms
8,092 KB |
testcase_01 | AC | 7 ms
8,080 KB |
testcase_02 | AC | 80 ms
8,116 KB |
testcase_03 | AC | 96 ms
8,092 KB |
ソースコード
#include <stdio.h> const int Mod = 1000000007; long long fact[400001], fact_inv[400001]; long long div_mod(long long x, long long y, long long z) { if (x % y == 0) return x / y; else return (div_mod((1 + x / y) * y - x, (z % y), y) * z + x) / y; } long long combination(int n, int k) { return fact[n] * fact_inv[k] % Mod * fact_inv[n-k] % Mod; } long long solve(int A, int B, int C, int D) { if (A != B) return 0; else if (A == 0) { if (D == 0) return 1; else return 0; } else return div_mod(combination(A + D - 1, D) * combination(A * 2 + D, A - 1) % Mod, A, Mod) * combination(A + B + C + D, C) % Mod; } int main() { int i, j, T, A, B, C, D; for (i = 1, fact[0] = 1; i <= 400000; i++) fact[i] = fact[i-1] * i % Mod; for (i = 399999, fact_inv[400000] = div_mod(1, fact[400000], Mod); i >= 0; i--) fact_inv[i] = fact_inv[i+1] * (i + 1) % Mod; scanf("%d", &T); for (i = 1; i <= T; i++) { scanf("%d %d %d %d", &A, &B, &C, &D); printf("%lld\n", solve(A, B, C, D)); } fflush(stdout); return 0; }