阶乘的约数和

import os
import sys
from collections import Counter

mod = 998244353

n = int(input())


# 线性筛
def is_(m):
    lis = []
    vis = [0] * (m + 1)
    vis[0] = vis[1] = 1
    for i in range(2, m + 1):
        if vis[i] == 0:
            lis.append(i)
        for x in lis:
            if i * x > m:
                break
            vis[i * x] = 1
            if i % x == 0:
                break
    return vis, lis


vis, lis = is_(n)


# 唯一分解定理
def f(x):
    li = []
    if vis[x] == 0:
        li.append(x)
        return li
    for i in range(2, x + 1):
        while x % i == 0:
            li.append(i)
            x //= i
        if x == 1:
            break
    return li


def ksm(a, b, k):
    return pow(a, b, k)


li = []
for i in range(1, n + 1):
    li += f(i)

dic = Counter(li)
ans = 1
for i in dic:
    ans = ans * (ksm(i, dic[i] + 1, mod) - 1) % mod * ksm(i - 1, mod - 2, mod) % mod
print(ans)

链接:阶乘的约数和

  • 3
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值