P1217 [USACO1.5] 回文质数 Prime Palindromes

[USACO1.5] 回文质数 Prime Palindromes

题目描述

因为 $151$ 既是一个质数又是一个回文数（从左到右和从右到左是看一样的），所以 $151$ 是回文质数。

写一个程序来找出范围 $[a,b](5\le a<b\le 100,000,000)$（一亿）间的所有回文质数。

输入格式

第一行输入两个正整数 $a$ 和 $b$。

输出格式

输出一个回文质数的列表，一行一个。

样例 #1

样例输入 #1

5 500

样例输出 #1

5
7
11
101
131
151
181
191
313
353
373
383

第一次，使用穷举法找质数：

#include <stdio.h>
int main(){
    int a,b,a1,flag,m;
    scanf("%d%d",&a,&b);
    for (int i = a; i <= b; ++i) {//遍历a到b的所有数
        flag = 0;//把所有数标记为质数
        a1 = i;//因为后面会涉及到i的运算，为了避免造成死循环，将i单独保存
        for (int j = 2; j <= i; ++j) {
            if (i % j == 0 && i!=j){//使用i，对2到i的数取余，若有第三个约数，标记为非质数
                flag = 1;
                break;
            }
        }
        if (flag == 0){//当i为质数
            m=0;//因为每次循环时，回文不同，所以将回文数设置为0
            while(flag == 0){//进入死循环
                //计算回文，例a1=123 ，m = (0*10)+(123%10)=3,下一次循环就是a1=123,m = (3*10)+(12%10)=32
                m = (m*10)+(a1%10);
                a1 /= 10;
                if (a1 == 0){//当除到最高位后，跳出循环
                    break;
                }
            }
            if (m == i){//如果是回文，则输出
                printf("%d\n",i);
            }
        }
    }

    return 0;
}

穷举法过于低效，超时了。

第二次，使用“试除法”找质数：

	试除法的优化在于，任何合数（非质数）都可以表示为两个较小数的乘积。而这两个较小数中，必然有一个小于或等于该数的平方根。所以我们只需要在2到sqrt(num)范围内寻找可能的因子即可。

#include <stdio.h>
#include <math.h>
int main(){
    int a,b,a1,flag,m;
    scanf("%d%d",&a,&b);
    for (int i = a; i <= b; ++i) {//遍历a到b的所有数
        flag = 0;//把所有数标记为质数
        a1 = i;//因为后面会涉及到i的运算，为了避免造成死循环，将i单独保存
        for (int j = 2; j <= sqrt(i); ++j) {
            if (i % j == 0 && i!=j){//使用i，对2到sqrt(i)的数取余，若有第三个约数，标记为非质数
                flag = 1;
                break;
            }
        }
        if (flag == 0){//当i为质数
            m=0;//因为每次循环时，回文不同，所以将回文数设置为0
            while(flag == 0){//进入死循环
                //计算回文，例a1=123 ，m = (0*10)+(123%10)=3,下一次循环就是a1=123,m = (3*10)+(12%10)=32
                m = (m*10)+(a1%10);
                a1 /= 10;
                if (a1 == 0){//当除到最高位后，跳出循环
                    break;
                }
            }
            if (m == i){//如果是回文，则输出
                printf("%d\n",i);
            }
        }
    }

    return 0;
}

可以看到时间效率有所提高，但是还是不够。

第三次，使用“试除法”找质数，并剔除偶数位数字，因为偶数位是不会出现回文质数的：

#include <stdio.h>
#include <math.h>

bool is_prime(int num) {
    if (num <= 1) {
        return false; // 0和1不是质数
    }

    if (num == 2) {
        return true; // 2是质数
    }

    if (num % 2 == 0) {
        return false; // 排除所有偶数（除了2）
    }


    int sqrt_num = (int) sqrt(num);
    for (int i = 3; i <= sqrt_num; i += 2) {//只要奇数
        if (num % i == 0) {
            return false; // 如果能整除，说明不是质数
        }
    }

    return true; // 到这里说明是质数
}

int reverse_number(int num) {
    int reversed = 0;
    while (num > 0) {
        reversed = (reversed * 10) + (num % 10);//取反序
        num /= 10;
    }
    return reversed;
}

int main() {
    int a, b;
    scanf("%d%d", &a, &b);

    for (int i = a; i <= b; ++i) {
        if ((i > 999 && i < 10000) || (i > 99999 && i < 1000000) ||
            (i > 9999999 && i < 100000000)) {//排除偶数位数字1000-9999，100000-999999
            continue;
        }
zhong
        int reversed = reverse_number(i);
        if (reversed == i) {//如果是回文
            if (is_prime(i)) {//如果是质数
                printf("%d\n", i);
            }
        }

    }

    return 0;
}

终于过了，这里要注意的是要先判断是否为回文，不然还是会超时。