Min_25筛学习小记

前言

听说大家都会了$Min\_25$或州阁筛了，虚的一批的我马上学了一下。

<script type="math/tex; mode=display" id="MathJax-Element-14"></script>

Min_25筛

首先这种筛法可以用来筛某种积性函数的前缀和，当然也不一定要积性函数，某些特殊的函数也能筛，听说常数和空间都比州阁筛优秀，以绝对优势踩爆州阁筛，时间复杂度和州阁筛一样，也是$O(\frac{n^{\frac{3}{4}}}{log(\sqrt n)})$，当然本蒟蒻并不会证。

<script type="math/tex; mode=display" id="MathJax-Element-16"></script>

筛质数的函数值

现有一函数$F(x)$，求

$$\sum_{i=2}^nF(i)[i是质数]$$
举个例子，假设 $F(x)=i^k$，那我们就是要求
$$\sum_{i=2}^ni^k[i是质数]$$
现在用 $Min\_25$筛求上式的值，首先先求出 $\sqrt n$内的所有质数并排好序， $p_i$表示第 $i$个质数，设一共有$v$个质数。

再设

$$g(n,j)=\sum_{i=2}^ni^k[i的最小质因子大于p_j或者i是质数]$$

不难推出转移式

$$g(n,j)=\left\{ \begin{array}{} g(n,j-1)-p_j^k[g(\lfloor\frac{n}{p_j}\rfloor,j-1)-g(p_{j-1},j-1)]& &{p_j^2\leq n}\\ g(n,j-1)& &{p_j^2\geq n} \end{array}\right.$$

其中$g(p_{j-1},j-1)$显然等于$\sum_{i=1}^jp_i^k$

那我们要求的显然就为$g(n,v)$，我们可以发现，上述转移式都是由$g(?,j-1)$转移到$g(?,j)$，并且第一维可能的取值最多只有$\sqrt n$个，有时我们可以先预处理出第一维所有可能的取值，然后按照第二维分层算即可。

Code

//ask(n)=1^k+2^k+3^k+....+n^k
//quick(n)=n^k
//zh[n]=p[1]^k+p[2]^k+....+p[n]^k
#define fo(i,j,l) for(int i=j;i<=l;++i)
#define fd(i,j,l) for(int i=j;i>=l;--i)
ll l=1,r;
while(l<=n){
    ll len=n/l; r=n/len;
    if(len<=P)id1[len]=++m;else id2[r]=++m;
    g[m]=ask(len)-1; w[m]=len;
    //h[m]=len-1;
    l=r+1;
}
fo(i,1,v)
fo(j,1,m)if((ll)p[i]*p[i]>w[j])break;
else{
    int op=w[j]/p[i]<=P?id1[w[j]/p[i]]:id2[n/(w[j]/p[i])];
    g[j]=g[j]-quick(p[i])*(g[op]-zh[i-1]);
    //h[j]=h[j]-(h[op]-i+1);
}

同样的我们也能筛一个范围内质数的个数，代码中的$h$数组求解的即为质数个数。

积性函数求前缀和

求积性函数的前缀和有两种求法（某些非积性的特殊函数也可以），一个是递归版，一个是非递归版，各有各的长处。
现在对积性函数$F(x)$求前缀和，即求$\sum_{i=1}^nF(i)$。

递归版

设

$$S(n,i)=\sum_{i=2}^nF(i)[i的最小质因子不小于p_i]$$

于是我们可以很容易的得到

$$S(n,i)=g(n,v)-g(p_{i-1},i-1)+\sum_{j=i}^v\ \ \ \sum_{e\geq 1且p_j^{e+1}\leq n}F(p_j^e)S(\lfloor\frac{n}{p_j^e}\rfloor,j+1)+F(p_j^{e+1})$$

那所求值显然为$S(n,1)+F(1)$，实现用递归。

Code

//zh[n]=F(p[1])+F(p[2])+.....+F(p[n])
//P=sqrt(n)
inline ll ask(ll x,int k)
{
    if(x<=1||p[k]>x)return 0;
    int wz=x<=P?id1[x]:id2[n/x];
    ll ans=g[wz]-zh[k-1];
    for(int i=k;i<=v&&(ll)p[i]*p[i]<=x;++i)
    for(ll l=p[i],e=1;l*p[i]<=x;l=l*p[i],++e)
    ans=ans+ask(x/l,i+1)*F(l)+F(p[i]*l);
    return ans;
}
...
ans=ask(n,1)+F(1);
...

非递归版

方法和筛质数时的方法类似。

设

$$S'(n,i)=\sum_{i=2}^nF(i)[i的最小质因子不小于p_i或i是质数]$$

于是我们一样可以很容易的得到

$$S'(n,i)=\left\{ \begin{array}{} S'(n,i+1)+\sum_{e\geq 1且\ p_i^{e+1}\ \leq n}\ \ F(p_i^e)[S(\lfloor\frac{n}{p_i^e}\rfloor,i+1)-g(p_i,i)]+F(p_i^{e+1})& &{p_i^2\leq n} \\ S'(n,i+1)& &{p_i^2\geq n} \end{array}\right.$$

那所求值显然为$S'(n,1)+F(1)$，关于实现一样先离散出第一维所有可能的取值，同时我们也能观察到转移都是从$S'(?,i+1)$转移到$S'(?,i)$，一样以第二维分层计算即可。

Code

//zh[n]=F(p[1])+F(p[2])+.....+F(p[n])
//P=sqrt(n)
//id1[],id2[]均为离散数组
#define fd(i,j,l) for(int i=j;i>=l;--i)
fd(i,u,1){
    ll ww=ksm(p[i],k),U=(ll)p[i]*p[i];
    fo(j,1,m)if(U<=w[j]){
        for(ll o=p[i],u=w[j]/p[i];o*p[i]<=uu;o=o*p[i],u=u/p[i])
        {
            ok=(u<=P)?id1[u]:id2[n/u];
            g[j]=g[j]+(g[ok]-zh[i])*F(o1);
            g[j]=g[j]+F(o*p[i]);
        }
    }else break;
}
ans=g[1]+F(1);

<script type="math/tex; mode=display" id="MathJax-Element-46"></script>

关于两种版本的比较

首先如果仅是对一个$n$求解，那显然是第一个版本更快，因为它比第二个版本少遍历了很多无用状态。
但我们再看第二个版本，（设第一维所有可能的取值为$a_1,a_2,a_3....a_u$），我们可以发现第二个版本中不仅求出了$S'(n,1)$，还求出了$S'(a_1,1),S'(a_2,1),......S'(a_u,1)$，而这是第一个版本做不到的，因此用哪个版本决定于题目。
事实证明当$n=10^9$时，第一个版本比第二个版本快$\frac{1}{3}$甚至快$\frac{1}{2}$。

例题

Libre OJ #6053.简单的函数 Solution传送门
51Nod 1847 奇怪的数学题