|
|
A275628
|
|
Pisot sequence E(31,51), a(n)=[a(n-1)^2/a(n-2)+1/2].
|
|
1
|
|
|
31, 51, 84, 138, 227, 373, 613, 1007, 1654, 2717, 4463, 7331, 12042, 19780, 32490, 53367, 87659, 143986, 236507, 388479, 638103, 1048127, 1721619, 2827875, 4644975, 7629684, 12532269, 20585095, 33812403, 55539146, 91226783, 149846127, 246132342, 404288926, 664071752, 1090782516
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,1
|
|
LINKS
|
|
|
FORMULA
|
It is known (Boyd, 1977) that this sequence does not satisfy a linear recurrence. - N. J. A. Sloane, Aug 07 2016
|
|
MAPLE
|
A[0]:= 31:
A[1]:= 51:
for n from 2 to 100 do
A[n]:= floor(A[n-1]^2/A[n-2]+1/2)
od:
|
|
MATHEMATICA
|
nxt[{a_, b_}]:={b, Floor[b^2/a+1/2]}; NestList[nxt, {31, 51}, 40][[All, 1]] (* Harvey P. Dale, Nov 02 2020 *)
|
|
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
}
(Python)
a, b = 31, 51
for i in range(1000):
c, d = divmod(b**2, a)
a, b = b, c + (0 if 2*d < a else 1)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|