OFFSET
1,2
COMMENTS
For n=100, a(n)/n = 0.93; for n=1000, a(n)/n = 0.799; for n=10000, a(n)/n = 0.7326.
EXAMPLE
For n=1, sigma(n)/n = 1/1, so the set of distinct values so far is {1}, and the number of terms of this set is 1.
For n=2, sigma(n)/n = 3/2, so the set of distinct values so far is {1, 2, 3}, and the number of terms of this set is 3.
For n=3, sigma(n)/n = 4/3, so the set of distinct values so far is {1, 2, 3, 4}, and the number of terms of this set is 4.
PROG
(PARI) nbnd(n, vtot) = {ab = sigma(n)/n; if (! vecsearch(vtot, num=numerator(ab)), vtot = vecsort(concat(vtot, num))); if (! vecsearch(vtot, den=denominator(ab)), vtot = vecsort(concat(vtot, den))); return (vtot); }
lista(nn) = {vtot = []; for (n=1, nn, vtot = nbnd(n, vtot); print1(#vtot, ", "); ); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Nov 20 2014
STATUS
approved
