結果

問題 No.1529 Constant Lcm
ユーザー shakayamishakayami
提出日時 2021-06-04 20:19:49
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 120 ms / 3,000 ms
コード長 1,244 bytes
コンパイル時間 182 ms
コンパイル使用メモリ 82,560 KB
実行使用メモリ 83,712 KB
最終ジャッジ日時 2024-04-29 23:35:34
合計ジャッジ時間 3,075 ms
ジャッジサーバーID
(参考情報)
judge2 / judge5
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 53 ms
55,296 KB
testcase_01 AC 55 ms
61,824 KB
testcase_02 AC 49 ms
55,296 KB
testcase_03 AC 50 ms
55,168 KB
testcase_04 AC 50 ms
55,296 KB
testcase_05 AC 49 ms
55,424 KB
testcase_06 AC 50 ms
55,296 KB
testcase_07 AC 49 ms
55,296 KB
testcase_08 AC 49 ms
55,296 KB
testcase_09 AC 49 ms
55,296 KB
testcase_10 AC 112 ms
81,792 KB
testcase_11 AC 81 ms
72,192 KB
testcase_12 AC 64 ms
65,536 KB
testcase_13 AC 103 ms
79,488 KB
testcase_14 AC 89 ms
74,752 KB
testcase_15 AC 91 ms
75,008 KB
testcase_16 AC 84 ms
73,088 KB
testcase_17 AC 76 ms
69,760 KB
testcase_18 AC 76 ms
69,760 KB
testcase_19 AC 91 ms
75,008 KB
testcase_20 AC 119 ms
83,328 KB
testcase_21 AC 118 ms
83,456 KB
testcase_22 AC 117 ms
83,712 KB
testcase_23 AC 118 ms
83,712 KB
testcase_24 AC 116 ms
83,584 KB
testcase_25 AC 120 ms
83,456 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