洛谷P2261 [CQOI2007] 余数求和 [数论分块]

  题目传送门

余数求和

题目背景

数学题,无背景

题目描述

给出正整数n和k,计算G(n, k)=k mod 1 + k mod 2 + k mod 3 + … + k mod n的值,其中k mod i表示k除以i的余数。例如G(10, 5)=5 mod 1 + 5 mod 2 + 5 mod 3 + 5 mod 4 + 5 mod 5 …… + 5 mod 10=0+1+2+1+0+5+5+5+5+5=29

输入输出格式

输入格式:

两个整数n k

输出格式:

答案

输入输出样例

输入样例#1:
10 5
输出样例#1:
29

说明

30%: n,k <= 1000

60%: n,k <= 10^6

100% n,k <= 10^9


  分析:

  之前没怎么写过数论分块(蒟蒻并不会莫比乌斯反演),于是做道题练下手。

  因为$a\%b=a-b*\lfloor \frac{a}{b}\rfloor$,所以我们所求的式子$\sum^n_{i=1}k\mod i$可以转化为$n*k-\sum^n_i{i*\lfloor\frac{k}{i}\rfloor}$。

  和式的部分就可以用整除分块来做,复杂度就是$O(\sqrt{n})$的。

  Code:

//It is made by HolseLee on 7th Nov 2018
//Luogu.org P2261
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std; typedef long long ll;
ll n,k,ans; int main()
{
cin>>n>>k;
ans=n*k;
for(ll l=,r; l<=n; l=r+) {
if( l<=k ) r=min(k/(k/l),n);
else r=n;
ans-=(r-l+)*(r+l)*(k/l)/;
}
cout<<ans<<'\n';
return ;
}
上一篇:PHP无限极分类,多种方法|很简单,这里说的很详细,其它地方说的很不好懂


下一篇:ExtJS 数据模型