|
|
A014004
|
|
Pisot sequence E(9,17), a(n) = floor( a(n-1)^2/a(n-2) + 1/2 ).
|
|
1
|
|
|
9, 17, 32, 60, 113, 213, 401, 755, 1422, 2678, 5043, 9497, 17885, 33682, 63432, 119459, 224972, 423680, 797898, 1502646, 2829867, 5329364, 10036557, 18901407, 35596190, 67036742, 126247353, 237756097, 447755619, 843238499, 1588034044, 2990674795, 5632206541
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,1
|
|
LINKS
|
|
|
FORMULA
|
Known not to satisfy any linear recurrence.
|
|
MATHEMATICA
|
nxt[{a_, b_}]:={b, Floor[b^2/a+1/2]}; NestList[nxt, {9, 17}, 40][[All, 1]] (* Harvey P. Dale, Sep 22 2017 *)
|
|
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
}
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|