A138702 - OEIS

The OEIS is supported by the many generous donors to the OEIS Foundation.

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Other ways to Give

A138702

a(n) = number of terms in the continued fraction of the absolute value of B_n, the n-th Bernoulli number.

1, 2, 2, 1, 2, 1, 2, 1, 2, 1, 3, 1, 6, 1, 2, 1, 7, 1, 7, 1, 4, 1, 4, 1, 6, 1, 2, 1, 6, 1, 7, 1, 7, 1, 2, 1, 10, 1, 2, 1, 8, 1, 2, 1, 3, 1, 5, 1, 10, 1, 3, 1, 7, 1, 7, 1, 6, 1, 6, 1, 17, 1, 2, 1, 7, 1, 10, 1, 2, 1, 7, 1, 23, 1, 2, 1, 2, 1, 5, 1, 18, 1, 5, 1, 16, 1, 2, 1, 10, 1, 14, 1, 6, 1, 2, 1, 18, 1, 2, 1

(list; graph; refs; listen; history; text; internal format)

OFFSET

0,2

COMMENTS

The continued fraction terms being counted include the initial 0, if there is one. (a(n), for all odd n >= 3, is 1.)

LINKS

Antti Karttunen, Table of n, a(n) for n = 0..16383

EXAMPLE

The 12th Bernoulli number is -691/2730. Now 691/2730 has the continued fraction 0 + 1/(3 + 1/(1 + 1/(19 + 1/(3 + 1/11)))), which has 6 terms (including the zero). So a(12) = 6.

PROG

(PARI)

lcf(x)=local(r); r=1; while(1, x-=floor(x); if(x==0, return(r)); x=1/x; r++)

a(n)=lcf(abs(bernfrac(n))) \\ Franklin T. Adams-Watters, May 14 2010

CROSSREFS

Cf. A138701, A138703.

Sequence in context: A052005 A377902 A346707 * A344339 A348364 A279620

Adjacent sequences: A138699 A138700 A138701 * A138703 A138704 A138705

KEYWORD

nonn

AUTHOR

Leroy Quet, Mar 26 2008

EXTENSIONS

More terms from Franklin T. Adams-Watters, May 14 2010

STATUS

approved