Search: a074081 -id:a074081
|
|
A073817
|
|
Tetranacci numbers with different initial conditions: a(n) = a(n-1) + a(n-2) + a(n-3) + a(n-4) starting with a(0)=4, a(1)=1, a(2)=3, a(3)=7.
|
|
+10
42
|
|
|
4, 1, 3, 7, 15, 26, 51, 99, 191, 367, 708, 1365, 2631, 5071, 9775, 18842, 36319, 70007, 134943, 260111, 501380, 966441, 1862875, 3590807, 6921503, 13341626, 25716811, 49570747, 95550687, 184179871, 355018116, 684319421, 1319068095, 2542585503
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,1
|
|
COMMENTS
|
These tetranacci numbers follow the same pattern as Lucas and generalized tribonacci(A001644) numbers: Binet's formula is a(n) = r1^n + r2^n + r3^n + r4^n, with r1, r2, r3, r4 roots of the characteristic polynomial.
For n >= 4, a(n) is the number of cyclic sequences consisting of n zeros and ones that do not contain four consecutive ones provided the positions of the zeros and ones are fixed on a circle. This is proved in Charalambides (1991) and Zhang and Hadjicostas (2015). For example, a(4)=15 because only the sequences 1110, 1101, 1011, 0111, 0011, 0101, 1001, 1010, 0110, 1100, 0001, 0010, 0100, 1000, 0000 avoid four consecutive ones on a circle. (For n=1,2,3 the statement is still true provided we allow the sequence to wrap around itself on a circle. For example, a(2)=3 because only the sequences 00, 01, 10 avoid four consecutive ones when wrapped around on a circle.) - Petros Hadjicostas, Dec 18 2016
|
|
LINKS
|
A. V. Zarelua, On Matrix Analogs of Fermat's Little Theorem, Mathematical Notes, vol. 79, no. 6, 2006, pp. 783-796. Translated from Matematicheskie Zametki, vol. 79, no. 6, 2006, pp. 840-855.
|
|
FORMULA
|
G.f.: (4 - 3*x - 2*x^2 - x^3)/(1 - x - x^2 - x^3 - x^4).
a(n) = 2*a(n-1) - a(n-5), with a(0)=4, a(1)=1, a(2)=3, a(3)=7, a(4)=15. - Vincenzo Librandi, Dec 20 2010
a(n) = Trace(M^n), where M is the 4 X 4 matrix [0, 0, 0, 1; 1, 0, 0, 1; 0, 1, 0, 1; 0, 0, 1, 1], the companion matrix to the monic polynomial x^4 - x^3 - x^2 - x - 1. It follows that the sequence satisfies the Gauss congruences: a(n*p^r) == a(n*p^(r-1)) (mod p^r) for positive integers n and r and all primes p. See Zarelua. - Peter Bala, Dec 31 2022
|
|
MATHEMATICA
|
a[0]=4; a[1]=1; a[2]=3; a[3]=7; a[4]=15; a[n_]:= 2*a[n-1] -a[n-5]; Array[a, 34, 0]
CoefficientList[Series[(4-3x-2x^2-x^3)/(1-x-x^2-x^3-x^4), {x, 0, 40}], x]
LinearRecurrence[{1, 1, 1, 1}, {4, 1, 3, 7}, 40] (* Harvey P. Dale, Jun 01 2015 *)
|
|
PROG
|
(PARI) Vec((4-3*x-2*x^2-x^3)/(1-x-x^2-x^3-x^4) + O(x^40)) \\ Michel Marcus, Jan 29 2016
(Magma) I:=[4, 1, 3, 7]; [n le 4 select I[n] else Self(n-1) +Self(n-2) +Self(n-3) +Self(n-4): n in [1..40]]; // G. C. Greubel, Feb 19 2019
(Sage) ((4-3*x-2*x^2-x^3)/(1-x-x^2-x^3-x^4)).series(x, 40).coefficients(x, sparse=False) # G. C. Greubel, Feb 19 2019
(GAP) a:=[4, 1, 3, 7];; for n in [5..40] do a[n]:=a[n-1]+a[n-2]+a[n-3] +a[n-4]; od; a; # G. C. Greubel, Feb 19 2019
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
Mario Catalani (mario.catalani(AT)unito.it), Aug 12 2002
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|
|
A074193
|
|
Sum of determinants of 2nd order principal minors of powers of the matrix ((1,1,0,0),(1,0,1,0),(1,0,0,1),(1,0,0,0)).
|
|
+10
2
|
|
|
6, -1, -3, -1, 17, -16, -15, 13, 81, -127, -58, 175, 329, -885, -31, 1424, 833, -5543, 2181, 9233, -2298, -31025, 27893, 49495, -54879, -150416, 245697, 204965, -526887, -570895, 1801670, 407711, -3882303, -946397, 11542929, -3442672, -24121039, 10317745, 64959629, -56727711, -127083514
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,1
|
|
COMMENTS
|
Let f(x) = x^4 - x^3 - x^2 - x - 1 be the characteristic polynomial for Tetranacci numbers (A000078). Let {x1,x2,x3,x4} be the roots of f(x). Then a(n) = (x1*x2)^n + (x1*x3)^n + (x1*x4)^n + (x2*x3)^n + (x2*x4)^n + (x3*x4)^n.
Let g(y) = y^6 + y^5 + 2*y^4 + 2*y^3 - 2*y^2 + y - 1 be the characteristic polynomial for a(n). Let {y1,y2,y3,y4,y5,y6} be the roots of g(y). Then a(n) = y1^n + y2^n + y3^n + y4^n + y5^n + y6^n. (End)
|
|
LINKS
|
|
|
FORMULA
|
a(n) = -a(n-1)-2a(n-2)-2a(n-3)+2a(n-4)-a(n-5)+a(n-6).
G.f.: (6+5x+8x^2+6x^3-4x^4+x^5)/(1+x+2x^2+2x^3-2x^4+x^5-x^6).
|
|
MATHEMATICA
|
CoefficientList[Series[(6+5*x+8*x^2+6*x^3-4*x^4+x^5)/(1+x+2*x^2+2*x^3-2*x^4+x^5-x^6), {x, 0, 50}], x]
|
|
PROG
|
(PARI) polsym(x^6 + x^5 + 2*x^4 + 2*x^3 - 2*x^2 + x - 1, 44) \\ Joerg Arndt, Oct 22 2020
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,sign
|
|
AUTHOR
|
Mario Catalani (mario.catalani(AT)unito.it), Aug 20 2002
|
|
STATUS
|
approved
|
|
|
|
|
A074453
|
|
Sum of determinants of 2nd order principal minors of powers of inverse of the matrix ((1,1,0,0),(1,0,1,0),(1,0,0,1),(1,0,0,0)).
|
|
+10
1
|
|
|
6, 1, -3, 1, 17, 16, -15, -13, 81, 127, -58, -175, 329, 885, -31, -1424, 833, 5543, 2181, -9233, -2298, 31025, 27893, -49495, -54879, 150416, 245697, -204965, -526887, 570895, 1801670, -407711, -3882303, 946397, 11542929, 3442672, -24121039, -10317745, 64959629, 56727711, -127083514
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,1
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
a(n)=a(n-1)-2a(n-2)+2a(n-3)+2a(n-4)+a(n-5)+a(n-6).
G.f.: (6-5x+8x^2-6x^3-4x^4-x^5)/(1-x+2x^2-2x^3-2x^4-x^5-x^6).
|
|
MATHEMATICA
|
CoefficientList[Series[(6-5*x+8*x^2-6*x^3-4*x^4-x^5)/(1-x+2*x^2-2*x^3-2*x^4-x^5-x^6), {x, 0, 40}], x]
LinearRecurrence[{1, -2, 2, 2, 1, 1}, {6, 1, -3, 1, 17, 16}, 50] (* Harvey P. Dale, Mar 16 2012 *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,sign
|
|
AUTHOR
|
Mario Catalani (mario.catalani(AT)unito.it), Aug 22 2002
|
|
STATUS
|
approved
|
|
|
Search completed in 0.010 seconds
|