OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
sigma(98) - 98 - 1 = 171 - 98 - 1 = 72 and 8*9 = 72 so 98 is in the sequence.
MAPLE
with(numtheory):
for n from 1 to 3000 do:
q:=convert(n, base, 10):n0:=nops(q):
pr:=product('q[i]', 'i'=1..n0):p:=sigma(n)-n-1:
if p=pr
then
printf(`%d, `, n):
else
fi:
od:
MATHEMATICA
Do[If[DivisorSigma[1, n]-n-1==Apply[Times, IntegerDigits[n]], Print[n]], {n, 2000}]
Select[Range[2, 2000], Total[Most[Rest[Divisors[#]]]]==Times@@ IntegerDigits[ #]&] (* Harvey P. Dale, Jul 20 2019 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Michel Lagneau, May 23 2016
STATUS
approved