結果
問題 | No.1396 Giri |
ユーザー |
|
提出日時 | 2021-02-14 21:59:13 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 245 ms / 2,000 ms |
コード長 | 644 bytes |
コンパイル時間 | 275 ms |
コンパイル使用メモリ | 12,800 KB |
実行使用メモリ | 22,016 KB |
最終ジャッジ日時 | 2024-07-22 09:37:54 |
合計ジャッジ時間 | 3,538 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 23 |
ソースコード
from math import gcd MOD = 998244353 lcm = lambda a, b: a * b // gcd(a, b) modinv = lambda a, mod=MOD: pow(a, mod - 2, mod) def primes(n): is_prime = [True] * (n + 1) is_prime[0] = False is_prime[1] = False for i in range(2, int(n ** 0.5) + 1): if not is_prime[i]: continue for j in range(i * 2, n + 1, i): is_prime[j] = False return (i for i in range(n + 1) if is_prime[i]) N = int(input()) P = tuple(primes(N)) L = 1 for p in P: tmp = N ct = 0 while tmp >= p: tmp //= p ct += 1 L *= pow(p, ct, MOD) L %= MOD print(L * modinv(max(P)) % MOD)