OFFSET
1,2
LINKS
H. P. Robinson, Letter to N. J. A. Sloane, Oct 1981
MAPLE
c:=0: H[0]:=0: for n from 1 to 10^4 do H[n]:=1/n+H[n-1]: if floor(H[n])-floor(H[n-1])=1 then c:=1+c: b[c]:=n-1: else c:=c: fi: od: seq(b[j], j=1..c); # Emeric Deutsch
MATHEMATICA
a[n_] := Ceiling[k /. FindRoot[HarmonicNumber[k] == n, {k, Exp[n]}, WorkingPrecision -> 100]] - 1;
Array[a, 26] (* Jean-François Alcover, Apr 10 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Jan 23 2006
STATUS
approved