結果

問題 No.1529 Constant Lcm
ユーザー tktk_snsn
提出日時 2021-06-04 20:17:06
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 862 ms / 3,000 ms
コード長 1,254 bytes
コンパイル時間 262 ms
コンパイル使用メモリ 82,460 KB
実行使用メモリ 97,792 KB
最終ジャッジ日時 2024-11-19 08:28:16
合計ジャッジ時間 10,439 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 24
権限があれば一括ダウンロードができます

ソースコード

diff # 

from collections import defaultdict, Counter
from itertools import chain


def prime_set(N):
    """
    Nまでの素数のsetを返す
    """
    if N < 4:
        return ({}, {}, {2}, {2, 3})[N]
    Nsq = int(N ** 0.5 + 0.5) + 1
    primes = {2, 3} | set(chain(range(5, N + 1, 6), range(7, N + 1, 6)))
    for i in range(5, Nsq, 2):
        if i in primes:
            primes -= set(range(i * i, N + 1, i * 2))
    return primes


def prime_sieve(N):
    #sieve[i] : iの最小の素因数
    sieve = [0] * (N + 1)
    prime = []
    for i in range(2, N + 1):
        if sieve[i] == 0:
            sieve[i] = i
            prime.append(i)
        for p in prime:
            if p > sieve[i] or i * p > N:
                break
            sieve[i * p] = p
    return sieve


N = int(input())
mod = 998244353
sieve = prime_sieve(N)
cnt = defaultdict(int)

for x in range(1, N):
    y = N - x
    if x > y:
        break

    P = []
    while sieve[x] != 0:
        s = sieve[x]
        P.append(s)
        x //= s
    while sieve[y] != 0:
        s = sieve[y]
        P.append(s)
        y //= s

    for k, v in Counter(P).items():
        cnt[k] = max(cnt[k], v)

ans = 1
for k, v in cnt.items():
    ans *= pow(k, v, mod)
    ans %= mod
print(ans)
0