結果

問題 No.1596 Distance Sum in 2D Plane
ユーザー NatsubiSogan
提出日時 2021-07-11 19:38:02
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 1,454 ms / 2,000 ms
コード長 1,142 bytes
コンパイル時間 141 ms
コンパイル使用メモリ 12,544 KB
実行使用メモリ 42,240 KB
最終ジャッジ日時 2024-07-02 03:11:59
合計ジャッジ時間 18,579 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

diff # 

#nCrなど
mod = 10 ** 9 + 7
class Combinatorics:
    def __init__(self, n):
        self.n = n
        self.fa = [1] * (self.n * 2 + 1)
        self.fi = [1] * (self.n * 2 + 1)

        for i in range(1, self.n * 2 + 1):
            self.fa[i] = self.fa[i - 1] * i % mod

        self.fi[-1] = pow(self.fa[-1], mod - 2, mod)

        for i in range(self.n * 2, 0, -1):
            self.fi[i - 1] = self.fi[i] * i % mod

    def comb(self, n, r):
        if n < r:return 0
        if n < 0 or r < 0:return 0
        return self.fa[n] * self.fi[r] % mod * self.fi[n - r] % mod

    def perm(self, n, r):
        if n < r:return 0
        if n < 0 or r < 0:return 0
        return self.fa[n] * self.fi[n - r] % mod
        
    def combr(self, n, r):
        if n == r == 0:return 1
        return self.comb(n + r - 1, r)

n, m = map(int, input().split())

C = Combinatorics(n)

ans = C.comb(2 * n, n) * n * 2
for _ in range(m):
    t, x, y = map(int, input().split())
    if t == 1:
        ans -= C.comb(x + y, x) * C.comb(n - x - 1 + n - y, n - y)
    else:
        ans -= C.comb(x + y, y) * C.comb(n - x + n - y - 1, n - x)
print(ans % mod)
0