[bzoj 1257][洛谷 2261] [CQOI2007]余数求和 {约数+除数分块}

最新推荐文章于 2019-10-17 20:41:29 发布

原创最新推荐文章于 2019-10-17 20:41:29 发布 · 300 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#bzoj #洛谷 #CQOI #2007

数论（数学）同时被 2 个专栏收录

81 篇文章

订阅专栏

分块

6 篇文章

订阅专栏

本文详细解析了洛谷平台上的P2261题目，通过数学推导和算法优化，提供了高效的解题思路。介绍了如何利用除数分块的方法减少计算复杂度，给出了解题代码。

题目

https://www.luogu.org/problemnew/show/P2261

解题思路

我们注意到 $k\ mod\ i=k-\left \lfloor \frac{k}{i} \right \rfloor*i$ ，可以转换成 $n*k-\sum_{i=1}^{n}\left \lfloor \frac{k}{i} \right \rfloor*i$
对于任意整数 $x\in \left [ 1,k \right ]$ ，设 $g(x)=\left \lfloor k\div \left \lfloor \frac{k}{x} \right \rfloor \right \rfloor$ 。所以 $g(x)\geq \left \lfloor k\div \frac{k}{x} \right \rfloor=x$ ，故 $\left \lfloor k\div g(x) \right \rfloor\leq \left \lfloor k/x \right \rfloor$ 。
另外 $\left \lfloor k/ g(x) \right \rfloor\geq \left \lfloor k/(k/\left \lfloor k/x \right \rfloor) \right \rfloor =\left \lfloor k/k*\left \lfloor k/x \right \rfloor \right \rfloor =\left \lfloor k/x \right \rfloor$ 故 $\left \lfloor k/g(x) \right \rfloor=\left \lfloor k/x \right \rfloor$ 可得 $\forall i \in \left [ x,\left \lfloor k \div \left \lfloor \frac{k}{x} \right \rfloor \right \rfloor \right ],\left \lfloor k/i \right \rfloor的值相等。$

如果还不懂，可以参考 https://www.luogu.org/problemnew/solution/P2261

代码`除数分块`

#include<cstdio>
#include<algorithm>
#define ll long long 
using namespace std;
ll n,k,ans; 
int main(){
	scanf("%lld%lld",&n,&k); ans=n*k;
	for(register int x=1,gx;x<=n;x=gx+1){
		gx=k/x?min(k/(k/x),n):n; 
		ans-=(k/x)*(x+gx)*(gx-x+1)/2; 
	}
	printf("%lld",ans); 
}