%我#34 2021年4月7日00:23:45

%S 1,2,3,4,5,6,7,8,9,10,11,12,14,15,16,18,20,21,22,24,25,27,28,30,32,33，

%电话：35,36,40,42,44,45,48,49,50,54,55,56,60,63,64,66,70,72,75,77,80,81,84，

%U 88,90,96,98,99100101105110112120121125126128131单位

%所有素因子都是回文的数字。

%C A002385的乘法闭包；A051038和A046368是子序列_Reinhard Zumkeller，2011年4月11日

%H Ivan Neretin，n表，n=1..10550的a（n）</a>

%H<a href=“/index/Pri#prime_factors”>与质数因子相关的索引条目序列</a>

%F总和{n>=1}1/a（n）=A002385}中的产品{p p/（p-1）=5.0949…-阿米拉姆·埃尔达尔，2020年9月27日

%e10=2*5是一个术语，因为2和5都是回文。

%e 110=2×5×11是一个术语，因为2、5和11是回文。

%p N:=5:#以获取最多N位数的所有项

%p digrev:=程序（t）局部L；L： =转换（t，基数，10）；

%p加（L[-i-1]*10^i，i=0..nops（L）-1）；

%p端程序：

%p p素数：=[2,3,5,7,11]：

%p代表d从3到N乘以2 do

%p m：=（d-1）/2；

%p p素数：=p素数，选择（i素数，[seq（seq（n*10^（m+1）+y*10^m+digrev（n），y=0..9），n=10^；

%日期：

%p p素数：=映射（op，[PP素数]）：

%p M：=10^N：

%p B:=矢量（M）；

%p B[1]：=1：

%p中的p代表p素数do

%p代表k从1到地板（log[p]（M））do

%p R：=[1美元地板（M/p^k）]；

%p B[p^k*R]：=B[p*k*R]+B[R]

%日期

%日期：

%p选择（t->B[t]>0，[1..M]美元）；#_罗伯特·伊斯雷尔，2015年7月5日

%p#备选方案

%p为A033620:=进程（n）

%numtheory[因子集]（n）do中d的p

%p如果不是A002113（op（1，d）），则

%p返回false；

%p end if；

%p端do；

%p为真；

%p端程序：

%p表示n从1到300 do

%p如果是A033620（n），则

%p打印f（“%d，”，n）；

%p end if；

%截止日期：#R.J.Mathar_，2015年9月9日

%t palQ[n_]：=反向[x=整数位数[n]]==x；选择[Range[131]，和@@palQ/@First/@FactorInteger[#]&]（*_Jayanta Basu_，2013年6月5日*）

%o（哈斯克尔）

%o a033620 n=a033620_列表！！（n-1）

%o a033620_list=过滤器chi[1..]，其中

%o chin=a136522 spf==1&&（n’==1|| chin’）其中

%o n’=n `div`spf

%o spf=a020639 n——参考a020639

%o--_Reinhard Zumkeller，2011年4月11日

%o（PARI）ispal（n）=n=数字（n）；对于（i=1，#n\2，如果（n[i]！=n[#n+1-i]，返回（0））；1

%o是（n）=如果（n<13，n>0，vecmin（应用（ispal，系数（n）[，1]））

%o（Python）

%o来自症状输入isprime，素数

%o定义pal（n）：s=str（n）；返回s==s[：：-1]

%o def ok（n）：为素数（n）中的f返回全部（pal（f））

%o打印（列表（过滤器（ok，范围（1132）））#_Michael S.Branicky_，2021年4月6日

%Y参考A002113、A002385、A046368、A051038。

%K nonn，基础，简单

%O 1、2

%A _N.J.A.Sloane，1998年5月17日