题面
这道题用到了数论分块
首先什么是数论分块 顾名思义就是分块在数论上的应用
不过非常神奇的一点是你如果不会分块但你还是可以会数论分块的
所求为
∑
i
=
1
n
k
m
o
d
i
\sum_{i=1}^{n}k mod i
i=1∑nkmodi
将这个式子改为:
∑
i
=
1
n
k
−
i
×
(
⌊
k
i
⌋
)
\sum_{i=1}^{n}k-i×(\left \lfloor \frac{k} {i} \right \rfloor)
i=1∑nk−i×(⌊ik⌋)
前面的
n
×
k
n×k
n×k直接先算出来
后面的
i
i
i是单增的
而
⌊
k
i
⌋
\left \lfloor \frac{k} {i} \right \rfloor
⌊ik⌋在一段范围内是相等的 总共只有
2
n
2\sqrt{n}
2n种取值
因此,每次求出这一段范围,直接计算等差数列即可
#include<bits/stdc++.h>
using namespace std;
inline int read()
{
int x=0,t=1;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=-1,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return x*t;
}
long long ans,n,k;
int main()
{
n=read();k=read();
ans=n*k;
for(long long l=1,r=0;l<=n;l=r+1)
{
if(k/l)r=min(n,k/(k/l));
else r=n;
ans-=(k/l)*(r-l+1)*(l+r)>>1;
}
printf("%lld\n",ans);
return 0;
}
数论分块推荐一道题:欧拉心算
(
b
z
o
j
4804
)
(bzoj4804)
(bzoj4804)
具体在这篇博客中哦!
数论
2333
2333
2333