cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A076661 Index of first term of the harmonic sequence having the same denominator as the partial harmonic sequence beginning with 1/n.

This page as a plain text file.
%I A076661 #4 Mar 31 2012 21:03:40
%S A076661 1,2,4,9,9,10,10,14,25,27,27,27,27,27,27,27,49,49,49,49,49,49,49,49,
%T A076661 49,50,50,125,125,125,125,125,125,125,143,143,143,143,143,136,136,136,
%U A076661 136,136,136,136,136,136,136,136,98,98,98,133,133,133,133,125,125,125,125
%N A076661 Index of first term of the harmonic sequence having the same denominator as the partial harmonic sequence beginning with 1/n.
%C A076661 Of more interest is the index of terms after which the denominators of the harmonic sequence always match the denominators of the partial harmonic sequence. Notice that 1/4+..1/21 has denominator 15519504, but 1/1+1/2+..1/21 has denominator 5173168.
%e A076661 a(4) = firstHM[4] = 9 because 1/4+1/5+1/6+1/7+1/8+1/9 has the same denominator (2520) as 1/1+1/2+..+1/8+1/9 (and the sums to 4,5,6,7 and 8 do not).
%t A076661 harmNumber[m_, n_] := HarmonicNumber[n] - HarmonicNumber[m - 1]; denH[n_] := Denominator[HarmonicNumber[n]]; denH[m_, n_] := Denominator[harmNumber[m, n]]; firstHM[m_] := Do[If[denH[k] == denH[m, k], Return[k], ], {k, m, 10^4}]
%Y A076661 Cf. A002805.
%K A076661 nonn
%O A076661 1,2
%A A076661 _Hollie L. Buchanan II_, Oct 24 2002