欧拉筛&&埃氏筛

简介: 欧拉筛&&埃氏筛

数论——欧拉筛、埃氏筛

记录一点关于数论的知识,该知识点本身不难,主要是学习一下思想~

埃氏筛

埃拉托斯特尼筛法,简称埃氏筛。

一个质数的倍数一定是合数,所以,假设P PP是质数,我们可以筛掉区间[ 1 , 1 e 7 ] 中所有P的倍数。

对于数列1~11:

先筛去2的倍数(4、6、8、10),然后筛去3的倍数(6、9),然后筛去5的倍数(10)。

至此,1~11内的所有合数都被筛完了, 2 3 5 7 11是数列中的质数。

为什么这样能筛去所有的合数呢,因为一个合数一定能被分解为几个质数的幂的乘积,并且这个数的质因子一定是小于它本身的,所以当我们从小到大将每个质数的倍数都筛去的话,当遍历到一个合数时,它一定已经被它的质因子给筛去了。

但筛除3的倍数时,我们还是从3的2倍开始筛,其实3 * 2 ,已经被2 * 3时筛过了。

所以我们每一次只需要从i*i开始筛!此时时间复杂度可以近似看成O ( n )了。

埃氏筛代码:

public class 数论_埃氏筛 {
    static final int N = (int) (1e7 + 5);
    static int[] st = new int[N];
    public static void E_sieve(int n) {
        for (int i = 2; i <= n / i; i++)
            if (st[i] == 0)
                for (int j = i * i; j <= n; j += i)
                    st[j] = 1; // j是i的一个倍数,j是合数,筛掉。
    }
}

欧拉筛

埃氏筛是筛去每个质数的倍数,但难免,会有合数会被其不同的质因子多次重复筛去。这就造成了时间浪费。

所以,我们不需要用一个for循环去筛除一个质数的所有倍数,我们将所有质数存储到primes[]中,然后枚举到第i个数时,就筛去所有的primes[j] * i。这样就在每一次遍历中,正好筛除了所有已知素数的i倍。

但是为了确保合数只被最小质因子筛掉,最小质因子要乘以最大的倍数,即要乘以最大的i , 所以不能提前筛,所以如果 i % primes [ j ] = = 0 , 我们就结束内层循环!

欧拉筛代码:

public class 数论_欧拉筛 {
    static final int N = (int) (1e7 + 5);
    static int cnt = 0;
    static int[] st = new int[N], primes = new int[N];
    static void ola(int n) {
        for (int i = 2; i <= n; i++) {
            if (st[i] == 0)
                primes[cnt++] = i;//将质数存到primes中
            for (int j = 0; primes[j] <= n / i; j++) {//要确保质数的第i倍是小于等于n的。
                st[primes[j] * i] = 1;
                if (i % primes[j] == 0)
                    break;//欧拉筛的核心思想就是确保每个合数只被最小质因数筛掉。或者说是被合数的最大因子筛掉。
            }
        }
    }
    public static void main(String[] args) {
        ola(8);
        for (int i = 0; i < cnt; i++) {
            System.out.println(primes[i]);
        }
    }
}

总结

素数筛法的思想仍然是通过已有条件筛出能够事先拿出的答案,从而达到优化的效果,其思想值得我们学习!

文章粗浅,希望对大家有帮助!

相关文章
|
6月前
|
Java C++
筛法求质数
筛法求质数
54 0
|
12月前
|
C++
筛质数、分解质因数和快速幂的应用
筛质数、分解质因数和快速幂的应用
61 0
|
12月前
数论——高斯消元
数论——高斯消元
59 0
|
5月前
|
机器学习/深度学习 资源调度 算法
杜教筛
【6月更文挑战第9天】
27 2
|
6月前
|
算法 测试技术 C++
【数学归纳法 组合数学】容斥原理
【数学归纳法 组合数学】容斥原理
|
6月前
D - 11(逆元好题)
D - 11(逆元好题)
|
12月前
一个求公约数和公倍数的有趣求法
一个求公约数和公倍数的有趣求法
56 0
|
算法
质数筛法:朴素素数筛,埃氏筛,欧式筛
质数筛法:朴素素数筛,埃氏筛,欧式筛
146 0