蓝桥杯-因数平方和

"""
题目来源:
https://www.lanqiao.cn/problems/2139/learning/?page=1&first_category_id=1&name=%E5%9B%A0%E6%95%B0%E5%B9%B3%E6%96%B9%E5%92%8C
"""
"""
在1~n中有n个数字是1的倍数, 因此因子1有n个
在1~n中有n//2个数字是2的倍数, 因此因子2有n//2个
在1~n中有n//3个数字是3的倍数, 因此因子3有n//3个
在1~n中有n//4个数字是4的倍数, 因此因子4有n//4个
在1~n中有n//5个数字是5的倍数, 因此因子5有n//5个
...

令n=100, 则:
因子1-1都是有100个, 因子2-2都是有50个 ... 因子12-12都是有8个, 因子13-14都是有7个, 
因子15-16都是有6个, 因子17-20都是有5个, 因子21-25都是有4个。
可以发现我们可以按照因子的个数进行区间划分:
在[1, 1]区间的因子的个数为100//1=100个
在[2, 2]区间的因子的个数为100//2=50个
...
在[12, 12]区间的因子的个数为100//12=8个
在[13, 14]区间的因子的个数为100//13=7个
在[15, 16]区间的因子的个数为100//15=6个
在[17, 20]区间的因子的个数为100//17=5个
在[21, 25]区间的因子的个数为100//21=4个
...
在[l1, r1]区间的因子的个数为100//l1个
在[r1 + 1, r2]区间的因子的个数为100//(r11 + 1)个
...
(数论关系: 区间右端点r1与左端点l1的关系: r1 = n // (n // l1))

首先求每个区间因子的平方和, 在乘以因子的个数,
然后将所有区间求得的结果求和即可得到答案(结果需要求模)
"""
MOD = 1000000007

# 求数字1-n的平方和
def S(x):
    return x * (x + 1) * (2 * x + 1) // 6

n = int(input())
ans = 0
l = 1
# 枚举左端点
while l <= n:
    r = min(n // (n // l), n)
    ans += (n // l) * (S(r) - S(l - 1))
    l = r + 1
print(ans % MOD)