OFFSET
1,1
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
EXAMPLE
digsum(10) = 1 + 0 = 1, PrimeDivisors(10) = PrimeDivisors(2 *5) = {2,5} and sopf(10) = 2 + 5 = 7 = 7*1.
digsum(154) = 1 + 5 + 4 = 10, PrimeDivisors(154) = PrimeDivisors(2 * 7 * 11) = {2,7,11} and sopf(154) = 2+7+11 = 20 = 2*10.
MATHEMATICA
Select[Range[2, 800], Divisible[Total[Select[Divisors[#], PrimeQ]], Total[ IntegerDigits[#]]]&] (* Harvey P. Dale, Sep 23 2012 *)
PROG
(Magma) [m:m in [2..730]| &+PrimeDivisors(m) mod &+Intseq(m) eq 0]; // Marius A. Burtea, Jul 11 2019
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Floor van Lamoen, Sep 23 2002 and Sep 30 2002
STATUS
approved