[Python]回文质数

题目描述

Oj Url：[USACO1.5]回文质数 Prime Palindromes
因为151既是一个质数又是一个回文数(从左到右和从右到左是看一样的)，所以 151 是回文质数。

写一个程序来找出范围[a,b](5 <= a < b <= 100,000,000)( 一亿)间的所有回文质数;

输入输出格式

输入格式：第 1 行: 二个整数 a 和 b .

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

输入输出样例

输入样例#1：

5 500

输出样例#1：

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

说明

Hint 1: Generate the palindromes and see if they are prime.

提示 1: 找出所有的回文数再判断它们是不是质数（素数）.

Hint 2: Generate palindromes by combining digits properly. You might need more than one of the loops like below.

提示 2: 要产生正确的回文数，你可能需要几个像下面这样的循环。

题目翻译来自NOCOW。

USACO Training Section 1.5

产生长度为5的回文数:

for (d1 = 1; d1 <= 9; d1+=2) {    // 只有奇数才会是素数
     for (d2 = 0; d2 <= 9; d2++) {
         for (d3 = 0; d3 <= 9; d3++) {
           palindrome = 10000*d1 + 1000*d2 +100*d3 + 10*d2 + d1;//(处理回文数...)
         }
     }
 }

程序代码

#打表法
def is_prime(radix):
    if radix == 1:
        return False
    if (radix == 2) or (radix ==3):
        return True
    if (radix % 6 != 1) and (radix % 6 != 5):
        return False
    for i in range(5, int(radix ** 0.5) + 1, 6):
        if (radix % i == 0) or (radix % (i + 2) == 0):
            return False
    return True
    
def is_palindrome(radix):
    for i in range(0, len(radix) // 2):
        if radix[i] != radix[len(radix) - i - 1]:
            return False
    return True

def Nums_of_Digits(radix):
    if (radix >= 10 and radix < 100 and radix != 11) or (radix >= 1000 and radix < 10000):
        return False
    if (radix >= 100000 and radix < 1000000) or (radix >= 10000000 and radix < 100000000):
        return False
    return True

a = 5
b = 100000000

with open("out.txt", "w") as fileobject: #将得到的结果输出至out.txt中
    for i in range(a, b + 1):
        if Nums_of_Digits(i):
            if is_palindrome(str(i)):
                if is_prime(i):
                    fileobject.write(str(i) + ',' + '\n')

form = [
5,
7,
11,
101,
131,
#此处省略约700个元素
9981899,
9989899
]

s = input().split()
a = int(s[0])
b = int(s[1])

for i in range(0, len(form)):
    if (form[i] >= a) and (form[i] <=b):
        print(form[i], end = '\n')