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

数论_埃氏筛法(求区间内多少素数)

时间:2018-12-08 13:17:08      阅读:196      评论:0      收藏:0      [点我收藏+]

标签:typedef   return   const   stream   cstring   bsp   朋友   筛法   古希腊   

埃拉托斯特尼(公元前276—公元前194)


埃拉托斯特尼是古希腊著名的数学家、地理学家、天文学家。他先在亚历山大港学习,后又转至雅典。公元前236年,托勒密三世指定他为亚历山大图书馆的图书管理员和馆长。他跟阿基米德是好朋友。埃拉托斯特尼的主要贡献包括:

埃拉托斯特尼筛法:寻找素数的方法。

地理常数测量:日地间距的测量(现在称一个这样的距离为一个天文单位)、地月间距的测量、测量赤道与黄道之间的偏角、地球半径测量等。

精确地图绘制:当时只有托勒密等级的人物能绘出同等级的地图。

算法数学原理:
埃拉托斯特尼筛法是快速筛选素数的算法,在处理大量整数是否是素数时有较高的效率。

例:筛选小于n的整数并记录结果。

解:计算sqrt(n),∵对于任意的z<n,若z为合数,不妨设z=a*b,则必有min(a,b)<sqrt(n)

∴所有小于n的合数均可被小于sqrt(n)的整数整除

从2开始,依次去除小于n的整数中能被其整除的数,最后剩下的就是素数

int prime[maxn];
bool is_prime[maxn];//is_prime[i]是true表示i是素数
//返回n以内素数的个数
int sieve(int n)
{
    int p=0;
    for(int i=0;i<=n;i++) is_prime[i]=true;
    is_prime[0]=is_prime[1]=false;
    for(int i=2;i<=n;i++){
        if(is_prime[i]){
            prime[p++]=i;
            for(int j=2*i;j<=n;j+=i) is_prime[j]=false;
        }
    }
    return p;
}

  例题:洛谷 P1865  ,

codevs3223 素数密度

[a,b]的素数,只需要

[2,b√]
的素数表即可。这样我们可以筛出这些素数,然后用这些素数筛[a,b]。

注意,j=max(2ll,(l+i-1)/i)*i的意思:

(l+i-1)/i表示大于等于a的i的倍数的最小值。
2则是普遍的i的倍数的最小值,一定不能小于2。
再乘i,j就可以枚举[a,b]之间的合数了。

 

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int SIZE=1000010;
bool vis[SIZE];
bool pri[SIZE];
typedef long long LL;
int main()
{
    LL l,r;
    scanf("%lld%lld",&l,&r);

    for(LL i=2;i<=sqrt(r);i++)
    {
        if(!vis[i])
        {
            for(LL j=i*i;j<=sqrt(r);j+=i) vis[j]=1;

            for(LL j=max(2ll,(l+i-1)/i)*i;j<=r;j+=i) pri[j-l]=1;  
        }
    }
    int ans=0;
    for(int i=0;i<=r-l;i++) if(!pri[i]) ans++;
    printf("%d",ans);

    return 0;
}

  

数论_埃氏筛法(求区间内多少素数)

标签:typedef   return   const   stream   cstring   bsp   朋友   筛法   古希腊   

原文地址:https://www.cnblogs.com/passion-sky/p/10087034.html

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