OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
MAPLE
filter:= proc(n) local t;
add(t[1]*t[2], t=ifactors(n)[2]) mod convert(convert(n, base, 10), `+`) = 0
end proc:
select(filter, [$1..1000]);
MATHEMATICA
Select[Range[650], Divisible[Plus @@ Times @@@ FactorInteger[#], Plus @@ IntegerDigits[#]] &] (* Amiram Eldar, Jun 21 2021 *)
PROG
(PARI) sopfr(n) = my(f=factor(n)); sum(k=1, #f~, f[k, 1]*f[k, 2]);
isok(k) = !(sopfr(k) % sumdigits(k)); \\ Michel Marcus, Jun 24 2021
CROSSREFS
KEYWORD
nonn,base
AUTHOR
J. M. Bergot and Robert Israel, Jun 21 2021
STATUS
approved