OFFSET
1,2
COMMENTS
The sequence also includes: 83860580242, 4578632504347, 4578632504465, 4578632504515. - Daniel Suteu, Nov 24 2020
LINKS
Eric Weisstein's World of Mathematics, Odd Divisor Function.
EXAMPLE
165 is in the sequence because the average number of odd divisors of {1..165} is an integer: A060831(165) / 165 = 495 / 165 = 3.
MATHEMATICA
s[n_] := Module[{c = 0, k = 1, sum = 0, seq = {}}, While[c < n, sum += DivisorSigma[0, k/2^IntegerExponent[k, 2]]; If[Divisible[sum, k], c++; AppendTo[seq, k]]; k++]; seq]; s[13] (* Amiram Eldar, Nov 18 2020 *)
PROG
(PARI) f(n) = my(n2=n\2); sum(k=1, sqrtint(n), n\k)*2-sqrtint(n)^2-sum(k=1, sqrtint(n2), n2\k)*2+sqrtint(n2)^2; \\ A060831
isok(k) = (f(k) % k) == 0; \\ Michel Marcus, Nov 25 2020
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Ilya Gutkovskiy, Nov 18 2020
STATUS
approved