login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A116601 a(0) = a(1) = 0; for n >= 2, a(n) = floor(sqrt(2^(n-2)-1)). 1
0, 0, 0, 1, 1, 2, 3, 5, 7, 11, 15, 22, 31, 45, 63, 90, 127, 181, 255, 362, 511, 724, 1023, 1448, 2047, 2896, 4095, 5792, 8191, 11585, 16383, 23170, 32767, 46340, 65535, 92681, 131071, 185363, 262143, 370727, 524287, 741455, 1048575, 1482910, 2097151, 2965820, 4194303, 5931641 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,6

COMMENTS

Produced by a Markov chain and designed to have a Sqrt[2] limiting ratio.

LINKS

G. C. Greubel, Table of n, a(n) for n = 0..440

FORMULA

M = {{3, 1}, {1, 3}}/N[Sqrt[8]]: Det[M]=1; w[0] = {0, 1}; w[n_] := w[n] = M.w[n - 1] a(n) = Floor[w[n][[1]]]

MATHEMATICA

M = {{3, 1}, {1, 3}}/N[Sqrt[8]] w[0] = {0, 1}; w[n_] := w[n] = M.w[n - 1] a0 = Table[Floor[w[n][[1]]], {n, 0, 25}] b0 = Table[N[a0[[n + 1]]/a0[[n]]], {n, 4, 24}]

Join[{0, 0}, Table[Floor[Sqrt[2^(n - 2) - 1]], {n, 2, 50}]] (* G. C. Greubel, Oct 28 2017 *)

CROSSREFS

Cf. A000225.

Sequence in context: A101977 A024793 A326590 * A218024 A236102 A241726

Adjacent sequences:  A116598 A116599 A116600 * A116602 A116603 A116604

KEYWORD

nonn

AUTHOR

Roger L. Bagula, Mar 28 2006

EXTENSIONS

Edited by N. J. A. Sloane, May 10 2007

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 24 23:45 EDT 2020. Contains 334581 sequences. (Running on oeis4.)