A108890 a(n) = floor(Product_{k=1..n} (Sum_{j=1..k} 1/j)).

1, 1, 2, 5, 13, 32, 83, 225, 638, 1871, 5651, 17537, 55772, 181348, 601754, 2034368, 6997316, 24456376, 86764857, 312157366, 1137925574, 4199870789, 15683541837, 59220398058, 225982562259, 871031643495, 3389581971377, 13311468152196

Offset

1,3

Links

Robert Israel, Table of n, a(n) for n = 1..1214

Maple

S:= ListTools:-PartialSums([seq(1/j, j=1..100)]):
P:= 1: R:= 1;
for k from 2 to 100 do
  P:= P * S[k];
  R:= R, floor(P);
od:
R; # Robert Israel, Nov 23 2022

Mathematica

f[n_] := Floor[ Product[ Sum[1/j, {j, k}], {k, n}]]; Table[ f[n], {n, 28}] (* Robert G. Wilson v, Jul 18 2005 *)

Prog

(PARI) for(n=1, 50, print1(floor(prod(k=1, n, sum(j=1, k, 1/j))), ", ")) \\ Lambert Klasen, Jul 18 2005

Crossrefs

Cf. A097423, A097424.

Floor(A097423/A097424).

Sequence in context: A298535 A110320 A219230 * A220739 A337282 A366117

Adjacent sequences: A108887 A108888 A108889 * A108891 A108892 A108893

Keyword

easy,nonn

Author

Leroy Quet, Jul 14 2005

Extensions

More terms from Lambert Klasen (lambert.klasen(AT)gmx.net), Hauke Worpel (hw1(AT)email.com), Ryan Propper and Robert G. Wilson v, Jul 18 2005

Status

approved