min_25 筛是由 min_25 大佬使用后普遍推广的一种新型算法,这个算法能在 \(O({n^{3\over 4}\over log~ n})\) 的复杂度内解决所有的积性函数前缀和求解问题(个人感觉套上素数定理证明的复杂度的话应该要把下面的 log 改成 ln ,不过也差不多啦~)

其实 min_25 筛的入门TXC 大佬的 blog 已经写的非常棒了QVQ

所以搬博客的话鉴于博主太懒了就不干了...直接帮 TXC 大佬安利博客完事

这篇博客主要的目的是证明网上大多没有的 min_25 筛的复杂度

所以你首先得学会 min_25 筛...

前置芝士

Min_25 筛

素数定理

proof

首先我们知道 min_25 筛做了两件事:

  1. 筛出一个所求函数的质数前缀和

  2. 递归/循环 得到所求函数的前缀和

对于第一个部分和第二个部分我们发现都有一个枚举质数的过程,那么我们需要先知道 n 范围内的质数个数

这就需要素数定理了,素数定理指:在 n 以内的素数个数是 \(O({n\over ln~ n})\) 的(下面为了方便起见就用 \(log ~n\) 代替这里的 \(ln~n\) 了 )

至于具体证明就不是本博客的任务了 QVQ

我们可以根据这一点结合埃氏筛法的复杂度得知:

我们在进行 min_25 中处理的 \(2\sqrt n\) 个数的函数前缀和时,一个数 x 被筛到过 \(O({\sqrt x\over log~ \sqrt x})\) 次

并且,无论是筛出所求函数关于 x 的质数前缀和还是全部的前缀和,都是这个次数(如果你写过循环代替递归版的 min_25 的话,就更加清楚这是为什么了)

这里可以贴出我写的一种 min_25 筛的版本(以筛 mu 的前缀和为例):

inline int ID(Rg int x){
return x<=sq?id1[x]:id2[n/x];
}
inline int calc(Rg int k,int j){
return w[k]>=p[j]?h[k]-j+1:0;
}
inline void Min_25(ll n){
sq=sqrt(n),m=0; int tot=0;
for(tot=1;1ll*p[tot]*p[tot]<=n;++tot); for(Rg ll l=1;l<=n;l=w[m]+1){
w[++m]=n/(n/l),
h[m]=w[m]-1,Mu[m]=0,
//这里的 h 指质数个数的前缀和 , Mu 指莫比乌斯函数的前缀和
if(w[m]<=sq) id1[w[m]]=m;
else id2[w[m]]=m;
}
for(Rg int j=1;j<=tot;++j){
for(Rg int i=m;i&&1ll*p[j]*p[j]<=w[i];--i)
h[i]-=h[ID(w[i]/p[j])]-j+1;
}
for(Rg int j=tot;j;--j){
for(Rg int i=m,k;i&&1ll*p[j]*p[j]<=w[i];--i)
Mu[i]-=Mu[k]-calc(k,j+1);
}
for(Rg int i=m;i;--i) Mu[i]-=h[i];
}

不压行真的好难受嘤嘤嘤...

我们发现上面筛出 mu 前缀和的时候两个处理部分都是用了极其类似的双重循环,(并且保证了答案是对的...),那么也就可以证明上面的说法了

所以我们就可以列出复杂度的式子:

\[O(Min\_25)=\sum_{x=1} ^{\sqrt n} O({\sqrt x\over log~ \sqrt x}) + \sum_{x=1} ^{\sqrt n} O({\sqrt {n\over x}\over log~ \sqrt {n\over x}})
\]

主定理一波就是:

\[\sum_{x=1} ^{\sqrt n} O({\sqrt {n\over x}\over log~ \sqrt {n\over x}})
\]

然后开始积分:

\[O\Big( \int_1^{\sqrt n} {\sqrt {n\over x}dx\over log~ \sqrt {n\over x}} \Big)
\]

然后正常做法就难搞了,那么我们假一点,反正求的是大致复杂度,那么我们让下面的表达式直接变为 \(log~ n\)(这里把常数直接略去) 这样我们的任务就只剩下求上面表达式的积分了

至于这样为什么可行? 我们发现下面的表达式中有一个 log ,虽然我们让下面的表达式都变大了,但在 log 后,不过是少个常数的问题

如果这里有疑问的话,我们也可以把下面的表达式改成最小的 \(log~\sqrt{n}={1\over2}log~\sqrt{n}\),然后我们再把 \(1\over 2\) 忽略,这样的结果和上面是一样的...

于是乎原来的式子就变成了:

\[O\Big( {\int_1^{\sqrt n} \sqrt {n\over x}dx\over log~ n} \Big)
\]

\[O\Big( {\sqrt n \int_1^{\sqrt n} x^{-{1\over2}} dx\over\ log~ n} \Big)
\]

上面的那个积分已经是个多项式了:

\[O\Big( {\sqrt n ·( 2{\sqrt n}^{1\over2} )\over\ log~ n} \Big)
\]

\[O\Big( { n^{1\over 2} · n^{1\over 4} \over\ log~ n} \Big)
\]

\[O\Big( {n^{3\over 4} \over\ log~ n} \Big)
\]

这样整个 min_25 筛的复杂度就证明完毕啦~

关于 min_25 筛的入门以及复杂度证明的更多相关文章

  1. min_25筛入门

    目录 1.什么是min_25筛 2.前置知识 2.1.数论函数 2.2.埃拉托色尼筛 2.3.欧拉筛 3.min_25筛 3.1.计算质数贡献 3.2.计算总贡献 3.3.实现 4.例题 4.1.[L ...

  2. 51Nod1222 最小公倍数计数 数论 Min_25 筛

    原文链接https://www.cnblogs.com/zhouzhendong/p/51Nod1222.html 题意 给定 $a,b$, 求 $$\sum_{n=a}^b \sum_{i=1}^n ...

  3. LOJ6053 简单的函数 【Min_25筛】【埃拉托斯特尼筛】

    先定义几个符号: []:若方括号内为一个值,则向下取整,否则为布尔判断 集合P:素数集合. 题目分析: 题目是一个积性函数.做法之一是洲阁筛,也可以采用Min_25筛. 对于一个可以进行Min_25筛 ...

  4. [复习]莫比乌斯反演,杜教筛,min_25筛

    [复习]莫比乌斯反演,杜教筛,min_25筛 莫比乌斯反演 做题的时候的常用形式: \[\begin{aligned}g(n)&=\sum_{n|d}f(d)\\f(n)&=\sum_ ...

  5. Min_25 筛

    Min_25 筛 yyb好神仙啊 干什么用的 可以在\(O(\frac{n^{\frac 34}}{\log n})\)的时间内求积性函数\(f(x)\)的前缀和. 别问我为什么是这个复杂度 要求\( ...

  6. 关于min_25筛的一些理解

    关于min_25筛的一些理解 如果想看如何筛个普通积性函数啥的,就别往下看了,下面没有的(QwQ). 下文中,所有的\(p\)都代表质数,\(P\)代表质数集合. 注意下文中定义的最小/最大质因子都是 ...

  7. LOJ #6202. 叶氏筛法(min_25 筛)

    题意 求 \([L, R]\) 之间的素数之和 . \(L≤10^{10},2×10^{10} \le R \le 10^{11}\) 题解 一个有点裸的 min_25筛 ? 现在我只会筛素数的前缀和 ...

  8. 【LOJ#572】Misaka Network 与求和(莫比乌斯反演,杜教筛,min_25筛)

    [LOJ#572]Misaka Network 与求和(莫比乌斯反演,杜教筛,min_25筛) 题面 LOJ \[ans=\sum_{i=1}^n\sum_{j=1}^n f(gcd(i,j))^k\ ...

  9. LOJ572. 「LibreOJ Round #11」Misaka Network 与求和 [莫比乌斯反演,杜教筛,min_25筛]

    传送门 思路 (以下令\(F(n)=f(n)^k\)) 首先肯定要莫比乌斯反演,那么可以推出: \[ ans=\sum_{T=1}^n \lfloor\frac n T\rfloor^2\sum_{d ...

随机推荐

  1. 家庭记账本小程序之查(java web基础版六)

    实现查询消费账单 1.main_left.jsp中该部分,调用search.jsp 2.search.jsp,提交到Servlet中的search方法 <%@ page language=&qu ...

  2. 前端获取checkbox复选框的值 通过数组形式传递

    html代码: <form role="form" class="select_people"> <div style="displ ...

  3. [SimplePlayer] 8. 音视频同步

    音频与视频在播放当中可能会由于种种原因(如:音视频并非在同一时间开始播放,或视频由于解码任务繁重导致输出图像延迟等)导致音频与视频的播放时间出现偏差,这种就是音视频的同步问题,本文会对音视频同步进行讨 ...

  4. 第六十五天 js操作

    1.闭包 // 函数的嵌套定义,定义在内部的函数都称之为 闭包 // 1.一个函数要使用另一个函数的局部变量 // 2.闭包会持久化包裹自身的函数的局部变量 // 3.解决循环绑定 function ...

  5. css经常使用的六种文本样式

    css当中经常使用的六种文本样式 css 文本样式是相对于内容进行的样式修饰,下面来说下几种常见的文本样式. 首行缩进 首行缩进是将段落的第一行缩进,这是常用的文本格式化效果.一般地,中文写作时开头空 ...

  6. 【BZOJ5503】[GXOI/GZOI2019]宝牌一大堆(动态规划)

    [BZOJ5503][GXOI/GZOI2019]宝牌一大堆(动态规划) 题面 BZOJ 洛谷 题解 首先特殊牌型直接特判. 然后剩下的部分可以直接\(dp\),直接把所有可以存的全部带进去大力\(d ...

  7. dajngo cache,throttling

    缓存 背景介绍: 动态网站的问题就在于它是动态的. 也就是说每次用户访问一个页面,服务器要执行数据库查询,启动模板,执行业务逻辑以及最终生成一个你所看到的网页,这一切都是动态即时生成的. 从处理器资源 ...

  8. 题解-洛谷P1981 表达式求值(模拟+处理优先级的递归)

    https://www.luogu.org/problemnew/show/P1981 (原题链接) 显然乘法的优先级高与加法,因此碰到乘号就要优先把一连串与乘号相连的数算出,很容易想到递归.可用普通 ...

  9. qt: 获取sql数据表的所有的字段;

    1. mysql 数据库: 转载: https://www.cnblogs.com/fuqia/p/8994080.html mysql安装成功后可以看到已经存在mysql.information_s ...

  10. 金融量化分析【day112】:量化平台的使用-第一个策略

    一.策略代码 # 导入函数库 import jqdata #初始化函数,设定基准等等 def initialize(context): set_benchmark('000300.XSHG') g.s ...