結果
問題 | No.1044 正直者大学 |
ユーザー |
|
提出日時 | 2020-05-06 09:08:30 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 23 ms / 2,000 ms |
コード長 | 921 bytes |
コンパイル時間 | 2,158 ms |
コンパイル使用メモリ | 193,608 KB |
最終ジャッジ日時 | 2025-01-10 07:09:50 |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 28 |
ソースコード
typedef long long ll; typedef long double ld; #include <bits/stdc++.h> using namespace std; ll const MOD = 1e9+7; const int MAX = 510000; long long fac[MAX], finv[MAX], inv[MAX]; // テーブルを作る前処理 void COMinit() { fac[0] = fac[1] = 1; finv[0] = finv[1] = 1; inv[1] = 1; for (int i = 2; i < MAX; i++){ fac[i] = fac[i - 1] * i % MOD; inv[i] = MOD - inv[MOD%i] * (MOD / i) % MOD; finv[i] = finv[i - 1] * inv[i] % MOD; } } // 二項係数計算 long long COM(int n, int k){ if (n < k) return 0; if (n < 0 || k < 0) return 0; return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD; } int main() { ll n,m,k; COMinit(); std::cin >> n>>m>>k; ll ans = 0; for (int i = 0; n+m-2*i>=k; i++) { ans += COM(n,i)*COM(m-1,m-i)%MOD*fac[n-1]%MOD*fac[m]%MOD; ans %= MOD; } std::cout << ans << std::endl; }