結果

問題 No.1529 Constant Lcm
ユーザー kozykozy
提出日時 2021-06-04 20:20:21
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,483 ms / 3,000 ms
コード長 995 bytes
コンパイル時間 173 ms
コンパイル使用メモリ 82,176 KB
実行使用メモリ 196,424 KB
最終ジャッジ日時 2024-04-29 23:39:00
合計ジャッジ時間 17,007 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 40 ms
52,224 KB
testcase_01 AC 80 ms
71,168 KB
testcase_02 AC 41 ms
52,224 KB
testcase_03 AC 42 ms
52,224 KB
testcase_04 AC 41 ms
51,840 KB
testcase_05 AC 41 ms
52,224 KB
testcase_06 AC 41 ms
52,352 KB
testcase_07 AC 41 ms
52,736 KB
testcase_08 AC 41 ms
51,968 KB
testcase_09 AC 41 ms
52,224 KB
testcase_10 AC 1,343 ms
185,944 KB
testcase_11 AC 532 ms
120,636 KB
testcase_12 AC 119 ms
79,360 KB
testcase_13 AC 1,166 ms
171,448 KB
testcase_14 AC 742 ms
138,496 KB
testcase_15 AC 778 ms
142,108 KB
testcase_16 AC 596 ms
127,128 KB
testcase_17 AC 343 ms
105,048 KB
testcase_18 AC 377 ms
106,520 KB
testcase_19 AC 780 ms
140,224 KB
testcase_20 AC 1,430 ms
196,172 KB
testcase_21 AC 1,449 ms
196,164 KB
testcase_22 AC 1,463 ms
196,304 KB
testcase_23 AC 1,450 ms
196,172 KB
testcase_24 AC 1,483 ms
196,164 KB
testcase_25 AC 1,451 ms
196,424 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def primes(n):
    is_prime = [True] * (n + 1)
    is_prime[0] = False
    is_prime[1] = False
    for i in range(2, int(n**0.5) + 1):
        if not is_prime[i]:
            continue
        for j in range(i * 2, n + 1, i):
            is_prime[j] = False
    return [i for i in range(n + 1) if is_prime[i]]
N=int(input())

mod=998244353
L=primes(N-1)
R=[0]*N
#コーナーに気をつけて
K=[list() for i in range(N)]
for i in range(len(L)):
  s=L[i]
  for j in range(s,N,s):
    K[j].append(s)
for i in range(1,N):
  a=i
  b=N-i
  if a>b:
    break
  x=dict()
  a2=a
  b2=b
  for j in range(len(K[a])):
    s=K[a][j]
    this=0
    while a2%s==0:
      a2//=s
      this+=1
    if s in x:
      x[s]+=this
    else:
      x[s]=this
  for j in range(len(K[b])):
    s=K[b][j]
    this=0
    while b2%s==0:
      b2//=s
      this+=1
    if s in x:
      x[s]+=this
    else:
      x[s]=this
  for i in x:
    R[i]=max(R[i],x[i])
ans=1
for i in L:
  ans*=pow(i,R[i],mod)
  ans%=mod
print(ans)
0