求质素(素数)

举个例子:求2020(可泛指n)内的所有质素(素数),方法如下:
存在于个人总结,算法多种多样,看个人的思维想法吧


方法一:O(n^2)

最最最最最普通的了,就两个for循环,不想多说

int p[2020],cnt=0;
void Prime(){
    for(int i=2;i<=2020;i++){
        bool fl=true;
        for(int j=2;j<i;j++)
            if(i%j==0){
                fl=false;
                break;
            }
        if(fl)p[cnt++]=i;
    }
}

想要更快一点,无非是将 j<i 改为 j<√i


方法二:O(nloglogn)埃拉特斯特尼筛法

从2到2020(2020可泛指n)开始,先筛去n内所有2的倍数,然后每次从下一个剩下的数(必然为质数)开始,筛去其n内所有的倍数,最后剩下的数都是质数。
设置p[]储存素数
vis[]为标记数组,即0为未标记的,1为标记的,也就是某个数的倍数

int p[2020],vis[2021],cnt=0;
void Prime(){
    memset(vis,0,sizeof vis);
    for(int i=2;i<=2020;i++){
        if(!vis[i])p[cnt++]=i;
        for(int j=i*i;j<2020;j+=i)
            vis[j]=1;
    }
}

方法三:O(n)欧拉线性筛法

在第二种方法里面,会发现有合数(ps:合数即为素数的倍数,如:18=2*9=3*6)被重复筛选过了,先上代码吧

int p[2020],vis[2021],cnt=0;
void Prime(){
    memset(vis,0,sizeof vis);
    for(int i=2;i<=2020;i++){
        if(!vis[i])p[cnt++]=i;
        //cout<<"i="<<i<<endl;
        for(int j=0;j<cnt&&i*p[j]<=2020;j++){
            //cout<<" j="<<j<<" prime["<<j<<"]="<<prime[j]<<" i*prime[j]="<<i*prime[j]<<endl;
            v[i*p[j]]=1;
            if(i%p[j]==0)break;
        }
    }
}

来解释一下里面的一些代码意思吧
首先 i*p[j] 表示地是质素的倍数了,然后标记
最核心的也就是那个跳出语句 if(i%p[j]==0)break
为什么要跳出呢,看举例子吧(本弱弱喜欢举例子,比较好理解,哈哈哈哈

i = 2
  j = 1 prime[1] = 2 i*prime[j] = 4
 i = 3
  j = 1 prime[1] = 2 i*prime[j] = 6
  j = 2 prime[2] = 3 i*prime[j] = 9
 i = 4
  j = 1 prime[1] = 2 i*prime[j] = 8
 i = 5
  j = 1 prime[1] = 2 i*prime[j] = 10
  j = 2 prime[2] = 3 i*prime[j] = 15
  j = 3 prime[3] = 5 i*prime[j] = 25
 i = 6
  j = 1 prime[1] = 2 i*prime[j] = 12
 i = 7
  j = 1 prime[1] = 2 i*prime[j] = 14
  j = 2 prime[2] = 3 i*prime[j] = 21
  j = 3 prime[3] = 5 i*prime[j] = 35
  j = 4 prime[4] = 7 i*prime[j] = 49
 i = 8
  j = 1 prime[1] = 2 i*prime[j] = 16
 i = 9
  j = 1 prime[1] = 2 i*prime[j] = 18
  j = 2 prime[2] = 3 i*prime[j] = 27
 i = 10
  j = 1 prime[1] = 2 i*prime[j] = 20
 ……
  ……

有没有看出点什么?,没有的话,我把它再单独拿出来

i = 2
  j = 1 prime[1] = 2 i*prime[j] = 4
 i = 4
  j = 1 prime[1] = 2 i*prime[j] = 8
 i = 8
  j = 1 prime[1] = 2 i*prime[j] = 16

很明显了吧,就是当 i 恰好是 p[] 的倍数的时候,就会标记 i*p[],并且跳出
i*p[] 被标记时,下一次就会标记它的倍数的,自然就只被标记了一次
倘若不跳出,就会出现合数被多次标记的情况,如:

 i = 2
  j = 1 prime[1] = 2 i*prime[j] = 4
 i = 3
  j = 1 prime[1] = 2 i*prime[j] = 6
  j = 2 prime[2] = 3 i*prime[j] = 9
 i = 4
  j = 1 prime[1] = 2 i*prime[j] = 8
  j = 2 prime[2] = 3 i*prime[j] = 12
 i = 5
  j = 1 prime[1] = 2 i*prime[j] = 10
  j = 2 prime[2] = 3 i*prime[j] = 15
  j = 3 prime[3] = 5 i*prime[j] = 25
 i = 6
  j = 1 prime[1] = 2 i*prime[j] = 12
  j = 2 prime[2] = 3 i*prime[j] = 18
  j = 3 prime[3] = 5 i*prime[j] = 30
 i = 7
  j = 1 prime[1] = 2 i*prime[j] = 14
  j = 2 prime[2] = 3 i*prime[j] = 21
  j = 3 prime[3] = 5 i*prime[j] = 35
  j = 4 prime[4] = 7 i*prime[j] = 49
 i = 8
  j = 1 prime[1] = 2 i*prime[j] = 16
  j = 2 prime[2] = 3 i*prime[j] = 24
  j = 3 prime[3] = 5 i*prime[j] = 40
  j = 4 prime[4] = 7 i*prime[j] = 56
 i = 9
  j = 1 prime[1] = 2 i*prime[j] = 18
  j = 2 prime[2] = 3 i*prime[j] = 27
  j = 3 prime[3] = 5 i*prime[j] = 45
  j = 4 prime[4] = 7 i*prime[j] = 63

再次提出来看吧

i = 4
  j = 1 prime[1] = 2 i*prime[j] = 8
  j = 2 prime[2] = 3 i*prime[j] = 12
 i = 5
  j = 1 prime[1] = 2 i*prime[j] = 10
  j = 2 prime[2] = 3 i*prime[j] = 15
  j = 3 prime[3] = 5 i*prime[j] = 25
 i = 6
  j = 1 prime[1] = 2 i*prime[j] = 12
  j = 2 prime[2] = 3 i*prime[j] = 18
  j = 3 prime[3] = 5 i*prime[j] = 30
现在应该懂了吧,over再者还有很多方法,就先到这吧
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 161,780评论 4 369
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 68,424评论 1 305
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 111,397评论 0 254
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 44,576评论 0 218
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 52,997评论 3 295
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 40,945评论 1 224
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 32,107评论 2 317
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 30,850评论 0 208
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 34,625评论 1 250
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 30,804评论 2 253
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 32,285评论 1 265
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 28,613评论 3 261
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 33,291评论 3 242
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 26,164评论 0 8
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 26,963评论 0 201
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 36,096评论 2 285
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 35,886评论 2 278