結果
| 問題 |
No.1044 正直者大学
|
| コンテスト | |
| ユーザー |
mamimume____mo
|
| 提出日時 | 2020-06-11 21:35:53 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 39 ms / 2,000 ms |
| コード長 | 1,523 bytes |
| コンパイル時間 | 841 ms |
| コンパイル使用メモリ | 85,760 KB |
| 実行使用メモリ | 26,880 KB |
| 最終ジャッジ日時 | 2024-06-24 03:21:44 |
| 合計ジャッジ時間 | 2,975 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 28 |
ソースコード
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <iostream>
#include <map>
#include <queue>
#include <deque>
#include <string>
#include <stack>
#include <vector>
#include <set>
#include <tuple>
#include <utility>
#include <functional>
using namespace std;
typedef long long ll;
typedef pair<int,int> P;
typedef tuple<int,int,int> T;
const int INF = 1000000000;
const int MOD = 1000000007;
const long long MAX = 1000100;
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(){
COMinit();
int n,m,k;
cin >> n >> m >> k;
//基本,分ける区間の数lを決めうって,nをl個に分ける場合の数×mをl個に分ける場合の数みたいなことをする
//実際には,oxの区間の組がl個あると/lしないといけない.(123=312=231みたいにローテーションしても同じ
//決め打ちするのは,幾つの区間に分けるか.
ll ans = 0;
ll fact = fac[n] * fac[m] % MOD;
for(int i = 1;i <= min(n,m);i++){
if(n + m - 2 * i >= k){
ll com = COM(n-1,i-1) * COM(m-1,i-1);
com %= MOD;
ans += (com * fact % MOD) * inv[i] % MOD;
ans %= MOD;
}
}
cout << ans << endl;
}
mamimume____mo