OFFSET
1,2
COMMENTS
Subsequence of A007602.
LINKS
David A. Corneth, Table of n, a(n) for n = 1..1228 (first 200 terms from Paola P. Lava, terms < 10^13)
MAPLE
R:=proc(w) local x, y, z; x:=w; y:=0; for z from 1 to ilog10(x)+1 do y:=10*y+(x mod 10); x:=trunc(x/10); od; y; end:
P:=proc(q) local a, n; for n from 1 to q do a:=convert(convert(n, base, 10), `*`);
if a>0 then if R(n) mod a=0 and n mod a=0 then print(n); fi; fi; od; end: P(10^12);
MATHEMATICA
pddQ[n_]:=Module[{pd=Times@@IntegerDigits[n]}, pd!=0&&Mod[n, pd] == Mod[ IntegerReverse[n], pd]==0]; Select[Range[22*10^4], pddQ] (* Harvey P. Dale, Jan 07 2023 *)
PROG
(PARI) is(n) = { my(d = digits(n), vp = vecprod(d)); if(vp != 0 && n%vp == 0 && fromdigits(Vecrev(d))%vp == 0, return(1) ); 0 } \\ David A. Corneth, Mar 30 2021
CROSSREFS
KEYWORD
nonn,base,easy
AUTHOR
Paolo P. Lava, Nov 02 2016
STATUS
approved