素数回文(高效判断素数法)

摘要:
例如,151既是素数又是回文。输出对于每组数据,从小到大输出a和b之间的所有合格素数回文,每组数据后一个空行。主题给出的最大回文素数是998899。步骤3:输出值为print[]的下标为true:forifprintf;结果是2357111317192329#includeboolis[9998990];//如果i是素数,theis[i]=fasle,否则[i]=true;整数[1000];//Prime用于保存文档素数表voiset()//有效的素数方法:所有和都等于N个素数的乘积{int i,j;/*foris[i]=0;由于布尔类型的默认值为false,*/i=2可以被注释掉;对于//偶数,true不是素数Bootest//判断a是否为回文数{inttemp=a;intb=0;while(temp!

Problem Description

xiaoou33对既是素数又是回文的数特别感兴趣。比如说151既是素数又是个回文。现在xiaoou333想要你帮助他找出某个范围内的素数回文数,请你写个程序找出 a 跟b 之间满足条件的数。(5 <= a < b <= 100,000,000);
 
Input
这里有许多组数据,每组包括两组数据a跟b。
 
Output
对每一组数据,按从小到大输出a,b之间所有满足条件的素数回文数(包括a跟b)每组数据之后空一行。
 
Sample Input
5 500
 
Sample Output
5
7
11
101
131
151
181
191
313
353
373
383
 
Author
xiaoou333
 
Source
zjut
 
Recommend
linle
 
 
分析:除了11外,任意偶数长度的回文都不是素数因为都会被11整除。题目给的范围里最大的回文素数就是9989899。所以数组只要开到这个数就够了,减少了很多数的筛选与数组的空间。所有的偶数不是素数,素数的倍数不是素数。

 素数筛法是这样的:

    1.开一个大的bool型数组prime[],大小就是n+1就可以了.先把所有的下标为奇数的标为true,下标为偶数(除了2)的标为false.

    2.然后:将素数的倍数筛掉

      for( i=3; i<=sqrt(n); i+=2 )

      {   if(prime[i])

          for( j=i+i; j<=n; j+=i ) prime[j]=false;

      }

    3.最后输出bool数组中的值为true的单元的下标,就是所求的n以内的素数了。

    原理很简单,就是当i是质(素)数的时候,i的所有的倍数必然是合数。如果i已经被判断不是质数了,那么再找到i后面的质数来把这个质

数的倍数筛掉。

一个简单的筛素数的过程:n=30。

    1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30

   

    第 1 步过后4 ... 28 30这14个单元被标成false,其余为true。

    第 2 步开始:

     i=3;  由于prime[3]=true, 把prime[6], [9], [12], [15], [18], [21], [24], [27], [30]标为false.

     i=4;  由于prime[4]=false,不在继续筛法步骤。

     i=5;  由于prime[5]=true, 把prime[10],[15],[20],[25],[30]标为false.

     i=6>sqrt(30)算法结束。

    第 3 步把prime[]值为true的下标输出来:

     for(i=2; i<=30; i++)

     if(prime[i]) printf("%d ",i);

    结果是 2 3 5 7 11 13 17 19 23 29

#include<stdio.h>    
bool is[9989900]; //如果i是素数,then is[i]=fasle, else is[i]=true; 
int prime[1000]; //prime用来存回文素数表  

void set()//高效判断素数法:所有和数都等于N个素数的乘积
{
    int i,j;   
    /*for(i=5;i<=3163;i++)   
    is[i]=0;  
    由于bool类型默认值是false,所以可以注释掉*/
    i=2;
    for(j=i*i;j<=9989899;j+=i)   //偶数,true的不是素数
        is[j]=true; 
    for(i=3;i<=3163;i=i+2)   
    {   
        if(is[i])
            continue;   
        for(j=i*i;j<=9989899;j+=i)    //素数的倍数不是素数
            is[j]=true;   
    }   
}   

bool test(int a)//判断a是不是回文数  
{   
    int temp=a;   
    int b=0;   
    while(temp!=0)   //每位倒置,看结果与原数是否一致
    {   
        b=b*10;   
        b+=temp%10;   
        temp/=10;   
    }   
    return a==b;   
}   

int main()   
{   
    int a,b;   
    int i,k=0;   
    set();   
    for(i=5;i<=9989899;i+=2)   //将最大范围的所有满足条件的数求出来
        if(!is[i]&&test(i))   
            prime[k++]=i;   
        while(~scanf("%d %d",&a,&b))   //输出给定范文的结果
        {   
            for(i=0;i<k;i++)   
            {   
                if(prime[i]<a)   
                    continue;   
                else if(prime[i]<=b)   
                    printf("%d
",prime[i]);   
                else  
                    break;   
            }   
            printf("
");   
        }   
        return 0;   
}
方法二:
#include <stdio.h>
#include <string.h>
int huiwen(int k)  
{
    int a[10],i=0,j;
    while (k>0)
    {
        a[i]=k%10;
        k/=10;
        i++;
    }
    for (j=0; j<i; j++)
        if (a[j]!=a[i-j-1])
            return 0;
    return 1;
}
int hwlength(int k)  
{
    int a[10],i=0;
    while (k>0)
    {
        a[i]=k%10;
        k/=10;
        i++;
    }
    return (i);
}
int prime(int k) 
{
    int i;
    for (i=3; i*i<=k; i+=2)
        if (k%i==0)
            return 0;
    return 1;
}
int main()
{
    int a,b,i,j;
    scanf("%d%d",&a,&b);
    for (i=a; i<=b; i++)
    {
        if (i%2==0&&i!=2)  
            continue;
        if (i%5==0&&i!=5)  
            continue;
        if (hwlength(i)%2==0&&i!=11)  
            continue;           
        if (!huiwen(i))    
            continue;
        if (prime(i))
            printf("%d
",i);
    }
    return 0;
}

免责声明:文章转载自《素数回文(高效判断素数法)》仅用于学习参考。如对内容有疑问,请及时联系本站处理。

上篇Git忽略提交规则(转)OpenLayers3基础教程——OL3 介绍interaction下篇

宿迁高防,2C2G15M,22元/月;香港BGP,2C5G5M,25元/月 雨云优惠码:MjYwNzM=

相关文章

素数算法大全

注意: 如果没有特殊说明, 以下讨论的都是针对n为素数时的时间复杂度 1. 根据概念判断: 如果一个正整数只有两个因子, 1和p,则称p为素数. 代码: bool isPrime(int n) { if(n < 2) return false; for(int i = 2; i < n; ++i) if...

素数筛法知识点整理

素数的定义:除了1和它本身之外,不能被其他整数整除。 一、判定一个正整数n是否为素数的方法: ①定义法:枚举2~n-1这n-2个正整数,如果它们均不能整除n,则可断定n为素数。代码如下:时间复杂度为O(n),如果n为10^9,就不能用此方法。 1 bool is_prime(int n){ 2 if(n==1)return false; 3...

MD5加密和RSA加密

1.MD5加密    MD5(单向散列算法)的全称是Message-Digest Algorithm 5(信息-摘要算法),MD5算法的使用不需要支付任何版权费用。      MD5的功能:       ①.输入任意长度的信息,经过处理,输出为128位的信息(数字指纹);      ②.不同的输入得到的不同的结果(唯一性);      ③.根据128位的输...

204. Count Primes 素数的个数

Description: Count the number of prime numbers less than a non-negative number, n. 求0-n的素数个数,使用埃拉托斯特尼筛法 要得到自然数n以内的全部素数,必须把不大于  的所有素数的倍数剔除,剩下的就是素数。 给出要筛数值的范围n,找出以内的素数。先用2去筛,即把2留下,...

使用c语言和GMP库实现伪随机算法笔记

步骤一:安装GMP库,两种方法选其一既可 1.使用源码安装方式: 2.使用在线安装的方式: 步骤二:使用GMP库随机生成一个大数,样本代码如下: 步骤三:使用GMP库随机生成一个大数,并判断生成的大数是不是素数,样本代码如下: 步骤四:根据生成的大素数,产生下一个大素数,样本代码如下: 步骤五:根据随机产生的大素数和随机数,使用BlumBlumShus算...

中国石油大学(华东)计算机复试C语言参考题库

目录 复试c语言 【研究创新型】8.1 谁能出线 【设计型】8.2 统计素数的个数 【设计型】8.3 数组逆序输出 【设计型】8.4 在屏幕上显示杨辉三角形 【设计型】8.5 求最大值 【设计型】8.6 二维数组 【设计型】8.11 存储并输出一个矩阵 【设计型】8.7 给数组中的元素按顺序编号 【设计型】8.8 求各位数字组成的最大数 【设计型】8...