欧拉筛&&埃氏筛

不管现实多么惨不忍睹,都要持之以恒地相信,这只是黎明前短暂的黑暗而已。不要惶恐眼前的难关迈不过去,不要担心此刻的付出没有回报,别再花时间等待天降好运。真诚做人,努力做事!你想要的,岁月都会给你。欧拉筛&&埃氏筛,希望对大家有帮助,欢迎收藏,转发!站点地址:www.bmabk.com,来源:原文

数论——欧拉筛、埃氏筛

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



埃氏筛

埃拉托斯特尼筛法,简称埃氏筛。
一个质数的倍数一定是合数,所以,假设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]);
        }
    }
}

总结

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

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

参考博客:素数筛法详解:埃氏筛和欧拉筛

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

文章由极客之音整理,本文链接:https://www.bmabk.com/index.php/post/159828.html

(0)
飞熊的头像飞熊bm

相关推荐

发表回复

登录后才能评论
极客之音——专业性很强的中文编程技术网站,欢迎收藏到浏览器,订阅我们!