OFFSET
1,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Harry J. Smith)
Vaclav Kotesovec, Plot of a(n)/n for n = 1..100000
FORMULA
a(n) = Sum_{i=1..n} floor(n/i^i). - Wesley Ivan Hurt, Sep 15 2017
G.f.: (1/(1 - x)) * Sum_{k>=1} x^(k^k)/(1 - x^(k^k)). - Seiichi Manyama, Aug 30 2021
Conjecture: a(n) ~ c * n, where c = A073009. - Vaclav Kotesovec, Aug 30 2021
EXAMPLE
a(7) = [7/1] + [7/4] + [7/27] + ... = 7 + 1 + 0 + 0 + ... = 8.
a(8) = [8/1] + [8/4] + [8/27] + [8/256] + ... = 8 + 2 + 0 + 0 + ... = 10.
MATHEMATICA
Flatten[{1, Table[Sum[Floor[n/k^k], {k, 1, Floor[N[Log[n]/LambertW[Log[n]]]] + 1}], {n, 2, 100}]}] (* Vaclav Kotesovec, Aug 30 2021 *)
PROG
(PARI) \p 10 v=[]; for(n=1, 120, v=concat(v, suminf(k=1, floor(n/k^k)))); v
(PARI) for (n=1, 1000, write("b062071.txt", n, " ", suminf(k=1, n\k^k)\1) ) \\ Harry J. Smith, Jul 31 2009
(PARI) a(n)=sum(k=1, exp(lambertw(log(n)))+1, n\k^k) \\ Charles R Greathouse IV, May 28 2015
(SageMath) [sum( floor(n/j^j) for j in (1..1+log(n)) ) for n in (1..100)] # G. C. Greubel, May 06 2022
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amarnath Murthy, Jun 13 2001
EXTENSIONS
More terms from Jason Earls, Jun 21 2001
STATUS
approved