login
A067054
a(n) = floor(s^s) where s = Sum_{k=1..n} 1/k.
2
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
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 * A069357 A039893 A369112
KEYWORD
easy,nonn
AUTHOR
Amarnath Murthy, Jan 02 2002
EXTENSIONS
Corrected and extended by Robert G. Wilson v, Jan 04 2002
STATUS
approved