結果
| 問題 | 
                            No.1035 Color Box
                             | 
                    
| コンテスト | |
| ユーザー | 
                             | 
                    
| 提出日時 | 2022-08-30 19:58:41 | 
| 言語 | PyPy3  (7.3.15)  | 
                    
| 結果 | 
                             
                                WA
                                 
                             
                            
                         | 
                    
| 実行時間 | - | 
| コード長 | 822 bytes | 
| コンパイル時間 | 163 ms | 
| コンパイル使用メモリ | 82,048 KB | 
| 実行使用メモリ | 84,096 KB | 
| 最終ジャッジ日時 | 2024-11-07 13:57:03 | 
| 合計ジャッジ時間 | 3,317 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge5 / judge4 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 2 | 
| other | AC * 6 WA * 30 | 
ソースコード
MOD = 10**9+7
mod = MOD
n,m = map(int,input().split())
N = n
def cmb(n, r, mod):
    if ( r<0 or r>n ):
        return 0
    r = min(r, n-r)
    return g1[n] * g2[r] * g2[n-r] % mod
g1 = [1, 1] # 元テーブル
g2 = [1, 1] #逆元テーブル
inverse = [0, 1] #逆元テーブル計算用テーブル
for i in range( 2, N + 1 ):
    g1.append( ( g1[-1] * i ) % mod )
    inverse.append( ( -inverse[mod % i] * (mod//i) ) % mod )
    g2.append( (g2[-1] * inverse[-1]) % mod )
fact = [1]*(N+1)
rfact = [1]*(N+1)
r = 1
for i in range(1, N+1):
    fact[i] = r = r * i % MOD
rfact[N] = r = pow(fact[N], MOD-2, MOD)
for i in range(N, 0, -1):
    rfact[i-1] = r = r * i % MOD
# nPk (mod MOD) を求める
def perm(n, k):
    return fact[n] * rfact[n-k] % MOD
ans = perm(m,m)
if n != m:
    ans *= cmb(n,n-m,mod)
print(ans)