結果

問題 No.1529 Constant Lcm
ユーザー kozykozy
提出日時 2021-06-04 20:20:21
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,412 ms / 3,000 ms
コード長 995 bytes
コンパイル時間 365 ms
コンパイル使用メモリ 82,660 KB
実行使用メモリ 196,360 KB
最終ジャッジ日時 2024-11-19 08:43:55
合計ジャッジ時間 15,619 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 39 ms
52,224 KB
testcase_01 AC 75 ms
70,912 KB
testcase_02 AC 37 ms
51,712 KB
testcase_03 AC 39 ms
52,096 KB
testcase_04 AC 38 ms
52,224 KB
testcase_05 AC 38 ms
52,352 KB
testcase_06 AC 38 ms
52,352 KB
testcase_07 AC 38 ms
52,224 KB
testcase_08 AC 38 ms
52,224 KB
testcase_09 AC 39 ms
52,352 KB
testcase_10 AC 1,240 ms
186,060 KB
testcase_11 AC 458 ms
120,656 KB
testcase_12 AC 108 ms
79,488 KB
testcase_13 AC 1,065 ms
171,328 KB
testcase_14 AC 687 ms
138,620 KB
testcase_15 AC 712 ms
142,108 KB
testcase_16 AC 541 ms
127,124 KB
testcase_17 AC 285 ms
105,064 KB
testcase_18 AC 299 ms
106,784 KB
testcase_19 AC 701 ms
140,096 KB
testcase_20 AC 1,365 ms
196,160 KB
testcase_21 AC 1,343 ms
196,236 KB
testcase_22 AC 1,369 ms
196,172 KB
testcase_23 AC 1,331 ms
196,044 KB
testcase_24 AC 1,412 ms
196,292 KB
testcase_25 AC 1,334 ms
196,360 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