【网络赛一：整除分块】C:Fear Factoring

难度

$3/10$
不学整除分块自己推很难的一题

题意

定义 $F(n)=\underset{i=1}{\sum ^n}\left[i\right|x]\times i$ 表示为 $n$ 的所有因子的和。

求 $S=\sum _{a\le n\le b}F(n)$

数据范围

$1\le a\le b\le 10^{12}$
$b-a\le 10^6$

思路

【1】首先考虑暴力枚举

求 $S$ 即求 $\sum _{1\le n\le b}F(n)-\sum _{1\le n\le a-1}F(n)$

记 $G(n)=\sum _{1\le k\le n}F(k)$

你每个值都是 $O(N)$ 枚举判断累加答案，时间复杂度 $O(b)=10^{12}$ ，不可行。

【2】考虑整除分块
普通的整除分块是 $g(n)=\underset{i=1}{\sum ^n}\lfloor \frac{n}{i}\rfloor$
但是易得我们所需求的分块为 $G(n)=\underset{i=1}{\sum ^n}\lfloor \frac{n}{i}\rfloor \times i$
于是我就不会了。。

我们枚举$G(12)$ 易得

i	1	2	3	4	5	6	7	8	9	10	11	12
$\frac{x}{i}$	12	6	4	3	2	2	1	1	1	1	1	1

很多连续的 $i$ 的情况下的系数 $\frac{x}{i}$ 也是相同的。
我们需要枚举每一个不同系数的左端点 $l$ 和右端点 $r$
首先，初始化 $l=1$
对于每一个 $l$ ，我们的新的区间的系数为 $\lfloor \frac{x}{l}\rfloor$
我们需要找到最大的 $r$ 满足 $\lfloor \frac{x}{l}\rfloor =\lfloor \frac{x}{r}\rfloor$
那么易得此时的 $r=\lfloor \frac{x}{\lfloor \frac{x}{l}\rfloor }\rfloor$
那么我们需要累加多少个$\frac{x}{i}$ 的答案呢？
显然为 $\underset{i=l}{\sum ^r}\lfloor \frac{x}{l}\rfloor \times i=\lfloor \frac{x}{l}\rfloor \underset{i=l}{\sum ^r}i=\lfloor \frac{x}{l}\rfloor \times \frac{(l+r)(r-l+1)}{2}$

算完这区间，明显易得 $l$ 更新为此时的 $r+1$

核心代码

时间复杂度 ：$O(\log b)$
Time(s)：125MS

/*
 _            __   __          _          _
| |           \ \ / /         | |        (_)
| |__  _   _   \ V /__ _ _ __ | |     ___ _
| '_ \| | | |   \ // _` | '_ \| |    / _ \ |
| |_) | |_| |   | | (_| | | | | |___|  __/ |
|_.__/ \__, |   \_/\__,_|_| |_\_____/\___|_|
        __/ |
       |___/
*/
#define ll unsigned long long			/// 不开 ull 会 WA 一个点
ll cal(ll x){
    ll l=1,r=0,ans=0;
    while(l<=x){
        r=x/(x/l);
        ans+=(x/l)*(r-l+1)*(l+r)/2;
        l=r+1;
    }
    return ans;
}
int main()
{
    ll a,b;
    cin >> a >> b;
    cout << cal(b) - cal(a-1) << endl;
    return 0;
}