2021-08-03SDUT 2021 Summer Individual Contest - 4(for 20)(补题)

20 篇文章 0 订阅
6 篇文章 0 订阅

C - Fear Factoring

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
题意: 求F(x)定义为,x的约数和。求 ∑ i = a b f ( i ) \sum_{i=a}^{b}f(i) \qquad i=abf(i)
思路: 枚举1到n的所有约数,会发现,所有因子的 i i i出现次数就是 n i \frac{n}{i} in,那么因子 i i i对答案的贡献就是⌊ n i \frac{n}{i} in ∗ i *i i
举例来说, 1到10的所有约数和,会发现有10个1,5个2,3个3,2个4,2个5,1个6,1个7,1个8,1个9,1个10,但这样枚举一定会超时,数据太大,要优化时间复杂度。
可以发现4到5的个数一样,6到10的个数一样,总会有区间 [ L , R ] [L,R] [L,R]的数的个数是一样的,所以我们可以直接利用等差数列求和求这段 [ L , R ] [L,R] [L,R]的和,并乘上对应的个数,这样就不用一一枚举因子。
区间的 L L L可以直接就是上一个区间的 R ′ R' R + 1 +1 +1,就是当前区间 L L L
区间的 R R R n / ( n / L ) n/(n/L) n/(n/L) ⌊ n / L ⌋ ⌊n/L⌋ n/L求的是因子L出现的次数,再用n/⌊n/L⌋,求得就是出现这个次数的因子当中最大的那个,因为⌊n/L⌋是向下取整,理论上 n / R n/R n/R n / L n/L n/L,不是同一个数,但因为向下取整变成了同一个数,我们定义了R为出现当前次数最右端的因子(最大),所以如果R如果加上1,变成 R + 1 R+1 R+1, n / ( R + 1 ) n/(R+1) n/(R+1),次数一定不是当前值,所以用n除以当前次数,算出来的就是 R R R
剩下的话,次数有了, [ L , R ] [L,R] [LR],有了,直接等差数列求和即可
举例的话,就像 ⌊ 10 / 6 ⌋ ⌊10/6⌋ 10/6 ⌊ 10 / 10 ⌋ ⌊10/10⌋ 10/10是一样的,向下取整相当于把分母凑了个最大数。
注意:要开unsigned long long,否则会wa
这是除法分块的内容,算是记个小笔记

#include <bits/stdc++.h>
using namespace std;
#define ll unsigned long long
ll F(ll n) {
    ll res = 0;
    for (ll L = 1, R; L <= n; L = R + 1) {
        R = n / (n / L);
        res += (n / L) * (L + R) * (R - L + 1) / 2;
    }
    return res;
}
int main() {
    ll a, b;
    cin >> a >> b;
    ll res = F(b) - F(a - 1);
    cout << res << endl;
}

To be continued
如果你有任何建议或者批评和补充,请留言指出,不胜感激

评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值