%I#15 2023年1月7日12:05:23
%S 1,2,3,4,5,6,7,8,9,11111 21221661211111113113112123111,
%电话:4224441614411111331131131131171113113131131131121112,
%电话:2113221312231121234243111331311311311313131142432426241111211112211116
%N对k进行编号,使其数字的乘积除以k和R(k),其中R(k)是k的反面数字。
%C A007602的子序列。
%H David A.Corneth,n的表,n的A(n)=1.1228(Paola P.Lava的前200个术语,术语<10^13)
%p R:=进程(w)局部x,y,z;x: =w;y: =0;对于z从1到ilog10(x)+1做y:=10*y+(xmod10);x: =trunc(x/10);od;年;结束时间:
%p p:=proc(q)局部a,n;对于从1到q的n,执行a:=转换(convert(n,base,10),`*`);
%p如果a>0,则如果R(n)mod a=0且n mod a=0,则打印(n);fi;fi;od;结束:P(10^12);
%t pddQ[n_]:=模块[{pd=Times@@IntegerDigits[n]},pd=0&&Mod[n,pd]==Mod[IntegerReverse[n],pd]==0];选择[范围[22*10^4],pddQ](*_哈维·P·戴尔,2023年1月7日*)
%o(PARI)是(n)={my(d=数字(n),vp=vecprod(d));如果(vp!=0&&n%vp==0&fromdigits(Vecrev(d),%vp==0,返回(1));0}\\_David A.Corneth_,2021年3月30日
%Y参考A004086、A007602。
%K nonn,基础,简单
%O 1,2号机组
%A_Paolo P.Lava,2016年11月2日
|