結果

問題 No.1529 Constant Lcm
ユーザー shakayamishakayami
提出日時 2021-06-04 20:19:49
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 100 ms / 3,000 ms
コード長 1,244 bytes
コンパイル時間 253 ms
コンパイル使用メモリ 82,184 KB
実行使用メモリ 83,972 KB
最終ジャッジ日時 2024-11-19 08:40:22
合計ジャッジ時間 3,045 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 38 ms
56,796 KB
testcase_01 AC 45 ms
61,884 KB
testcase_02 AC 40 ms
56,384 KB
testcase_03 AC 40 ms
55,940 KB
testcase_04 AC 42 ms
56,208 KB
testcase_05 AC 41 ms
55,832 KB
testcase_06 AC 41 ms
56,104 KB
testcase_07 AC 40 ms
57,100 KB
testcase_08 AC 41 ms
55,872 KB
testcase_09 AC 40 ms
55,908 KB
testcase_10 AC 94 ms
81,964 KB
testcase_11 AC 69 ms
72,428 KB
testcase_12 AC 53 ms
65,664 KB
testcase_13 AC 87 ms
79,884 KB
testcase_14 AC 81 ms
75,096 KB
testcase_15 AC 76 ms
75,420 KB
testcase_16 AC 70 ms
73,296 KB
testcase_17 AC 63 ms
69,668 KB
testcase_18 AC 61 ms
70,260 KB
testcase_19 AC 77 ms
75,316 KB
testcase_20 AC 98 ms
83,972 KB
testcase_21 AC 98 ms
83,668 KB
testcase_22 AC 97 ms
83,796 KB
testcase_23 AC 100 ms
83,568 KB
testcase_24 AC 97 ms
83,708 KB
testcase_25 AC 95 ms
83,876 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from collections import defaultdict
mod=998244353
import random
#N=int(input())
def solve(N):
    isPrime=[True for i in range(N+1)]
    isPrime[0]=False;isPrime[1]=False
    Prime=[]
    for i in range(2,N+1):
        if isPrime[i]:
            for j in range(2*i,N+1,i):
                isPrime[j]=False
            Prime.append(i)
    def primefact(M):
        res=defaultdict(int)
        p=2
        while(p*p<=M):
            while(M%p==0):
                res[p]+=1
                M//=p
            p+=1
        if M>1:res[M]+=1
        return res
    PF=primefact(N)
    ans=1
    for p in Prime:
        if N%p==0:
            continue
        k=0
        while(p**(k+1)<N):k+=1
        ans*=pow(p,k,mod)
        ans%=mod
    for p in PF:
        A=N//(p**PF[p])
        if A>1:
            cnt=PF[p]*2
            pr=p
            while(A>pr):
                cnt+=1
                pr*=p
            ans*=pow(p,cnt,mod)
            ans%=mod
        else:
            cnt=(PF[p]-1)*2
            ans*=pow(p,cnt,mod)
            ans%=mod
    return ans

from math import gcd
def lcm(x,y):return (x*y)//gcd(x,y)
def naive(N):
    L=1
    for i in range(1,N):
        L=lcm(L,i*(N-i))
    L%=mod
    return L
print(solve(int(input())))
0