ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

数论分块

2022-02-24 11:00:09  阅读:128  来源: 互联网

标签:lfloor le 分块 数论 ll rfloor ab ans


前置知识

  • 对于\(N^*=\{\lfloor n/i \rfloor|i\in[1,n]\}\).
    (1) 满足 \(\lfloor n/i \rfloor,x\in N^*\) 的最大的 \(i\) 为 \(\lfloor n/x \rfloor\).

(1)证明:
显然有 \(x\le n/i < x+1\)

\(\Leftrightarrow ix \le n < i(x+1)\)

\(\Leftrightarrow n/(x+1) < i \le n/x\)

又由于 \(i\in Z\) ,所以 \(\lfloor n/(x+1) \rfloor +1 \le i \le \lfloor n/x \rfloor\)
Q.E.D.

稍微变一下就是 \(\lfloor n/\lfloor n/i \rfloor \rfloor\)

而且我们证明最后得到的结论也很好用:

\[\lfloor n/(x+1) \rfloor +1 \le i \le \lfloor n/x \rfloor \]

然后给出一个关于向下取整的性质:

  • \(\lfloor n/ab \rfloor=\lfloor \lfloor n/a \rfloor/b \rfloor\)

证明:
令 \(a/b=\lfloor n/a \rfloor +s\ (0 \le s < 1)\)
那么 \(\lfloor n/ab \rfloor=\lfloor 1/b(\lfloor n/a \rfloor+s) \rfloor =\lfloor \lfloor n/a \rfloor b \rfloor\)

对此我们可以计算 \(\sum_{i=1}^{n}\lfloor n/i \rfloor\)

OI wiki上的图明了的解释了这一切

1

发现函数值是连续一段一段的,直接有了写方的代码。

ll l=1,r=0,ans=0;
while(l<=n){
	r=n/(n/l);
	ans=ans+(r-l+1)*(n/l);
	l=r+1;
}
printf("%d\n",ans);

整数分块

用来计算形如 \(\sum_{i=1}^{n}f(i)g(\lfloor \frac{n}{i} \rfloor)\) 的和式。

实际上是上面的那个东西的推广。

其实没多大差别,我们只需要能够求出一段的连续区间的权值和即可,这个我们用前缀和即可。

ex:

ll l=1,r=0,ans=0;
while(l<=n){
	r=n/(n/l);
	ans=ans+(sum[r]-sum[l-1])*(n/l);
	l=r+1;
}
printf("%d\n",ans);

没了,结束。

标签:lfloor,le,分块,数论,ll,rfloor,ab,ans
来源: https://www.cnblogs.com/cbdsopa/p/15930746.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有