2018.07.17 CQOI2017 余数求和(整除分块)

洛谷传送门
bzoj传送门
这道题要用到学习莫比乌斯反演时掌握的整除分块算法,也就是对于一个数 n n n n 除以 1 1 n n 的数商的取值只有 sqrt(n) s q r t ( n ) 种,然后这道题对于商相同的一段余数,它们会构成一个等差数列,于是直接上整除分块统计答案即可。

代码如下:

#include
using namespace std;
int main(){
    long long n,k;
    long long ans=0;
    scanf("%lld%lld",&n,&k); 
    for(long long l=1,r;l<=n;l=r+1){
        if(k/l==0)r=n;
        else r=min(k/(k/l),n);
        long long a1=k%r,an=a1+(r-l)*(k/l);
        ans+=(a1+an)*(r-l+1)>>1;
    }
    printf("%lld",ans);
    return 0;
}

你可能感兴趣的:(数学,#,小技巧,#,整除分块)