OFFSET
1,5
COMMENTS
This sequence is the case k=3 of a family of sequences with recurrences a(2*n+1) = a(2*n) + a(2*n-1), a(2*n+2) = k*a(2*n-1) + a(2*n), a(1)=0, a(2)=1. Values of k, for k >= 0, are given by A057979 (k=0), A158780 (k=1), A002965 (k=2), this sequence (k=3). See "Family of sequences for k" link for other connected sequences.
It seems that the ratio of two successive numbers with even, or two successive numbers with odd, indices approaches sqrt(k) for these sequences as n-> infinity.
This algorithm can be found in a historical figure named "Villardsche Figur" of the 13th century. There you can see a geometrical interpretation.
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..1000
W. Beinert, Villardscher Teilungskanon, Lexikon der Typographie
W. Limbrunner, Das Quadrat, ein Wunder der Geometrie. (in German)
Willibald Limbrunner, Family of sequences for k
M-T. Zenner, Villard de Honnecourt and Euclidean Geoometry, Nexus Network Journal 4 (2002) 65-78.
Index entries for linear recurrences with constant coefficients, signature (0,2,0,2).
MATHEMATICA
LinearRecurrence[{0, 2, 0, 2}, {0, 1, 1, 1}, 40] (* G. C. Greubel, Feb 18 2023 *)
PROG
(Magma) I:=[0, 1, 1, 1]; [n le 4 select I[n] else 2*(Self(n-2) +Self(n-4)): n in [1..40]]; // G. C. Greubel, Feb 18 2023
(SageMath)
@CachedFunction
def a(n): # a = A160444
if (n<5): return ((n+1)//3)
else: return 2*(a(n-2) + a(n-4))
[a(n) for n in range(1, 41)] # G. C. Greubel, Feb 18 2023
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Willibald Limbrunner (w.limbrunner(AT)gmx.de), May 14 2009
EXTENSIONS
Edited by R. J. Mathar, May 14 2009
STATUS
approved