A014006 Pisot sequence E(10,18), a(n)=[ a(n-1)^2/a(n-2)+1/2 ].

10, 18, 32, 57, 102, 183, 328, 588, 1054, 1889, 3386, 6069, 10878, 19498, 34949, 62644, 112286, 201267, 360761, 646646, 1159081, 2077595, 3723986, 6675060, 11964714, 21446156, 38441170, 68903889, 123506801, 221379811, 396812324, 711266397, 1274909717

Offset

0,1

Links

Colin Barker, Table of n, a(n) for n = 0..1000

D. W. Boyd, Some integer sequences related to the Pisot sequences, Acta Arithmetica, 34 (1979), 295-305.

D. W. Boyd, Linear recurrence relations for some generalized Pisot sequences, Advances in Number Theory ( Kingston ON, 1991) 333-340, Oxford Sci. Publ., Oxford Univ. Press, New York, 1993.

Formula

Known not to satisfy any linear recurrence.

Prog

(PARI) pisotE(nmax, a1, a2) = {
  a=vector(nmax); a[1]=a1; a[2]=a2;
  for(n=3, nmax, a[n] = floor(a[n-1]^2/a[n-2]+1/2));
  a
}
pisotE(50, 10, 18) \\ Colin Barker, Jul 28 2016

Crossrefs

Sequence in context: A162828 A190038 A233695 * A090995 A363769 A153360

Adjacent sequences: A014003 A014004 A014005 * A014007 A014008 A014009

Keyword

nonn

Author

Simon Plouffe

Status

approved