%I#42 2023年2月28日06:03:09
%编号:1,2,5,12,163623411961872570856714987392762990591
%对k进行编号,使k*prime(k)是回文。
%C a(12)>3.7*10^12_Giovanni Resta_,2013年6月28日
%e 4119是序列中的,因为4119素数是39119,4199×39119=161131161是回文。
%p ispal:=proc(n)局部L;
%p L:=换算(n,基数,10);
%p L=列表工具:-反向(L);
%p端程序:
%p R:=NULL:计数:=0:p:=1:
%当计数<11 do时,k从1开始为p
%p:=下一素数(p);
%p如果ispal(k*p),则R:=R,k;计数:=计数+1 fi
%日期:
%p R;编号_罗伯特·伊斯雷尔,2023年2月22日
%t palQ[n_]:=起始数字[Reverse[IntegerDigits[n]]==n;t={};Do[If[palQ[Prime[n]*n],AppendTo[t,n]],{n,15*10^5}];t(*Jayanta Basu,2013年5月11日*)
%o(PARI)ispal(n)=我的(d=数字(n));d==Vecrev(d);\\A002113号
%o isok(k)=ispal(k*prime(k))\\阿莱克桑德鲁·彼得雷斯库,2023年2月22日
%o(Python)
%o来自sympy进口筛
%o def ok(n):返回n和(s:=str(n*筛[n]))==s[::-1]
%o打印([k表示范围内的k(10**6),如果ok(k)])#_Michael S.Branicky_,2023年2月22日
%Y参考A002113、A033286、A084121、A084123。
%K基,nonn,更多
%O 1,2号机组
%A Giovanni Resta,2003年5月14日
|