OFFSET
1,2
COMMENTS
Numbers that are divisible by the sum of their digits and whose trajectory under iteration of sum of squares of digits map includes 1.
LINKS
EXAMPLE
133 is a term because 133/7 = 19 and its trajectory under iteration of sum of squares of digits map is 133 -> 19 -> 82 -> 68 -> 100 -> 1.
MAPLE
q:= proc(n) local m, s; m, s:= n, {};
if irem(n, add(i, i=convert(n, base, 10)))>0 then return false fi;
do if m=1 then return true
elif m in s then return false
else s, m:= s union {m}, add(i^2, i=convert(m, base, 10))
fi
od
end:
select(q, [$1..2000])[];
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Alois P. Heinz, Jan 03 2022
STATUS
approved