|
%I #4 Jun 24 2013 23:04:56
%S 3,3,3,2,3,3,3,2,3,3,3,3,2,3,3,3,2,3,3,3,3,2,3,3,3,2,3,3,3,3,2,3,3,3,
%T 3,2,3,3,3,2,3,3,3,3,2,3,3,3,2,3,3,3,3,2,3,3,3,2,3,3,3,3,2,3,3,3,3,2,
%U 3,3,3,2,3,3,3,3,2,3,3,3,2,3,3,3,3,2
%N Difference sequence for A226894.
%C The number of numbers log k between consecutive harmonic numbers is 1 or 2, so that the difference sequence for A226894 consists of 2's and 3's.
%H Clark Kimberling, <a href="/A226895/b226895.txt">Table of n, a(n) for n = 1..1000</a>
%e log 1 < log 2 < h(1) < log 3 < log 4 < h(2) < ..., so that the positions of h(1) and h(2) are 3 and 6, whence a(1) = 6 - 3 = 3.
%t z = 300; h[n_] := N[HarmonicNumber[n], 100]; t1 = Table[h[n], {n, 1, z}]; t2 = Table[N[Log[n], 100], {n, 1, 3 z}]; t3 = Union[t1, t2]; p[n_] := Position[t3, h[n]]; Flatten[Table[p[n], {n, 1, 3 z}]] (* A226894 *)
%t Differences[%] (* A226895 *)
%t Complement[Range[z], %%] (* A226896 *)
%Y Cf. A226894 , A226896 (complement).
%K nonn,easy
%O 1,1
%A _Clark Kimberling_, Jun 21 2013
|
