2017 J 安徽省赛 数论 + 分块

题意:

定以

f(n)=(i=1nϕ(ni))%(n+1)

g(n)=ni=1f(i)
n <= 1000000

分析:
设n的质因子分解为
n=ki=1paii
那么n、n^2…..n^n的表示:
ni=ki=1paiii
所以 ϕ(ni)=ϕ(n)ni1

f(n)=(i=1nϕ(ni))=(ϕ(n)i=0n1ni)%(n+1)

我们考虑 n1i=0ni%(n+1)

如果n为偶数:
S(n)=n1i=0ni=n/21i=0(n2i+n2i+1)=n/21i=0n2i(n+1)
所以S(n) % (n + 1) = 0, f(n) = 0;
若n为奇数
S(n)=1+n1i=1ni=1+n1/2i=1(n2i1+n2i)=1+n1/2i=1n2i1(n+1)
此时S(n) % (n + 1) = 1, f(n) = phi(n)

预处理欧拉数, 记个sum就可…

这题卡样例, 第三组给错了也是蛮骚的

### Python 实现数论分块 数论分块是一种优化技巧,在处理某些涉及整除问题的区间查询时非常有效。其核心思想在于利用整除性质减少重复计算次数。 对于给定范围 `[l, r]` 和常量 `d`,当遍历该范围内每一个元素并执行模运算或整除操作时,可能会遇到大量冗余计算。数论分块通过预先计算若干个连续区间的贡献来加速这一过程。 #### 数论分块原理说明 假设要解决形如 $\sum_{i=l}^{r}\lfloor\frac{n}{i}\rfloor$ 的求和问题: - 对于任意正整数 $n$ 及其因子 $k=\lfloor \frac{n}{m} \rfloor$ - 存在一个最大值 $j=min(n/(n//i),r)$ 使得对于所有的 $t∈[i,j]$ 都满足 $(n/t)=(n/i)$ - 这意味着我们可以一次性处理整个区间内具有相同商的结果,而不是逐项累加 #### 示例代码实现 下面是一个简单的 Python 函数用来展示如何应用数论分块方法来进行上述类型的求和: ```python def num_theory_block_sum(n): i = 1 total = 0 while i <= n: block_end = min(n // (n // i), n) # 计算当前块结束位置 value = n // i # 当前块内的公共值 count = block_end - i + 1 # 块大小 total += value * count # 更新总和 i = block_end + 1 # 移动到下一个未处理的位置 return total if __name__ == "__main__": test_value = int(input("Enter an integer to calculate sum of floor divisions up to this number: ")) result = num_theory_block_sum(test_value) print(f"The calculated sum is {result}.") ``` 此函数接收一个参数 `n` 并返回从 1 到 n 所有整数与其对应的最大不超过自身的倍数值之间的向下取整后的商之和[^1]。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值