|
|
A006138
|
|
a(n) = a(n-1) + 3*a(n-2).
(Formerly M1399)
|
|
15
|
|
|
1, 2, 5, 11, 26, 59, 137, 314, 725, 1667, 3842, 8843, 20369, 46898, 108005, 248699, 572714, 1318811, 3036953, 6993386, 16104245, 37084403, 85397138, 196650347, 452841761, 1042792802, 2401318085, 5529696491, 12733650746, 29322740219, 67523692457, 155491913114
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
The binomial transform of a(n) is b(n) = A006190(n+1), which satisfies b(n) = 3*b(n-1) + b(n-2). - Paul Barry, May 21 2006
An elephant sequence, see A175654. For the corner squares four A[5] vectors, with decimal values 69, 261, 321 and 324, lead to this sequence. For the central square these vectors lead to the companion sequence A105476 (without the first leading 1). - Johannes W. Meijer, Aug 15 2010
Pisano period lengths: 1, 3, 1, 6, 24, 3, 24, 6, 3, 24, 120, 6, 156, 24, 24, 12, 16, 3, 90, 24, ... - R. J. Mathar, Aug 10 2012
The sequence is the INVERT transform of A016116: (1, 1, 2, 2, 4, 4, 8, 8, ...). - Gary W. Adamson, Jul 17 2015
|
|
REFERENCES
|
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
|
|
LINKS
|
|
|
FORMULA
|
G.f.: (1+x)/(1-x-3*x^2). - Paul Barry, May 21 2006
a(n) = Sum_{k=0..n} C(floor((2n-k)/2),n-k)*3^floor(k/2). - Paul Barry, Feb 02 2007
a(n) = (i*sqrt(3))^(n-1)*(i*sqrt(3)*ChebyshevU(n, 1/(2*i*sqrt(3))) + ChebyshevU(n-1, 1/(2*i*sqrt(3)))), where i=sqrt(-1). - G. C. Greubel, Nov 19 2019
|
|
MAPLE
|
|
|
MATHEMATICA
|
Table[(I*Sqrt[3])^(n-1)*(I*Sqrt[3]*ChebyshevU[n, 1/(2*I*Sqrt[3])] + ChebyshevU[n-1, 1/(2*I*Sqrt[3])]), {n, 0, 40}]//Simplify (* G. C. Greubel, Nov 19 2019 *)
|
|
PROG
|
(PARI) main(size)={my(v=vector(size), i); v[1]=1; v[2]=2; for(i=3, size, v[i]=v[i-1]+3*v[i-2]); return(v); } /* Anders Hellström, Jul 17 2015 */
(Magma) [n le 2 select n else Self(n-1)+3*Self(n-2): n in [1..40]]; // Vincenzo Librandi, Sep 15 2016
(Sage)
P.<x> = PowerSeriesRing(ZZ, prec)
return P((1+x)/(1-x-3*x^2)).list()
(GAP) a:=[1, 2];; for n in [3..40] do a[n]:=a[n-1]+3*a[n-2]; od; a; # G. C. Greubel, Nov 19 2019
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|