码迷,mamicode.com
首页 > 编程语言 > 详细

【算法杂谈】埃氏素数筛

时间:2016-12-02 21:38:29      阅读:164      评论:0      收藏:0      [点我收藏+]

标签:应用   com   技术分享   turn   space   素数   cout   来讲   方法   

【今天我们来讲讲筛子】

【埃氏筛的基本思想】

简单来说就是把不大于技术分享(n为数据范围)以内的素数的倍数全都去掉,那么剩下的就是2~n之间的素数了。

【举个例子】

我们假设现在n是25。

第一步:先把2作为筛子,那么所有2的倍数都被筛掉了。

则当前序列为:2 3 5 7 9 11 13 15 17 19 21 23 25

第二步:剩下的序列中第一个素数是3,将序列中3的倍数划掉。

则当前序列为:2 3 5 7 11 13 17 19 23 25

25仍然大于3的平方,所以我们还要继续筛

现在序列中第一个素数是5,同样将序列中5的倍数划掉。

则当前序列为:2 3 5 7 11 13 17 19 23

因为23小于5的平方,跳出循环

So 最后的答案就是:2 3 5 7 11 13 17 19 23 啦

是不是很简单?

哦,忘说了,这种方法的时间复杂度是O(n log n)。

想要更快的吗,还有线性筛呢!!!

【参考代码】

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<vector>
#include<algorithm>
using namespace std;
const long long maxn=10000007+10;
const long long maxp=700000;
int vis[maxn];
long long prime[maxp];
long long n;

long long gen()
{
    long long m=(long long)sqrt(n+0.5);
    memset(vis,0,sizeof(vis));
    vis[2]=0;
    for(long long i=3;i<=m;i=i+2)
	{
        if(!vis[i])for(long long j=i*i;j<=n;j+=i) vis[j]=1;
        if(i*i>n)break;
    }
    long long c=1;
    prime[0]=2;
    for(long long i=3;i<=n;i=i+2) if(!vis[i]) prime[c++]=i;
    return c;
}
int main()
{
    long long c;
    cin>>n;
    c=gen();
    for(long long i=0;i<c;i++) printf("%lld ",prime[i]);
    cout<<endl<<c;
    return 0;
}

 【实际应用】

找素数啊,这还用说吗???

【算法杂谈】埃氏素数筛

标签:应用   com   技术分享   turn   space   素数   cout   来讲   方法   

原文地址:http://www.cnblogs.com/lijiaxin-blog-cpp/p/6126827.html

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