結果
問題 | No.2902 ZERO!! |
ユーザー |
![]() |
提出日時 | 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 |
ソースコード
mod = 998244353class Prime():def __init__(self,N):self.N = Nself.Is_prime = [1 for i in range(N+1)]self.Is_prime[0] = self.Is_prime[1] = 0self.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:continueself.P.append(p)#self.Fact[p].append((p,1)) ##########################for d in range(2,N//p + 1):self.Is_prime[p*d] = 0q = 1r = pwhile (p*d) % (p*r) == 0:q += 1r *= p#self.Fact[p*d].append((p,q)) ######################def is_prime(self,p):return self.Is_prime[p]def primelist(self):return self.Pdef fact(self,n): #########のところを外して使う.return self.Fact[n]N = int(input())P = Prime(N)Q = []for p in range(2,N + 1):q = 0pp = pif P.is_prime(p):while pp <= N:q += N // pppp *= pQ.append(q)M = max(Q)C = [0 for i in range(M + 1)]for q in Q:C[q] += 1X = [1 for i in range(M + 1)]res = 0for 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 % modans = -(M + 1)ans += sum(X) % modans %= modprint(ans)