P2424 约数和（整除分块）

P2424 约数和（整除分块）

题目入口
题目描述
对于一个数 XX，函数 f(X)f(X) 表示 XX 所有约数的和。例如：f(6)=1+2+3+6=12f(6)=1+2+3+6=12。对于一个 XX，Smart 可以很快的算出 f(X)f(X)。现在的问题是，给定两个正整数 X,Y(X<Y)X,Y(X<Y)，Smart 希望尽快地算出 f(X)+f(X+1)+……+f(Y)f(X)+f(X+1)+……+f(Y)的值，你能帮助 Smart 算出这个值吗？

输入格式
输入文件仅一行，两个正整数 XX 和 Y(X<Y)Y(X<Y)，表示需要计算 f(X)+f(X+1)+\dots +f(Y)f(X)+f(X+1)+⋯+f(Y)。

输出格式
输出只有一行，为 f(X)+f(X+1)+\dots+f(Y)f(X)+f(X+1)+⋯+f(Y) 的值。

输入#1
2 4

输出#1
14

输入#2
123 321

输出#2
72543

思路
我们知道，在1-n中，若a<n，则a的倍数有n/a个，因此一种暴力的做法便是枚举约数1-y在x-y的倍数，复杂度O（y），超时，于是我们得想新方法。我们注意到我们需要知道的只是以1-y为约数的数有多少个，而大部分数的个数相同，如1-100内以51，52，53……100为约数的个数都只有一个，于是我们想到整除分块。易得答案为∑( y / i - x-1 / i ) * i ( 1 <=i<=y )。

代码

#include<bits/stdc++.h>
#define int long long
using namespace std;
int x,y;
signed main()
{
    cin>>x>>y;
    int res=0;
    for(int l=1,r;l<=y;l=r+1)
    {
        if(y/l)
            r=y/(y/l);
        else break;
        res+=(y/l)*(r-l+1)*(l+r)/2;
    }
    for(int l=1,r;l<=y;l=r+1)
    {
        if((x-1)/l)
             r=(x-1)/((x-1)/l);
        else break;
        res-=(x-1)/l*(r-l+1)*(l+r)/2;
    }
    cout<<res<<endl;
    return 0;
}