結果

問題 No.2902 ZERO!!
ユーザー PNJ
提出日時 2024-09-27 20:27:41
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 344 ms / 2,000 ms
コード長 1,334 bytes
コンパイル時間 226 ms
コンパイル使用メモリ 82,428 KB
実行使用メモリ 108,960 KB
最終ジャッジ日時 2024-09-27 20:27:50
合計ジャッジ時間 7,024 ms
ジャッジサーバーID
(参考情報)
judge1 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff #

mod = 998244353

class Prime():
  def __init__(self,N):
    self.N = N
    self.Is_prime = [1 for i in range(N+1)]
    self.Is_prime[0] = self.Is_prime[1] = 0
    self.P = []
    #self.Fact = [[] for i in range(N+1)]   #################
    #self.Fact[1] = [(1,1)]  ################################
    for p in range(2,N + 1):
      if self.Is_prime[p] == 0:
        continue
      self.P.append(p)
      #self.Fact[p].append((p,1))  ##########################
      for d in range(2,N//p + 1):
        self.Is_prime[p*d] = 0
        q = 1
        r = p
        while (p*d) % (p*r) == 0:
          q += 1
          r *= p
        #self.Fact[p*d].append((p,q))  ######################
  
  def is_prime(self,p):
    return self.Is_prime[p]
  
  def primelist(self):
    return self.P
  
  def fact(self,n): #########のところを外して使う.
    return self.Fact[n]
  
N = int(input())
P = Prime(N)
Q = []
for p in range(2,N + 1):
  q = 0
  pp = p
  if P.is_prime(p):
    while pp <= N:
      q += N // pp
      pp *= p
    Q.append(q)

M = max(Q)
C = [0 for i in range(M + 1)]
for q in Q:
  C[q] += 1

X = [1 for i in range(M + 1)]
res = 0
for q in range(1,M + 1):
  if C[q]:
    for r in range(1,q + 1):
      c = pow((q // r) + 1,C[q],mod)
      X[r] = X[r] * c % mod
ans = -(M + 1)
ans += sum(X) % mod
ans %= mod
print(ans)
0