电子商务网站规划,网站建设基础条件,网站开发兼职合同,网络推广服务费素数筛
素数筛的作用是筛选出[2,N]范围内的所有素数#xff0c;本次主要讲解两种方法#xff0c;分别是埃氏筛和欧拉筛。证明时会提到唯一分解定理#xff0c;如果不知道的小伙伴可以先去学一学#xff0c;那我们开始啦#xff01;
1.埃氏筛
主要思想#xff1a;当找到…素数筛
素数筛的作用是筛选出[2,N]范围内的所有素数本次主要讲解两种方法分别是埃氏筛和欧拉筛。证明时会提到唯一分解定理如果不知道的小伙伴可以先去学一学那我们开始啦
1.埃氏筛
主要思想当找到一个素数时利用该素数把该素数的所有倍数筛掉。
时间复杂度 O ( n l o g ( l o g ( n ) ) ) O(nlog(log(n))) O(nlog(log(n)))
上代码 //每个数的最小质因子//pre[i]表示i的最小质因子book[1] 1;//记录是否为素数1表示不是素数book[0] 1;for(int i 2;ibook.length;i) {if(book[i]0) {// i是素数筛掉素数的倍数 i2 6 222for(int j ii;jbook.lengthj0;ji) {book[j] 1;} }}问题 为什么遍历到i时若i没有被标记为合数也就是没有被i前面的数筛掉则一定是素数 为什么for循环遍历到sqrt(N)就可以了
先自己想一想哦提示是唯一分解定理。
答案
还记得唯一分解定理吗一个正整数可以用若干个质数表示假设当前正整数是n它可以用质数 p 1 , p 2 . . . p k p_1,p_2...p_k p1,p2...pk表示 p 1 , p 2 . . . p k p_1,p_2...p_k p1,p2...pk一定比q小。假设q是合数那么遍历到qq一定会被 p 1 , p 2 . . . p k p_1,p_2...p_k p1,p2...pk筛掉。如果q是质数呢他只能写出1*q的形式它会被自己筛掉。其实也就是证明sqrt(N)后面的合数一定会被小于sqrt(N)的数筛掉。设 N n N \sqrt{N}nN N nN且 a ∗ b n a*bn a∗bn,若ab则 a n N a\sqrt{n}\sqrt{N} an N 若a是素数则n会被a筛掉若a是合数则a可以继续分解为更小的素数而a和n都会被这个更小的素数筛掉所以即便 N n \sqrt{N}n N n但是仅用小于 N \sqrt{N} N 的数就可以把n筛掉所以可以遍历到sqrt(N)。
2.欧拉筛
主要思想埃氏筛的一部分时间耗在了重复的筛某些合数比如18会被2和3筛掉。欧拉筛保证每个合数只被筛一次因此也保证了 O ( n ) O(n) O(n)的时间复杂度。
时间复杂度 O ( n ) O(n) O(n)
上代码 int count 0;for (int i 2; i 20000005; i) {//线性if (!visit[i]) {//如果i是一个质数prime[count] i;//记录当前已经找出来的所有的质数}for (int j 0; j count i * prime[j] 20000005; j) {visit[i * prime[j]] true;//用prime[j]筛掉了i * prime[j]。if (i % prime[j] 0) break;//保证每个合数只被最小的质因子筛掉}}问题
为什么if语句满足后可以提前退出循环两个for循环嵌套如何实现的线性复杂度
先自己想一想哦提示是prime[j]是i的因子你可以把式子写出来看看。
再讲答案之前先来捋一捋欧拉筛的结构因为它不像埃氏筛那么直接。
首先一个for循环接着如果当前的i是素数则用另一个数组prime存一下这个数组只存素数。
再来一个for循环这个for循环就是用来筛合数的遍历之前找到的所有素数然后筛掉 p r i m e [ j ] ∗ i prime[j]*i prime[j]∗i。当满足if语句时这一轮的筛合数可以提前退出了。
答案
若此时if语句条件满足了则prime[j]是i的因子因此有 i k ∗ p r i m e [ j ] ik*prime[j] ik∗prime[j]。如果此时没有退出for循环会有 p r i m e [ j 1 ] ∗ i prime[j1]*i prime[j1]∗i被prime[j1]筛掉。 p r i m e [ j 1 ] ∗ i p r i m e [ j 1 ] ∗ k ∗ p r i m e [ j ] k ‘ ∗ p r i m e [ j ] prime[j1]*iprime[j1]*k*prime[j]k^*prime[j] prime[j1]∗iprime[j1]∗k∗prime[j]k‘∗prime[j]这说明了什么说明被prime[j1]筛掉的 p r i m e [ j 1 ] ∗ i prime[j1]*i prime[j1]∗i也会被prime[j]筛掉这就重复筛了怎么办我们让每个数都被其最小的质因子筛掉那么这里prime[j]就是 p r i m e [ j 1 ] ∗ i prime[j1]*i prime[j1]∗i最小的质因子因此j就不继续增大了直接退出该循环。因为保证了每个数只被筛一次第二个for循环总共被执行n次所有的数被筛完代码也就结束了。
例题
埃氏筛——最小质因子之和
参考代码
import java.io.BufferedReader;
import java.io.IOException;
import java.io.InputStreamReader;
import java.io.StreamTokenizer;
import java.util.Scanner;public class 最小质因子之和Easy {
public static void main(String[] args) throws IOException{//进行预处理f();//求2-n每个数对应的最小质因子sum();//求前缀和数组StreamTokenizer sc new StreamTokenizer(new BufferedReader(new InputStreamReader(System.in))); Scanner scanner new Scanner(System.in);sc.nextToken();int t (int)sc.nval;while(t-- 0) {sc.nextToken();int n (int)sc.nval;System.out.println(res[n]);}
}
static int book[] new int[4000000];
static int pre[] new int[4000000];
private static void f() {//埃氏筛模板//每个数的最小质因子//pre[i]表示i的最小质因子book[1] 1;//记录是否为素数1表示不是素数book[0] 1;for(int i 2;ibook.length;i) {if(book[i]0) {// i是素数筛掉素数的倍数 i2 6 222pre[i] i;//求的是质数的最小质因子for(int j ii;jbook.lengthj0;ji) {if(book[j]0) {pre[j] i;}book[j] 1;}}}}
static long res[] new long[4000000];
private static void sum() {//一次求出i 2- n// 2-i的最小质因子之和前缀和数组可以在Onfor(int i2;ires.length;i) {res[i] res[i-1]pre[i];}
}
}欧拉筛——最小质因子之和困难版
参考代码
import java.io.BufferedReader;
import java.io.IOException;
import java.io.InputStreamReader;public class 最小质因子之和Hard {public static void main(String[] args) throws IOException {BufferedReader br new BufferedReader(new InputStreamReader(System.in));int[] prime new int[20000005];int[] f new int[20000005];boolean[] visit new boolean[20000005];int count 0;for (int i 2; i 20000005; i) {//线性if (!visit[i]) {//如果i是一个质数prime[count] i;f[i] i;}for (int j 0; j count i * prime[j] 20000005; j) {visit[i * prime[j]] true;f[i * prime[j]] prime[j];if (i % prime[j] 0) break;}}long[] sum new long[20000005];//前缀和数组for (int i 2; i f.length; i) {//System.out.println(f[i]);sum[i] sum[i - 1] f[i];}int t Integer.parseInt(br.readLine());while (t-- 0) {int n Integer.parseInt(br.readLine());System.out.println(sum[n]);}}
}
文章转载自： http://www.morning.hqmfn.cn.gov.cn.hqmfn.cn http://www.morning.zymgs.cn.gov.cn.zymgs.cn http://www.morning.jbkcs.cn.gov.cn.jbkcs.cn http://www.morning.webpapua.com.gov.cn.webpapua.com http://www.morning.ndmbz.cn.gov.cn.ndmbz.cn http://www.morning.ytrbq.cn.gov.cn.ytrbq.cn http://www.morning.nfbkz.cn.gov.cn.nfbkz.cn http://www.morning.fosfox.com.gov.cn.fosfox.com http://www.morning.nlzpj.cn.gov.cn.nlzpj.cn http://www.morning.srgyj.cn.gov.cn.srgyj.cn http://www.morning.qlwfz.cn.gov.cn.qlwfz.cn http://www.morning.dndjx.cn.gov.cn.dndjx.cn http://www.morning.pngph.cn.gov.cn.pngph.cn http://www.morning.sbjhm.cn.gov.cn.sbjhm.cn http://www.morning.wphfl.cn.gov.cn.wphfl.cn http://www.morning.kqzrt.cn.gov.cn.kqzrt.cn http://www.morning.hdscx.cn.gov.cn.hdscx.cn http://www.morning.mdjzydr.com.gov.cn.mdjzydr.com http://www.morning.dxrbp.cn.gov.cn.dxrbp.cn http://www.morning.kmbgl.cn.gov.cn.kmbgl.cn http://www.morning.ppghc.cn.gov.cn.ppghc.cn http://www.morning.qmpbs.cn.gov.cn.qmpbs.cn http://www.morning.divocn.com.gov.cn.divocn.com http://www.morning.kqpsj.cn.gov.cn.kqpsj.cn http://www.morning.gstg.cn.gov.cn.gstg.cn http://www.morning.wnjbn.cn.gov.cn.wnjbn.cn http://www.morning.kwqt.cn.gov.cn.kwqt.cn http://www.morning.jcffp.cn.gov.cn.jcffp.cn http://www.morning.bfgbz.cn.gov.cn.bfgbz.cn http://www.morning.wqjpl.cn.gov.cn.wqjpl.cn http://www.morning.kwdfn.cn.gov.cn.kwdfn.cn http://www.morning.gblrn.cn.gov.cn.gblrn.cn http://www.morning.bpwfr.cn.gov.cn.bpwfr.cn http://www.morning.rqsnl.cn.gov.cn.rqsnl.cn http://www.morning.rgnq.cn.gov.cn.rgnq.cn http://www.morning.nlglm.cn.gov.cn.nlglm.cn http://www.morning.ai-wang.cn.gov.cn.ai-wang.cn http://www.morning.tmbfz.cn.gov.cn.tmbfz.cn http://www.morning.qnbck.cn.gov.cn.qnbck.cn http://www.morning.rqpgk.cn.gov.cn.rqpgk.cn http://www.morning.gklxm.cn.gov.cn.gklxm.cn http://www.morning.rjmg.cn.gov.cn.rjmg.cn http://www.morning.cpwmj.cn.gov.cn.cpwmj.cn http://www.morning.wlxfj.cn.gov.cn.wlxfj.cn http://www.morning.xjwtq.cn.gov.cn.xjwtq.cn http://www.morning.sqskm.cn.gov.cn.sqskm.cn http://www.morning.wfkbk.cn.gov.cn.wfkbk.cn http://www.morning.qtryb.cn.gov.cn.qtryb.cn http://www.morning.jnvivi.com.gov.cn.jnvivi.com http://www.morning.plwfx.cn.gov.cn.plwfx.cn http://www.morning.dyhlm.cn.gov.cn.dyhlm.cn http://www.morning.ymfzd.cn.gov.cn.ymfzd.cn http://www.morning.dhrbj.cn.gov.cn.dhrbj.cn http://www.morning.haolipu.com.gov.cn.haolipu.com http://www.morning.svtxeu.com.gov.cn.svtxeu.com http://www.morning.ljxps.cn.gov.cn.ljxps.cn http://www.morning.kflpf.cn.gov.cn.kflpf.cn http://www.morning.xglgm.cn.gov.cn.xglgm.cn http://www.morning.fthqc.cn.gov.cn.fthqc.cn http://www.morning.bysey.com.gov.cn.bysey.com http://www.morning.zkgpg.cn.gov.cn.zkgpg.cn http://www.morning.xqnzn.cn.gov.cn.xqnzn.cn http://www.morning.tbnn.cn.gov.cn.tbnn.cn http://www.morning.fhrt.cn.gov.cn.fhrt.cn http://www.morning.rccbt.cn.gov.cn.rccbt.cn http://www.morning.tqdqc.cn.gov.cn.tqdqc.cn http://www.morning.zyytn.cn.gov.cn.zyytn.cn http://www.morning.gfprf.cn.gov.cn.gfprf.cn http://www.morning.brwei.com.gov.cn.brwei.com http://www.morning.hcsqznn.cn.gov.cn.hcsqznn.cn http://www.morning.rnzjc.cn.gov.cn.rnzjc.cn http://www.morning.lcwhn.cn.gov.cn.lcwhn.cn http://www.morning.rfwrn.cn.gov.cn.rfwrn.cn http://www.morning.zqcdl.cn.gov.cn.zqcdl.cn http://www.morning.kzcfp.cn.gov.cn.kzcfp.cn http://www.morning.gxhqt.cn.gov.cn.gxhqt.cn http://www.morning.cbnxq.cn.gov.cn.cbnxq.cn http://www.morning.knlyl.cn.gov.cn.knlyl.cn http://www.morning.rdkqt.cn.gov.cn.rdkqt.cn http://www.morning.kmqms.cn.gov.cn.kmqms.cn
