結果

問題 No.1529 Constant Lcm
ユーザー 👑 KazunKazun
提出日時 2021-04-29 01:59:17
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 849 ms / 3,000 ms
コード長 1,318 bytes
コンパイル時間 309 ms
コンパイル使用メモリ 81,940 KB
実行使用メモリ 97,620 KB
最終ジャッジ日時 2024-04-18 17:00:45
合計ジャッジ時間 10,624 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 40 ms
54,144 KB
testcase_01 AC 74 ms
73,344 KB
testcase_02 AC 41 ms
53,888 KB
testcase_03 AC 46 ms
54,016 KB
testcase_04 AC 46 ms
53,760 KB
testcase_05 AC 40 ms
54,272 KB
testcase_06 AC 42 ms
53,632 KB
testcase_07 AC 45 ms
53,376 KB
testcase_08 AC 42 ms
53,632 KB
testcase_09 AC 42 ms
53,504 KB
testcase_10 AC 794 ms
96,980 KB
testcase_11 AC 327 ms
83,620 KB
testcase_12 AC 95 ms
77,228 KB
testcase_13 AC 616 ms
93,164 KB
testcase_14 AC 414 ms
87,056 KB
testcase_15 AC 486 ms
88,012 KB
testcase_16 AC 367 ms
84,952 KB
testcase_17 AC 245 ms
81,144 KB
testcase_18 AC 267 ms
81,244 KB
testcase_19 AC 418 ms
86,752 KB
testcase_20 AC 849 ms
96,988 KB
testcase_21 AC 846 ms
96,992 KB
testcase_22 AC 842 ms
97,620 KB
testcase_23 AC 813 ms
97,488 KB
testcase_24 AC 800 ms
97,144 KB
testcase_25 AC 804 ms
97,396 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def Smallest_Prime_Factor(N):
    """0,1,2,...,Nの最小の素因数のリスト(0,1については1にしている)
    """

    if N==0:
        return [1]

    N=abs(N)
    L=list(range(N+1))
    L[0]=L[1]=1

    x=4
    while x<=N:
        L[x]=2
        x+=2

    x=9
    while x<=N:
        if L[x]==x:
            L[x]=3
        x+=6

    x=5
    Flag=0
    while x*x<=N:
        if L[x]==x:
            y=x*x
            while y<=N:
                if L[y]==y:
                    L[y]=x
                y+=x<<1
        x+=2+2*Flag
        Flag^=1

    return L

def Faster_Prime_Factorization(N,L):
    """

    L:Smallest_Prime_Factors(N)で求めたリスト
    """
    N=abs(N)

    D=[]
    while N>1:
        a=L[N]
        k=0
        while L[N]==a:
            k+=1
            N//=a
        D.append([a,k])
    return D
#================================================
from collections import defaultdict
N=int(input())
Mod=998244353
D=defaultdict(int)

P=Smallest_Prime_Factor(N)

D=defaultdict(int)

for i in range(1,(N+1)//2+1):
    A=defaultdict(int)
    for p,e in Faster_Prime_Factorization(i,P):
        A[p]+=e

    for p,e in Faster_Prime_Factorization(N-i,P):
        A[p]+=e

    for p in A:
        D[p]=max(D[p],A[p])

X=1
for p,e in D.items():
    X*=pow(p,e,Mod)
    X%=Mod

print(X)
0