【BZOJ 1257】[CQOI2007]余数之和sum

分析发现所有a/b=x....y 如果x相同的话y是一个等差数列于是可以利用等差求和，枚举x最多只枚举到根号n，然后再暴力枚举根号n 以下的数字，复杂度根号n

今天研究莫比乌斯反演的时候突然觉得自己以前的做法很傻

求解

∑k%i

=∑k-floor(k/i)*i

=n*k-∑(k/i)*i

然后发现其实floor(k/i)在一定的范围内是定值变化可以用等差数列求和也就是(l+r)*(r-l+1)*floor(k/i)/2;

然后每一次只用求出l和r就好了显然k/(k/i)就是k/i向下取整的最后一个位置（即r）然后。。。就没有然后了

以前的

#include<cstdio>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
LL n,k,ans,l;


int main(){
	scanf("%lld%lld",&n,&k);
	if(n>=k)ans+=(n-k)*k;
	LL last=min(k,n);
	l=min(k,n);
	for(LL i=2;i*i<=k;i++){
		LL x=k/i+1;
		if(x>l)continue;
		ans+=(k%last+k%x)*(last-x+1)/2;
		last=x-1;
	}
	for(LL i=2;i<=last;i++){
		if(i>l)continue;
		ans+=k%i;
	}
	printf("%lld",ans);
	return 0;
}//10 12

now

#include<cstdio>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
int n,k;
LL ans;

int main(){
	scanf("%d%d",&n,&k);
	ans+=(LL)n*k;
	n=min(n,k);
	for(int l,r,j,i=1;i<=n;i=r+1){
		j=k/i;l=k/(j+1)+1,r=k/j;
		r=min(r,n);
		ans-=(LL)(l+r)*(r-l+1)*j/2;
	}
	cout<<ans;
	return 0;
}