結果

問題 No.391 CODING WAR
ユーザー _polarbear08_polarbear08
提出日時 2019-06-10 08:09:57
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 679 ms / 2,000 ms
コード長 1,460 bytes
コンパイル時間 141 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 22,656 KB
最終ジャッジ日時 2024-10-06 00:30:55
合計ジャッジ時間 6,037 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 31 ms
10,752 KB
testcase_01 AC 29 ms
10,752 KB
testcase_02 AC 30 ms
10,752 KB
testcase_03 AC 30 ms
10,752 KB
testcase_04 AC 29 ms
10,752 KB
testcase_05 AC 30 ms
10,752 KB
testcase_06 AC 29 ms
10,752 KB
testcase_07 AC 30 ms
10,752 KB
testcase_08 AC 29 ms
10,880 KB
testcase_09 AC 679 ms
22,656 KB
testcase_10 AC 549 ms
22,528 KB
testcase_11 AC 448 ms
22,528 KB
testcase_12 AC 28 ms
10,752 KB
testcase_13 AC 607 ms
22,656 KB
testcase_14 AC 529 ms
19,584 KB
testcase_15 AC 596 ms
20,992 KB
testcase_16 AC 370 ms
16,640 KB
testcase_17 AC 407 ms
17,920 KB
testcase_18 AC 306 ms
15,488 KB
testcase_19 AC 302 ms
15,616 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
stdin = sys.stdin

sys.setrecursionlimit(10 ** 7)

def li(): return map(int, stdin.readline().split())
def li_(): return map(lambda x: int(x) - 1, stdin.readline().split())
def lf(): return map(float, stdin.readline().split())
def ls(): return stdin.readline().split()
def ns(): return stdin.readline().rstrip()
def lc(): return list(ns())
def ni(): return int(stdin.readline())
def nf(): return float(stdin.readline())


# nの逆元のリスト
def inv_mod(n: int, mod: int) -> list:
    inv = [0, 1]
    for i in range(2, n + 1):
        inv.append(mod - ((mod // i) * inv[mod % i]) % mod)
    return inv


# nの階乗のリスト
def fact(n: int, mod: int) -> list:
    fac = [1, 1]
    res = 1
    for i in range(2, n + 1):
        res = res * i % mod
        fac.append(res)
    return fac


# nの階乗の逆元のリスト
def fact_inv(n: int, inv: list, mod: int) -> list:
    facInv = [1, 1]
    for i in range(2, n + 1):
        facInv.append(facInv[i - 1] * inv[i] % mod)
    return facInv


# 二項係数
def nCr(n: int, r: int, mod: int, fac: list, facInv: list) -> int:
    if not (0 <= r and r <= n):
        return 0

    return ((fac[n] * facInv[r]) % mod) * facInv[n - r] % mod


n,m = li()
MOD = 10**9 + 7

inv = inv_mod(m, MOD)
fac = fact(m, MOD)
fac_inv = fact_inv(m, inv, MOD)

ans = 0

for i in range(m+1):
    ans += (-1)**(i%2) * nCr(m, i, MOD, fac, fac_inv) * pow(m-i, n, MOD)

    ans += MOD
    ans %= MOD

print(ans)
0