数论分块

面对这样一个问题
∑ i = 1 n ⌊ n i ⌋ \sum_{i = 1}^{n}{\lfloor \frac{n}{{i}} \rfloor} i=1nin
可以选择用O(n)的复杂度去暴力求解,但我们可以发现n/i的值在某一段中是相等的,比方说5/ 3, 5/ 4 ,5/5的大小都为1,所以可以用这一段的长度Len 乘以这一段的值就是可以快速求和。

代码:

int ans = 0;
for(int l = 1, r = 0; l <= n; l = r + 1) {
    r = n / (n / l);  			// 求区间的右端,这是一个数学规律 
    ans += (r - l + 1) * (n / l);
} 	

代码很简短,但是数论分块往往是隐藏在一些数学公式中的。

就比方说下面这个式子:
∑ i = 1 n k % i \sum_{i = 1}^{n}{k \% i} i=1nk%i
看上去和数论分块的那个形式不大一样,但是知道
k % i = k − ( k / i ) i k\%i = k - (k / i)i k%i=k(k/i)i
就可以把原来的式子化成
n k − ∑ i = 1 n ( k / i ) i nk - \sum_{i = 1}^{n} {(k/i)i} nki=1n(k/i)i
后面就用等差数列求和就行了。

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

int main()
{
    // freopen("i.txt", "r", stdin);
    ll n ,k;
    cin >> n >> k;
    ll ans = n * k;
    for(ll l = 1,r = 0; l <= n; l=r+1) {
        if(k / l != 0) r =s min(k / (k / l), n);
        else r = n;    
        ans -= (k / l) * (r - l + 1) * (l + r) / 2;
    }
    cout << ans << endl;
    return 0;
}
  • 9
    点赞
  • 21
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值