A067054 a(n) = floor(s^s) where s = Sum_{k=1..n} 1/k.

1, 1, 3, 4, 6, 8, 11, 15, 18, 23, 28, 33, 39, 46, 53, 61, 70, 79, 89, 100, 111, 123, 137, 150, 165, 181, 197, 215, 233, 252, 273, 294, 316, 340, 364, 389, 416, 443, 472, 502, 533, 565, 598, 633, 669, 706, 745, 784, 825, 868, 911, 956, 1003, 1051, 1100, 1151

Offset

1,3

Links

Harry J. Smith, Table of n, a(n) for n = 1..1000

Example

a(6) = floor(1 + 1/2 + 1/3 + 1/4 + 1/5 + 1/6)^(1 + 1/2 + 1/3 + 1/4 + 1/5 + 1/6)) = floor(2.45^2.45) = floor(8.983734461...) = 8.

Maple

f:= proc(n) local v;
v:= floor((Psi(n+1)+gamma)^(Psi(n+1)+gamma)):
while not v::integer do
    Digits:= 2*Digits;
    v:= floor((Psi(n+1)+gamma)^(Psi(n+1)+gamma)):
  od;
  v
end proc:
seq(f(n), n=1..1000); # Robert Israel, Dec 21 2015

Mathematica

Table[ Floor[ Sum[ (1/i), {i, 1, n} ]^Sum[ (1/i), {i, 1, n} ]], {n, 1, 60} ]

Prog

(PARI) { s=0; for (n=1, 1000, s+=1/n; write("b067054.txt", n, " ", floor(s^s)) ) } \\ Harry J. Smith, May 01 2010

Crossrefs

Sequence in context: A365271 A242254 A107770 * A039893 A050115 A167711

Adjacent sequences: A067051 A067052 A067053 * A067055 A067056 A067057

Keyword

easy,nonn

Author

Amarnath Murthy, Jan 02 2002

Extensions

Corrected and extended by Robert G. Wilson v, Jan 04 2002

Status

approved