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.
%I A091656 #6 May 13 2013 01:54:07 %S A091656 1,2,5,9,9,13,26,63,68,68,68,87,121,121,165,207,207,221,221,287,289, %T A091656 325,428,440,483,544,544,544,544,544,558,558,558,966,1035,1035,1146, %U A091656 1146,1332,1332,1332,1665,1665,1665,1665,1665,1727,1727,2052,2157,2331,2331 %N A091656 Least number k such that the continued fraction expansion of H(k) contains the numbers 1, 2, ..., n, where H(k) is the k-th Harmonic number. %H A091656 Charles R Greathouse IV, <a href="/A091656/b091656.txt">Table of n, a(n) for n = 1..250</a> %e A091656 a(6) = 13 because CF( H(13)) = 3 + [5, 1, 1, 4, 2, 1, 3, 2, 1, 3, 1, 4, 1, 6], the first six integers are present. %t A091656 f[n_] := Block[{k = 1}, While[ StringPosition[ ToString[ Union[ ContinuedFraction[ Sum[1/i, {i, 1, k}]]]], StringDrop[ ToString[ Table[i, {i, n}]], -1]] == {}, k++ ]; k]; Table[ f[n], {n, 1, 52}] %o A091656 (PARI) list(lim)=my(v=vector(lim\1),n,t,H,i=1);while(1,H+=1/n++;t=vecsort(contfrac(H),,8);if(#t>=i&&t[i]==i,v[i]=n;print1(n":"i", ");if(i++>#v,return(v));H-=1/n;n--)) \\ _Charles R Greathouse IV_, Jan 25 2012 %Y A091656 Cf. A055573, A091655. %K A091656 nonn %O A091656 1,2 %A A091656 _Robert G. Wilson v_, Jan 26 2004