Luogu P3935 Calculating

题目大意

给定 l , r l,r l,r,求 ∑ i = l r ∑ k ∣ i 1 \sum_{i=l}^r\sum_{k|i}1 i=lrki1

数据范围  1 ⩽ l , r ⩽ 1.6 × 1 0 14 3  second 1\leqslant l,r\leqslant 1.6\times 10^{14}\quad3\ \text{second} 1l,r1.6×10143 second

题解

可以发现, 原式等价于:
Sum d ( r ) − Sum d ( l − 1 ) \text{Sum}_{d}(r)-\text{Sum}_{d}(l-1) Sumd(r)Sumd(l1)其中 d d d 是约数个数函数。
某些学过杜教筛的(比如我)会开始考虑,然后发现打死都想不出来(至少我想不出来)……可能有的大佬想出来了,结果时间复杂度是 Θ ( n 2 3 ) \Theta(n^{\frac 23}) Θ(n32),发现死都卡不过去……好吧。
Sum d ( n ) = ∑ i = 1 n ∑ k ∣ i 1 = ∑ k = 1 n ∑ k ∣ i i ⩽ n 1 = ∑ k = 1 n ⌊ n k ⌋ \begin{aligned} \text{Sum}_{d}(n)&=\sum_{i=1}^n\sum_{k|i}1\\ &=\sum_{k=1}^n\sum_{k|i}^{i\leqslant n}1\\ &=\sum_{k=1}^n\Big\lfloor\frac nk \Big\rfloor\\ \end{aligned} Sumd(n)=i=1nki1=k=1nkiin1=k=1nkn

心情舒畅。除法分块,时间复杂度 Θ ( n ) \Theta(\sqrt n) Θ(n )

#include<bits/stdc++.h>
using namespace std;
#define mod 998244353
long long ds(long long n){
	long long ans=0;
	for(long long l=1,r;l<=n;l=r+1){
		r=n/(n/l); if(r>n) r=n;
		ans+=(n/l)*(r-l+1);
		ans=ans%mod;
	} return ans;
}
int main(void){
	long long l,r;
	scanf("%lld%lld",&l,&r);
	printf("%lld",(ds(r)-ds(l-1)+mod)%mod);
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值