結果

問題 No.1044 正直者大学
ユーザー yuly3
提出日時 2020-05-02 18:44:52
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 92 ms / 2,000 ms
コード長 1,437 bytes
コンパイル時間 323 ms
コンパイル使用メモリ 82,472 KB
実行使用メモリ 76,168 KB
最終ジャッジ日時 2024-12-31 00:22:36
合計ジャッジ時間 3,479 ms
ジャッジサーバーID
(参考情報)
judge2 / judge4
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys

sys.setrecursionlimit(10 ** 7)
rl = sys.stdin.readline


class Combination:
    def __init__(self, n: int, mod: int):
        self.mod = mod
        self.fact = [0] * (n + 1)
        self.factinv = [0] * (n + 1)
        self.inv = [0] * (n + 1)
        
        self.fact[0] = self.fact[1] = 1
        self.factinv[0] = self.factinv[1] = 1
        self.inv[1] = 1
        for i in range(2, n + 1):
            self.fact[i] = (self.fact[i - 1] * i) % mod
            self.inv[i] = (-self.inv[mod % i] * (mod // i)) % mod
            self.factinv[i] = (self.factinv[i - 1] * self.inv[i]) % mod
    
    def ncr(self, n: int, r: int):
        if r < 0 or n < r:
            return 0
        r = min(r, n - r)
        return self.fact[n] * self.factinv[r] % self.mod * self.factinv[n - r] % self.mod
    
    def nhr(self, n: int, r: int):
        return self.ncr(n + r - 1, r)
    
    def npr(self, n: int, r: int):
        if r < 0 or n < r:
            return 0
        return self.fact[n] * self.factinv[n - r] % self.mod


def solve():
    N, M, K = map(int, rl().split())
    MOD = 10 ** 9 + 7
    com = Combination(200010, MOD)
    
    ans = 0
    for i in range(1, N + 1):
        if N - i + M - i < K:
            continue
        ans = (ans + com.ncr(N, i) * com.nhr(i, M - i)) % MOD
    
    ans = ans * com.fact[N - 1] % MOD
    ans = ans * com.fact[M] % MOD
    print(ans)


if __name__ == '__main__':
    solve()
0