結果

問題 No.1529 Constant Lcm
ユーザー tktk_snsntktk_snsn
提出日時 2021-06-04 20:17:06
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 896 ms / 3,000 ms
コード長 1,254 bytes
コンパイル時間 601 ms
コンパイル使用メモリ 82,304 KB
実行使用メモリ 97,920 KB
最終ジャッジ日時 2024-04-29 23:28:28
合計ジャッジ時間 11,297 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 46 ms
53,504 KB
testcase_01 AC 101 ms
76,928 KB
testcase_02 AC 46 ms
53,632 KB
testcase_03 AC 47 ms
53,888 KB
testcase_04 AC 47 ms
53,760 KB
testcase_05 AC 47 ms
53,376 KB
testcase_06 AC 46 ms
53,760 KB
testcase_07 AC 46 ms
53,888 KB
testcase_08 AC 47 ms
53,504 KB
testcase_09 AC 46 ms
53,888 KB
testcase_10 AC 804 ms
97,024 KB
testcase_11 AC 347 ms
83,968 KB
testcase_12 AC 110 ms
77,056 KB
testcase_13 AC 692 ms
93,312 KB
testcase_14 AC 447 ms
86,528 KB
testcase_15 AC 495 ms
88,192 KB
testcase_16 AC 381 ms
84,352 KB
testcase_17 AC 242 ms
80,896 KB
testcase_18 AC 268 ms
80,896 KB
testcase_19 AC 451 ms
86,912 KB
testcase_20 AC 880 ms
97,536 KB
testcase_21 AC 896 ms
97,408 KB
testcase_22 AC 830 ms
97,792 KB
testcase_23 AC 847 ms
97,792 KB
testcase_24 AC 821 ms
97,536 KB
testcase_25 AC 815 ms
97,920 KB
権限があれば一括ダウンロードができます

ソースコード

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