結果

問題 No.1529 Constant Lcm
ユーザー tktk_snsntktk_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
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 43 ms
53,888 KB
testcase_01 AC 96 ms
76,928 KB
testcase_02 AC 43 ms
53,760 KB
testcase_03 AC 42 ms
53,760 KB
testcase_04 AC 43 ms
53,504 KB
testcase_05 AC 42 ms
54,144 KB
testcase_06 AC 41 ms
53,888 KB
testcase_07 AC 48 ms
54,144 KB
testcase_08 AC 42 ms
53,760 KB
testcase_09 AC 43 ms
53,504 KB
testcase_10 AC 755 ms
96,896 KB
testcase_11 AC 330 ms
83,712 KB
testcase_12 AC 96 ms
76,800 KB
testcase_13 AC 643 ms
93,312 KB
testcase_14 AC 427 ms
86,528 KB
testcase_15 AC 457 ms
88,192 KB
testcase_16 AC 349 ms
84,480 KB
testcase_17 AC 220 ms
81,152 KB
testcase_18 AC 245 ms
81,152 KB
testcase_19 AC 421 ms
87,144 KB
testcase_20 AC 806 ms
97,792 KB
testcase_21 AC 862 ms
97,664 KB
testcase_22 AC 807 ms
97,408 KB
testcase_23 AC 796 ms
97,792 KB
testcase_24 AC 797 ms
97,792 KB
testcase_25 AC 787 ms
97,536 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