洛谷 P2424:约数和 ← 整数分块算法 + 约数
【题目来源】
https://www.luogu.com.cn/problem/P2424
【题目描述】
对于一个数 X,函数 f(X) 表示 X 所有约数的和。例如:f(6)=1+2+3+6=12。对于一个 X,Smart 可以很快的算出 f(X)。现在的问题是,给定两个正整数 X,Y(X<Y),Smart 希望尽快地算出 f(X)+f(X+1)+……+f(Y)的值,你能帮助 Smart 算出这个值吗?
【输入格式】
输入文件仅一行,两个正整数 X 和 Y(X<Y),表示需要计算 f(X)+f(X+1)+⋯+f(Y)。
【输出格式】
输出只有一行,为 f(X)+f(X+1)+⋯+f(Y) 的值。
【输入样例】
123 321
【输出样例】
72543
【数据范围】
对于 20% 的数据有 1≤X<Y≤10^5。
对于 60% 的数据有 1≤X<Y≤1×10^7。
对于 100% 的数据有 1≤X<Y≤2×10^9。
【算法分析】
● 洛谷 P2424 要求计算:∑f(i),i=1~n。其中,f(i) 表示 i 的所有约数之和。直接计算每个数的约数之和再累加,复杂度太高。我们用交换求和顺序的技巧:
(1)枚举每个可能的约数 d,统计它在 1∼n 中作为约数出现的次数。
(2)对于约数 d,它在 1∼n 中作为约数出现的次数是 ⌊n/d⌋,每次贡献 d。因此:∑f(i)=d⋅⌊n/d⌋,d=1~n。
例如:若 i=1~6,则 ∑f(i)=f(1)+f(2)+f(3)+f(4)+f(5)+f(6)=1+(1+2)+(1+3)+(1+2+4)+(1+5)+(1+2+3+6)
=1×⌊6/1⌋+2×⌊6/2⌋+3×⌊6/3⌋+4×⌊6/4⌋+5×⌊6/5⌋+6×⌊6/6⌋。
● 对于块 [le,ri],⌊n/d⌋=k 为常数,需要计算:∑d⋅k=k⋅∑d,d=le~ri。区间 [le,ri] 内所有 d 的和是一个等差数列:∑d=(le+ri)⋅(ri−le+1)/2,d=le~ri。
● 注意:这道题交换了求和顺序,从“枚举每个数 i,求它的所有约数之和”变成了“枚举每个约数 d,统计它在多少个数中出现过”。这个转换改变了枚举的对象(从 i 变成了 d),但 d 本身的顺序依然是 1, 2, 3, ... 递增的,没有被打乱。
● 本题代码与“洛谷 P3935:Calculating:https://blog.csdn.net/hnjzsyjyj/article/details/162990202”及其类似。
【算法代码】
#include <bits/stdc++.h> using namespace std; typedef long long LL; LL cal(LL n) { LL t=0; for(LL le=1,ri=0; le<=n; le=ri+1) { LL k=n/le; ri=n/k; t=t+k*(ri+le)*(ri-le+1)/2; } return t; } int main() { ios::sync_with_stdio(0); cin.tie(0); LL le,ri; cin>>le>>ri; LL ans=cal(ri)-cal(le-1); cout<<ans<<"\n"; return 0; } /* in:123 321 out:72543 */
【参考文献】
https://blog.csdn.net/hnjzsyjyj/article/details/162990202
https://blog.csdn.net/hnjzsyjyj/article/details/163011369
https://blog.csdn.net/hnjzsyjyj/article/details/162819219