結果
問題 | No.1520 Zigzag Sum |
ユーザー |
|
提出日時 | 2021-02-03 19:24:54 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 1,021 ms / 2,000 ms |
コード長 | 1,137 bytes |
コンパイル時間 | 84 ms |
コンパイル使用メモリ | 12,800 KB |
実行使用メモリ | 42,220 KB |
最終ジャッジ日時 | 2024-09-25 11:33:03 |
合計ジャッジ時間 | 7,212 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 7 |
ソースコード
mod = 10 ** 9 + 7#順列class Permutation:#前計算def __init__(self, n):self.n = nself.fact = [1] * (self.n + 1)self.fact_inv = [1] * (self.n + 1)for i in range(1, self.n + 1):self.fact[i] = self.fact[i - 1] * i % modself.fact_inv[-1] = invmod(self.fact[-1], mod)for i in range(self.n, 0, -1):self.fact_inv[i - 1] = self.fact_inv[i] * i % moddef perm(self, n, r):if n < r:return 0if n < 0 or r < 0:return 0return self.fact[n] * self.fact_inv[n - r] % mod#拡張ユークリッドの互除法def extgcd(a, b, d = 0):g = aif b == 0:x, y = 1, 0else:x, y, g = extgcd(b, a % b)x, y = y, x - a // b * yreturn x, y, g#mod pにおける逆元def invmod(a, p):x, y, g = extgcd(a, p)x %= preturn xhwmax = 4 * 10 ** 5P = Permutation(hwmax)t = int(input())for _ in range(t):h, w = map(int, input().split())if h == 1 or w == 1:print(0)continueprint(2 * P.perm(h + w - 3, h - 1) % mod * P.fact_inv[h - 2] % mod)