結果
| 問題 |
No.1529 Constant Lcm
|
| コンテスト | |
| ユーザー |
AEn
|
| 提出日時 | 2022-09-02 01:26:00 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 593 ms / 3,000 ms |
| コード長 | 938 bytes |
| コンパイル時間 | 351 ms |
| コンパイル使用メモリ | 82,164 KB |
| 実行使用メモリ | 92,384 KB |
| 最終ジャッジ日時 | 2024-11-14 17:53:45 |
| 合計ジャッジ時間 | 8,470 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 24 |
ソースコード
from collections import defaultdict
N = int(input())
mod = 998244353
def smallest_prime_factors(n):
"""各数の最小素因数を求める"""
spf = [0]*(n+1)
for i in range(n+1):
spf[i] = i
for i in range(2,n+1):
if i*i>n:
break
if spf[i] == i:
for j in range(i*i, n+1, i):
if spf[j] == j:
spf[j] = i
return spf
def fac(x, spf, ret) -> list:
while x!=1:
if spf[x] in ret:
ret[spf[x]] += 1
else:
ret[spf[x]] = 1
x //= spf[x]
return ret
res = [0]*N
spf = smallest_prime_factors(N-1)
for i in range(N//2):
ret = {}
a, b = i+1, N-i-1
ret = fac(a, spf, ret)
ret = fac(b, spf, ret)
for key in ret.keys():
res[key] = max(ret[key], res[key])
ans = 1
for i in range(N):
if res[i] != 0:
ans *= pow(i, res[i], mod)
ans %= mod
print(ans)
AEn