{算法}乱谈线性筛法

最新推荐文章于 2023-10-13 22:15:24 发布

FarmerJohnLYH

最新推荐文章于 2023-10-13 22:15:24 发布

阅读量2.7k

点赞数 2

文章标签：算法 x 函数

本文链接：https://blog.csdn.net/FarmerJohnOfZS/article/details/51702456

版权

Concept

指在线性时间 $O(n)$ 用筛选的方法求一些东西的算法。
伪目录
1. 求质数表
2. 筛出欧拉函数φ

Ideas

As we all know, $x={p_1}^{a_1}·{p_2}^{a_2}·...·{p_n}^{a_n}$ 即 $x=\prod_{i=1}^n {p_i}^{a_i}$
其中 $p_i为第i个质数,a_i为其对应的质数$
我们还可以通过乘法定理得知， $x的约数共有\prod_{i=1}^n(a_i+1)$ 个
当然，这是后话了

For prime Numbers table

最简单的，我们使用线性筛法筛质数表
先来看看普通的筛法

fillchar(x,sizeof(x),true);//假设所有数都是质数(Tips x[1]=false)
for i:=2 to n do
if x[i] then  //x[i]表示i是否为质数
    for j:=2 to n div i do
    x[i*j]:=false;  //筛掉i*j(非质数)

显然复杂度为 $O(n^2)$
我们发现，有的数字被筛了多次，如24
当然最好的情况是每个数字只被筛1次

看Code之间，先回忆一句话

使每个数字只被筛1次

也许我们需要换个思路(线性筛法)
那线性筛法又是如何做到使每个数字只被筛1次的呢？
这就是我们所需要证明的

Code

fillchar(bz,sizeof(bz),true);
for i:=2 to n do//枚举数i
begin
    if(bz[i])then//i为素数
    begin
        inc(tot);
        s[tot]:=i;//素数数组
    end;
    for j:=1 to tot do
    begin 
        bz[i*s[j]]:=false;//筛去i*s[j]
        if(i mod s[j]=0)then break;//什么鬼？怎么就Break了？
    end;
end;

相信都对if(i mod s[j]=0)then break;一脸懵逼吧~
这也正是线性筛法的神奇所在

好了不说废话
用上 $x=\prod_{i=1}^n {p_i}^{a_i}$

$\because$ 为使令循环break掉的就是 $p_1$ (最小的素数)
$\therefore$ 不能让 $x被p_l(l>1) 筛掉$
$设x=p_1*a=p_l*b$

$p_l∗b=\prod_{i=1}^n {p_i}^{a_i}*b$
$=(u∗p_2∗p_3∗...∗p_{m−1}∗p_m)∗p_1$
$=a∗p_1$
显然 $a<b$

$\because$ a没有比 $p_1$ 要小的质因数
$\therefore$ 一定会枚举到 $b$ 前用 $a$ 将 $x=a*p_1=b*p_l$ 筛掉。
$\therefore$ 每个 $x$ 只会被 $p_1*a$ 筛掉

至于某些人问会不会筛掉质数……
我推荐给你一个小学奥数班！

筛出欧拉函数φ

读作φ(fai)
φ(n)为小于n并与n互质的数的个数,是积性函数
即 $如果gcd(a,b)=1那么φ(a*b)=φ(a)*φ(b)$
另 $φ(1)=1$

性质

$如果x为素数，那么φ(x)=x−1$
小于 $x$ 的数不都与 $x$ 互质吗？
$如果x为素数，那么φ(x^k)=x^k−\frac {x^k}{x} =x^{k−1}(x−1)$
显然只有 $x$ 的倍数才不互质 $x^k$
$φ(x)=\prod_{j=1}^k φ({p_j}^{a_j})$
$φ(x)$ 为积性函数

于是乎~

设 $x=p_1$
由性质2得
$①φ(x^k)=x^{k−1}(x−1)$
$②φ(x^{k+1})=x^k(x−1)$
$①代入②$
$φ(x^{k+1})=φ(x^k)·x$
$φ(x∗i)$
$=φ(x^{a1+1})∗φ({p_2}^{a_2})∗....∗φ({p_k}^{a_k})$
$=x∗φ(x^{a_1})∗φ({p_2}^{a_2})∗....∗φ({p_k}^{a_k})$
$=x∗φ(i)$
$若x|i那么 φ(i∗x)=φ(i)∗x$
$否则φ(i∗x)=φ(i)∗φ(x)=φ(i)∗(x−1)$

Code

恕暂不提供Pascal语言

int prime[M/3],phi[M];  
bool flag[M];  
void get_prime()  
{  
    int i,j,k;  
    memset(flag,false,sizeof(flag));  
    k=0;  
    for(i=2;i<M;i++){  
        if(!flag[i]){                              
            prime[k++]=i;  
            phi[i]=i-1;  
        }  
        for(j=0;j<k&&i*prime[j]<M;j++){  
            flag[i*prime[j]]=true;              
            if(i%prime[j]==0){  
                phi[i*prime[j]]=phi[i]*prime[j];  
                    break;  
            }  
            else  
                phi[i*prime[j]]=phi[i]*(prime[j]-1);  
        }  
    }  
}