結果

問題 No.1785 Inequality Signs
ユーザー NatsubiSoganNatsubiSogan
提出日時 2021-12-12 22:46:38
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 610 ms / 2,000 ms
コード長 1,501 bytes
コンパイル時間 162 ms
コンパイル使用メモリ 81,976 KB
実行使用メモリ 85,172 KB
最終ジャッジ日時 2024-07-22 23:28:29
合計ジャッジ時間 19,719 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 62 ms
67,624 KB
testcase_01 AC 235 ms
80,400 KB
testcase_02 AC 419 ms
84,420 KB
testcase_03 AC 444 ms
81,876 KB
testcase_04 AC 142 ms
78,712 KB
testcase_05 AC 317 ms
82,436 KB
testcase_06 AC 441 ms
82,172 KB
testcase_07 AC 160 ms
78,684 KB
testcase_08 AC 418 ms
81,884 KB
testcase_09 AC 207 ms
79,952 KB
testcase_10 AC 397 ms
84,216 KB
testcase_11 AC 267 ms
81,236 KB
testcase_12 AC 268 ms
81,480 KB
testcase_13 AC 204 ms
79,624 KB
testcase_14 AC 350 ms
83,492 KB
testcase_15 AC 285 ms
81,688 KB
testcase_16 AC 260 ms
80,808 KB
testcase_17 AC 378 ms
83,768 KB
testcase_18 AC 361 ms
83,564 KB
testcase_19 AC 324 ms
82,448 KB
testcase_20 AC 331 ms
83,352 KB
testcase_21 AC 326 ms
82,572 KB
testcase_22 AC 263 ms
80,716 KB
testcase_23 AC 251 ms
80,776 KB
testcase_24 AC 475 ms
81,948 KB
testcase_25 AC 288 ms
81,700 KB
testcase_26 AC 327 ms
81,552 KB
testcase_27 AC 610 ms
82,584 KB
testcase_28 AC 427 ms
82,452 KB
testcase_29 AC 418 ms
82,052 KB
testcase_30 AC 456 ms
82,616 KB
testcase_31 AC 436 ms
82,272 KB
testcase_32 AC 440 ms
82,156 KB
testcase_33 AC 63 ms
69,448 KB
testcase_34 AC 227 ms
79,424 KB
testcase_35 AC 294 ms
80,836 KB
testcase_36 AC 469 ms
82,300 KB
testcase_37 AC 459 ms
82,828 KB
testcase_38 AC 449 ms
83,136 KB
testcase_39 AC 278 ms
82,060 KB
testcase_40 AC 307 ms
82,376 KB
testcase_41 AC 363 ms
85,064 KB
testcase_42 AC 452 ms
82,704 KB
testcase_43 AC 491 ms
83,620 KB
testcase_44 AC 417 ms
82,432 KB
testcase_45 AC 306 ms
82,800 KB
testcase_46 AC 424 ms
81,952 KB
testcase_47 AC 384 ms
84,276 KB
testcase_48 AC 383 ms
85,172 KB
testcase_49 AC 191 ms
79,296 KB
testcase_50 AC 428 ms
82,320 KB
testcase_51 AC 217 ms
79,384 KB
testcase_52 AC 178 ms
79,252 KB
testcase_53 AC 207 ms
80,508 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

class Combinatorics:
	def __init__(self, n: int, mod: int) -> None:
		self.n = n
		self.mod = mod
		self.fa = [1] * (self.n + 1)
		self.fi = [1] * (self.n + 1)

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

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

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

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

	def perm(self, n: int, r: int) -> int:
		if n < r:return 0
		if n < 0 or r < 0:return 0
		return self.fa[n] * self.fi[n - r] % self.mod
		
	def combr(self, n: int, r: int) -> int:
		if n == r == 0:return 1
		return self.comb(n + r - 1, r)
import typing
# 拡張Euclidの互除法
def extgcd(a: int, b: int, d: int = 0) -> typing.Tuple[int, int, int]:
	g = a
	if b == 0:
		x, y = 1, 0
	else:
		x, y, g = extgcd(b, a % b)
		x, y = y, x - a // b * y
	return x, y, g
 
# mod p における逆元
def invmod(a: int, p: int) -> int:
	x, y, g = extgcd(a, p)
	x %= p
	return x
n, k = map(int, input().split())
mod = 10 ** 9 + 7
C = Combinatorics(n, mod)
cnt = 1
for i in range(1, n + 1):
	cnt *= (k - i + 1)
	cnt *= invmod(i, mod)
	cnt %= mod
if k < n: cnt = 0
ans = 0
ans += cnt * C.comb(n - 1, 0)
for i in range(1, n + 1):
	cnt *= k + i
	cnt *= invmod(k + i - n, mod)
	cnt %= mod
	if k + i == n: cnt = 1
	ans += cnt * C.comb(n - 1, i)
	ans %= mod
print(ans % (mod))
0