結果
| 問題 |
No.2527 H and W
|
| コンテスト | |
| ユーザー |
shobonvip
|
| 提出日時 | 2023-11-03 21:34:57 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 128 ms / 2,000 ms |
| コード長 | 2,428 bytes |
| コンパイル時間 | 352 ms |
| コンパイル使用メモリ | 82,476 KB |
| 実行使用メモリ | 108,612 KB |
| 最終ジャッジ日時 | 2024-09-25 19:30:21 |
| 合計ジャッジ時間 | 3,920 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 23 |
ソースコード
import os
import sys
from io import BytesIO, IOBase
BUFSIZE = 8192
class FastIO(IOBase):
newlines = 0
def __init__(self, file):
self._fd = file.fileno()
self.buffer = BytesIO()
self.writable = "x" in file.mode or "r" not in file.mode
self.write = self.buffer.write if self.writable else None
def read(self):
while True:
b = os.read(self._fd, max(os.fstat(self._fd).st_size, BUFSIZE))
if not b:
break
ptr = self.buffer.tell()
self.buffer.seek(0, 2), self.buffer.write(b), self.buffer.seek(ptr)
self.newlines = 0
return self.buffer.read()
def readline(self):
while self.newlines == 0:
b = os.read(self._fd, max(os.fstat(self._fd).st_size, BUFSIZE))
self.newlines = b.count(b"\n")+(not b)
ptr = self.buffer.tell()
self.buffer.seek(0, 2), self.buffer.write(b), self.buffer.seek(ptr)
self.newlines -= 1
return self.buffer.readline()
def flush(self):
if self.writable:
os.write(self._fd, self.buffer.getvalue())
self.buffer.truncate(0), self.buffer.seek(0)
class IOWrapper(IOBase):
def __init__(self, file):
self.buffer = FastIO(file)
self.flush = self.buffer.flush
self.writable = self.buffer.writable
self.write = lambda s:self.buffer.write(s.encode("ascii"))
self.read = lambda:self.buffer.read().decode("ascii")
self.readline = lambda:self.buffer.readline().decode("ascii")
sys.stdin, sys.stdout = IOWrapper(sys.stdin), IOWrapper(sys.stdout)
input = lambda:sys.stdin.readline().rstrip("\r\n")
mod = 998244353
N = 3 * 10**6 + 5
fact = [1]*(N+1)
factinv = [1]*(N+1)
for i in range(2, N+1):
fact[i] = fact[i-1] * i % mod
factinv[-1] = pow(fact[-1], mod-2, mod)
for i in range(N-1, 1, -1):
factinv[i] = factinv[i+1] * (i+1) % mod
def cmb(a, b):
if (a < b) or (b < 0):
return 0
return fact[a] * factinv[b] % mod * factinv[a-b] % mod
def makediv(n):
lower_divisors , upper_divisors = [], []
i = 1
while i*i <= n:
if n % i == 0:
lower_divisors.append(i)
if i != n // i:
upper_divisors.append(n//i)
i += 1
return lower_divisors + upper_divisors[::-1]
# 0 <= x <= h
# 0 <= y <= w
# x*w + y*h - x*y = k
# x*y - x*w - y*h + k = 0
# (x-h)(y-w) = hw - k
# (h-x)(w-y) = hw-k
h,w,k = map(int,input().split())
if k == 0:
print(1)
else:
ans = 0
for d in makediv(k):
x = h - d
y = w - k//d
if x < 0 or y < 0 or x > h or y > w:
continue
ans += cmb(h, x) * cmb(w, y) % mod
# x-h < 0 and y-w < 0
print(ans % mod)
shobonvip