OFFSET
1,3
COMMENTS
Terms of A061205, sorted in increasing order, with duplicates removed.
LINKS
Jon E. Schoenfield, Table of n, a(n) for n = 1..10000 (first 1000 terms from Alois P. Heinz)
EXAMPLE
12*21 = 252, so 252 is a term.
156*651 = 101556, so 101556 is a term. (It can also be written as 273*372; see A203924.)
MAPLE
a:= proc(n) option remember; local k, d; for k from 1+a(n-1) do
for d in numtheory[divisors](k) do if k = d*(s-> parse(cat(
seq(s[-i], i=1..length(s)))))(""||d) then return k fi od od
end: a(1):=0:
seq(a(n), n=1..60); # Alois P. Heinz, May 27 2018
MATHEMATICA
a={0}; h=-1; For[k=0, k<=2701, k++, For[m=1, m<=DivisorSigma[0, k], m++, d=Divisors[k]; If[k/Part[d, m] == FromDigits[Reverse[IntegerDigits[Part[d, m]]]] && k>h , AppendTo[a, k]; h=k]]]; a (* Stefano Spezia, Jan 28 2023 *)
PROG
(PARI) isok(n) = if (n==0, return (1), fordiv(n, d, if (n/d == fromdigits(Vecrev(digits(d))), return (1))); return (0)); \\ Michel Marcus, May 28 2018
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Jon E. Schoenfield, May 27 2018
STATUS
approved