結果
問題 | 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 = 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)