[go: up one dir, main page]

login
Search: a074081 -id:a074081
     Sort: relevance | references | number | modified | created      Format: long | short | data
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.
+0
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
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
Martin Burtscher, Igor Szczyrba, and Rafał Szczyrba, Analytic Representations of the n-anacci Constants and Generalizations Thereof, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.5.
C. A. Charalambides, Lucas numbers and polynomials of order k and the length of the longest circular success run, The Fibonacci Quarterly, 29 (1991), 290-297.
Spiros D. Dafnis, Andreas N. Philippou, and Ioannis E. Livieris, An Alternating Sum of Fibonacci and Lucas Numbers of Order k, Mathematics (2020) Vol. 9, 1487.
Petros Hadjicostas, Cyclic Compositions of a Positive Integer with Parts Avoiding an Arithmetic Sequence, Journal of Integer Sequences, 19 (2016), #16.8.2.
Tony D. Noe and Jonathan Vos Post, Primes in Fibonacci n-step and Lucas n-step Sequences, J. of Integer Sequences, Vol. 8 (2005), Article 05.4.4
J. L. Ramírez and V. F. Sirvent, A Generalization of the k-Bonacci Sequence from Riordan Arrays, The Electronic Journal of Combinatorics, 22(1) (2015), #P1.38.
Yüksel Soykan, Gaussian Generalized Tetranacci Numbers, arXiv:1902.03936 [math.NT], 2019.
Yüksel Soykan, Tetranacci and Tetranacci-Lucas Quaternions, arXiv:1902.05868 [math.RA], 2019.
Yüksel Soykan, Summation Formulas for Generalized Tetranacci Numbers, Asian Journal of Advanced Research and Reports (2019) Vol. 7, No. 2, Article No. AJARR.52434, 1-12.
Yüksel Soykan, Properties of Generalized (r, s, t, u)-Numbers, Earthline J. of Math. Sci. (2021) Vol. 5, No. 2, 297-327.
Kai Wang, Identities for generalized enneanacci numbers, Generalized Fibonacci Sequences (2020).
Eric Weisstein's World of Mathematics, Fibonacci n-Step Number.
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.
L. Zhang and P. Hadjicostas, On sequences of independent Bernoulli trials avoiding the pattern '11..1', Math. Scientist, 40 (2015), 89-96.
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) = A000078(n+2) + 2*A000078(n+1) + 3*A000078(n) + 4*A000078(n-1). - Advika Srivastava, Aug 22 2019
a(n) = 8*a(n-3) - a(n-5) - 2*a(n-6) - 4*a(n-7). - Advika Srivastava, Aug 25 2019
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
Cf. A000078, A001630, A001644, A000032, A106295 (Pisano periods). Two other versions: A001648, A074081.
KEYWORD
nonn,easy
AUTHOR
Mario Catalani (mario.catalani(AT)unito.it), Aug 12 2002
EXTENSIONS
Typo in definition corrected by Vincenzo Librandi, Dec 20 2010
STATUS
approved
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)).
+0
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
OFFSET
0,1
COMMENTS
From Kai Wang, Oct 21 2020: (Start)
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)
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).
abs(a(n)) = abs(A074453(n)). - Joerg Arndt, Oct 22 2020
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
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)).
+0
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
OFFSET
0,1
COMMENTS
a(n) is the reflected (A074058) sequence of sequence A074193.
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).
abs(a(n)) = abs(A074193(n)). - Joerg Arndt, Oct 22 2020
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.008 seconds