A084058 -id:A084058 - cp's OEIS frontend

A001541 a(0) = 1, a(1) = 3; for n > 1, a(n) = 6*a(n-1) - a(n-2).

1, 3, 17, 99, 577, 3363, 19601, 114243, 665857, 3880899, 22619537, 131836323, 768398401, 4478554083, 26102926097, 152139002499, 886731088897, 5168247530883, 30122754096401, 175568277047523, 1023286908188737, 5964153172084899, 34761632124320657

Offset: 0

N. J. A. Sloane

Chebyshev polynomials of the first kind evaluated at 3.
This sequence gives the values of x in solutions of the Diophantine equation x^2 - 8*y^2 = 1, the corresponding values of y are in A001109. For n > 0, the ratios a(n)/A001090(n) may be obtained as convergents to sqrt(8): either successive convergents of [3; -6] or odd convergents of [2; 1, 4]. - Lekraj Beedassy, Sep 09 2003 [edited by Jon E. Schoenfield, May 04 2014]
Also gives solutions to the equation x^2 - 1 = floor(x*r*floor(x/r)) where r = sqrt(8). - Benoit Cloitre, Feb 14 2004
Appears to give all solutions greater than 1 to the equation: x^2 = ceiling(x*r*floor(x/r)) where r = sqrt(2). - Benoit Cloitre, Feb 24 2004
This sequence give numbers n such that (n-1)*(n+1)/2 is a perfect square. Remark: (i-1)*(i+1)/2 = (i^2-1)/2 = -1 = i^2 with i = sqrt(-1) so i is also in the sequence. - Pierre CAMI, Apr 20 2005
a(n) is prime for n = {1, 2, 4, 8}. Prime a(n) are {3, 17, 577, 665857}, which belong to A001601(n). a(2k-1) is divisible by a(1) = 3. a(4k-2) is divisible by a(2) = 17. a(8k-4) is divisible by a(4) = 577. a(16k-8) is divisible by a(8) = 665857. - Alexander Adamchuk, Nov 24 2006
The upper principal convergents to 2^(1/2), beginning with 3/2, 17/12, 99/70, 577/408, comprise a strictly decreasing sequence; essentially, numerators=A001541 and denominators=A001542. - Clark Kimberling, Aug 26 2008
Also index of sequence A082532 for which A082532(n) = 1. - Carmine Suriano, Sep 07 2010
Numbers n such that sigma(n-1) and sigma(n+1) are both odd numbers. - Juri-Stepan Gerasimov, Mar 28 2011
Also, numbers such that floor(a(n)^2/2) is a square: base 2 analog of A031149, A204502, A204514, A204516, A204518, A204520, A004275, A001075. - M. F. Hasler, Jan 15 2012
Numbers such that 2n^2 - 2 is a square. Also integer square roots of the expression 2*n^2 + 1, at values of n given by A001542. Also see A228405 regarding 2n^2 -+ 2^k generally for k >= 0. - Richard R. Forberg, Aug 20 2013
Values of x (or y) in the solutions to x^2 - 6xy + y^2 + 8 = 0. - Colin Barker, Feb 04 2014
Panda and Ray call the numbers in this sequence the Lucas-balancing numbers C_n (see references and links).
Partial sums of X or X+1 of Pythagorean triples (X,X+1,Z). - Peter M. Chema, Feb 03 2017
a(n)/A001542(n) is the closest rational approximation to sqrt(2) with a numerator not larger than a(n), and 2*A001542(n)/a(n) is the closest rational approximation to sqrt(2) with a denominator not larger than a(n). These rational approximations together with those obtained from the sequences A001653 and A002315 give a complete set of closest rational approximations to sqrt(2) with restricted numerator or denominator. a(n)/A001542(n) > sqrt(2) > 2*A001542(n)/a(n). - A.H.M. Smeets, May 28 2017
x = a(n), y = A001542(n) are solutions of the Diophantine equation x^2 - 2y^2 = 1 (Pell equation). x = 2*A001542(n), y = a(n) are solutions of the Diophantine equation x^2 - 2y^2 = -2. Both together give the set of fractional approximations for sqrt(2) obtained from limited fractions obtained from continued fraction representation to sqrt(2). - A.H.M. Smeets, Jun 22 2017
a(n) is the radius of the n-th circle among the sequence of circles generated as follows: Starting with a unit circle centered at the origin, every subsequent circle touches the previous circle as well as the two limbs of hyperbola x^2 - y^2 = 1, and lies in the region y > 0. - Kaushal Agrawal, Nov 10 2018
All of the positive integer solutions of a*b+1=x^2, a*c+1=y^2, b*c+1=z^2, x+z=2*y, 0A001542(n), b=A005319(n), c=A001542(n+1), x=A001541(n), y=A001653(n+1), z=A002315(n) with 0Michael Somos, Jun 26 2022

			99^2 + 99^2 = 140^2 + 2. - _Carmine Suriano_, Jan 05 2015
G.f. = 1 + 3*x + 17*x^2 + 99*x^3 + 577*x^4 + 3363*x^5 + 19601*x^6 + 114243*x^7 + ...

Julio R. Bastida, Quadratic properties of a linearly recurrent sequence. Proceedings of the Tenth Southeastern Conference on Combinatorics, Graph Theory and Computing (Florida Atlantic Univ., Boca Raton, Fla., 1979), pp. 163--166, Congress. Numer., XXIII-XXIV, Utilitas Math., Winnipeg, Man., 1979. MR0561042 (81e:10009)
J. W. L. Glaisher, On Eulerian numbers (formulas, residues, end-figures), with the values of the first twenty-seven, Quarterly Journal of Mathematics, vol. 45, 1914, pp. 1-51.
G. K. Panda, Some fascinating properties of balancing numbers, In Proc. of Eleventh Internat. Conference on Fibonacci Numbers and Their Applications, Cong. Numerantium 194 (2009), 185-189.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 257-258.
P.-F. Teilhet, Query 2376, L'Intermédiaire des Mathématiciens, 11 (1904), 138-139. - N. J. A. Sloane, Mar 08 2022

T. D. Noe, Table of n, a(n) for n = 0..200
I. Adler, Three Diophantine equations - Part II, Fib. Quart., 7 (1969), 181-193.
Christian Aebi and Grant Cairns, Lattice Equable Parallelograms, arXiv:2006.07566 [math.NT], 2020.
Jean-Paul Allouche, Zeta-regularization of arithmetic sequences, EPJ Web of Conferences (2020) Vol. 244, 01008.
Hacène Belbachir, Soumeya Merwa Tebtoub, and László Németh, Ellipse Chains and Associated Sequences, J. Int. Seq., Vol. 23 (2020), Article 20.8.5.
H. Brocard, Notes élémentaires sur le problème de Peel, Nouvelle Correspondance Mathématique, 4 (1878), 161-169.
John M. Campbell, An Integral Representation of Kekulé Numbers, and Double Integrals Related to Smarandache Sequences, arXiv preprint arXiv:1105.3399 [math.GM], 2011.
P. Catarino, H. Campos, and P. Vasco, On some identities for balancing and cobalancing numbers, Annales Mathematicae et Informaticae, 45 (2015) pp. 11-24.
Kwang-Wu Chen and Yu-Ren Pan, Greatest Common Divisors of Shifted Horadam Sequences, J. Int. Seq., Vol. 23 (2020), Article 20.5.8.
S. Falcon, Relationships between Some k-Fibonacci Sequences, Applied Mathematics, 2014, 5, 2226-2234.
Morgan Fiebig, aBa Mbirika, and Jürgen Spilker, Period patterns, entry points, and orders in the Lucas sequences: theory and applications, arXiv:2408.14632 [math.NT], 2024. See p. 5.
Robert Frontczak, A Note on Hybrid Convolutions Involving Balancing and Lucas-Balancing Numbers, Applied Mathematical Sciences, Vol. 12, 2018, No. 25, 1201-1208.
Robert Frontczak, Sums of Balancing and Lucas-Balancing Numbers with Binomial Coefficients, International Journal of Mathematical Analysis (2018) Vol. 12, No. 12, 585-594.
Robert Frontczak, Powers of Balancing Polynomials and Some Consequences for Fibonacci Sums, International Journal of Mathematical Analysis (2019) Vol. 13, No. 3, 109-115.
Robert Frontczak, Identities for generalized balancing numbers, Notes on Number Theory and Discrete Mathematics (2019) Vol. 25, No. 2, 169-180.
Robert Frontczak and Taras Goy, Additional close links between balancing and Lucas-balancing polynomials, arXiv:2007.14048 [math.NT], 2020.
Robert Frontczak and Taras Goy, More Fibonacci-Bernoulli relations with and without balancing polynomials, arXiv:2007.14618 [math.NT], 2020.
Robert Frontczak and Taras Goy, Lucas-Euler relations using balancing and Lucas-balancing polynomials, arXiv:2009.09409 [math.NT], 2020.
O. Khadir, K. Liptai, and L. Szalay, On the Shifted Product of Binary Recurrences, J. Int. Seq. 13 (2010), 10.6.1.
J. M. Katri and D. R. Byrkit, Problem E1976, Amer. Math. Monthly, 75 (1968), 683-684.
Tanya Khovanova, Recursive Sequences
D. H. Lehmer, A cotangent analogue of continued fractions, Duke Math. J., 4 (1935), 323-340.
D. H. Lehmer, A cotangent analogue of continued fractions, Duke Math. J., 4 (1935), 323-340. [Annotated scanned copy]
D. H. Lehmer, Lacunary recurrence formulas for the numbers of Bernoulli and Euler, Annals Math., 36 (1935), 637-649.
Dino Lorenzini and Z. Xiang, Integral points on variable separated curves, Preprint 2016.
aBa Mbirika, Janeè Schrader, and Jürgen Spilker, Pell and associated Pell braid sequences as GCDs of sums of k consecutive Pell, balancing, and related numbers, arXiv:2301.05758 [math.NT], 2023. See also J. Int. Seq. (2023) Vol. 26, Art. 23.6.4.
A. Patra, G. K. Panda, and T. Khemaratchatakumthorn, Exact divisibility by powers of the balancing and Lucas-balancing numbers, Fibonacci Quart., 59:1 (2021), 57-64; see C(n).
Robert Phillips, Polynomials of the form 1+4ke+4ke^2, 2008.
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
Prasanta K. Ray, Curious congruences for balancing numbers, Int. J. Contemp. Math. Sci. 7 (2012), 881-889.
Soumeya M. Tebtoub, Hacène Belbachir, and László Németh, Integer sequences and ellipse chains inside a hyperbola, Proceedings of the 1st International Conference on Algebras, Graphs and Ordered Sets (ALGOS 2020), hal-02918958 [math.cs], 17-18.
Bibhu Prasad Tripathy and Bijan Kumar Patel, On balancing and Lucas-balancing numbers expressible as product of two k-Fibonacci numbers, arXiv:2508.12238 [math.NT], 2025.
N. J. Wildberger, Pell's equation without irrational numbers, J. Int. Seq. 13 (2010), 10.4.3, Section 4.
Index entries for two-way infinite sequences
Index entries for linear recurrences with constant coefficients, signature (6,-1).
Index entries for sequences related to Chebyshev polynomials.

Bisection of A001333. A003499(n) = 2a(n).
Cf. A055997 = numbers n such that n(n-1)/2 is a square.
Row 1 of array A188645.
Cf. A046090, A001109, A053142, A084130, A001601, A056771, A077420, A005319, A082532, A001542.
Cf. A055792 (terms squared), A132592.

a001541 n = a001541_list !! (n-1)
a001541_list =
1 : 3 : zipWith (-) (map (* 6) $ tail a001541_list) a001541_list
-- Reinhard Zumkeller, Oct 06 2011
(Scheme, with memoization-macro definec)
(definec (A001541 n) (cond ((zero? n) 1) ((= 1 n) 3) (else (- (* 6 (A001541 (- n 1))) (A001541 (- n 2))))))
;; Antti Karttunen, Oct 04 2016

[n: n in [1..10000000] |IsSquare(8*(n^2-1))]; // Vincenzo Librandi, Nov 18 2010

a[0]:=1: a[1]:=3: for n from 2 to 26 do a[n]:=6*a[n-1]-a[n-2] od: seq(a[n], n=0..20); # Zerinvary Lajos, Jul 26 2006
A001541:=-(-1+3*z)/(1-6*z+z**2); # Simon Plouffe in his 1992 dissertation

Table[Simplify[(1/2) (3 + 2 Sqrt[2])^n + (1/2) (3 - 2 Sqrt[2])^n], {n, 0, 20}] (* Artur Jasinski, Feb 10 2010 *)
a[ n_] := If[n == 0, 1, With[{m = Abs @ n}, m Sum[4^i Binomial[m + i, 2 i]/(m + i), {i, 0, m}]]]; (* Michael Somos, Jul 11 2011 *)
a[ n_] := ChebyshevT[ n, 3]; (* Michael Somos, Jul 11 2011 *)
LinearRecurrence[{6, -1}, {1, 3}, 50] (* Vladimir Joseph Stephan Orlovsky, Feb 12 2012 *)

{a(n) = real((3 + quadgen(32))^n)}; /* Michael Somos, Apr 07 2003 */

{a(n) = subst( poltchebi( abs(n)), x, 3)}; /* Michael Somos, Apr 07 2003 */

{a(n) = if( n<0, a(-n), polsym(1 - 6*x + x^2, n) [n+1] / 2)}; /* Michael Somos, Apr 07 2003 */

{a(n) = polchebyshev( n, 1, 3)}; /* Michael Somos, Jul 11 2011 */

a(n)=([1,2,2;2,1,2;2,2,3]^n)[3,3] \\ Vim Wenders, Mar 28 2007

G.f.: (1-3*x)/(1-6*x+x^2). - Barry E. Williams and Wolfdieter Lang, May 05 2000
E.g.f.: exp(3*x)*cosh(2*sqrt(2)*x). Binomial transform of A084128. - Paul Barry, May 16 2003
From N. J. A. Sloane, May 16 2003: (Start)
a(n) = sqrt(8*((A001109(n))^2) + 1).
a(n) = T(n, 3), with Chebyshev's T-polynomials A053120. (End)
a(n) = ((3+2*sqrt(2))^n + (3-2*sqrt(2))^n)/2.
a(n) = cosh(2*n*arcsinh(1)). - Herbert Kociemba, Apr 24 2008
a(n) ~ (1/2)*(sqrt(2) + 1)^(2*n). - Joe Keane (jgk(AT)jgk.org), May 15 2002
For all elements x of the sequence, 2*x^2 - 2 is a square. Limit_{n -> infinity} a(n)/a(n-1) = 3 + 2*sqrt(2). - Gregory V. Richardson, Oct 10 2002 [corrected by Peter Pein, Mar 09 2009]
a(n) = 3*A001109(n) - A001109(n-1), n >= 1. - Barry E. Williams and Wolfdieter Lang, May 05 2000
For n >= 1, a(n) = A001652(n) - A001652(n-1). - Charlie Marion, Jul 01 2003
From Paul Barry, Sep 18 2003: (Start)
a(n) = ((-1+sqrt(2))^n + (1+sqrt(2))^n + (1-sqrt(2))^n + (-1-sqrt(2))^n)/4 (with interpolated zeros).
E.g.f.: cosh(x)*cosh(sqrt(2)x) (with interpolated zeros). (End)
For n > 0, a(n)^2 + 1 = 2*A001653(n-1)*A001653(n). - Charlie Marion, Dec 21 2003
a(n)^2 + a(n+1)^2 = 2*(A001653(2*n+1) - A001652(2*n)). - Charlie Marion, Mar 17 2003
a(n) = Sum_{k >= 0} binomial(2*n, 2*k)*2^k = Sum_{k >= 0} A086645(n, k)*2^k. - Philippe Deléham, Feb 29 2004
a(n)*A002315(n+k) = A001652(2*n+k) + A001652(k) + 1; for k > 0, a(n+k)*A002315(n) = A001652(2*n+k) - A001652(k-1). - Charlie Marion, Mar 17 2003
For n > k, a(n)*A001653(k) = A011900(n+k) + A053141(n-k-1). For n <= k, a(n)*A001653(k) = A011900(n+k) + A053141(k-n). - Charlie Marion, Oct 18 2004
A053141(n+1) + A055997(n+1) = a(n+1) + A001109(n+1). - Creighton Dement, Sep 16 2004
a(n+1) - A001542(n+1) = A090390(n+1) - A046729(n) = A001653(n); a(n+1) - 4*A079291(n+1) = (-1)^(n+1). Formula generated by the floretion - .5'i + .5'j - .5i' + .5j' - 'ii' + 'jj' - 2'kk' + 'ij' + .5'ik' + 'ji' + .5'jk' + .5'ki' + .5'kj' + e. - Creighton Dement, Nov 16 2004
a(n) = sqrt( A055997(2*n) ). - Alexander Adamchuk, Nov 24 2006
a(2n) = A056771(n). a(2*n+1) = 3*A077420(n). - Alexander Adamchuk, Feb 01 2007
a(n) = (A000129(n)^2)*4 + (-1)^n. - Vim Wenders, Mar 28 2007
2*a(k)*A001653(n)*A001653(n+k) = A001653(n)^2 + A001653(n+k)^2 + A001542(k)^2. - Charlie Marion, Oct 12 2007
a(n) = A001333(2*n). - Ctibor O. Zizka, Aug 13 2008
A028982(a(n)-1) + 2 = A028982(a(n)+1). - Juri-Stepan Gerasimov, Mar 28 2011
a(n) = 2*A001108(n) + 1. - Paul Weisenhorn, Dec 17 2011
a(n) = sqrt(2*x^2 + 1) with x being A001542(n). - Zak Seidov, Jan 30 2013
a(2n) = 2*a(n)^2 - 1 = a(n)^2 + 2*A001542(n)^2. a(2*n+1) = 1 + 2*A002315(n)^2. - Steven J. Haker, Dec 04 2013
a(n) = 3*a(n-1) + 4*A001542(n-1); e.g., a(4) = 99 = 3*17 + 4*12. - Zak Seidov, Dec 19 2013
a(n) = cos(n * arccos(3)) = cosh(n * log(3 + 2*sqrt(2))). - Daniel Suteu, Jul 28 2016
From Ilya Gutkovskiy, Jul 28 2016: (Start)
Inverse binomial transform of A084130.
Exponential convolution of A000079 and A084058.
Sum_{n>=0} (-1)^n*a(n)/n! = cosh(2*sqrt(2))/exp(3) = 0.4226407909842764637... (End)
a(2*n+1) = 2*a(n)*a(n+1) - 3. - Timothy L. Tiffin, Oct 12 2016
a(n) = a(-n) for all n in Z. - Michael Somos, Jan 20 2017
a(2^n) = A001601(n+1). - A.H.M. Smeets, May 28 2017
a(A298210(n)) = A002350(2*n^2). - A.H.M. Smeets, Jan 25 2018
a(n) = S(n, 6) - 3*S(n-1, 6), for n >= 0, with S(n, 6) = A001109(n+1), (Chebyshev S of A049310). See the first comment and the formula a(n) = T(n, 3). - Wolfdieter Lang, Nov 22 2020
From Peter Bala, Dec 31 2021: (Start)
a(n) = [x^n] (3*x + sqrt(1 + 8*x^2))^n.
The Gauss congruences a(n*p^k) == a(n*p^(k-1)) hold for all prime p and positive integers n and k.
O.g.f. A(x) = 1 + x*d/dx(log(B(x))), where B(x) = 1/sqrt(1 - 6*x + x^2) is the o.g.f. of A001850. (End)
From Peter Bala, Aug 17 2022: (Start)
Sum_{n >= 1} 1/(a(n) - 2/a(n)) = 1/2.
Sum_{n >= 1} (-1)^(n+1)/(a(n) + 1/a(n)) = 1/4.
Sum_{n >= 1} 1/(a(n)^2 - 2) = 1/2 - 1/sqrt(8). (End)
From  Peter Bala, Jun 23 2025: (Start)
Product_{n >= 0} (1 + 1/a(2^n)) = sqrt(2).
Product_{n >= 0} (1 - 1/(2*a(2^n))) = (4/7)*sqrt(2). See A002812. (End)

A098158 Triangle T(n,k) with diagonals T(n,n-k) = binomial(n, 2*k).

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 0, 3, 1, 0, 0, 1, 6, 1, 0, 0, 0, 5, 10, 1, 0, 0, 0, 1, 15, 15, 1, 0, 0, 0, 0, 7, 35, 21, 1, 0, 0, 0, 0, 1, 28, 70, 28, 1, 0, 0, 0, 0, 0, 9, 84, 126, 36, 1, 0, 0, 0, 0, 0, 1, 45, 210, 210, 45, 1, 0, 0, 0, 0, 0, 0, 11, 165, 462, 330, 55, 1, 0, 0, 0, 0, 0, 0, 1, 66, 495, 924

Offset: 0

Paul Barry, Aug 29 2004

Row sums are A011782. Inverse is A065547.
Triangle T(n,k), 0 <= k <= n, read by rows, given by [0, 1, -1, 0, 0, 0, 0, 0, ...] DELTA [1, 0, 1, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Jul 29 2006
Sum of entries in column k is A001519(k+1) (the odd-indexed Fibonacci numbers). - Philippe Deléham, Dec 02 2008
Number of permutations of length n avoiding simultaneously the patterns 123 and 132 with k left-to-right minima. A left-to-right minimum in a permutation a(1)a(2)...a(n) is position i such that a(j) > a(i) for all j < i. - Tian Han, Nov 16 2023

			Rows begin
  1;
  0, 1;
  0, 1, 1;
  0, 0, 3, 1;
  0, 0, 1, 6, 1;

G. C. Greubel, Rows n = 0..49 of triangle, flattened
D. Dumont and J. Zeng, Polynomes d'Euler et les fractions continues de Stieltjes-Rogers, Ramanujan J. 2 (1998) 3, 387-410.
Tian Han and Sergey Kitaev, Joint distributions of statistics over permutations avoiding two patterns of length 3, arXiv:2311.02974 [math.CO], 2023.

Cf. A098157, A034839.

Cf. A119900. - Philippe Deléham, Dec 02 2008

Flat(List([0..12], n-> List([0..n], k-> Binomial(n, 2*(n-k)) ))); # G. C. Greubel, Aug 01 2019

[Binomial(n, 2*(n-k)): k in [0..n], n in [0..12]]; // G. C. Greubel, Aug 01 2019

Table[Binomial[n, 2*(n-k)], {n,0,12}, {k,0,n}]//Flatten (* Michael De Vlieger, Oct 12 2016 *)

{T(n,k)=polcoeff(polcoeff((1-x*y)/((1-x*y)^2-x^2*y)+x*O(x^n), n, x) + y*O(y^k),k,y)} (Hanna)

T(n,k) = binomial(n, 2*(n-k));
for(n=0,12, for(k=0,n, print1(T(n,k), ", "))) \\ G. C. Greubel, Aug 01 2019

[[binomial(n, 2*(n-k)) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Aug 01 2019

T(n,k) = binomial(n,2*(n-k)).
From Tom Copeland, Oct 10 2016: (Start)
E.g.f.: exp(t*x) * cosh(t*sqrt(x)).
O.g.f.: (1/2) * ( 1 / (1 - (1 + sqrt(1/x))*x*t) + 1 / (1 - (1 - sqrt(1/x))*x*t) ).
Row polynomial: x^n * ((1 + sqrt(1/x))^n + (1 - sqrt(1/x))^n) / 2. (End)
Column k is generated by the polynomial Sum_{j=0..floor(k/2)} C(k, 2j) * x^(k-j). - Paul Barry, Jan 22 2005
G.f.: (1-x*y)/((1-x*y)^2 - x^2*y). - Paul D. Hanna, Feb 25 2005
Sum_{k=0..n} x^k*T(n,k)= A009116(n), A000007(n), A011782(n), A006012(n), A083881(n), A081335(n), A090139(n), A145301(n), A145302(n), A145303(n), A143079(n) for x = -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively. - Philippe Deléham, Dec 04 2006, Oct 15 2008, Oct 19 2008
T(n,k) = T(n-1,k-1) + Sum_{i=0..k-1} T(n-2-i,k-1-i); T(0,0)=1; T(n,k)=0 if n < 0 or k < 0 or n < k. E.g.: T(8,5) = T(7,4) + T(6,4) + T(5,3) + T(4,2) + T(3,1) + T(2,0) = 7+15+5+1+0+0 = 28. - Philippe Deléham, Dec 04 2006
Sum_{k=0..n} T(n,k)*x^(n-k) = A000012(n), A011782(n), A001333(n), A026150(n), A046717(n), A084057(n), A002533(n), A083098(n), A084058(n), A003665(n), A002535(n), A133294(n), A090042(n), A125816(n), A133343(n), A133345(n), A120612(n), A133356(n), A125818(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18 respectively. - Philippe Deléham, Dec 24 2007
Sum_{k=0..n} T(n,k)*(-x)^(n-k) = A000012(n), A146559(n), A087455(n), A138230(n), A006495(n), A138229(n) for x = 0,1,2,3,4,5 respectively. - Philippe Deléham, Nov 14 2008
T(n,k) = A085478(k,n-k). - Philippe Deléham, Dec 02 2008
T(n,k) = 2*T(n-1,k-1) + T(n-2,k-1) - T(n-2,k-2), T(0,0) = T(1,1) = 1, T(1,0) = 0 and T(n,k) = 0 if k < 0 or if k > n. - Philippe Deléham, Mar 15 2012

A083100 a(n) = 2a(n-1) + 7a(n-2).

Original entry on oeis.org

1, 9, 25, 113, 401, 1593, 5993, 23137, 88225, 338409, 1294393, 4957649, 18976049, 72655641, 278143625, 1064876737, 4076758849, 15607654857, 59752621657, 228758827313, 875786006225, 3352883803641, 12836269650857, 49142725927201

Offset: 0

Mario Catalani (mario.catalani(AT)unito.it), Apr 23 2003

a(n) = a(n-1) + 8*A015519(n). a(n)/A015519(n+1) converges to sqrt(8).

a(n-1) is the number of compositions of n when there is 1 type of 1 and 8 types of other natural numbers. - Milan Janjic, Aug 13 2010

G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (2, 7).

Essentially a duplicate of A084058.

The following sequences (and others) belong to the same family: A001333, A000129, A026150, A002605, A046717, A015518, A084057, A063727, A002533, A002532, A083098, A083099, A083100, A015519.

I:=[1,1]; [n le 2 select I[n] else 2*Self(n-1) + 7*Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 08 2018

CoefficientList[Series[(1 + 7 x)/(1 - 2 x - 7 x^2), {x, 0, 25}], x] (* Or *) a[n_] := Simplify[((1 + Sqrt[8])^n + (1 - Sqrt[8])^n)/2]; Array[a, 25, 0] (* Or *) LinearRecurrence[{2, 7}, {1, 1}, 28] (* Or *) Table[ MatrixPower[{{1, 2}, {4, 1}}, n][[1, 1]], {n, 0, 25}] (* Robert G. Wilson v, Sep 18 2013 *)

a(n)=([0,1; 7,2]^n*[1;9])[1,1] \\ Charles R Greathouse IV, Apr 06 2016

x='x+O('x^30); Vec((1+7*x)/(1-2*x-7*x^2)) \\ G. C. Greubel, Jan 08 2018

G.f.: (1+7*x)/(1-2*x-7*x^2).
a(n) = binomial transform of 1,8,8,64,64,512. - Al Hakanson (hawkuu(AT)gmail.com), Aug 17 2009
If p[1]=1, and p[i]=8,(i>1), and if A is Hessenberg matrix of order n defined by: A[i,j]=p[j-i+1], (i<=j), A[i,j]=-1, (i=j+1), and A[i,j]=0 otherwise. Then, for n>=1, a(n-1)=det A. - Milan Janjic, Apr 29 2010
G.f.: G(0)/(2*x) - 1/x, where G(k)= 1 + 1/(1 - x*(8*k-1)/(x*(8*k+7) - 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 28 2013

A084128 a(n) = 4a(n-1) + 4a(n-2), a(0)=1, a(1)=2.

Original entry on oeis.org

1, 2, 12, 56, 272, 1312, 6336, 30592, 147712, 713216, 3443712, 16627712, 80285696, 387653632, 1871757312, 9037643776, 43637604352, 210700992512, 1017354387456, 4912221519872, 23718303629312, 114522100596736, 552961616904192, 2669934870003712

Offset: 0

Paul Barry, May 16 2003

Original name was: Generalized Fibonacci sequence.

Binomial transform of A084058.

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
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.
Tanya Khovanova, Recursive Sequences
Index entries for linear recurrences with constant coefficients, signature (4,4).

Cf. A001333 A001541, A002203, A057087, A084058, A084128, A201730.

Appears in A086346, A086347 and A086348. - Johannes W. Meijer, Aug 01 2010

[2^(n-1)*Evaluate(DicksonFirst(n,-1), 2): n in [0..40]]; // G. C. Greubel, Oct 13 2022

a:=proc(n) option remember; if n=0 then 1 elif n=1 then 2 else
4*a(n-1)+4*a(n-2); fi; end: seq(a(n), n=0..40); # Wesley Ivan Hurt, Jan 31 2017
a := n -> (2*I)^n*ChebyshevT(n, -I):
seq(simplify(a(n)), n = 0..23); # Peter Luschny, Dec 03 2023

CoefficientList[Series[(2 z - 1)/(4 z^2 + 4 z - 1), {z, 0, 100}], z] (* Vladimir Joseph Stephan Orlovsky, Jul 01 2011 *)
Table[2^(n-1) LucasL[n, 2], {n, 0, 20}] (* Vladimir Reshetnikov, Oct 07 2016 *)
LinearRecurrence[{4,4},{1,2},30] (* Harvey P. Dale, Mar 01 2018 *)

a(n)=if(n<0,0,polsym(4+4*x-x^2,n)[n+1]/2)

[lucas_number2(n,4,-4)/2 for n in range(0, 23)] # Zerinvary Lajos, May 14 2009

a(n) = 2^n * A001333(n).
G.f.: (1-2*x)/(1-4*x-4*x^2).
a(n) = 4*a(n-1) + 4*a(n-2), a(0)=1, a(1)=2.
a(n) = (2 + 2*sqrt(2))^n/2 + (2 - 2*sqrt(2))^n/2.
E.g.f.: exp(2*x)*cosh(2*x*sqrt(2)).
From Johannes W. Meijer, Aug 01 2010: (Start)
Lim_{k->infinity} a(n+k)/a(k) = A084128(n) + 2*A057087(n-1)*sqrt(2).
Lim_{n->infinity} A084128(n)/A057087(n-1) = sqrt(2). (End)
a(n) = Sum_{k=0..n} A201730(n,k)*7^k. - Philippe Deléham, Dec 06 2011
G.f.: G(0)/2, where G(k)= 1 + 1/(1 - x*(4*k-2)/(x*(4*k+2) - 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 27 2013
a(n) = 2^(n-1)*A002203(n). - Vladimir Reshetnikov, Oct 07 2016

A124182 A skewed version of triangular array A081277.

Original entry on oeis.org

1, 0, 1, 0, 1, 2, 0, 0, 3, 4, 0, 0, 1, 8, 8, 0, 0, 0, 5, 20, 16, 0, 0, 0, 1, 18, 48, 32, 0, 0, 0, 0, 7, 56, 112, 64, 0, 0, 0, 0, 1, 32, 160, 256, 128, 0, 0, 0, 0, 0, 9, 120, 432, 576, 256, 0, 0, 0, 0, 0, 1, 50, 400, 1120, 1280, 512

Offset: 0

Philippe Deléham, Dec 05 2006

Triangle T(n,k), 0 <= k <= n, read by rows given by [0, 1, -1, 0, 0, 0, 0, 0, 0, ...] DELTA [1, 1, 0, 0, 0, 0, 0, 0, 0,...] where DELTA is the operator defined in A084938. Falling diagonal sums in A052980.

			Triangle begins:
  1;
  0, 1;
  0, 1, 2;
  0, 0, 3, 4;
  0, 0, 1, 8,  8;
  0, 0, 0, 5, 20, 16;
  0, 0, 0, 1, 18, 48,  32;
  0, 0, 0, 0,  7, 56, 112,  64;
  0, 0, 0, 0,  1, 32, 160, 256,  128;
  0, 0, 0, 0,  0,  9, 120, 432,  576,  256;
  0, 0, 0, 0,  0,  1,  50, 400, 1120, 1280, 512;

Cf. A025192 (column sums). Diagonals include A011782, A001792, A001793, A001794, A006974, A006975, A006976.

T(0,0)=T(1,1)=1, T(n,k)=0 if n < k or if k < 0, T(n,k) = T(n-2,k-1) + 2*T(n-1,k-1).
Sum_{k=0..n} x^k*T(n,k) = (-1)^n*A090965(n), (-1)^n*A084120(n), (-1)^n*A006012(n), A033999(n), A000007(n), A001333(n), A084059(n) for x = -4, -3, -2, -1, 0, 1, 2 respectively.
Sum_{k=0..floor(n/2)} T(n-k,k) = Fibonacci(n-1) = A000045(n-1).
Sum_{k=0..n} T(n,k)*x^(n-k) = A000012(n), A011782(n), A001333(n), A026150(n), A046717(n), A084057(n), A002533(n), A083098(n), A084058(n), A003665(n), A002535(n), A133294(n), A090042(n), A125816(n), A133343(n), A133345(n), A120612(n), A133356(n), A125818(n) for x = -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17 respectively. - Philippe Deléham, Dec 26 2007
Sum_{k=0..n} T(n,k)*(-x)^(n-k) = A011782(n), A000012(n), A146559(n), A087455(n), A138230(n), A006495(n), A138229(n) for x= 0,1,2,3,4,5,6 respectively. - Philippe Deléham, Nov 14 2008
G.f.: (1-y*x)/(1-2y*x-y*x^2). - Philippe Deléham, Dec 04 2011
Sum_{k=0..n} T(n,k)^2 = A002002(n) for n > 0. - Philippe Deléham, Dec 04 2011

A201701 Riordan triangle ((1-x)/(1-2x), x^2/(1-2x)).

Original entry on oeis.org

1, 1, 0, 2, 1, 0, 4, 3, 0, 0, 8, 8, 1, 0, 0, 16, 20, 5, 0, 0, 0, 32, 48, 18, 1, 0, 0, 0, 64, 112, 56, 7, 0, 0, 0, 0, 128, 256, 160, 32, 1, 0, 0, 0, 0, 256, 576, 432, 120, 9, 0, 0, 0, 0, 0, 512, 1280, 1120, 400, 50, 1, 0, 0, 0, 0, 0

Offset: 0

Philippe Deléham, Dec 03 2011

Triangle T(n,k), read by rows, given by (1,1,0,0,0,0,0,0,0,...) DELTA (0,1,-1,0,0,0,0,0,0,0,...) where DELTA is the operator defined in A084938.
Skewed version of triangle in A200139.
Triangle without zeros: A207537.
For the version with negative odd numbered columns, which is Riordan ((1-x)/(1-2*x), -x^2/(1-2*x)) see comments on A028297 and A039991. - Wolfdieter Lang, Aug 06 2014
This is an example of a stretched Riordan array in the terminology of Section 2 of Corsani et al. - Peter Bala, Jul 14 2015

			The triangle T(n,k) begins:
  n\k      0     1     2     3     4    5   6  7 8 9 10 11 ...
  0:       1
  1:       1     0
  2:       2     1     0
  3:       4     3     0     0
  4:       8     8     1     0     0
  5:      16    20     5     0     0    0
  6:      32    48    18     1     0    0   0
  7:      64   112    56     7     0    0   0  0
  8:     128   256   160    32     1    0   0  0 0
  9:     256   576   432   120     9    0   0  0 0 0
  10:    512  1280  1120   400    50    1   0  0 0 0  0
  11:   1024  2816  2816  1232   220   11   0  0 0 0  0  0
  ...  reformatted and extended. - _Wolfdieter Lang_, Aug 06 2014

C. Corsani, D. Merlini, and R. Sprugnoli, Left-inversion of combinatorial sums, Discrete Mathematics, 180 (1998) 107-122.

Columns include A011782, A001792, A001793, A001794, A006974, A006975, A006976.
Diagonals sums are in A052980.
Cf. A028297, A081265, A124182, A131577, A039991 (zero-columns deleted, unsigned and zeros appended).
Cf. A098158, A200139, A207537.
Cf. A028297 (signed version, zeros deleted). Cf. A034839.

(* The function RiordanArray is defined in A256893. *)
RiordanArray[(1 - #)/(1 - 2 #)&, #^2/(1 - 2 #)&, 11] // Flatten (* Jean-François Alcover, Jul 16 2019 *)

T(n,k) = 2*T(n-1,k) + T(n-2,k-1) with T(0,0) = T(1,0) = 1, T(1,1) = 0 and T(n,k) = 0 for k<0 or for n

Sum_{k=0..n} T(n,k)^2 = A002002(n) for n>0.

Sum_{k=0..n} T(n,k)*x^k = A138229(n), A006495(n), A138230(n), A087455(n), A146559(n), A000012(n), A011782(n), A001333(n), A026150(n), A046717(n), A084057(n), A002533(n), A083098(n), A084058(n), A003665(n), A002535(n), A133294(n), A090042(n), A125816(n), A133343(n), A133345(n), A120612(n), A133356(n), A125818(n) for x = -6,-5,-4,-3,-2,-1,0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17 respectively.

G.f.: (1-x)/(1-2*x-y*x^2). - Philippe Deléham, Mar 03 2012

From Peter Bala, Jul 14 2015: (Start)

Factorizes as A034839 * A007318 = (1/(1 - x), x^2/(1 - x)^2) * (1/(1 - x), x/(1 - x)) as a product of Riordan arrays.

T(n,k) = Sum_{i = k..floor(n/2)} binomial(n,2*i) *binomial(i,k). (End)

Name changed, keyword:easy added, crossrefs A028297 and A039991 added, and g.f. corrected by Wolfdieter Lang, Aug 06 2014

A202023 Triangle T(n,k), read by rows, given by (1, 0, 1, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (0, 1, -1, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 3, 0, 0, 1, 6, 1, 0, 0, 1, 10, 5, 0, 0, 0, 1, 15, 15, 1, 0, 0, 0, 1, 21, 35, 7, 0, 0, 0, 0, 1, 28, 70, 28, 1, 0, 0, 0, 0, 1, 36, 126, 84, 9, 0, 0, 0, 0, 0, 1, 45, 210, 210, 45, 1, 0, 0, 0, 0, 0

Offset: 0

Philippe Deléham, Dec 10 2011

Riordan array (1/(1-x), x^2/(1-x)^2).
A skewed version of triangular array A085478.
Mirror image of triangle in A098158.
Sum_{k, 0<=k<=n} T(n,k)*x^k = A138229(n), A006495(n), A138230(n),A087455(n), A146559(n), A000012(n), A011782(n), A001333(n),A026150(n), A046717(n), A084057(n), A002533(n), A083098(n),A084058(n), A003665(n), A002535(n), A133294(n), A090042(n),A125816(n), A133343(n), A133345(n), A120612(n), A133356(n), A125818(n) for x = -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18 respectively.
Sum_{k, 0<=k<=n} T(n,k)*x^(n-k) = A009116(n), A000007(n), A011782(n), A006012(n), A083881(n), A081335(n), A090139(n), A145301(n), A145302(n), A145303(n), A143079(n) for x = -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively.
From Gus Wiseman, Jul 08 2025: (Start)
After the first row this is also the number of subsets of {1..n-1} with k maximal runs (sequences of consecutive elements increasing by 1) for k = 0..n. For example, row n = 5 counts the following subsets:
  {}  {1}        {1,3}    .  .  .
      {2}        {1,4}
      {3}        {2,4}
      {4}        {1,2,4}
      {1,2}      {1,3,4}
      {2,3}
      {3,4}
      {1,2,3}
      {2,3,4}
      {1,2,3,4}
Requiring n-1 gives A202064.
For anti-runs instead of runs we have A384893.
(End)

			Triangle begins :
1
1, 0
1, 1, 0
1, 3, 0, 0
1, 6, 1, 0, 0
1, 10, 5, 0, 0, 0
1, 15, 15, 1, 0, 0, 0
1, 21, 35, 7, 0, 0, 0, 0
1, 28, 70, 28, 1, 0, 0, 0, 0

Column k = 1 is A000217.
Column k = 2 is A000332.
Row sums are A011782 (or A000079 shifted right).
Removing all zeros gives A034839 (requiring n-1 A034867).
Last nonzero term in each row appears to be A093178, requiring n-1 A124625.
Reversing rows gives A098158, without zeros A109446.
Without the k = 0 column we get A210039.
Row maxima appear to be A214282.
A116674 counts strict partitions by number of maximal runs, for anti-runs A384905.
A268193 counts integer partitions by number of maximal runs, for anti-runs A384881.
Cf. A000045, A000071, A007318, A010027, A053538, A084938, A202064, A208342, A210034, A384175, A384893.

Table[Length[Select[Subsets[Range[n-1]],Length[Split[#,#2==#1+1&]]==k&]],{n,0,10},{k,0,n}] (* Gus Wiseman, Jul 08 2025 *)

T(n,k) = binomial(n,2k).
G.f.: (1-x)/((1-x)^2-y*x^2).
T(n,k)= Sum_{j, j>=0} T(n-1-j,k-1)*j with T(n,0)=1 and T(n,k)= 0 if k<0 or if nT(n,k) = 2*T(n-1,k) + T(n-2,k-1) - T(n-2,k) for n>1, T(0,0) = T(1,0) = 1, T(1,1) = 0, T(n,k) = 0 if k>n or if k<0. - Philippe Deléham, Nov 10 2013

A084097 Square array whose rows have e.g.f. exp(x)cosh(sqrt(k)x), k>=0, read by ascending antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 4, 1, 1, 1, 4, 7, 8, 1, 1, 1, 5, 10, 17, 16, 1, 1, 1, 6, 13, 28, 41, 32, 1, 1, 1, 7, 16, 41, 76, 99, 64, 1, 1, 1, 8, 19, 56, 121, 208, 239, 128, 1, 1, 1, 9, 22, 73, 176, 365, 568, 577, 256, 1, 1, 1, 10, 25, 92, 241, 576, 1093, 1552, 1393, 512, 1

Offset: 0

Paul Barry, May 11 2003

Rows are the binomial transforms of expansions of cosh(sqrt(k)*x), k >= 0.

			Array, A(n,k), begins:
.n\k.........0..1...2...3....4.....5......6......7.......8........9.......10
.0: A000012..1..1...1...1....1.....1......1......1.......1........1........1
.1: A000079..1..1...2...4....8....16.....32.....64.....128......256......512
.2: A001333..1..1...3...7...17....41.....99....239.....577.....1393.....3363
.3: A026150..1..1...4..10...28....76....208....568....1552.....4240....11584
.4: A046717..1..1...5..13...41...121....365...1093....3281.....9841....29525
.5: A084057..1..1...6..16...56...176....576...1856....6016....19456....62976
.6: A002533..1..1...7..19...73...241....847...2899...10033....34561...119287
.7: A083098..1..1...8..22...92...316...1184...4264...15632....56848...207488
.8: A084058..1..1...9..25..113...401...1593...5993...23137....88225...338409
.9: A003665..1..1..10..28..136...496...2080...8128...32896...130816...524800
10: A002535..1..1..11..31..161...601...2651..10711...45281...186961...781451
11: A133294..1..1..12..34..188...716...3312..13784...60688...259216..1125312
12: A090042..1..1..13..37..217...841...4069..17389...79537...350353..1575613
13: A125816..1..1..14..40..248...976...4928..21568..102272...463360..2153984
14: A133343..1..1..15..43..281..1121...5895..26363..129361...601441..2884575
15: A133345..1..1..16..46..316..1276...6976..31816..161296...768016..3794176
16: A120612..1..1..17..49..353..1441...8177..37969..198593...966721..4912337
17: A133356..1..1..18..52..392..1616...9504..44864..241792..1201408..6271488
18: A125818..1..1..19..55..433..1801..10963..52543..291457..1476145..7907059
25: A083578
- _Robert G. Wilson v_, Jan 02 2013
Antidiagonal triangle, T(n,k), begins:
  1;
  1,  1;
  1,  1,  1;
  1,  1,  2,  1;
  1,  1,  3,  4,  1;
  1,  1,  4,  7,  8,   1;
  1,  1,  5, 10, 17,  16,   1;
  1,  1,  6, 13, 28,  41,  32,    1;
  1,  1,  7, 16, 41,  76,  99,   64,    1;
  1,  1,  8, 19, 56, 121, 208,  239,  128,    1;
  1,  1,  9, 22, 73, 176, 365,  568,  577,  256,   1;
  1,  1, 10, 25, 92, 241, 576, 1093, 1552, 1393, 512,  1;

G. C. Greubel, Antidiagonals n = 0..50, flattened

Cf  A140895, A221131.
Rows: A000012, A000079, A001333, A026150, A046717, A084057, A002533, A083098, A084058, A003665,
Rows: A002535, A133294, A090042, A125816, A133343, A133345, A120612, A133356, A125818, A083578.
Columns: A000012, A000012, A000027, A016777, A028884, A134593.
Rows include A011782, A001333, A026150, A046717, A002533.

function A084097(n,k)
  if k eq 0 then return 1;
  else return k*2^(k-1)*(&+[ Binomial(k-j,j)*((n-k-1)/4)^j/(k-j): j in [0..Floor(k/2)]]);
  end if; return A084097; end function;
[A084097(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Oct 15 2022

T[j_, k_] := Expand[((1 + Sqrt[j])^k + (1 - Sqrt[j])^k)/2]; T[1, 0] = 1; Table[ T[j - k, k], {j, 0, 11}, {k, 0, j}] // Flatten (* Robert G. Wilson v, Jan 02 2013 *)

def A084097(n,k):
    if (k==0): return 1
    else: return k*2^(k-1)*sum( binomial(k-j,j)*((n-k-1)/4)^j/(k-j) for j in range( (k+2)//2 ) )
flatten([[A084097(n,k) for k in range(n+1)] for n in range(15)]) # G. C. Greubel, Oct 15 2022

From Robert G. Wilson v, Jan 02 2013: (Start)
A(n, k) = (1/2)*( (1 + sqrt(n))^k + (1 - sqrt(n))^k ) (array).
T(n, k) = A(n-k, k). (End)
T(n, k) = Sum_{j=0..floor(k/2)} binomial(k-j, j)*((n-k-1)/4)^j/(k-j), with T(n, 0) = 1 (antidiagonal triangle T(n,k)). - G. C. Greubel, Oct 15 2022

Edited by N. J. A. Sloane, Jul 14 2010

A084150 A Pell related sequence.

Original entry on oeis.org

0, 0, 1, 3, 14, 50, 199, 749, 2892, 11028, 42301, 161799, 619706, 2372006, 9081955, 34767953, 133109592, 509594856, 1950956857, 7469077707, 28594853414, 109473250778, 419110475455, 1604533706357, 6142840740900, 23517417426300

Offset: 0

Paul Barry, May 16 2003

Binomial transform of expansion of (sinh(sqrt(2)x))^2/4 = (0, 0, 1, 0, 8, 0, 64, ...). Inverse binomial transform of A006668.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (3,5,-7).

Cf. A006646, A006668, A015519, A084058.

[n le 3 select Floor((n-1)/2) else 3*Self(n-1) +5*Self(n-2) -7*Self(n-3): n in [1..41]]; // G. C. Greubel, Oct 11 2022

LinearRecurrence[{3,5,-7}, {0,0,1}, 41] (* G. C. Greubel, Oct 11 2022 *)

A084058 = BinaryRecurrenceSequence(2,7,1,1)
def A084150(n): return (A084058(n) - 1)/8
[A084150(n) for n in range(41)] # G. C. Greubel, Oct 11 2022

a(n) = ( (1+sqrt(8))^n + (1-sqrt(8))^n - 2 )/16.
E.g.f.: (1/4)*exp(x)*( sinh(sqrt(2)*x) )^2.
G.f.: x^2 / ( (1-x)*(1-2*x-7*x^2) ). - R. J. Mathar, Feb 05 2011
a(n) = (A015519(n) - A015519(n-1) - 1)/8 = (A084058(n) - 1)/8. - G. C. Greubel, Oct 11 2022

Showing 1-10 of 11 results. Next

cp's OEIS Frontend

A111013 Prime numbers in A084058.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A001541 a(0) = 1, a(1) = 3; for n > 1, a(n) = 6*a(n-1) - a(n-2).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

PARI

Formula

A098158 Triangle T(n,k) with diagonals T(n,n-k) = binomial(n, 2*k).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

GAP

Magma

Mathematica

PARI

PARI

Sage

Formula

A083100 a(n) = 2*a(n-1) + 7*a(n-2).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Formula

A084128 a(n) = 4*a(n-1) + 4*a(n-2), a(0)=1, a(1)=2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Sage

Formula

A124182 A skewed version of triangular array A081277.

Views

Author

Keywords

A083100 a(n) = 2a(n-1) + 7a(n-2).

A084128 a(n) = 4a(n-1) + 4a(n-2), a(0)=1, a(1)=2.

A201701 Riordan triangle ((1-x)/(1-2x), x^2/(1-2x)).

A084097 Square array whose rows have e.g.f. exp(x)cosh(sqrt(k)x), k>=0, read by ascending antidiagonals.