結果
問題 | No.2902 ZERO!! |
ユーザー |
|
提出日時 | 2024-09-27 20:54:32 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 206 ms / 2,000 ms |
コード長 | 1,792 bytes |
コンパイル時間 | 153 ms |
コンパイル使用メモリ | 82,364 KB |
実行使用メモリ | 92,760 KB |
最終ジャッジ日時 | 2024-09-27 20:54:38 |
合計ジャッジ時間 | 5,638 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 41 |
ソースコード
import sysinput = lambda :sys.stdin.readline()[:-1]ni = lambda :int(input())na = lambda :list(map(int,input().split()))yes = lambda :print("yes");Yes = lambda :print("Yes");YES = lambda : print("YES")no = lambda :print("no");No = lambda :print("No");NO = lambda : print("NO")#######################################################################def factorial(n):if n == 0:return 1return n * factorial(n-1)def g(n, p):res = 0while n:res += n // pn //= preturn resdef factorize(n):res = []for i in range(2, int(n ** 0.5) + 1):if n % i == 0:cnt = 0while n % i == 0:n //= icnt += 1res.append((i, cnt))if n > 1:res.append((n, 1))return resdef gg(n, a):res = 10 ** 18for p, k in factorize(a):res = min(res, g(n, p) // k)return resdef greedy(n):res = 0for b in range(2, factorial(n) + 1):res += gg(n, b)print(b, gg(n, b))return resmod = 998244353def solve(n):prime = [True] * (n + 1)prime[0] = prime[1] = FalseP = []F = []for i in range(2, n + 1):if prime[i]:P.append(i)x = nres = 0while x:res += x // ix //= iF.append(res)for j in range(i * 2, n + 1, i):prime[j] = Falseans = 0for m in range(1, n + 1):res = 1for i, p in enumerate(P):e = F[i] // mif e == 0:breakres *= (1 + e)res %= mod#print(m, res)ans += res - 1ans %= modreturn ansN = ni()print(solve(N))