码迷,mamicode.com
首页 > 其他好文 > 详细

【模板】素数筛法

时间:2019-09-25 15:45:30      阅读:84      评论:0      收藏:0      [点我收藏+]

标签:int   停止   筛选   div   namespace   怎么   表示   有一个   stream   

(好-----快------)

素数是什么就不用介绍了吧。。。先介绍判断素数的方法

判断素数

先看朴素算法

for (i=2;i<n;i++) 
if (n%i==0) break;
if (i==n) printf(“n is prime”);
else printf(“n is not prime”);

(真的好朴素。。)

用时O(n)

(肯定不行啊,吃枣药丸的。。)

怎么优化呢?

不难发现,如果a是n的约数,那么n/a也是n的约数

所以就有以下代码:

n2=int(sqrt(n)+0.5);
for (i=2;i<=n2;i++) if (n%i==0) break;
if (i>n2) printf(“n is prime”);
else printf(“n is not prime”);

这个用时O(n^0.5)

一般使用的话这样其实差不多了。。(其实有一个更快的只是我不会。。)

筛选素数

进入正题啦!

目标:

尽量快的找出1---n之间的所有素数

朴素算法:

for(int i=2;i<=n;i++)
    {
        for(int j=2;j<=sqrt(i);j++)
        {
            if(i%j==0)
            {
               flag=1;
               break;
            }
        }
        if(!flag)
        cout<<i<<" ";
        flag=0;
    } 

其实也不是很朴素啦。。用了上面介绍的优化,但是速度还是很慢

怎么优化呢?

对于任意数a,a所有的整数倍都是合数

int flag[2333]; int prime[2333]; int primes;
for (int i=2;i<=n;i++) flag[i]=0;
for (int i=2;i<=n;i++) for (int j=i*i;j<=n;j+=i) flag[j]=1; //表示j被标记为合数
for (int i=2;i<=n;i++) if (flag[i]==0) prime[primes++]=i;

有了这个想法,还可以继续优化

对于质数a,a的所有整数倍都是合数

int flag[2333]; int prime[2333]; int primes;
for (int i=2;i<=n;i++) 
if (flag[i]==0) 
{
    prime[primes++]=i;
    for (int j=i*i;j<=n;j+=i) flag[j]=1;
}

还可以继续优化吗?

其实还可以进一步优化,只要每个合数都只被自己能整除的最小的质数标记一次,速度就可以达到线性级别(好快!)

说白了就是每个质数尽可能多的标记合数(也不是很好懂的样子?)

看代码吧:

#include<iostream>
using namespace std;
int prime[100001];
int f[100001],primes;
int n;
int main()
{
    cin>>n;
    for(int i=2;i<=n;i++)
    {
        if(!f[i])//如果i还没有被标记,就说明i是新出现的素数 
        {
            prime[primes]=i;
            f[i]=-1;
            primes++;//标记 
        }
        for(int j=0;j<primes;j++)//枚举当前找到的所有素数 
        {
            int ans=prime[j]*i;
            if(ans>n)
            break;//越界就停止 
            f[ans]=j+1;//标记每一个含有prime[j] 这个因数的合数 
            if(i%prime[j]==0)//如果当前 i可以整除prime[j],就说明后面的prime[j]都不会是ans的最小质因子了 
            break;
        }
    }
    cout<<primes<<endl;//素数个数 
    for(int i=0;i<primes;i++)
    {
        cout<<prime[i]<<" ";//所有素数 
    }
    return 0;
} 

(RP++!)

【模板】素数筛法

标签:int   停止   筛选   div   namespace   怎么   表示   有一个   stream   

原文地址:https://www.cnblogs.com/Daz-Os0619/p/11469775.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!