結果
| 問題 |
No.1396 Giri
|
| コンテスト | |
| ユーザー |
👑 SPD_9X2
|
| 提出日時 | 2021-02-14 21:38:02 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 103 ms / 2,000 ms |
| コード長 | 1,220 bytes |
| コンパイル時間 | 308 ms |
| コンパイル使用メモリ | 82,116 KB |
| 実行使用メモリ | 85,376 KB |
| 最終ジャッジ日時 | 2024-07-22 09:02:51 |
| 合計ジャッジ時間 | 2,709 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 23 |
ソースコード
"""
すべての倍数
ある数字を除いた全体のlcmがわかればおk
算数
一番大きい素数…?
全部のlcmを取った時、
素数のN乗を考えたとき、1回しか出てこないやつの内最大
"""
import math
import sys
from sys import stdin
mod = 998244353
def inverse(a,mod): #aのmodを法にした逆元を返す
return pow(a,mod-2,mod)
def Sieve(n): #n以下の素数全列挙(O(nloglogn)) retは素数が入ってる。divlisはその数字の素因数が一つ入ってる
ret = []
divlis = [-1] * (n+1) #何で割ったかのリスト(初期値は-1)
flag = [True] * (n+1)
flag[0] = False
flag[1] = False
ind = 2
while ind <= n:
if flag[ind]:
ret.append(ind)
ind2 = ind ** 2
while ind2 <= n:
flag[ind2] = False
divlis[ind2] = ind
ind2 += ind
ind += 1
return ret,divlis
N = int(stdin.readline())
prime,tmp = Sieve(N)
nmax = 1
ans = 1
for p in prime:
pm = p
while pm * p <= N:
pm *= p
if N // pm == 1:
nmax = max(nmax,p)
ans *= pm
ans %= mod
ans *= inverse(nmax,mod)
print (ans % mod)
SPD_9X2