結果

問題 No.1396 Giri
ユーザー HAPPA
提出日時 2021-02-14 23:26:55
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 636 ms / 2,000 ms
コード長 1,571 bytes
コンパイル時間 179 ms
コンパイル使用メモリ 82,464 KB
実行使用メモリ 98,580 KB
最終ジャッジ日時 2024-07-22 11:27:09
合計ジャッジ時間 7,176 ms
ジャッジサーバーID
(参考情報)
judge5 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys
from collections import defaultdict

sys.setrecursionlimit(10 ** 7)
input = sys.stdin.readline
f_inf = float('inf')
mod = 998244353


class Eratosthenes:
    def __init__(self, n):
        self.n = n
        self.min_factor = [-1] * (n + 1)
        self.min_factor[0], self.min_factor[1] = 0, 1
        self.primes = []
        self.get_primes()

    def get_primes(self):
        is_prime = [True] * (self.n + 1)
        is_prime[0] = is_prime[1] = False
        for i in range(2, self.n + 1):
            if not is_prime[i]:
                continue
            self.primes.append(i)
            self.min_factor[i] = i
            for j in range(i * 2, self.n + 1, i):
                is_prime[j] = False
                if self.min_factor[j] == -1:
                    self.min_factor[j] = i

    def return_primes(self):
        return self.primes

    def prime_factorization(self, n):
        res = []
        while n != 1:
            prime = self.min_factor[n]
            exp = 0
            while self.min_factor[n] == prime:
                exp += 1
                n //= prime
            res.append((prime, exp))
        return res


def resolve():
    n = int(input())

    er = Eratosthenes(n)
    tr = er.primes[-1]
    pf = defaultdict(int)
    for i in range(1, n + 1):
        if i == tr:
            continue
        for p, ex in er.prime_factorization(i):
            pf[p] = max(pf[p], ex)

    res = 1
    for k, v in pf.items():
        res *= pow(k, v, mod)
        res %= mod
    print(res)


if __name__ == '__main__':
    resolve()
0