結果
問題 |
No.2527 H and W
|
ユーザー |
![]() |
提出日時 | 2023-11-23 21:56:29 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 278 ms / 2,000 ms |
コード長 | 923 bytes |
コンパイル時間 | 171 ms |
コンパイル使用メモリ | 81,920 KB |
実行使用メモリ | 247,936 KB |
最終ジャッジ日時 | 2024-09-26 08:20:56 |
合計ジャッジ時間 | 6,212 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 23 |
ソースコード
def make_divisors(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] p = 998244353 fact = [1, 1] # fact[n] = (n! mod p) factinv = [1, 1] # factinv[n] = ((n!)^(-1) mod p) inv = [0, 1] # factinv 計算用 H,W,K = map(int,input().split()) mod = 998244353 D = make_divisors(K) def cmb(n, r, p): if (r < 0) or (n < r): return 0 r = min(r, n - r) return fact[n] * factinv[r] * factinv[n-r] % p for i in range(2, max(W,H) + 1): fact.append((fact[-1] * i) % mod) inv.append((-inv[p % i] * (p // i)) % mod) factinv.append((factinv[-1] * inv[-1]) % mod) ans = 0 for d in D: if d<=H and K//d<=W: ans = (ans + cmb(H,d,mod)*cmb(W,K//d,mod) ) % mod print(ans)