OFFSET
1,2
COMMENTS
For most values (except 9,2025 and 99225) number of divisors of n = sum of digits of n, see A057531.
The above comment is wrong: for 16 out of the first 34 terms of the sequence, the number of divisors of n does not equal the sum of the digits of n. - Harvey P. Dale, Dec 31 2015
Since A000005(n) is odd, n must be a square. - Robert Israel, Oct 31 2019
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= proc(m) local n;
n:= m^2;
n mod numtheory:-tau(n) = 0 and n mod convert(convert(n, base, 10), `+`) = 0
end proc:
map(`^`, select(filter, [seq(i, i=1..10000, 2)]), 2); # Robert Israel, Oct 31 2019
MATHEMATICA
Select[Range[1, 5*10^6, 2], Divisible[#, DivisorSigma[0, #]] && Divisible[ #, Total[ IntegerDigits[#]]]&] (* Harvey P. Dale, Dec 31 2015 *)
PROG
(Magma) [k:k in [1..6000001 by 2]| IsIntegral(k/NumberOfDivisors(k)) and IsIntegral(k/&+Intseq(k))]; // Marius A. Burtea, Oct 31 2019
CROSSREFS
KEYWORD
nonn,easy,base
AUTHOR
Asher Auel, Sep 03 2000
EXTENSIONS
More terms from Harvey P. Dale, Dec 31 2015
STATUS
approved
