积性函数与莫比乌斯反演
文章平均质量分 85
积性函数,筛法,推公式。
KXL5180
这个作者很懒,什么都没留下…
展开
-
HDU6834 Yukikaze and Smooth numbers
题目链接 题意:给出一个n和k,要求输出1到n中有多少个数的最大质因子都小于等于k T组输入,n和k的范围都是1e9. 做法:如果n≤kn \leq kn≤k 那么答案是nnn. 如果 k<nk< nk<n: 首先考虑一个容斥很容易得到答案的表达如下: ans=n−∑i1pi+∑i<j1pipj−...(−1)k1p1p2...pk(min(p)>k)ans =n -\sum_{i}\frac{1}{p_i}+\sum_{i<j}\frac{1}{p_ip_j} -..原创 2020-08-07 20:15:36 · 703 阅读 · 0 评论 -
UOJ188【UR #13】Sanrd (Min_25筛)
这里 题意:说白了就是求这个∑i=1npsmax(i)\sum_{i=1}^{n}p_{smax}(i)∑i=1npsmax(i) 其中psmaxp_{smax}psmax表示次小质因子,规定1和质数的次小质因子为0 做法:这道题是一道min25筛法的题,但不过不是普通的积性函数求和,需要对该算法有一定的理解。 首先我们还是以s(n,j)s(n,j)s(n,j)表示在1−n1-n1−n范围内...原创 2020-02-18 15:28:22 · 371 阅读 · 0 评论 -
Min_25筛 代码及模板
LOJ简单的函数 以这道题为例子,讲一下min_25筛如何用代码实现。 首先min_25是一种亚线性筛,可以处理1e9以上的数据,老版min_25筛的复杂度为O(n0.75logn)O(\frac{n^{0.75}}{logn})O(lognn0.75),新版已经优化到O(n23)O(n^{\frac{2}{3}})O(n32),复杂度,我也不会证明。 据说常数比较大,确实是跑得快,但不过差距...原创 2020-01-19 15:54:17 · 603 阅读 · 0 评论 -
HDU6683 Rikka with Geometric Sequence 多校九(推导+杜教筛+分块)
http://acm.hdu.edu.cn/showproblem.php?pid=6683 题意:问1-n这些数字中有多少子序列是等比数列。 做法:这道题有点悬。。。。。 我们接下来来一下,数学推导(瞎JB乱搞): 我们设等比数列的公比,等比数列的长度为,首项末项 这个必定是一个整数,所以可以得到,这个是显而易见的。但不过这个好像不好弄。 因为如果你统计p的数目,这只是n范围内的...原创 2019-08-20 17:11:47 · 410 阅读 · 0 评论 -
BZOJ 4407: 于神之怒加强版 积性函数线性筛+反演
https://www.lydsy.com/JudgeOnline/problem.php?id=4407 题意:很简单不说了。 做法:废话不多说我们开始开心的推公式吧 首先为了方便我把表示,x=min(n,m),首先我们可以得到: 我们枚举gcd可以得到如下式子: 我们发现后面的那部分可以直接套用莫比乌斯反演的结论: 我们在枚举t可以得到如下式子: 我们发现td不好...原创 2019-07-27 09:16:55 · 202 阅读 · 0 评论 -
P1829 [国家集训队]Crash的数字表格 / JZPTAB
https://www.luogu.org/problem/P1829 题意:下面一起说 做法: 我们要求的就是下面这个: 首先讲一个东西,关于gcd的题多做几道后不要再套公式了,直接记住这个就可以了 然后我们开心的推公式吧. 我们把枚举d提到前面去 然后我们看见后面的一部分我们直接除以d可以得到 为什么会多一个d出来呢,因为i,j都除以了d所以后面乘以一个d^2 接下...原创 2019-07-26 14:48:47 · 191 阅读 · 0 评论 -
P4980 【模板】Polya定理
https://www.luogu.org/problem/P4980 题意:很简单的,就一个最简单的polya定理。 做法:网上随便看看吧,写得都不错。 我们要求的就是上面的式子,我们容易知道|G|=n,|B|=n; 但不过重点就是c(g); 但不动的时候c(g)=n这点没有毛病吧,因为都没有动。 但动k个点时候,这个时候就要看k了,随便找一找规律就知道这个就是 gcd(n,i...原创 2019-10-10 21:08:32 · 276 阅读 · 0 评论 -
Grisaia (推公式) 2018 第十届四川省程序设计竞赛
https://www.oj.swust.edu.cn/problem/show/2810 题意:很简单就是求题目上面那个式子的和。 做法:直接开始推公式。 然后,将前后两部分分开 然后对于前面的一部分直接可以由公式计算得到: 我们现在考虑后面一部分怎么计算。我刚刚拿到的时候一点头绪都没有。 但不过我们做题做多了,经常可以发现这种式子一般可以化成,两个上界为n的和...原创 2019-08-10 15:59:35 · 377 阅读 · 5 评论 -
P3768 简单的数学题(杜教筛+欧拉函数反演或者莫比乌斯反演)
https://www.luogu.org/problem/P3768 题意很简单就是求这个: 为了方便我们用(i,j)表示gcd(i,j),然后开始快乐的推公式吧: 看见后面那一坨直接莫比乌斯反演: 然后根据套路枚举t: 后面一部分就是等差数列通项公式,所以后面就是变形为 然后再根据套路令T=td,枚举T: 我们观察这部分刚好是狄利克雷卷积的标准形式 因此这部...原创 2019-08-08 20:09:46 · 313 阅读 · 2 评论 -
P2257 YY的GCD 莫比乌斯反演
https://www.luogu.org/problem/P2257 感觉这道题还是很难,公式推到一半觉得推不动了,而且觉得那复杂度还是很高。 首先设函数如下: 我们开始推公式: 我们的答案ans 我们令d/n=t我们可以得到 得到这样的公式已经很可以了,但不过还不能得到答案,因为这样复杂度还很高 我们另tp=T 然后我们换一下枚举的变量,枚举变量T ...原创 2019-07-25 20:37:12 · 246 阅读 · 0 评论 -
莫比乌斯反演简单证明+莫比乌斯反演入门题Sky Code
首先莫比乌斯反演 定理:和都是算术函数,并且满足如下条件: 则有如下结论: 首先我们说一下莫比乌斯函数 1.首先当n=1时,并且表示当n=1时才为1 2.当时,为互异素数时, 3.只要当n含有任何素因子的幂次大于等于2,则函数值为0; 4. 这里从定义上面去证明: 根据定义以及基本运算法则: 稍微思考一下可以得到(d,e)就是(d*e)|n的一个二元组合,因此可以交换两个sum...原创 2019-07-25 16:56:05 · 302 阅读 · 0 评论