%I#16 2019年9月12日01:55:18
%S 3,4,5,6,11,335558785775611691991958633685129707921,
%电话:16958285961867275572768501667227661053104397934013,
%电话:43844958859448345817988338897185
%N回文p,这样pi(p)是回文素数。
%C来自Zak Seidov的建议,2005年2月2日。
%C a(16)>32*10^12_多诺万·约翰逊,2009年12月3日
%t NextPalindrome[n_]:=块[{l=Floor[Log[10,n]+1],idn=IntegerDigits[n]},如果[Union[idn]=={9},Return[2],如果[l<2,Return[n+1],如果[FromDigits[Reverse[Take[idn,Ceiling[l/2]]]]]FromDiges[Take[idn,-Ciling[1/2]],FromDigs[Join[Take[2,Ceiling[l/2]],Reverse[Take[id,Floor[l/2]2]]]],idfhn=FromDigits[Take[idn,Ceiling[l/2]]]+1;idp=FromDigits[Join[IntegerDigits[idfhn],Drop[Reverse[Integer Digits[idfhn]],Mod[l,2]]]]];
%t p=0;a={};Do[p=下一个回文[p];q=PrimePi[p];如果[PrimeQ[q],r=整数位数[q];如果[Reverse[r]==r,则打印[{p,q}];附录[a,p]]],{n,10^6}];一
%t palQ[n_]:=反向[x=整数位数[n]]==x;t={};Do[If[palQ[n]&&PrimeQ[x=PrimePi[n]]&&palQ[x],AppendTo[t,n]],{n,10^6}];t(*Jayanta Basu,2013年6月24日*)
%Y参见A046941、A046942、A103357、A103338、A103403。
%K nonn,基础,更多
%O 1,1号机组
%A _Robert G.Wilson v _,2005年2月3日
%E a(15),来自Donovan Johnson,2009年12月3日
%E a(16)-a(17),来自恰瓦乌,2019年9月4日
%E a(18)-a(19)摘自Giovanni Resta,2019年9月12日
|