A048579 Pisot sequence L(3,8).

3, 8, 22, 61, 170, 474, 1322, 3688, 10289, 28705, 80084, 223427, 623341, 1739065, 4851835, 13536184, 37764738, 105360229, 293945582, 820081790, 2287954586, 6383187935, 17808521403, 49684176276, 138614392311, 386721712946, 1078918868168, 3010086801751

Offset

0,1

Links

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

Mathematica

RecurrenceTable[{a[0] == 3, a[1] == 8, a[n] == Ceiling[a[n - 1]^2/a[n - 2]]}, a, {n, 0, 30}] (* Bruno Berselli, Feb 05 2016 *)

Prog

(Magma) Lxy:=[3, 8]; [n le 2 select Lxy[n] else Ceiling(Self(n-1)^2/Self(n-2)): n in [1..30]]; // Bruno Berselli, Feb 05 2016
(PARI) pisotL(nmax, a1, a2) = {
  a=vector(nmax); a[1]=a1; a[2]=a2;
  for(n=3, nmax, a[n] = ceil(a[n-1]^2/a[n-2]));
  a
}
pisotL(50, 3, 8) \\ Colin Barker, Aug 07 2016

Crossrefs

See A008776 for definitions of Pisot sequences.

Sequence in context: A318820 A200752 A279378 * A121449 A025566 A027036

Adjacent sequences: A048576 A048577 A048578 * A048580 A048581 A048582

Keyword

nonn

Author

David W. Wilson

Extensions

Incorrect formula removed by Jinyuan Wang, Mar 10 2020

Status

approved