A020741 Pisot sequence T(6,10), a(n) = floor(a(n-1)^2/a(n-2)).

6, 10, 16, 25, 39, 60, 92, 141, 216, 330, 504, 769, 1173, 1789, 2728, 4159, 6340, 9664, 14730, 22451, 34219, 52155, 79492, 121157, 184660, 281447, 428963, 653797, 996474, 1518759, 2314790, 3528046, 5377208, 8195575, 12491138, 19038142, 29016639, 44225184

Offset

0,1

Links

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

Mathematica

nxt[{a_, b_}]:={b, Floor[b^2/a]}; NestList[nxt, {6, 10}, 40][[All, 1]] (* Harvey P. Dale, Dec 27 2021 *)

Prog

(PARI) pisotT(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]));
  a
}
pisotT(50, 6, 10) \\ Colin Barker, Jul 29 2016

Crossrefs

See A008776 for definitions of Pisot sequences.

Sequence in context: A230276 A079329 A302748 * A274287 A090991 A019533

Adjacent sequences: A020738 A020739 A020740 * A020742 A020743 A020744

Keyword

nonn

Author

David W. Wilson

Status

approved