login
Search: a323182 -id:a323182
     Sort: relevance | references | number | modified | created      Format: long | short | data
Antidiagonal sums of A323182.
+20
2
1, 1, 2, 8, 27, 105, 492, 2584, 14893, 93625, 637342, 4663856, 36455959, 302825585, 2661650680, 24662914640, 240141823417, 2450053360913, 26125165902810, 290487741343352, 3361177509359859, 40396577112990745, 503447944487902244, 6496090993661295784, 86660903426459105701
OFFSET
0,3
EXAMPLE
a(1) = 0 + 1 = 1.
a(2) = -1 + 2 + 1 = 2.
a(3) = 0 + 3 + 4 + 1 = 8.
MATHEMATICA
Table[Sum[ChebyshevU[k, n - k], {k, 0, n}], {n, 0, 30}] (* Vaclav Kotesovec, Jan 20 2019 *)
PROG
(PARI) {a(n) = sum(k=0, n, polchebyshev(k, 2, n-k))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jan 20 2019
STATUS
approved
a(n)^2 is a triangular number: a(n) = 6*a(n-1) - a(n-2) with a(0)=0, a(1)=1.
(Formerly M4217 N1760)
+10
193
0, 1, 6, 35, 204, 1189, 6930, 40391, 235416, 1372105, 7997214, 46611179, 271669860, 1583407981, 9228778026, 53789260175, 313506783024, 1827251437969, 10650001844790, 62072759630771, 361786555939836, 2108646576008245, 12290092900109634, 71631910824649559, 417501372047787720
OFFSET
0,3
COMMENTS
8*a(n)^2 + 1 = 8*A001110(n) + 1 = A055792(n+1) is a perfect square. - Gregory V. Richardson, Oct 05 2002
For n >= 2, A001108(n) gives exactly the positive integers m such that 1,2,...,m has a perfect median. The sequence of associated perfect medians is the present sequence. Let a_1,...,a_m be an (ordered) sequence of real numbers, then a term a_k is a perfect median if Sum_{j=1..k-1} a_j = Sum_{j=k+1..m} a_j. See Puzzle 1 in MSRI Emissary, Fall 2005. - Asher Auel, Jan 12 2006
(a(n), b(n)) where b(n) = A082291(n) are the integer solutions of the equation 2*binomial(b,a) = binomial(b+2,a). - Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de); comment revised by Michael Somos, Apr 07 2003
This sequence gives the values of y in solutions of the Diophantine equation x^2 - 8y^2 = 1. It also gives the values of the product xy where (x,y) satisfies x^2 - 2y^2 = +-1, i.e., a(n) = A001333(n)*A000129(n). a(n) also gives the inradius r of primitive Pythagorean triangles having legs whose lengths are consecutive integers, with corresponding semiperimeter s = a(n+1) = {A001652(n) + A046090(n) + A001653(n)}/2 and area rs = A029549(n) = 6*A029546(n). - Lekraj Beedassy, Apr 23 2003 [edited by Jon E. Schoenfield, May 04 2014]
n such that 8*n^2 = floor(sqrt(8)*n*ceiling(sqrt(8)*n)). - Benoit Cloitre, May 10 2003
For n > 0, ratios a(n+1)/a(n) may be obtained as convergents to continued fraction expansion of 3+sqrt(8): either successive convergents of [6;-6] or odd convergents of [5;1, 4]. - Lekraj Beedassy, Sep 09 2003
a(n+1) + A053141(n) = A001108(n+1). Generating floretion: - 2'i + 2'j - 'k + i' + j' - k' + 2'ii' - 'jj' - 2'kk' + 'ij' + 'ik' + 'ji' + 'jk' - 2'kj' + 2e ("jes" series). - Creighton Dement, Dec 16 2004
Kekulé numbers for certain benzenoids (see the Cyvin-Gutman reference). - Emeric Deutsch, Jun 19 2005
Number of D steps on the line y=x in all Delannoy paths of length n (a Delannoy path of length n is a path from (0,0) to (n,n), consisting of steps E=(1,0), N=(0,1) and D=(1,1)). Example: a(2)=6 because in the 13 (=A001850(2)) Delannoy paths of length 2, namely (DD), (D)NE, (D)EN, NE(D), NENE, NEEN, NDE, NNEE, EN(D), ENNE, ENEN, EDN and EENN, we have altogether six D steps on the line y=x (shown between parentheses). - Emeric Deutsch, Jul 07 2005
Define a T-circle to be a first-quadrant circle with integral radius that is tangent to the x- and y-axes. Such a circle has coordinates equal to its radius. Let C(0) be the T-circle with radius 1. Then for n > 0, define C(n) to be the smallest T-circle that does not intersect C(n-1). C(n) has radius a(n+1). Cf. A001653. - Charlie Marion, Sep 14 2005
Numbers such that there is an m with t(n+m)=2t(m), where t(n) are the triangular numbers A000217. For instance, t(20)=2*t(14)=210, so 6 is in the sequence. - Floor van Lamoen, Oct 13 2005
One half the bisection of the Pell numbers (A000129). - Franklin T. Adams-Watters, Jan 08 2006
Pell trapezoids: for n > 0, a(n) = (A000129(n-1)+A000129(n+1))*A000129(n)/2; see also A084158. - Charlie Marion, Apr 01 2006
Tested for 2 < p < 27: If and only if 2^p - 1 (the Mersenne number M(p)) is prime then M(p) divides a(2^(p-1)). - Kenneth J Ramsey, May 16 2006
If 2^p - 1 is prime then M(p) divides a(2^(p-1)-1). - Kenneth J Ramsey, Jun 08 2006; comment corrected by Robert Israel, Mar 18 2007
If 8*n+5 and 8*n+7 are twin primes then their product divides a(4*n+3). - Kenneth J Ramsey, Jun 08 2006
If p is an odd prime, then if p == 1 or 7 (mod 8), then a((p-1)/2) == 0 (mod p) and a((p+1)/2) == 1 (mod p); if p == 3 or 5 (mod 8), then a((p-1)/2) == 1 (mod p) and a((p+1)/2) == 0 (mod p). Kenneth J Ramsey's comment about twin primes follows from this. - Robert Israel, Mar 18 2007
a(n)*(a(n+b) - a(b-2)) = (a(n+1)+1)*(a(n+b-1) - a(b-1)). This identity also applies to any series a(0) = 0 a(1) = 1 a(n) = b*a(n-1) - a(n-2). - Kenneth J Ramsey, Oct 17 2007
For n < 0, let a(n) = -a(-n). Then (a(n+j) + a(k+j)) * (a(n+b+k+j) - a(b-j-2)) = (a(n+j+1) + a(k+j+1)) * (a(n+b+k+j-1) - a(b-j-1)). - Charlie Marion, Mar 04 2011
Sequence gives y values of the Diophantine equation: 0+1+2+...+x = y^2. If (a,b) and (c,d) are two consecutive solutions of the Diophantine equation: 0+1+2+...+x = y^2 with a<c then a+b = c-d and ((d+b)^2, d^2-b^2) is a solution too. If (a,b), (c,d) and (e,f) are three consecutive solutions of the Diophantine equation 0+1+2+...+x = y^2 with a < c < e then (8*d^2, d*(f-b)) is a solution too. - Mohamed Bouhamida, Aug 29 2009
If (p,q) and (r,s) are two consecutive solutions of the Diophantine equation: 0+1+2+...+x = y^2 with p < r then r = 3*p+4*q+1 and s = 2*p+3*q+1. - Mohamed Bouhamida, Sep 02 2009
a(n)/A002315(n) converges to cos^2(Pi/8) (see A201488). - Gary Detlefs, Nov 25 2009
Binomial transform of A086347. - Johannes W. Meijer, Aug 01 2010
If x=a(n), y=A055997(n+1) and z = x^2+y, then x^4 + y^3 = z^2. - Bruno Berselli, Aug 24 2010
In general, if b(0)=1, b(1)=k and for n > 1, b(n) = 6*b(n-1) - b(n-2), then
for n > 0, b(n) = a(n)*k-a(n-1); e.g.,
for k=2, when b(n) = A038725(n), 2 = 1*2 - 0, 11 = 6*2 - 1, 64 = 35*2 - 6, 373 = 204*2 - 35;
for k=3, when b(n) = A001541(n), 3 = 1*3 - 0, 17 = 6*3 - 1; 99 = 35*3 - 6; 577 = 204*3 - 35;
for k=4, when b(n) = A038723(n), 4 = 1*4 - 0, 23 = 6*4 - 1; 134 = 35*4 - 6; 781 = 204*4 - 35;
for k=5, when b(n) = A001653(n), 5 = 1*5 - 0, 29 = 6*5 - 1; 169 = 35*5 - 6; 985 = 204*5 - 35.
- Charlie Marion, Dec 08 2010
See a Wolfdieter Lang comment on A001653 on a sequence of (u,v) values for Pythagorean triples (x,y,z) with x=|u^2-v^2|, y=2*u*v and z=u^2+v^2, with u odd and v even, generated from (u(0)=1,v(0)=2), the triple (3,4,5), by a substitution rule given there. The present a(n) appears there as b(n). The corresponding generated triangles have catheti differing by one length unit. - Wolfdieter Lang, Mar 06 2012
a(n)*a(n+2k) + a(k)^2 and a(n)*a(n+2k+1) + a(k)*a(k+1) are triangular numbers. Generalizes description of sequence. - Charlie Marion, Dec 03 2012
a(n)*a(n+2k) + a(k)^2 is the triangular square A001110(n+k). a(n)*a(n+2k+1) + a(k)*a(k+1) is the triangular oblong A029549(n+k). - Charlie Marion, Dec 05 2012
From Richard R. Forberg, Aug 30 2013: (Start)
The squares of a(n) are the result of applying triangular arithmetic to the squares, using A001333 as the "guide" on what integers to square, as follows:
a(2n)^2 = A001333(2n)^2 * (A001333(2n)^2 - 1)/2;
a(2n+1)^2 = A001333(2n+1)^2 * (A001333(2n+1)^2 + 1)/2. (End)
For n >= 1, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,...,5}. - Milan Janjic, Jan 25 2015
Panda and Rout call these "balancing numbers" and note that the period of the sequence modulo a prime p is the same as that modulo p^2 when p = 13, 31, 1546463. But these are precisely the p in A238736 such that p^2 divides A000129(p - (2/p)), where (2/p) is a Jacobi symbol. In light of the above observation by Franklin T. Adams-Watters that the present sequence is one half the bisection of the Pell numbers, i.e., a(n) = A000129(2*n)/2, it follows immediately that modulo a fixed prime p, or any power thereof, the period of a(n) is half that of A000129(n). - John Blythe Dobson, Mar 06 2015
The triangular number = square number identity Tri((T(n, 3) - 1)/2) = S(n-1, 6)^2 with Tri, T, and S given in A000217, A053120 and A049310, is the special case k = 1 of the k-family of identities Tri((T(n, 2*k+1) - 1)/2) = Tri(k)*S(n-1, 2*(2*k+1))^2, k >= 0, n >= 0, with S(-1, x) = 0. For k=2 see A108741(n) for S(n-1, 10)^2. This identity boils down to the identities S(n-1, 2*x)^2 = (T(2*n, x) - 1)/(2*(x^2-1)) and 2*T(n, x)^2 - 1 = T(2*n, x) with x = 2*k+1. - Wolfdieter Lang, Feb 01 2016
a(2)=6 is perfect. For n=2*k, k > 0, k not equal to 1, a(n) is a multiple of a(2) and since every multiple (beyond 1) of a perfect number is abundant, then a(n) is abundant. sigma(a(4)) = 504 > 408 = 2*a(4). For n=2*k+1, k > 0, a(n) mod 10 = A000012(n), so a(n) is odd. If a(n) is a prime number, it is deficient; otherwise a(n) has one or two distinct prime factors and is therefore deficient again. So for n=2k+1, k > 0, a(n) is deficient. sigma(a(5)) = 1260 < 2378 = 2*a(5). - Muniru A Asiru, Apr 14 2016
Behera & Panda call these the balancing numbers, and A001541 are the balancers. - Michel Marcus, Nov 07 2017
In general, a second-order linear recurrence with constant coefficients having a signature of (c,d) will be duplicated by a third-order recurrence having a signature of (x,c^2-c*x+d,-d*x+c*d). The formulas of Olivares and Bouhamida in the formula section which have signatures of (7,-7,1) and (5,5,-1), respectively, are specific instances of this general rule for x = 7 and x = 5. - Gary Detlefs, Jan 29 2021
Note that 6 is the largest triangular number in the sequence, because it is proved that 8 and 9 are the largest perfect powers which are consecutive (Catalan's conjecture). 0 and 1 are also in the sequence because they are also perfect powers and 0*1/2 = 0^2 and 8*9/2 = (2*3)^2. - Metin Sariyar, Jul 15 2021
REFERENCES
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) - From N. J. A. Sloane, May 30 2012
A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, pp. 193, 197.
D. M. Burton, The History of Mathematics, McGraw Hill, (1991), p. 213.
L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 10.
P. Franklin, E. F. Beckenbach, H. S. M Coxeter, N. H. McCoy, K. Menger, and J. L. Synge, Rings And Ideals, No 8, The Carus Mathematical Monographs, The Mathematical Association of America, (1967), pp. 144-146.
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 B(n).
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).
P.-F. Teilhet, Query 2376, L'Intermédiaire des Mathématiciens, 11 (1904), 138-139. - N. J. A. Sloane, Mar 08 2022
LINKS
Indranil Ghosh, Table of n, a(n) for n = 0..1304 (terms 0..200 from T. D. Noe)
Marco Abrate, Stefano Barbero, Umberto Cerruti, and Nadir Murru, Polynomial sequences on quadratic curves, Integers, Vol. 15, 2015, #A38.
Irving Adler, Three Diophantine equations - Part II, Fib. Quart., 7 (1969), pp. 181-193.
Seyed Hassan Alavi, Ashraf Daneshkhah, and Cheryl E. Praeger, Symmetries of biplanes, arXiv:2004.04535 [math.GR], 2020. See v_n in Lemma 7.9 p. 21.
Jean-Paul Allouche, Zeta-regularization of arithmetic sequences, EPJ Web of Conferences (2020) Vol. 244, 01008.
Dario Alpern for Diophantine equation a^4+b^3=c^2.
Kasper Andersen, Lisa Carbone, and D. Penta, Kac-Moody Fibonacci sequences, hyperbolic golden ratios, and real quadratic fields, Journal of Number Theory and Combinatorics, Vol 2, No. 3 pp 245-278, 2011. See Section 9.
Francesca Arici and Jens Kaad, Gysin sequences and SU(2)-symmetries of C*-algebras, arXiv:2012.11186 [math.OA], 2020.
Muniru A. Asiru, All square chiliagonal numbers, International Journal of Mathematical Education in Science and Technology, Volume 47, 2016 - Issue 7.
Jeremiah Bartz, Bruce Dearden, and Joel Iiams, Classes of Gap Balancing Numbers, arXiv:1810.07895 [math.NT], 2018.
Jeremiah Bartz, Bruce Dearden, Joel Iiams, and Julia Peterson, Powers of Two as Sums of Two Balancing Numbers, Combinatorics, Graph Theory and Computing (SEICCGTC 2021) Springer Proc. Math. Stat., Vol 448, pp. 383-392. See p. 384.
Raymond A. Beauregard and Vladimir A. Dobrushkin, Powers of a Class of Generating Functions, Mathematics Magazine, Volume 89, Number 5, December 2016, pp. 359-363.
A. Behera and G. K. Panda, On the Square Roots of Triangular Numbers, Fib. Quart., 37 (1999), pp. 98-105.
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.
Elwyn Berlekamp and Joe P. Buhler, Puzzle Column, Emissary, MSRI Newsletter, Fall 2005. Problem 1, (6 MB).
Kisan Bhoi and Prasanta Kumar Ray, On the Diophantine equation Bn1+Bn2=2^a1+2^a2+2^a3, arXiv:2212.06372 [math.NT], 2022.
Daniel Birmajer, Juan B. Gil, and Michael D. Weiner, On the Enumeration of Restricted Words over a Finite Alphabet, J. Int. Seq. 19 (2016) # 16.1.3, example 12.
Paula Catarino, Helena Campos, and Paulo Vasco, On some identities for balancing and cobalancing numbers, Annales Mathematicae et Informaticae, 45 (2015) pp. 11-24.
E. K. Çetinalp, N. Yilmaz, and Ö. Deveci, The balancing-like sequences in groups, Acta Univ. Apulensis Math. (2023) No. 73, 139-153. See p. 144.
S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (pp. 301, 302, P_{13}).
Mahadi Ddamulira, Repdigits as sums of three balancing numbers, Mathematica Slovaca, (2019) hal-02405969.
Tomislav Doslic, Planar polycyclic graphs and their Tutte polynomials, Journal of Mathematical Chemistry, Volume 51, Issue 6, 2013, pp. 1599-1607.
D. B. Eperson, Triangular numbers, Math. Gaz., 47 (1963), 236-237.
Sergio Falcon, Relationships between Some k-Fibonacci Sequences, Applied Mathematics, 2014, 5, 2226-2234.
Bernadette Faye, Florian Luca, and Pieter Moree, On the discriminator of Lucas sequences, arXiv:1708.03563 [math.NT], 2017.
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.
Rigoberto Flórez, Robinson A. Higuita, and Antara Mukherjee, Alternating Sums in the Hosoya Polynomial Triangle, Article 14.9.5 Journal of Integer Sequences, Vol. 17 (2014).
Aviezri S. Fraenkel, On the recurrence f(m+1)= b(m)*f(m)-f(m-1) and applications, Discrete Mathematics 224 (2000), pp. 273-279.
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 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.
Robert Frontczak and Kalika Prasad, Balancing polynomials, Fibonacci numbers and some new series for $\pi$, Mediterranean Journal of Mathematics (2023) Vol. 20, Article number: 207.
Bill Gosper, The Triangular Squares, 2014.
H. Harborth, Fermat-like binomial equations, Applications of Fibonacci numbers, Proc. 2nd Int. Conf., San Jose/Ca., August 1986, 1-5 (1988).
Brian Hayes, Calculemus!, American Scientist, 96 (Sep-Oct 2008), 362-366.
Milan Janjic, On Linear Recurrence Equations Arising from Compositions of Positive Integers, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.7.
Michael A. Jones, Proof Without Words: The Square of a Balancing Number Is a Triangular Number, The College Mathematics Journal, Vol. 43, No. 3 (May 2012), p. 212.
Refik Keskin and Olcay Karaatli, Some New Properties of Balancing Numbers and Square Triangular Numbers, Journal of Integer Sequences, Vol. 15 (2012), Article #12.1.4.
Omar Khadir, Kalman Liptai, and Laszlo Szalay, On the Shifted Product of Binary Recurrences, J. Int. Seq. 13 (2010), 10.6.1.
Tanya Khovanova, Recursive Sequences
Phil Lafer, Discovering the square-triangular numbers, Fib. Quart., 9 (1971), 93-105.
Ioana-Claudia Lazăr, Lucas sequences in t-uniform simplicial complexes, arXiv:1904.06555 [math.GR], 2019.
Kalman Liptai, Fibonacci Balancing Numbers, Fib. Quart. 42 (4) (2004) 330-340.
Madras College, St Andrews, Square Triangular Numbers
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.
Roger B. Nelson, Multi-Polygonal Numbers, Mathematics Magazine, Vol. 89, No. 3 (June 2016), pp. 159-164.
G. K. Panda, Sequence balancing and cobalancing numbers, Fib. Q., Vol. 45, No. 3 (2007), 265-271. See p. 266.
G. K. Panda and S. S. Rout, Periodicity of Balancing Numbers, Acta Mathematica Hungarica 143 (2014), 274-286.
G. K. Panda and Ravi Kumar Davala, Perfect Balancing Numbers, Fibonacci Quart. 53 (2015), no. 3, 261-264.
Ashish Kumar Pandey and B. K. Sharma, On Inequalities Related to a Generalized Euler Totient Function and Lucas Sequences, J. Int. Seq. (2023) Vol. 26, Art. 23.8.6.
Poo-Sung Park, Ramanujan's Continued Fraction for a Puzzle, College Mathematics Journal, 2005, 363-365.
Michael Penn, Balancing Numbers, Youtube video, 2020.
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
B. Polster and M. Ross, Marching in squares, arXiv preprint arXiv:1503.04658 [math.HO], 2015.
Kalika Prasad, Munesh Kumari, Ritanjali Mohanty, and Hrishikesh Mahato, Spinor Algebra of k-Balancing and k-Lucas-Balancing Numbers, Journal of Algebra and Its Applications (2024).
Kalika Prasad, Munesh Kumari, and Jagmohan Tanti, Octonions and hyperbolic octonions with the k-balancing and k-Lucas balancing numbers, The Journal of Analysis, (2024), Vol. 32, No. 3, 1281-1296.
Kalika Prasad, Munesh Kumari, and Jagmohan Tanti, Generalized k-balancing and k-Lucas balancing numbers and associated polynomials, Kyungpook Mathematical Journal (2023), Vol. 63, No. 4, 539-550.
Kalika Prasad, Hrishikesh Mahato, and Munesh Kumari, Some properties of r-circulant matrices with k-balancing and k-Lucas balancing numbers, Boletín de la Sociedad Matemática Mexicana (2023), Vol. 29, No. 2, Article no. 44.
Helmut Prodinger, How to sum powers of balancing numbers efficiently, arXiv:2008.03916 [math.NT], 2020.
K. J. Ramsey, Relation of Mersenne Primes To Square Triangular Numbers [edited by K. J. Ramsey, May 14 2011]
Kenneth Ramsay and Andras Erszegi, Relation of Square Triangular Numbers To Mersenne Primes, digest of 4 messages in Triangular_and_Fibonacci_Numbers Yahoo Group, May 15 - Jun 28, 2006.
Kenneth Ramsey, Generalized Proof re Square Triangular Numbers, digest of 2 messages in Triangular_and_Fibonacci_Numbers Yahoo group, May 27, 2005 - Oct 10, 2011.
Salah E. Rihane, Bernadette Faye, Florian Luca, and Alain Togbe, An exponential Diophantine equation related to the difference between powers of two consecutive Balancing numbers, arXiv:1811.03015 [math.NT], 2018.
Sci.math Newsgroup, Square numbers which are triangular [Cached copy]
R. A. Sulanke, Bijective recurrences concerning Schroeder paths, Electron. J. Combin. 5 (1998), Research Paper 47, 11 pp.
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.
Ahmet Tekcan, Merve Tayat, and Meltem E. Ozbek, The diophantine equation 8x^2-y^2+8x(1+t)+(2t+1)^2=0 and t-balancing numbers, ISRN Combinatorics, Volume 2014, Article ID 897834, 5 pages.
Eric Weisstein's World of Mathematics, Binomial coefficient.
Eric Weisstein's World of Mathematics, Square Triangular Number.
Eric Weisstein's World of Mathematics, Triangular Number.
FORMULA
G.f.: x / (1 - 6*x + x^2). - Simon Plouffe in his 1992 dissertation.
a(n) = S(n-1, 6) = U(n-1, 3) with U(n, x) Chebyshev's polynomials of the second kind. S(-1, x) := 0. Cf. triangle A049310 for S(n, x).
a(n) = sqrt(A001110(n)).
a(n) = A001542(n)/2.
a(n) = sqrt((A001541(n)^2-1)/8) (cf. Richardson comment).
a(n) = 3*a(n-1) + sqrt(8*a(n-1)^2+1). - R. J. Mathar, Oct 09 2000
a(n) = A000129(n)*A001333(n) = A000129(n)*(A000129(n)+A000129(n-1)) = ceiling(A001108(n)/sqrt(2)). - Henry Bottomley, Apr 19 2000
a(n) ~ (1/8)*sqrt(2)*(sqrt(2) + 1)^(2*n). - Joe Keane (jgk(AT)jgk.org), May 15 2002
Limit_{n->infinity} a(n)/a(n-1) = 3 + 2*sqrt(2). - Gregory V. Richardson, Oct 05 2002
a(n) = ((3 + 2*sqrt(2))^n - (3 - 2*sqrt(2))^n) / (4*sqrt(2)). - Gregory V. Richardson, Oct 13 2002. Corrected for offset 0, and rewritten. - Wolfdieter Lang, Feb 10 2015
a(2*n) = a(n)*A003499(n). 4*a(n) = A005319(n). - Mario Catalani (mario.catalani(AT)unito.it), Mar 21 2003
a(n) = floor((3+2*sqrt(2))^n/(4*sqrt(2))). - Lekraj Beedassy, Apr 23 2003
a(-n) = -a(n). - Michael Somos, Apr 07 2003
For n >= 1, a(n) = Sum_{k=0..n-1} A001653(k). - Charlie Marion, Jul 01 2003
For n > 0, 4*a(2*n) = A001653(n)^2 - A001653(n-1)^2. - Charlie Marion, Jul 16 2003
For n > 0, a(n) = Sum_{k = 0..n-1}((2*k+1)*A001652(n-1-k)) + A000217(n). - Charlie Marion, Jul 18 2003
a(2*n+1) = a(n+1)^2 - a(n)^2. - Charlie Marion, Jan 12 2004
a(k)*a(2*n+k) = a(n+k)^2 - a(n)^2; e.g., 204*7997214 = 40391^2 - 35^2. - Charlie Marion, Jan 15 2004
For j < n+1, a(k+j)*a(2*n+k-j) - Sum_{i = 0..j-1} a(2*n-(2*i+1)) = a(n+k)^2 - a(n)^2. - Charlie Marion, Jan 18 2004
From Paul Barry, Feb 06 2004: (Start)
a(n) = A000129(2*n)/2;
a(n) = ((1+sqrt(2))^(2*n) - (1-sqrt(2))^(2*n))*sqrt(2)/8;
a(n) = Sum_{i=0..n} Sum_{j=0..n} A000129(i+j)*n!/(i!*j!*(n-i-j)!)/2. (End)
E.g.f.: exp(3*x)*sinh(2*sqrt(2)*x)/(2*sqrt(2)). - Paul Barry, Apr 21 2004
A053141(n+1) + A055997(n+1) = A001541(n+1) + a(n+1). - Creighton Dement, Sep 16 2004
a(n) = Sum_{k=0..n} binomial(2*n, 2*k+1)*2^(k-1). - Paul Barry, Oct 01 2004
a(n) = A001653(n+1) - A038723(n); (a(n)) = chuseq[J]( 'ii' + 'jj' + .5'kk' + 'ij' - 'ji' + 2.5e ), apart from initial term. - Creighton Dement, Nov 19 2004, modified by Davide Colazingari, Jun 24 2016
a(n+1) = Sum_{k=0..n} A001850(k)*A001850(n-k), self convolution of central Delannoy numbers. - Benoit Cloitre, Sep 28 2005
a(n) = 7*(a(n-1) - a(n-2)) + a(n-3), a(1) = 0, a(2) = 1, a(3) = 6, n > 3. Also a(n) = ( (1 + sqrt(2) )^(2*n) - (1 - sqrt(2) )^(2*n) ) / (4*sqrt(2)). - Antonio Alberto Olivares, Oct 23 2003
a(n) = 5*(a(n-1) + a(n-2)) - a(n-3). - Mohamed Bouhamida, Sep 20 2006
Define f(x,s) = s*x + sqrt((s^2-1)*x^2+1); f(0,s)=0. a(n) = f(a(n-1),3), see second formula. - Marcos Carreira, Dec 27 2006
The perfect median m(n) can be expressed in terms of the Pell numbers P() = A000129() by m(n) = P(n + 2) * (P(n + 2) + (P(n + 1)) for n >= 0. - Winston A. Richards (ugu(AT)psu.edu), Jun 11 2007
For k = 0..n, a(2*n-k) - a(k) = 2*a(n-k)*A001541(n). Also, a(2*n+1-k) - a(k) = A002315(n-k)*A001653(n). - Charlie Marion, Jul 18 2007
[A001653(n), a(n)] = [1,4; 1,5]^n * [1,0]. - Gary W. Adamson, Mar 21 2008
a(n) = Sum_{k=0..n-1} 4^k*binomial(n+k,2*k+1). - Paul Barry, Apr 20 2009
a(n+1)^2 - 6*a(n+1)*a(n) + a(n)^2 = 1. - Charlie Marion, Dec 14 2010
a(n) = A002315(m)*A011900(n-m-1) + A001653(m)*A001652(n-m-1) - a(m) = A002315(m)*A053141(n-m-1) + A001653(m)*A046090(n-m-1) + a(m) with m < n; otherwise a(n) = A002315(m)*A053141(m-n) - A001653(m)*A011900(m-n) + a(m) = A002315(m)*A053141(m-n) - A001653(m)*A046090(m-n) - a(m) = (A002315(n) - A001653(n))/2. - Kenneth J Ramsey, Oct 12 2011
16*a(n)^2 + 1 = A056771(n). - James R. Buddenhagen, Dec 09 2011
A010054(A000290(a(n))) = 1. - Reinhard Zumkeller, Dec 17 2011
In general, a(n+k)^2 - A003499(k)*a(n+k)*a(n) + a(n)^2 = a(k)^2. - Charlie Marion, Jan 11 2012
a(n+1) = Sum_{k=0..n} A101950(n,k)*5^k. - Philippe Deléham, Feb 10 2012
PSUM transform of a(n+1) is A053142. PSUMSIGN transform of a(n+1) is A084158. BINOMIAL transform of a(n+1) is A164591. BINOMIAL transform of A086347 is a(n+1). BINOMIAL transform of A057087(n-1). - Michael Somos, May 11 2012
a(n+k) = A001541(k)*a(n) + sqrt(A132592(k)*a(n)^2 + a(k)^2). Generalizes formula dated Oct 09 2000. - Charlie Marion, Nov 27 2012
a(n) + a(n+2*k) = A003499(k)*a(n+k); a(n) + a(n+2*k+1) = A001653(k+1)*A002315(n+k). - Charlie Marion, Nov 29 2012
From Peter Bala, Dec 23 2012: (Start)
Product_{n >= 1} (1 + 1/a(n)) = 1 + sqrt(2).
Product_{n >= 2} (1 - 1/a(n)) = (1/3)*(1 + sqrt(2)). (End)
G.f.: G(0)*x/(2-6*x), where G(k) = 1 + 1/(1 - x*(8*k-9)/( x*(8*k-1) - 3/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Aug 12 2013
G.f.: H(0)*x/2, where H(k) = 1 + 1/( 1 - x*(6-x)/(x*(6-x) + 1/H(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Feb 18 2014
a(n) = (a(n-1)^2 - a(n-3)^2)/a(n-2) + a(n-4) for n > 3. - Patrick J. McNab, Jul 24 2015
a(n-k)*a(n+k) + a(k)^2 = a(n)^2, a(n+k) + a(n-k) = A003499(k)*a(n), for n >= k >= 0. - Alexander Samokrutov, Sep 30 2015
Dirichlet g.f.: (PolyLog(s,3+2*sqrt(2)) - PolyLog(s,3-2*sqrt(2)))/(4*sqrt(2)). - Ilya Gutkovskiy, Jun 27 2016
4*a(n)^2 - 1 = A278310(n) for n > 0. - Bruno Berselli, Nov 24 2016
From Klaus Purath, Jan 18 2020: (Start)
a(n) = (a(n-3) + a(n+3))/198.
a(n) = Sum_{i=1..n} A001653(i), n>=1.
a(n) = sinh( 2 * n * arccsch(1) ) / ( 2 * sqrt(2) ). - Federico Provvedi, Feb 01 2021
(End)
a(n) = A002965(2*n)*A002965(2*n+1). - Jon E. Schoenfield, Jan 08 2022
a(n) = A002965(4*n)/2. - Gerry Martens, Jul 14 2023
a(n) = Sum_{k = 0..n-1} (-1)^(n+k+1)*binomial(n+k, 2*k+1)*8^k. - Peter Bala, Jul 17 2023
EXAMPLE
G.f. = x + 6*x^2 + 35*x^3 + 204*x^4 + 1189*x^5 + 6930*x^6 + 40391*x^7 + ...
6 is in the sequence since 6^2 = 36 is a triangular number: 36 = 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8. - Michael B. Porter, Jul 02 2016
MAPLE
a[0]:=1: a[1]:=6: for n from 2 to 26 do a[n]:=6*a[n-1]-a[n-2] od: seq(a[n], n=0..26); # Emeric Deutsch
with (combinat):seq(fibonacci(2*n, 2)/2, n=0..20); # Zerinvary Lajos, Apr 20 2008
MATHEMATICA
Transpose[NestList[Flatten[{Rest[#], ListCorrelate[{-1, 6}, #]}]&, {0, 1}, 30]][[1]] (* Harvey P. Dale, Mar 23 2011 *)
CoefficientList[Series[x/(1-6x+x^2), {x, 0, 30}], x] (* Harvey P. Dale, Mar 23 2011 *)
LinearRecurrence[{6, -1}, {0, 1}, 50] (* Vladimir Joseph Stephan Orlovsky, Feb 12 2012 *)
a[ n_]:= ChebyshevU[n-1, 3]; (* Michael Somos, Sep 02 2012 *)
Table[Fibonacci[2n, 2]/2, {n, 0, 20}] (* Vladimir Reshetnikov, Sep 16 2016 *)
TrigExpand@Table[Sinh[2 n ArcCsch[1]]/(2 Sqrt[2]), {n, 0, 10}] (* Federico Provvedi, Feb 01 2021 *)
PROG
(PARI) {a(n) = imag((3 + quadgen(32))^n)}; /* Michael Somos, Apr 07 2003 */
(PARI) {a(n) = subst( poltchebi( abs(n+1)) - 3 * poltchebi( abs(n)), x, 3) / 8}; /* Michael Somos, Apr 07 2003 */
(PARI) {a(n) = polchebyshev( n-1, 2, 3)}; /* Michael Somos, Sep 02 2012 */
(PARI) is(n)=ispolygonal(n^2, 3) \\ Charles R Greathouse IV, Nov 03 2016
(Sage) [lucas_number1(n, 6, 1) for n in range(27)] # Zerinvary Lajos, Jun 25 2008
(Sage) [chebyshev_U(n-1, 3) for n in (0..20)] # G. C. Greubel, Dec 23 2019
(Haskell)
a001109 n = a001109_list !! n :: Integer
a001109_list = 0 : 1 : zipWith (-)
(map (* 6) $ tail a001109_list) a001109_list
-- Reinhard Zumkeller, Dec 17 2011
(Magma) [n le 2 select n-1 else 6*Self(n-1)-Self(n-2): n in [1..30]]; // Vincenzo Librandi, Jul 25 2015
(GAP) a:=[0, 1];; for n in [3..25] do a[n]:=6*a[n-1]-a[n-2]; od; a; # Muniru A Asiru, Dec 18 2018
CROSSREFS
Chebyshev sequence U(n, m): A000027 (m=1), A001353 (m=2), this sequence (m=3), A001090 (m=4), A004189 (m=5), A004191 (m=6), A007655 (m=7), A077412 (m=8), A049660 (m=9), A075843 (m=10), A077421 (m=11), A077423 (m=12), A097309 (m=13), A097311 (m=14), A097313 (m=15), A029548 (m=16), A029547 (m=17), A144128 (m=18), A078987 (m=19), A097316 (m=33).
Cf. A323182.
KEYWORD
nonn,easy,nice
EXTENSIONS
Additional comments from Wolfdieter Lang, Feb 10 2000
Duplication of a formula removed by Wolfdieter Lang, Feb 10 2015
STATUS
approved
a(n) = 4*a(n-1) - a(n-2) with a(0) = 0, a(1) = 1.
(Formerly M3499 N1420)
+10
188
0, 1, 4, 15, 56, 209, 780, 2911, 10864, 40545, 151316, 564719, 2107560, 7865521, 29354524, 109552575, 408855776, 1525870529, 5694626340, 21252634831, 79315912984, 296011017105, 1104728155436, 4122901604639, 15386878263120, 57424611447841, 214311567528244
OFFSET
0,3
COMMENTS
3*a(n)^2 + 1 is a square. Moreover, 3*a(n)^2 + 1 = (2*a(n) - a(n-1))^2.
Consecutive terms give nonnegative solutions to x^2 - 4*x*y + y^2 = 1. - Max Alekseyev, Dec 12 2012
Values y solving the Pellian x^2 - 3*y^2 = 1; corresponding x values given by A001075(n). Moreover, we have a(n) = 2*a(n-1) + A001075(n-1). - Lekraj Beedassy, Jul 13 2006
Number of spanning trees in 2 X n grid: by examining what happens at the right-hand end we see that a(n) = 3*a(n-1) + 2*a(n-2) + 2*a(n-3) + ... + 2*a(1) + 1, where the final 1 corresponds to the tree ==...=| !. Solving this we get a(n) = 4*a(n-1) - a(n-2).
Complexity of 2 X n grid.
A016064 also describes triangles whose sides are consecutive integers and in which an inscribed circle has an integer radius. A001353 is exactly and precisely mapped to the integer radii of such inscribed circles, i.e., for each term of A016064, the corresponding term of A001353 gives the radius of the inscribed circle. - Harvey P. Dale, Dec 28 2000
n such that 3*n^2 = floor(sqrt(3)*n*ceiling(sqrt(3)*n)). - Benoit Cloitre, May 10 2003
For n>0, ratios a(n+1)/a(n) may be obtained as convergents of the continued fraction expansion of 2+sqrt(3): either as successive convergents of [4;-4] or as odd convergents of [3;1, 2]. - Lekraj Beedassy, Sep 19 2003
Ways of packing a 3 X (2*n-1) rectangle with dominoes, after attaching an extra square to the end of one of the sides of length 3. With reference to A001835, therefore: a(n) = a(n-1) + A001835(n-1) and A001835(n) = 3*A001835(n-1) + 2*a(n-1). - Joshua Zucker and the Castilleja School Math Club, Oct 28 2003
a(n+1) is a Chebyshev transform of 4^n, where the sequence with g.f. G(x) is sent to the sequence with g.f. (1/(1+x^2))G(x/(1+x^2)). - Paul Barry, Oct 25 2004
This sequence is prime-free, because a(2n) = a(n) * (a(n+1)-a(n-1)) and a(2n+1) = a(n+1)^2 - a(n)^2 = (a(n+1)+a(n)) * (a(n+1)-a(n)). - Jianing Song, Jul 06 2019
Numbers such that there is an m with t(n+m) = 3*t(m), where t(n) are the triangular numbers A000217. For instance, t(35) = 3*t(20) = 630, so 35 - 20 = 15 is in the sequence. - Floor van Lamoen, Oct 13 2005
a(n) = number of distinct matrix products in (A + B + C + D)^n where commutator [A,B] = 0 but neither A nor B commutes with C or D. - Paul D. Hanna and Max Alekseyev, Feb 01 2006
For n > 1, middle side (or long leg) of primitive Pythagorean triangles having an angle nearing Pi/3 with larger values of sides. [Complete triple (X, Y, Z), X < Y < Z, is given by X = A120892(n), Y = a(n), Z = A120893(n), with recurrence relations X(i+1) = 2*{X(i) - (-1)^i} + a(i); Z(i+1) = 2*{Z(i) + a(i)} - (-1)^i.] - Lekraj Beedassy, Jul 13 2006
From Dennis P. Walsh, Oct 04 2006: (Start)
Number of 2 X n simple rectangular mazes. A simple rectangular m X n maze is a graph G with vertex set {0, 1, ..., m} X {0, 1, ..., n} that satisfies the following two properties: (i) G consists of two orthogonal trees; (ii) one tree has a path that sequentially connects (0,0),(0,1), ..., (0,n), (1,n), ...,(m-1,n) and the other tree has a path that sequentially connects (1,0), (2,0), ..., (m,0), (m,1), ..., (m,n). For example, a(2) = 4 because there are four 2 X 2 simple rectangular mazes:
__ __ __ __
| | | |__ | | | | __|
| __| | __| | |__| | __|
(End)
[1, 4, 15, 56, 209, ...] is the Hankel transform of [1, 1, 5, 26, 139, 758, ...](see A005573). - Philippe Deléham, Apr 14 2007
The upper principal convergents to 3^(1/2), beginning with 2/1, 7/4, 26/15, 97/56, comprise a strictly decreasing sequence; numerators=A001075, denominators=A001353. - Clark Kimberling, Aug 27 2008
From Gary W. Adamson, Jun 21 2009: (Start)
A001353 and A001835 = bisection of continued fraction [1, 2, 1, 2, 1, 2, ...], i.e., of [1, 3, 4, 11, 15, 41, ...].
For n>0, a(n) equals the determinant of an (n-1) X (n-1) tridiagonal matrix with ones in the super and subdiagonals and (4, 4, 4, ...) as the main diagonal. [Corrected by Johannes Boot, Sep 04 2011]
A001835 and A001353 = right and next to right borders of triangle A125077. (End)
a(n) is equal to the permanent of the (n-1) X (n-1) Hessenberg matrix with 4's along the main diagonal, i's along the superdiagonal and the subdiagonal (i is the imaginary unit), and 0's everywhere else. - John M. Campbell, Jun 09 2011
2a(n) is the number of n-color compositions of 2n consisting of only even parts; see Guo in references. - Brian Hopkins, Jul 19 2011
Pisano period lengths: 1, 2, 6, 4, 3, 6, 8, 4, 18, 6, 10, 12, 12, 8, 6, 8, 18, 18, 5, 12, ... - R. J. Mathar, Aug 10 2012
From Michel Lagneau, Jul 08 2014: (Start)
a(n) is defined also by the recurrence a(1)=1; for n>1, a(n+1) = 2*a(n) + sqrt(3*a(n)^2 + 1) where a(n) is an integer for every n. This sequence is generalizable by the sequence b(n,m) of parameter m with the initial condition b(1,m) = 1, and for n > 1 b(n+1,m) = m*b(n,m) + sqrt((m^2 - 1)*b(n,m)^2 + 1) for m = 2, 3, 4, ... where b(n,m) is an integer for every n.
The first corresponding sequences are
b(n,2) = a(n) = A001353(n);
b(n,3) = A001109(n);
b(n,4) = A001090(n);
b(n,5) = A004189(n);
b(n,6) = A004191(n);
b(n,7) = A007655(n);
b(n,8) = A077412(n);
b(n,9) = A049660(n);
b(n,10) = A075843(n);
b(n,11) = A077421(n);
....................
We obtain a general sequence of polynomials {b(n,x)} = {1, 2*x, 4*x^2 - 1, 8*x^3 - 4*x, 16*x^4 - 12*x^2 + 1, 32*x^5 - 32*x^3 + 6*x, ...} with x = m where each b(n,x) is a Gegenbauer polynomial defined by the recurrence b(n,x)- 2*x*b(n-1,x) + b(n-2,x) = 0, the same relation as the Chebyshev recurrence, but with the initial conditions b(x,0) = 1 and b(x,1) = 2*x instead b(x,0) = 1 and b(x,1) = x for the Chebyshev polynomials. (End)
If a(n) denotes the n-th term of the above sequence and we construct a triangle whose sides are a(n) - 1, a(n) + 1 and sqrt(3a(n)^2 + 1), then, for every n the measure of one of the angles of the triangle so constructed will always be 120 degrees. This result of ours was published in Mathematics Spectrum (2012/2013), Vol. 45, No. 3, pp. 126-128. - K. S. Bhanu and Dr. M. N. Deshpande, Professor (Retd), Department of Statistics, Institute of Science, Nagpur (India).
For n >= 1, a(n) equals the number of 01-avoiding words of length n - 1 on alphabet {0, 1, 2, 3}. - Milan Janjic, Jan 25 2015
For n > 0, 10*a(n) is the number of vertices and roots on level n of the {4, 5} mosaic (see L. Németh Table 1 p. 6). - Michel Marcus, Oct 30 2015
(2 + sqrt(3))^n = A001075(n) + a(n)*sqrt(3), n >= 0; integers in the quadratic number field Q(sqrt(3)). - Wolfdieter Lang, Feb 16 2018
A strong divisibility sequence, that is, gcd(a(n), a(m)) = a(gcd(n, m)) for all positive integers n and m. - Michael Somos, Dec 12 2019
The Cholesky decomposition A = C C* for tridiagonal A with A[i,i] = 4 and A[i+1,i] = A[i,i+1] = -1, as it arises in the discretized 2D Laplace operator (Poisson equation...), has nonzero elements C[i,i] = sqrt(a(i+1)/a(i)) = -1/C[i+1,i], i = 1, 2, 3, ... - M. F. Hasler, Mar 12 2021
The triples (a(n-1), 2a(n), a(n+1)), n=2,3,..., are exactly the triples (a,b,c) of positive integers a < b < c in arithmetic progression such that a*b+1, b*c+1, and c*a+1 are perfect squares. - Bernd Mulansky, Jul 10 2021
REFERENCES
J. Austin and L. Schneider, Generalized Fibonacci sequences in Pythagorean triple preserving sequences, Fib. Q., 58:1 (2020), 340-350.
Bastida, Julio R., 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)
G. Everest, A. van der Poorten, I. Shparlinski and T. Ward, Recurrence Sequences, Amer. Math. Soc., 2003; p. 163.
F. Faase, On the number of specific spanning subgraphs of the graphs G X P_n, Ars Combin. 49 (1998), 129-154.
R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 329.
J. D. E. Konhauser et al., Which Way Did the Bicycle Go?, MAA 1996, p. 104.
Serge Lang, Introduction to Diophantine Approximations, Addison-Wesley, New York, 1966.
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).
LINKS
Marco Abrate, Stefano Barbero, Umberto Cerruti, and Nadir Murru, Polynomial sequences on quadratic curves, Integers, Vol. 15, 2015, #A38.
Christian Aebi and Grant Cairns, Lattice Equable Parallelograms, arXiv:2006.07566 [math.NT], 2020.
Christian Aebi and Grant Cairns, Equable Parallelograms on the Eisenstein Lattice, arXiv:2401.08827 [math.CO], 2024. See p. 14.
W. K. Alt, Enumeration of Domino Tilings on the Projective Grid Graph, A Thesis Presented to The Division of Mathematics and Natural Sciences, Reed College, May 2013.
K. Andersen, L. Carbone, and D. Penta, Kac-Moody Fibonacci sequences, hyperbolic golden ratios, and real quadratic fields, Journal of Number Theory and Combinatorics, Vol 2, No. 3 pp 245-278, 2011. See Section 9.
Francesca Arici and Jens Kaad, Gysin sequences and SU(2)-symmetries of C*-algebras, arXiv:2012.11186 [math.OA], 2020.
Krassimir T. Atanassov and Anthony G. Shannon, On intercalated Fibonacci sequences, Notes on Number Theory and Discrete Mathematics (2020) Vol. 26, No. 3, 218-223.
Daniel Birmajer, Juan B. Gil, and Michael D. Weiner, Linear recurrence sequences with indices in arithmetic progression and their sums, arXiv preprint arXiv:1505.06339 [math.NT], 2015.
D. Birmajer, J. B. Gil, and M. D. Weiner, On the Enumeration of Restricted Words over a Finite Alphabet, J. Int. Seq. 19 (2016) # 16.1.3, example 12.
K. S. Bhanu and M. N. Deshpande, Integral triangles with 120° angle Mathematics Spectrum, 45 (3) (2012/2013), 126-128.
Latham Boyle and Paul J. Steinhardt, Self-Similar One-Dimensional Quasilattices, arXiv preprint arXiv:1608.08220 [math-ph], 2016.
Fabrizio Canfora, Maxim Kurkov, Luigi Rosa, and Patrizia Vitale, The Gribov problem in Noncommutative QED, arXiv preprint arXiv:1505.06342 [hep-th], 2015.
Niccolò Castronuovo, On the number of fixed points of the map gamma, arXiv:2102.02739 [math.NT], 2021. Mentions this sequence.
Z. Cinkir, Effective Resistances, Kirchhoff index and Admissible Invariants of Ladder Graphs, arXiv preprint arXiv:1503.06353 [math.CO], 2015.
J. B. Cosgrave and K. Dilcher, A role for generalized Fermat numbers, Math. Comp. 86 (2017), 899-933 (see also paper #10).
M. N. Deshpande, One Interesting Family of Diophantine Triplets, International Journal of Mathematical Education In Science and Technology, Vol. 33 (No. 2, Mar-Apr), 2002.
M. N. Deshpande, Hansruedi Widmer and Zachary Franco, Simultaneous Squares from Arithmetic Progressions: 10622, The American Mathematical Monthly Vol. 106, No. 9 (Nov., 1999), 867-868.
Tomislav Doslic, Planar polycyclic graphs and their Tutte polynomials, Journal of Mathematical Chemistry, Volume 51, Issue 6, 2013, pp. 1599-1607.
G. Dresden and Y. Li, Periodic Weighted Sums of Binomial Coefficients, arXiv:2210.04322 [math.NT], 2022.
E. I. Emerson, Recurrent Sequences in the Equation DQ^2=R^2+N, Fib. Quart., 7 (1969), pp. 231-242.
F. Faase, On the number of specific spanning subgraphs of the graphs G X P_n, Preliminary version of paper that appeared in Ars Combin. 49 (1998), 129-154.
Felix Flicker, Time quasilattices in dissipative dynamical systems, arXiv:1707.09371 [nlin.CD], 2017. Also SciPost Phys. 5, 001 (2018).
D. Fortin, B-spline Toeplitz Inverse Under Corner Perturbations, International Journal of Pure and Applied Mathematics, Volume 77, No. 1, 2012, 107-118. - From N. J. A. Sloane, Oct 22 2012
Dale Gerdemann, Fractal images from (4, -1) recursion, YouTube, Oct 27 2014.
Juan B. Gil and Jessica A. Tomasko, Fibonacci colored compositions and applications, arXiv:2108.06462 [math.CO], 2021.
Andrew Granville and Zhi-Wei Sun, Values of Bernoulli polynomials, Pacific J. Math. 172 (1996), 117-137, at p. 119.
T. N. E. Greville, Table for third-degree spline interpolations with equally spaced arguments, Math. Comp., 24 (1970), 179-183.
Y.-H. Guo, n-Colour even self-inverse compositions, Proc. Indian Acad. Sci. (Math. Sci.), 120 (2010), 27-33.
B. Hopkins, Spotted tilings and n-color compositions, INTEGERS 12B (2012/2013), #A6.
A. F. Horadam, Special properties of the sequence W_n(a,b; p,q), Fib. Quart., 5.5 (1967), 424-434. Case a=0,b=1; p=4, q=-1.
W. D. Hoskins, Table for third-degree spline interpolation using equi-spaced knots, Math. Comp., 25 (1971), 797-801.
Milan Janjic, On Linear Recurrence Equations Arising from Compositions of Positive Integers, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.7.
Tanya Khovanova, Recursive Sequences
Seong Ju Kim, R. Stees, and L. Taalman, Sequences of Spiral Knot Determinants, Journal of Integer Sequences, Vol. 19 (2016), #16.1.4
Clark Kimberling, Best lower and upper approximates to irrational numbers, Elemente der Mathematik, 52 (1997) 122-126.
Germain Kreweras, Complexite et circuits Euleriens dans les sommes tensorielles de graphes, J. Combin. Theory, B 24 (1978), 202-212.
Pablo Lam-Estrada, Myriam Rosalía Maldonado-Ramírez, José Luis López-Bonilla, Fausto Jarquín-Zárate, The sequences of Fibonacci and Lucas for each real quadratic fields Q(Sqrt(d)), arXiv:1904.13002 [math.NT], 2019.
Wolfdieter Lang, On polynomials related to powers of the generating function of Catalan's numbers, Fib. Quart. 38,5 (2000) 408-419; Eq.(44), lhs, m=6.
Ioana-Claudia Lazăr, Lucas sequences in t-uniform simplicial complexes, arXiv:1904.06555 [math.GR], 2019.
Hojoo Lee, Problems in elementary number theory Problem I 18.
E. Keith Lloyd, The Standard Deviation of 1, 2,..., n: Pell's Equation and Rational Triangles, Math. Gaz. vol 81 (1997), 231-243.
Dino Lorenzini and Z. Xiang, Integral points on variable separated curves, Preprint 2016.
Valcho Milchev and Tsvetelina Karamfilova, Domino tiling in grid - new dependence, arXiv:1707.09741 [math.HO], 2017.
László Németh, Trees on hyperbolic honeycombs, arXiv:1510.08311 [math.CO], 2015.
Hideyuki Ohtskua, proposer, Problem B-1351, Elementary Problems and Solutions, The Fibonacci Quarterly, Vol. 62, No. 3 (2024), p. 258.
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.
Ariel D. Procaccia and Jamie Tucker-Foltz, Compact Redistricting Plans Have Many Spanning Trees, Harvard Univ. (2021).
P. Raff, Analysis of the Number of Spanning Trees of K_2 x P_n. Contains sequence, recurrence, generating function, and more. [From Paul Raff, Mar 06 2009]
Ryan Stees, Sequences of Spiral Knot Determinants, Senior Honors Projects, Paper 84, James Madison Univ., May 2016.
F. V. Waugh and M. W. Maxfield, Side-and-diagonal numbers, Math. Mag., 40 (1967), 74-83.
Eric Weisstein's World of Mathematics, Ladder Graph
Eric Weisstein's World of Mathematics, Spanning Tree
Jianqiang Zhao, Finite Multiple zeta Values and Finite Euler Sums, arXiv preprint arXiv:1507.04917 [math.NT], 2015.
FORMULA
G.f.: x/(1-4*x+x^2).
a(n) = ((2 + sqrt(3))^n - (2 - sqrt(3))^n)/(2*sqrt(3)).
a(n) = sqrt((A001075(n)^2 - 1)/3).
a(n) = 2*a(n-1) + sqrt(3*a(n-1)^2 + 1). - Lekraj Beedassy, Feb 18 2002
a(n) = -a(-n) for all integer n. - Michael Somos, Sep 19 2008
Limit_{n->infinity} a(n)/a(n-1) = 2 + sqrt(3). - Gregory V. Richardson, Oct 06 2002
Binomial transform of A002605.
E.g.f.: exp(2*x)*sinh(sqrt(3)*x)/sqrt(3).
a(n) = S(n-1, 4) = U(n-1, 2); S(-1, x) := 0, Chebyshev's polynomials of the second kind A049310.
a(n+1) = Sum_{k=0..floor(n/2)} binomial(n-k, k)(-1)^k*4^(n - 2*k). - Paul Barry, Oct 25 2004
a(n) = Sum_{k=0..n-1} binomial(n+k,2*k+1)*2^k. - Paul Barry, Nov 30 2004
a(n) = 3*a(n-1) + 3*a(n-2) - a(n-3), n>=3. - Lekraj Beedassy, Jul 13 2006
a(n) = -A106707(n). - R. J. Mathar, Jul 07 2006
M^n * [1,0] = [A001075(n), A001353(n)], where M = the 2 X 2 matrix [2,3; 1,2]; e.g., a(4) = 56 since M^4 * [1,0] = [97, 56] = [A001075(4), A001353(4)]. - Gary W. Adamson, Dec 27 2006
Sequence satisfies 1 = f(a(n), a(n+1)) where f(u, v) = u^2 + v^2 - 4*u*v. - Michael Somos, Sep 19 2008
Rational recurrence: a(n) = (17*a(n-1)*a(n-2) - 4*(a(n-1)^2 + a(n-2)^2))/a(n-3) for n > 3. - Jaume Oliver Lafont, Dec 05 2009
If p[i] = Fibonacci(2i) and if A is the 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) = det A. - Milan Janjic, May 08 2010
a(n) = C_{n-1}^{(1)}(2), where C_n^{(m)}(x) is the Gegenbauer polynomial. - Eric W. Weisstein, Jul 16 2011
a(n) = -i*sin(n*arccos(2))/sqrt(3). - Eric W. Weisstein, Jul 16 2011
a(n) = sinh(n*arccosh(2))/sqrt(3). - Eric W. Weisstein, Jul 16 2011
a(n) = b such that Integral_{x=0..Pi/2} (sin(n*x))/(2-cos(x)) dx = c + b*log(2). - Francesco Daddi, Aug 02 2011
a(n) = sqrt(A098301(n)) = sqrt([A055793 / 3]), base 3 analog of A031150. - M. F. Hasler, Jan 16 2012
a(n+1) = Sum_{k=0..n} A101950(n,k)*3^k. - Philippe Deléham, Feb 10 2012
1, 4, 15, 56, 209, ... = INVERT(INVERT(1, 2, 3, 4, 5, ...)). - David Callan, Oct 13 2012
Product_{n >= 1} (1 + 1/a(n)) = 1 + sqrt(3). - Peter Bala, Dec 23 2012
Product_{n >= 2} (1 - 1/a(n)) = 1/4*(1 + sqrt(3)). - Peter Bala, Dec 23 2012
a(n+1) = (A001834(n) + A001835(n))/2. a(n+1) + a(n) = A001834(n). a(n+1) - a(n) = A001835(n). - Richard R. Forberg, Sep 04 2013
a(n) = -(-i)^(n+1)*Fibonacci(n, 4*i), i = sqrt(-1). - G. C. Greubel, Jun 06 2019
a(n)^2 - a(m)^2 = a(n+m) * a(n-m), a(n+2)*a(n-2) = 16*a(n+1)*a(n-1) - 15*a(n)^2, a(n+3)*a(n-2) = 15*a(n+2)*a(n-1) - 14*a(n+1)*a(n) for all integer n, m. - Michael Somos, Dec 12 2019
a(n) = 2^n*Sum_{k >= n} binomial(2*k,2*n-1)*(1/3)^(k+1). Cf. A102591. - Peter Bala, Nov 29 2021
a(n) = Sum_{k > 0} (-1)^((k-1)/2)*binomial(2*n, n+k)*(k|12), where (k|12) is the Kronecker symbol. - Greg Dresden, Oct 11 2022
Sum_{k=0..n} a(k) = (a(n+1) - a(n) - 1)/2. - Prabha Sivaramannair, Sep 22 2023
a(2n+1) = A001835(n+1) * A001834(n). - M. Farrokhi D. G., Oct 15 2023
Sum_{n>=1} arctan(1/(4*a(n)^2)) = Pi/12 (A019679) (Ohtskua, 2024). - Amiram Eldar, Aug 29 2024
EXAMPLE
For example, when n = 3:
****
.***
.***
can be packed with dominoes in 4 different ways: 3 in which the top row is tiled with two horizontal dominoes and 1 in which the top row has two vertical and one horizontal domino, as shown below, so a(2) = 4.
---- ---- ---- ||--
.||| .--| .|-- .|||
.||| .--| .|-- .|||
G.f. = x + 4*x^2 + 15*x^3 + 56*x^4 + 209*x^5 + 780*x^6 + 2911*x^7 + 10864*x^8 + ...
MAPLE
A001353 := proc(n) option remember; if n <= 1 then n else 4*A001353(n-1)-A001353(n-2); fi; end;
A001353:=z/(1-4*z+z**2); # Simon Plouffe in his 1992 dissertation.
seq( simplify(ChebyshevU(n-1, 2)), n=0..20); # G. C. Greubel, Dec 23 2019
MATHEMATICA
a[n_] := (MatrixPower[{{1, 2}, {1, 3}}, n].{{1}, {1}})[[2, 1]]; Table[ a[n], {n, 0, 30}] (* Robert G. Wilson v, Jan 13 2005 *)
Table[GegenbauerC[n-1, 1, 2]], {n, 0, 30}] (* Zerinvary Lajos, Jul 14 2009 *)
Table[-((I Sin[n ArcCos[2]])/Sqrt[3]), {n, 0, 30}] // FunctionExpand (* Eric W. Weisstein, Jul 16 2011 *)
Table[Sinh[n ArcCosh[2]]/Sqrt[3], {n, 0, 30}] // FunctionExpand (* Eric W. Weisstein, Jul 16 2011 *)
Table[ChebyshevU[n-1, 2], {n, 0, 30}] (* Eric W. Weisstein, Jul 16 2011 *)
a[0]:=0; a[1]:=1; a[n_]:= a[n]= 4a[n-1] - a[n-2]; Table[a[n], {n, 0, 30}] (* Alonso del Arte, Jul 19 2011 *)
LinearRecurrence[{4, -1}, {0, 1}, 30] (* Sture Sjöstedt, Dec 06 2011 *)
Round@Table[Fibonacci[2n, Sqrt[2]]/Sqrt[2], {n, 0, 30}] (* Vladimir Reshetnikov, Sep 15 2016 *)
PROG
(PARI) M = [ 1, 1, 0; 1, 3, 1; 0, 1, 1]; for(i=0, 30, print1(([1, 0, 0]*M^i)[2], ", ")) \\ Lambert Klasen (Lambert.Klasen(AT)gmx.net), Jan 25 2005
(PARI) {a(n) = real( (2 + quadgen(12))^n / quadgen(12) )}; /* Michael Somos, Sep 19 2008 */
(PARI) {a(n) = polchebyshev(n-1, 2, 2)}; /* Michael Somos, Sep 19 2008 */
(PARI) concat(0, Vec(x/(1-4*x+x^2) + O(x^30))) \\ Altug Alkan, Oct 30 2015
(Sage) [lucas_number1(n, 4, 1) for n in range(30)] # Zerinvary Lajos, Apr 22 2009
(Sage) [chebyshev_U(n-1, 2) for n in (0..20)] # G. C. Greubel, Dec 23 2019
(Haskell)
a001353 n = a001353_list !! n
a001353_list =
0 : 1 : zipWith (-) (map (4 *) $ tail a001353_list) a001353_list
-- Reinhard Zumkeller, Aug 14 2011
(GAP) a:=[0, 1];; for n in [3..30] do a[n]:=4*a[n-1]-a[n-2]; od; a; # Muniru A Asiru, Feb 16 2018
(Magma) I:=[0, 1]; [n le 2 select I[n] else 4*Self(n-1)-Self(n-2): n in [1..30]]; // G. C. Greubel, Jun 06 2019
(Python)
a001353 = [0, 1]
for n in range(30): a001353.append(4*a001353[-1] - a001353[-2])
print(a001353) # Gennady Eremin, Feb 05 2022
CROSSREFS
A bisection of A002530.
Cf. A125077.
A row of A116469.
Chebyshev sequence U(n, m): A000027 (m=1), this sequence (m=2), A001109 (m=3), A001090 (m=4), A004189 (m=5), A004191 (m=6), A007655 (m=7), A077412 (m=8), A049660 (m=9), A075843 (m=10), A077421 (m=11), A077423 (m=12), A097309 (m=13), A097311 (m=14), A097313 (m=15), A029548 (m=16), A029547 (m=17), A144128 (m=18), A078987 (m=19), A097316 (m=33).
Cf. A323182.
KEYWORD
nonn,easy,nice
STATUS
approved
a(n) = 10*a(n-1) - a(n-2); a(0) = 0, a(1) = 1.
+10
61
0, 1, 10, 99, 980, 9701, 96030, 950599, 9409960, 93149001, 922080050, 9127651499, 90354434940, 894416697901, 8853812544070, 87643708742799, 867583274883920, 8588189040096401, 85014307126080090, 841554882220704499, 8330534515080964900, 82463790268588944501, 816307368170808480110
OFFSET
0,3
COMMENTS
Indices of square numbers which are also generalized pentagonal numbers.
If t(n) denotes the n-th triangular number, t(A105038(n))=a(n)*a(n+1). - Robert Phillips (bobanne(AT)bellsouth.net), May 25 2008
The n-th term is a(n) = ((5+sqrt(24))^n - (5-sqrt(24))^n)/(2*sqrt(24)). - Sture Sjöstedt, May 31 2009
For n >= 2, a(n) equals the permanent of the (n-1) X (n-1) tridiagonal matrix with 10's along the main diagonal, and i's along the superdiagonal and the subdiagonal (i is the imaginary unit). - John M. Campbell, Jul 08 2011
a(n) and b(n) (A001079) are the nonnegative proper solutions of the Pell equation b(n)^2 - 6*(2*a(n))^2 = +1. See the cross reference to A001079 below. - Wolfdieter Lang, Jun 26 2013
For n >= 1, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,...,9}. - Milan Janjic, Jan 25 2015
For n > 1, this also gives the number of (n-1)-decimal-digit numbers which avoid a particular two-digit number with distinct digits. For example, there are a(5) = 9701 4-digit numbers which do not include "39" as a substring; see Wikipedia link. - Charles R Greathouse IV, Jan 14 2016
All possible solutions for y in Pell equation x^2 - 24*y^2 = 1. The values for x are given in A001079. - Herbert Kociemba, Jun 05 2022
Dickson on page 384 gives the Diophantine equation "(20) 24x^2 + 1 = y^2" and later states "... three consecutive sets (x_i, y_i) of solutions of (20) or 2x^2 + 1 = 3y^2 satisfy x_{n+1} = 10x_n - x_{n-1}, y_{n+1} = 10y_n - y_{n-1} with (x_1, y_1) = (0, 1) or (1, 1), (x_2, y_2) = (1, 5) or (11, 9), respectively." The first set of values (x_n, y_n) = (A001079(n-1), a(n-1)). - Michael Somos, Jun 19 2023
REFERENCES
L. E. Dickson, History of the Theory of Numbers, Vol. II, Diophantine Analysis. AMS Chelsea Publishing, Providence, Rhode Island, 1999, p. 384.
LINKS
K. Andersen, L. Carbone, and D. Penta, Kac-Moody Fibonacci sequences, hyperbolic golden ratios, and real quadratic fields, Journal of Number Theory and Combinatorics, Vol 2, No. 3 pp 245-278, 2011. See Section 9.
D. Birmajer, J. B. Gil, and M. D. Weiner, On the Enumeration of Restricted Words over a Finite Alphabet, J. Int. Seq. 19 (2016) # 16.1.3, Example 12.
E. I. Emerson, Recurrent Sequences in the Equation DQ^2=R^2+N, Fib. Quart., 7 (1969), pp. 231-242.
D. Fortin, B-spline Toeplitz inverse under corner perturbations, International Journal of Pure and Applied Mathematics, Volume 77, No. 1, 2012, 107-118. - From N. J. A. Sloane, Oct 22 2012
A. F. Horadam, Special properties of the sequence W_n(a,b; p,q), Fib. Quart., 5.5 (1967), 424-434. Case a=0,b=1; p=10, q=-1.
Milan Janjic, On Linear Recurrence Equations Arising from Compositions of Positive Integers, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.7.
Tanya Khovanova, Recursive Sequences
Pablo Lam-Estrada, Myriam Rosalía Maldonado-Ramírez, José Luis López-Bonilla, and Fausto Jarquín-Zárate, The sequences of Fibonacci and Lucas for each real quadratic fields Q(Sqrt(d)), arXiv:1904.13002 [math.NT], 2019.
Wolfdieter Lang, On polynomials related to powers of the generating function of Catalan's numbers, Fib. Quart. 38,5 (2000) 408-419; Eq.(44), lhs, m=12.
Roger B. Nelson, Multi-Polygonal Numbers, Mathematics Magazine, Vol. 89, No. 3 (June 2016), pp. 159-164.
Robert Phillips, Polynomials of the form 1+4ke+4ke^2, 2008.
Wikipedia, Curse of 39
Jianqiang Zhao, Finite Multiple zeta Values and Finite Euler Sums, arXiv:1507.04917 [math.NT], 2015.
FORMULA
a(n) = S(2*n-1, sqrt(12))/sqrt(12) = S(n-1, 10); S(n, x) := U(n, x/2), Chebyshev polynomials of 2nd kind, A049310. S(-1, x) := 0.
A001079(n) = sqrt(24*(a(n)^2)+1), that is a(n) = sqrt((A001079(n)^2-1)/24).
From Barry E. Williams, Aug 18 2000: (Start)
a(n) = ( (5+2*sqrt(6))^n - (5-2*sqrt(6))^n )/(4*sqrt(6)).
G.f.: x/(1-10*x+x^2). (End)
a(-n) = -a(n). - Michael Somos, Sep 05 2006
From Mohamed Bouhamida, May 26 2007: (Start)
a(n) = 9*(a(n-1) + a(n-2)) - a(n-3).
a(n) = 11*(a(n-1) - a(n-2)) + a(n-3).
a(n) = 10*a(n-1) - a(n-2). (End)
a(n+1) = Sum_{k=0..n} A101950(n,k)*9^k. - Philippe Deléham, Feb 10 2012
From Peter Bala, Dec 23 2012: (Start)
Product {n >= 1} (1 + 1/a(n)) = 1/2*(2 + sqrt(6)).
Product {n >= 2} (1 - 1/a(n)) = 1/5*(2 + sqrt(6)). (End)
a(n) = (A054320(n-1) + A072256(n))/2. - Richard R. Forberg, Nov 21 2013
a(2*n - 1) = A046173(n).
E.g.f.: exp(5*x)*sinh(2*sqrt(6)*x)/(2*sqrt(6)). - Stefano Spezia, Dec 12 2022
a(n) = Sum_{k = 0..n-1} binomial(n+k, 2*k+1)*8^k = Sum_{k = 0..n-1} (-1)^(n+k+1)* binomial(n+k, 2*k+1)*12^k. - Peter Bala, Jul 18 2023
EXAMPLE
a(2)=10 and (3(-8)^2-(-8))/2=10^2, a(3)=99 and (3(81)^2-(81))/2=99^2. - Michael Somos, Sep 05 2006
G.f. = x + 10*x^2 + 99*x^3 + 980*x^4 + 9701*x^5 + 96030*x^6 + ...
MAPLE
A004189 := proc(n)
option remember;
if n <= 1 then
n ;
else
10*procname(n-1)-procname(n-2) ;
end if;
end proc:
seq(A004189(n), n=0..20) ; # R. J. Mathar, Apr 30 2017
seq( simplify(ChebyshevU(n-1, 5)), n=0..20); # G. C. Greubel, Dec 23 2019
MATHEMATICA
Table[GegenbauerC[n-1, 1, 5], {n, 0, 30}] (* Vladimir Joseph Stephan Orlovsky, Sep 11 2008; modified by G. C. Greubel, Jun 06 2019 *)
LinearRecurrence[{10, -1}, {0, 1}, 20] (* Jean-François Alcover, Nov 15 2017 *)
ChebyshevU[Range[21] -2, 5] (* G. C. Greubel, Dec 23 2019 *)
PROG
(PARI) {a(n) = subst(poltchebi(n+1) - 5*poltchebi(n), 'x, 5) / 24}; /* Michael Somos, Sep 05 2006 */
(PARI) a(n)=([9, 1; 8, 1]^(n-1)*[1; 1])[1, 1] \\ Charles R Greathouse IV, Jan 14 2016
(PARI) vector(21, n, n--; polchebyshev(n-1, 2, 5) ) \\ G. C. Greubel, Dec 23 2019
(Sage) [lucas_number1(n, 10, 1) for n in range(22)] # Zerinvary Lajos, Jun 25 2008
(Sage) [chebyshev_U(n-1, 5) for n in (0..20)] # G. C. Greubel, Dec 23 2019
(Magma) [ n eq 1 select 0 else n eq 2 select 1 else 10*Self(n-1)-Self(n-2): n in [1..20] ]; // Vincenzo Librandi, Aug 19 2011
(GAP) m:=5;; a:=[0, 1];; for n in [3..20] do a[n]:=2*m*a[n-1]-a[n-2]; od; a; # G. C. Greubel, Dec 23 2019
CROSSREFS
Chebyshev sequence U(n, m): A000027 (m=1), A001353 (m=2), A001109 (m=3), A001090 (m=4), this sequence (m=5), A004191 (m=6), A007655 (m=7), A077412 (m=8), A049660 (m=9), A075843 (m=10), A077421 (m=11), A077423 (m=12), A097309 (m=13), A097311 (m=14), A097313 (m=15), A029548 (m=16), A029547 (m=17), A144128 (m=18), A078987 (m=19), A097316 (m=33).
Cf. A323182.
KEYWORD
easy,nonn
STATUS
approved
a(n) = 8*a(n-1) - a(n-2); a(0) = 0, a(1) = 1.
(Formerly M4554 N1936)
+10
52
0, 1, 8, 63, 496, 3905, 30744, 242047, 1905632, 15003009, 118118440, 929944511, 7321437648, 57641556673, 453811015736, 3572846569215, 28128961537984, 221458845734657, 1743541804339272, 13726875588979519, 108071462907496880, 850844827670995521, 6698687158460467288
OFFSET
0,3
COMMENTS
This sequence gives the values of y in solutions of the Diophantine equation x^2 - 15*y^2 = 1; the corresponding values of x are in A001091. - Vincenzo Librandi, Nov 12 2010 [edited by Jon E. Schoenfield, May 02 2014]
For n >= 2, a(n) equals the permanent of the (n-1) X (n-1) tridiagonal matrix with 8's along the main diagonal, and i's along the superdiagonal and the subdiagonal (i is the imaginary unit). - John M. Campbell, Jul 08 2011
For n >= 1, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,...,7}. - Milan Janjic, Jan 25 2015
From Klaus Purath, Jul 25 2024: (Start)
For any three consecutive terms (x, y, z) y^2 - x*z = 1 always applies.
a(n) = (t(i+2n) - t(i))/(t(i+n+1) - t(i+n-1)) where (t) is any recurrence t(k) = 9t(k-1) - 9t(k-2) + t(k-3) or t(k) = 8t(k-1) - t(k-2) without regard to initial values.
In particular, if the recurrence (t) of the form (9,-9,1) has the initial values t(0) = 1, t(1) = 2, t(2) = 9, a(n) = t(n) - 1 applies. (End)
REFERENCES
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) - From N. J. A. Sloane, May 30 2012
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).
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000 (terms 0..100 from T. D. Noe)
Marco Abrate, Stefano Barbero, Umberto Cerruti, and Nadir Murru, Polynomial sequences on quadratic curves, Integers, Vol. 15, 2015, #A38.
K. Andersen, L. Carbone, and D. Penta, Kac-Moody Fibonacci sequences, hyperbolic golden ratios, and real quadratic fields, Journal of Number Theory and Combinatorics, Vol 2, No. 3 pp 245-278, 2011. See Section 9.
H. Brocard, Notes élémentaires sur le problème de Peel, Nouvelle Correspondance Mathématique, 4 (1878), 337-343.
E. I. Emerson, Recurrent Sequences in the Equation DQ^2=R^2+N, Fib. Quart., 7 (1969), pp. 231-242.
A. F. Horadam, Special properties of the sequence W_n(a,b; p,q), Fib. Quart., 5.5 (1967), 424-434. Case a=0,b=1; p=8, q=-1.
Milan Janjic, On Linear Recurrence Equations Arising from Compositions of Positive Integers, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.7.
Tanya Khovanova, Recursive Sequences
Wolfdieter Lang, On polynomials related to powers of the generating function of Catalan's numbers, Fib. Quart. 38,5 (2000) 408-419; Eq.(44), lhs, m=10.
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.
Shobhna Singh and Felix Flicker, Exact Solution to the Quantum and Classical Dimer Models on the Spectre Aperiodic Monotiling, arXiv:2309.14447 [cond-mat.str-el], 2023.
FORMULA
15*a(n)^2 - A001091(n)^2 = -1.
a(n) = sqrt((A001091(n)^2 - 1)/15).
a(n) = S(2*n-1, sqrt(10))/sqrt(10) = S(n-1, 8); S(n, x) := U(n, x/2), Chebyshev polynomials of 2nd kind, A049310, with S(-1, x) := 0.
From Barry E. Williams, Aug 18 2000: (Start)
a(n) = ((4+sqrt(15))^n - (4-sqrt(15))^n)/(2*sqrt(15)).
G.f.: x/(1-8*x+x^2). (End)
Limit_{n->infinity} a(n)/a(n-1) = 4 + sqrt(15). - Gregory V. Richardson, Oct 13 2002
[A070997(n-1), a(n)] = [1,6; 1,7]^n * [1,0]. - Gary W. Adamson, Mar 21 2008
a(-n) = -a(n). - Michael Somos, Apr 05 2008
a(n+1) = Sum_{k=0..n} A101950(n,k)*7^k. - Philippe Deléham, Feb 10 2012
From Peter Bala, Dec 23 2012: (Start)
Product_{n >= 1} (1 + 1/a(n)) = (1/3)*(3 + sqrt(15)).
Product_{n >= 2} (1 - 1/a(n)) = (1/8)*(3 + sqrt(15)).
(End)
a(n) = A136325(n)/3. - Greg Dresden, Sep 12 2019
E.g.f.: exp(4*x)*sinh(sqrt(15)*x)/sqrt(15). - Stefano Spezia, Dec 12 2022
a(n) = Sum_{k = 0..n-1} binomial(n+k, 2*k+1)*6^k = Sum_{k = 0..n-1} (-1)^(n+k+1)* binomial(n+k, 2*k+1)*10^k. - Peter Bala, Jul 17 2023
EXAMPLE
G.f. = x + 8*x^2 + 63*x^3 + 496*x^4 + 3905*x^5 + 30744*x^6 + 242047*x^7 + ...
MAPLE
A001090:=1/(1-8*z+z**2); # Simon Plouffe in his 1992 dissertation
seq( simplify(ChebyshevU(n-1, 4)), n=0..20); # G. C. Greubel, Dec 23 2019
MATHEMATICA
Table[GegenbauerC[n-1, 1, 4]], {n, 0, 20}] (* Vladimir Joseph Stephan Orlovsky, Sep 11 2008 *)
LinearRecurrence[{8, -1}, {0, 1}, 30] (* Harvey P. Dale, Aug 29 2012 *)
a[n_]:= ChebyshevU[n-1, 4]; (* Michael Somos, May 28 2014 *)
CoefficientList[Series[x/(1-8*x+x^2), {x, 0, 20}], x] (* G. C. Greubel, Dec 20 2017 *)
PROG
(PARI) {a(n) = subst(poltchebi(n+1) - 4 * poltchebi(n), x, 4) / 15}; /* Michael Somos, Apr 05 2008 */
(PARI) {a(n) = polchebyshev(n-1, 2, 4)}; /* Michael Somos, May 28 2014 */
(PARI) my(x='x+O('x^30)); concat([0], Vec(x/(1-8*x-x^2))) \\ G. C. Greubel, Dec 20 2017
(SageMath) [lucas_number1(n, 8, 1) for n in range(22)] # Zerinvary Lajos, Jun 25 2008
(SageMath) [chebyshev_U(n-1, 4) for n in (0..20)] # G. C. Greubel, Dec 23 2019
(Magma) I:=[0, 1]; [n le 2 select I[n] else 8*Self(n-1) - Self(n-2): n in [1..30]]; // G. C. Greubel, Dec 20 2017
(GAP) m:=4;; a:=[0, 1];; for n in [3..20] do a[n]:=2*m*a[n-1]-a[n-2]; od; a; # G. C. Greubel, Dec 23 2019
CROSSREFS
Equals one-third A136325.
Chebyshev sequence U(n, m): A000027 (m=1), A001353 (m=2), A001109 (m=3), this sequence (m=4), A004189 (m=5), A004191 (m=6), A007655 (m=7), A077412 (m=8), A049660 (m=9), A075843 (m=10), A077421 (m=11), A077423 (m=12), A097309 (m=13), A097311 (m=14), A097313 (m=15), A029548 (m=16), A029547 (m=17), A144128 (m=18), A078987 (m=19), A097316 (m=33).
Cf. A323182.
KEYWORD
nonn,easy
EXTENSIONS
More terms from Wolfdieter Lang, Aug 02 2000
STATUS
approved
Standard deviation of A007654.
(Formerly M4948)
+10
45
0, 1, 14, 195, 2716, 37829, 526890, 7338631, 102213944, 1423656585, 19828978246, 276182038859, 3846719565780, 53577891882061, 746243766783074, 10393834843080975, 144767444036350576, 2016350381665827089, 28084137899285228670, 391161580208327374291, 5448177985017298011404
OFFSET
1,3
COMMENTS
a(n) corresponds also to one-sixth the area of Fleenor-Heronian triangle with middle side A003500(n). - Lekraj Beedassy, Jul 15 2002
a(n) give all (nontrivial, integer) solutions of Pell equation b(n+1)^2 - 48*a(n+1)^2 = +1 with b(n+1)=A011943(n), n>=0.
For n>=3, a(n) equals the permanent of the (n-2) X (n-2) tridiagonal matrix with 14's along the main diagonal, and i's along the superdiagonal and the subdiagonal (i is the imaginary unit). - John M. Campbell, Jul 08 2011
For n>1, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,...,13}. - Milan Janjic, Jan 25 2015
6*a(n)^2 = 6*S(n-1, 14)^2 is the triangular number Tri((T(n, 7) - 1)/2) with Tri = A000217 and T = A053120. This is instance k = 3 of the general k-identity given in a comment to A001109. - Wolfdieter Lang, Feb 01 2016
REFERENCES
D. A. Benaron, personal communication.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
Indranil Ghosh, Table of n, a(n) for n = 1..874 (terms 1..100 from T. D. Noe)
R. Flórez, R. A. Higuita and A. Mukherjee, Alternating Sums in the Hosoya Polynomial Triangle, Article 14.9.5 Journal of Integer Sequences, Vol. 17 (2014).
D. S. Hale, 3165. Perfect Squares of the Form 48n^2+1, Math. Gaz., Oct. 1966, page 307.
Milan Janjic, On Linear Recurrence Equations Arising from Compositions of Positive Integers, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.7.
Tanya Khovanova, Recursive Sequences
Murray S. Klamkin, Perfect Squares of the Form (m^2 - 1)a_n^2 + t, Math. Mag., 1969, page 111.
E. K. Lloyd, The standard deviation of 1, 2, ..., n, Pell's equation and rational triangles, The Mathematical Gazette, Vol. 81, No. 491 (Jul., 1997), pp. 231-243.
Dino Lorenzini and Z. Xiang, Integral points on variable separated curves, Preprint 2016.
FORMULA
a(n) = 14*a(n-1) - a(n-2).
G.f.: x^2/(1-14*x+x^2).
a(n+1) ~ 1/24*sqrt(3)*(2 + sqrt(3))^(2*n). - Joe Keane (jgk(AT)jgk.org), May 15 2002
a(n+1) = S(n-1, 14), n>=0, with S(n, x) := U(n, x/2) Chebyshev's polynomials of the second kind. S(-1, x) := 0. See A049310.
a(n+1) = ( (7+4*sqrt(3))^n - (7-4*sqrt(3))^n )/(8*sqrt(3)).
a(n+1) = sqrt((A011943(n)^2 - 1)/48), n>=0.
Chebyshev's polynomials U(n-2, x) evaluated at x=7.
a(n) = A001353(2n)/4. - Lekraj Beedassy, Jul 15 2002
4*a(n+1) + A046184(n) = A055793(n+2) + A098301(n+1) 4*a(n+1) + A098301(n+1) + A055793(n+2) = A046184(n+1) (4*a(n+1))^2 = A098301(2n+1) (conjectures). - Creighton Dement, Nov 02 2004
(4*a(n))^2 = A103974(n)^2 - A011922(n-1)^2. - Paul D. Hanna, Mar 06 2005
From Mohamed Bouhamida, May 26 2007: (Start)
a(n) = 13*( a(n-1) + a(n-2) ) - a(n-3).
a(n) = 15*( a(n-1) - a(n-2) ) + a(n-3). (End)
a(n) = b such that (-1)^n/4*Integral_{x=-Pi/2..Pi/2} (sin((2*n-2)*x))/(2-sin(x)) dx = c+b*log(3). - Francesco Daddi, Aug 02 2011
a(n+2) = Sum_{k=0..n} A101950(n,k)*13^k. - Philippe Deléham, Feb 10 2012
Product {n >= 1} (1 + 1/a(n)) = 1/3*(3 + 2*sqrt(3)). - Peter Bala, Dec 23 2012
Product {n >= 2} (1 - 1/a(n)) = 1/7*(3 + 2*sqrt(3)). - Peter Bala, Dec 23 2012
a(n) = (A028230(n) - A001570(n))/2. - Richard R. Forberg, Nov 14 2013
E.g.f.: 1 - exp(7*x)*(12*cosh(4*sqrt(3)*x) - 7*sqrt(3)*sinh(4*sqrt(3)*x))/12. - Stefano Spezia, Dec 11 2022
EXAMPLE
G.f. = x^2 + 14*x^3 + 195*x^4 + 2716*x^5 + 37829*x^6 + 526890*x^7 + ...
MAPLE
0, seq(orthopoly[U](n, 7), n=0..30); # Robert Israel, Feb 04 2016
MATHEMATICA
Table[GegenbauerC[n, 1, 7], {n, 0, 20}] (* Vladimir Joseph Stephan Orlovsky, Sep 11 2008 *)
LinearRecurrence[{14, -1}, {0, 1}, 20] (* Vincenzo Librandi, Feb 02 2016 *)
ChebyshevU[Range[21] -2, 7] (* G. C. Greubel, Dec 23 2019 *)
Table[Sum[Binomial[n, 2 k - 1]*7^(n - 2 k + 1)*48^(k - 1), {k, 1, n}], {n, 0, 15}] (* Horst H. Manninger, Jan 16 2022 *)
PROG
(Sage) [lucas_number1(n, 14, 1) for n in range(0, 20)] # Zerinvary Lajos, Jun 25 2008
(Sage) [chebyshev_U(n, 7) for n in (-1..20)] # G. C. Greubel, Dec 23 2019
(Magma) [n le 2 select n-1 else 14*Self(n-1)-Self(n-2): n in [1..70]]; // Vincenzo Librandi, Feb 02 2016
(PARI) concat(0, Vec((x^2/(1-14*x+x^2) + O(x^30)))) \\ Michel Marcus, Feb 02 2016
(PARI) vector(21, n, polchebyshev(n-2, 2, 7) ) \\ G. C. Greubel, Dec 23 2019
(GAP) m:=7;; a:=[0, 1];; for n in [3..20] do a[n]:=2*m*a[n-1]-a[n-2]; od; a; # G. C. Greubel, Dec 23 2019
CROSSREFS
Chebyshev sequence U(n, m): A000027 (m=1), A001353 (m=2), A001109 (m=3), A001090 (m=4), A004189 (m=5), A004191 (m=6), this sequence (m=7), A077412 (m=8), A049660 (m=9), A075843 (m=10), A077421 (m=11), A077423 (m=12), A097309 (m=13), A097311 (m=14), A097313 (m=15), A029548 (m=16), A029547 (m=17), A144128 (m=18), A078987 (m=19), A097316 (m=33).
Cf. A323182.
KEYWORD
nonn,easy
EXTENSIONS
Chebyshev comments from Wolfdieter Lang, Nov 08 2002
STATUS
approved
Chebyshev U(n,x) polynomial evaluated at x=8.
+10
42
1, 16, 255, 4064, 64769, 1032240, 16451071, 262184896, 4178507265, 66593931344, 1061324394239, 16914596376480, 269572217629441, 4296240885694576, 68470281953483775, 1091228270370045824, 17391182043967249409
OFFSET
0,2
COMMENTS
For positive n, a(n) equals the permanent of the n X n tridiagonal matrix with 16's along the main diagonal, and i's along the superdiagonal and the subdiagonal (i is the imaginary unit). - John M. Campbell, Jul 08 2011
For n>=2, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,...,15}. - Milan Janjic, Jan 23 2015
LINKS
Tanya Khovanova, Recursive Sequences
Pablo Lam-Estrada, Myriam Rosalía Maldonado-Ramírez, José Luis López-Bonilla, and Fausto Jarquín-Zárate, The sequences of Fibonacci and Lucas for each real quadratic fields Q(Sqrt(d)), arXiv:1904.13002 [math.NT], 2019.
FORMULA
a(n) = 16*a(n-1) - a(n-2), n>=1, a(-1)=0, a(0)=1.
a(n) = S(n, 16) with S(n, x) := U(n, x/2), Chebyshev's polynomials of the second kind. See A049310.
G.f.: 1/(1 - 16*x + x^2).
a(n) = (((8 + 3*sqrt(7))^(n+1) - (8 - 3*sqrt(7))^(n+1)))/(6*sqrt(7)).
a(n) = sqrt((A001081(n+1)^2-1)/63).
a(n) = Sum_{k=0..n} A101950(n,k)*15^k. - Philippe Deléham, Feb 10 2012
Product {n >= 0} (1 + 1/a(n)) = 1/7*(7 + 3*sqrt(7)). - Peter Bala, Dec 23 2012
Product {n >= 1} (1 - 1/a(n)) = 1/16*(7 + 3*sqrt(7)). - Peter Bala, Dec 23 2012
MAPLE
seq( simplify(ChebyshevU(n, 8)), n=0..20); # G. C. Greubel, Dec 22 2019
MATHEMATICA
Table[GegenbauerC[n, 1, 8], {n, 0, 20}] (* Vladimir Joseph Stephan Orlovsky, Sep 11 2008 *)
CoefficientList[Series[1/(1-16x+x^2), {x, 0, 20}], x] (* Vincenzo Librandi, Dec 24 2012 *)
LinearRecurrence[{16, -1}, {1, 16}, 30] (* G. C. Greubel, Jan 18 2018 *)
ChebyshevU[Range[21] -1, 8] (* G. C. Greubel, Dec 22 2019 *)
PROG
(Sage) [lucas_number1(n, 16, 1) for n in range(1, 20)] # Zerinvary Lajos, Jun 25 2008
(Sage) [chebyshev_U(n, 8) for n in (0..20)] # G. C. Greubel, Dec 22 2019
(Magma) I:=[1, 16, 255]; [n le 3 select I[n] else 16*Self(n-1)-Self(n-2): n in [1..20]]; // Vincenzo Librandi, Dec 24 2012
(PARI) vector( 21, n, polchebyshev(n-1, 2, 8) ) \\ G. C. Greubel, Jan 18 2018
(GAP) m:=8;; a:=[1, 2*m];; for n in [3..20] do a[n]:=2*m*a[n-1]-a[n-2]; od; a; # G. C. Greubel, Dec 22 2019
CROSSREFS
Chebyshev sequence U(n, m): A000027 (m=1), A001353 (m=2), A001109 (m=3), A001090 (m=4), A004189 (m=5), A004191 (m=6), A007655 (m=7), this sequence (m=8), A049660 (m=9), A075843 (m=10), A077421 (m=11), A077423 (m=12), A097309 (m=13), A097311 (m=14), A097313 (m=15), A029548 (m=16), A029547 (m=17), A144128 (m=18), A078987 (m=19), A097316 (m=33).
Cf. A323182.
KEYWORD
nonn,easy
AUTHOR
Wolfdieter Lang, Nov 08 2002
STATUS
approved
Expansion of 1/(1 - 12*x + x^2).
+10
32
1, 12, 143, 1704, 20305, 241956, 2883167, 34356048, 409389409, 4878316860, 58130412911, 692686638072, 8254109243953, 98356624289364, 1172025382228415, 13965947962451616, 166419350167190977, 1983066254043840108, 23630375698358890319, 281581442126262843720
OFFSET
0,2
COMMENTS
Chebyshev's polynomials U(n,x) evaluated at x=6.
a(n) give all (nontrivial, integer) solutions of Pell equation b(n)^2 - 35*a(n)^2 = +1 with b(n)=A023038(n+1), n>=0.
For positive n, a(n) equals the permanent of the tridiagonal matrix of order n with 12's along the main diagonal, and i's along the superdiagonal and the subdiagonal (i is the imaginary unit). - John M. Campbell, Jul 08 2011
For n>=1, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,...,11}. - Milan Janjic, Jan 26 2015
a(n) = -a(-2-n) for all n in Z. - Michael Somos, Jun 29 2019
LINKS
K. Andersen, L. Carbone, and D. Penta, Kac-Moody Fibonacci sequences, hyperbolic golden ratios, and real quadratic fields, Journal of Number Theory and Combinatorics, Vol 2, No. 3 pp 245-278, 2011. See Section 9.
D. Birmajer, J. B. Gil, and M. D. Weiner, On the Enumeration of Restricted Words over a Finite Alphabet, J. Int. Seq. 19 (2016) # 16.1.3, example 12.
Milan Janjic, On Linear Recurrence Equations Arising from Compositions of Positive Integers, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.7.
Tanya Khovanova, Recursive Sequences
FORMULA
a(n) = S(n, 12) with S(n, x) := U(n, x/2) Chebyshev's polynomials of the second kind. See A049310.
a(n) = ((6+sqrt(35))^(n+1) - (6-sqrt(35))^(n+1))/(2*sqrt(35)).
a(n) = sqrt((A023038(n)^2 - 1)/35).
[A077417(n), a(n)] = the 2 X 2 matrix [1,10; 1,11]^(n+1) * [1,0]. - Gary W. Adamson, Mar 19 2008
a(n) = 12*a(n-1) - a(n-2) for n>1, a(0)=1, a(1)=12. - Philippe Deléham, Nov 17 2008
a(n) = b such that (-1)^(n+1)*Integral_{x=0..Pi/2} (sin((n+1)*x))/(6+cos(x)) dx = c + b*(log(2)+log(3)-log(7)). - Francesco Daddi, Aug 01 2011
a(n) = Sum_{k=0..n} A101950(n,k)*11^k. - Philippe Deléham, Feb 10 2012
From Peter Bala, Dec 23 2012 (Start):
Product_{n>=0} (1 + 1/a(n)) = 1/5*(5 + sqrt(35)).
Product_{n>=1} (1 - 1/a(n)) = 1/12*(5 + sqrt(35)). (End)
E.g.f.: exp(6*x)*(35*cosh(sqrt(35)*x) + 6*sqrt(35)*sinh(sqrt(35)*x))/35. - Stefano Spezia, Dec 14 2022
EXAMPLE
G.f. = 1 + 12*x + 143*x^2 + 1704*x^3 + 20305*x^4 + 241956*x^5 + ...
MAPLE
seq( simplify(ChebyshevU(n, 6)), n=0..20); # G. C. Greubel, Dec 23 2019
MATHEMATICA
Table[GegenbauerC[n, 1, 6], {n, 0, 20}] (* Vladimir Joseph Stephan Orlovsky, Sep 11 2008 *)
CoefficientList[Series[1/(1-12*x+x^2), {x, 0, 30}], x] (* T. D. Noe, Aug 01 2011 *)
LinearRecurrence[{12, -1}, {1, 12}, 30] (* Harvey P. Dale, Feb 17 2016 *)
a[n_]:= ChebyshevU[n, 6]; (* Michael Somos, Jun 29 2019 *)
PROG
(Sage) [lucas_number1(n, 12, 1) for n in range(1, 20)] # Zerinvary Lajos, Jun 25 2008
(Sage) [chebyshev_U(n, 6) for n in (0..20)] # G. C. Greubel, Dec 23 2019
(Magma) I:=[1, 12]; [n le 2 select I[n] else 12*Self(n-1)-Self(n-2): n in [1..20]]; // Vincenzo Librandi, Jun 13 2012
(PARI) Vec(1/(1-12*x+x^2)+O(x^99)) \\ Charles R Greathouse IV, Sep 23 2012
(PARI) {a(n) = polchebyshev(n, 2, 6)}; \\ Michael Somos, Jun 29 2019
(GAP) m:=8;; a:=[1, 2*m];; for n in [3..20] do a[n]:=2*m*a[n-1]-a[n-2]; od; a; # G. C. Greubel, Dec 23 2019
CROSSREFS
Chebyshev sequence U(n, m): A000027 (m=1), A001353 (m=2), A001109 (m=3), A001090 (m=4), A004189 (m=5), this sequence (m=6), A007655 (m=7), A077412 (m=8), A049660 (m=9), A075843 (m=10), A077421 (m=11), A077423 (m=12), A097309 (m=13), A097311 (m=14), A097313 (m=15), A029548 (m=16), A029547 (m=17), A144128 (m=18), A078987 (m=19), A097316 (m=33).
Cf. A323182.
KEYWORD
nonn,easy
EXTENSIONS
Chebyshev comments and a(n) formulas from Wolfdieter Lang, Nov 08 2002
STATUS
approved
Chebyshev sequence U(n,11)=S(n,22) with Diophantine property.
+10
29
1, 22, 483, 10604, 232805, 5111106, 112211527, 2463542488, 54085723209, 1187422368110, 26069206375211, 572335117886532, 12565303387128493, 275864339398940314, 6056450163389558415, 132966039255171344816
OFFSET
0,2
COMMENTS
b(n)^2 - 30*(2*a(n))^2 = 1 with the companion sequence b(n)=A077422(n+1).
For positive n, a(n) equals the permanent of the n X n tridiagonal matrix with 22's along the main diagonal, and i's along the subdiagonal and the superdiagonal (i is the imaginary unit). - John M. Campbell, Jul 08 2011
For n>=2, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,...,21}. - Milan Janjic, Jan 25 2015
FORMULA
a(n) = 22*a(n-1) - a(n-1), a(-1)=0, a(0)=1.
a(n) = S(n, 22) with S(n, x) := U(n, x/2) Chebyshev's polynomials of the 2nd kind. See A049310.
a(n) = (ap^(n+1) - am^(n+1))/(ap - am) with ap := 11+2*sqrt(30) and am := 11-2*sqrt(30).
a(n) = Sum_{k=0..floor(n/2)} (-1)^k*binomial(n-k, k)*22^(n-2*k).
a(n) = sqrt((A077422(n+1)^2-1)/30)/2.
G.f.: 1/(1-22*x+x^2). - Philippe Deléham, Nov 18 2008
a(n) = Sum_{k, 0<=k<=n} A101950(n,k)*21^k. - Philippe Deléham, Feb 10 2012
Product {n >= 0} (1 + 1/a(n)) = 1/5*(5 + sqrt(30)). - Peter Bala, Dec 23 2012
Product {n >= 1} (1 - 1/a(n)) = 1/11*(5 + sqrt(30)). - Peter Bala, Dec 23 2012
MAPLE
seq( simplify(ChebyshevU(n, 11)), n=0..20); # G. C. Greubel, Dec 23 2019
MATHEMATICA
Table[GegenbauerC[n, 1, 11], {n, 0, 20}] (* Vladimir Joseph Stephan Orlovsky, Sep 11 2008 *)
CoefficientList[Series[1/(1-22x+x^2), {x, 0, 20}], x] (* Vincenzo Librandi, Dec 24 2012 *)
ChebyshevU[Range[21] -1, 11] (* G. C. Greubel, Dec 23 2019 *)
PROG
(Sage) [lucas_number1(n, 22, 1) for n in range(1, 20)] # Zerinvary Lajos, Jun 25 2008
(Sage) [chebyshev_U(n, 11) for n in (0..20)] # G. C. Greubel, Dec 23 2019
(Magma) I:=[1, 22]; [n le 2 select I[n] else 22*Self(n-1)-Self(n-2): n in [1..20]]; // Vincenzo Librandi, Dec 24 2012
(PARI) vector( 21, n, polchebyshev(n-1, 2, 11) ) \\ G. C. Greubel, Dec 23 2019
(GAP) m:=11;; a:=[1, 2*m];; for n in [3..20] do a[n]:=2*m*a[n-1]-a[n-2]; od; a; # G. C. Greubel, Dec 23 2019
CROSSREFS
Chebyshev sequence U(n, m): A000027 (m=1), A001353 (m=2), A001109 (m=3), A001090 (m=4), A004189 (m=5), A004191 (m=6), A007655 (m=7), A077412 (m=8), A049660 (m=9), A075843 (m=10), this sequence (m=11), A077423 (m=12), A097309 (m=13), A097311 (m=14), A097313 (m=15), A029548 (m=16), A029547 (m=17), A144128 (m=18), A078987 (m=19), A097316 (m=33).
Cf. A323182.
KEYWORD
nonn,easy
AUTHOR
Wolfdieter Lang, Nov 29 2002
STATUS
approved
Numbers k such that 99*k^2 + 1 is a square.
+10
26
0, 1, 20, 399, 7960, 158801, 3168060, 63202399, 1260879920, 25154396001, 501827040100, 10011386405999, 199725901079880, 3984506635191601, 79490406802752140, 1585823629419851199, 31636982181594271840
OFFSET
0,3
COMMENTS
From Wolfdieter Lang, Nov 08 2002: (Start)
Chebyshev's polynomials U(n,x) evaluated at x=10.
The a(n) give all (unsigned, integer) solutions of Pell equation b(n)^2 - 99*a(n)^2 = +1 with b(n)= A001085(n). (End)
For n>=2, a(n) equals the permanent of the (n-1) X (n-1) tridiagonal matrix with 20's along the main diagonal, and i's along the superdiagonal and the subdiagonal (i is the imagianry unit). - John M. Campbell, Jul 08 2011
For n>=1, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,...,19}. - Milan Janjic, Jan 25 2015
REFERENCES
A. H. Beiler, "The Pellian", ch. 22 in Recreations in the Theory of Numbers: The Queen of Mathematics Entertains. Dover, New York, New York, pp. 248-268, 1966.
L. E. Dickson, History of the Theory of Numbers, Vol. II, Diophantine Analysis. AMS Chelsea Publishing, Providence, Rhode Island, 1999, pp. 341-400.
Peter G. L. Dirichlet, Lectures on Number Theory (History of Mathematics Source Series, V. 16); American Mathematical Society, Providence, Rhode Island, 1999, pp. 139-147.
LINKS
Tanya Khovanova, Recursive Sequences
J. J. O'Connor and E. F. Robertson, Pell's Equation
Eric Weisstein's World of Mathematics, Pell Equation.
FORMULA
a(n) = ((10+3*sqrt(11))^n - (10-3*sqrt(11))^n) / (6*sqrt(11)).
a(n) = 20*a(n-1) - a(n-2), n>=1, a(0)=0, a(1)=1.
a(n) = S(n-1, 20), with S(n, x) := U(n, x/2), Chebyshev's polynomials of the second kind. S(-1, x) := 0. See A049310.
G.f.: x/(1 - 20*x + x^2).
a(n) = sqrt((A001085(n)^2 - 1)/99).
Lim_{n->inf.} a(n)/a(n-1) = 10 + 3*sqrt(11).
a(n+1) = Sum_{k=0..n} A101950(n,k)*19^k. - Philippe Deléham, Feb 10 2012
Product_{n>=1} (1 + 1/a(n)) = 1/3*(3 + sqrt(11)). - Peter Bala, Dec 23 2012
Product_{n>=2} (1 - 1/a(n)) = 3/20*(3 + sqrt(11)). - Peter Bala, Dec 23 2012
MAPLE
seq( simplify(ChebyshevU(n-1, 10)), n=0..20); # G. C. Greubel, Dec 22 2019
MATHEMATICA
Table[GegenbauerC[n-1, 1, 10], {n, 0, 20}] (* Vladimir Joseph Stephan Orlovsky, Sep 11 2008 *)
CoefficientList[Series[x/(1-20x+x^2), {x, 0, 20}], x] (* Vincenzo Librandi, Dec 24 2012 *)
ChebyshevU[Range[22] -2, 10] (* G. C. Greubel, Dec 22 2019 *)
LinearRecurrence[{20, -1}, {0, 1}, 20] (* Harvey P. Dale, Dec 03 2023 *)
PROG
(Sage) [lucas_number1(n, 20, 1) for n in range(0, 20)] # Zerinvary Lajos, Jun 25 2008
(Sage) [chebyshev_U(n-1, 10) for n in (0..20)] # G. C. Greubel, Dec 22 2019
(Magma) I:=[0, 1]; [n le 2 select I[n] else 20*Self(n-1)-Self(n-2): n in [1..20]]; // Vincenzo Librandi, Dec 24 2012
(PARI) vector( 22, n, polchebyshev(n-2, 2, 10) ) \\ G. C. Greubel, Dec 22 2019
(GAP) m:=10;; a:=[0, 1];; for n in [3..20] do a[n]:=2*m*a[n-1]-a[n-2]; od; a; # G. C. Greubel, Dec 22 2019
CROSSREFS
Cf. A001084.
Chebyshev sequence U(n, m): A000027 (m=1), A001353 (m=2), A001109 (m=3), A001090 (m=4), A004189 (m=5), A004191 (m=6), A007655 (m=7), A077412 (m=8), A049660 (m=9), this sequence (m=10), A077421 (m=11), A077423 (m=12), A097309 (m=13), A097311 (m=14), A097313 (m=15), A029548 (m=16), A029547 (m=17), A144128 (m=18), A078987 (m=19), A097316 (m=33).
Cf. A323182.
KEYWORD
nonn,easy
AUTHOR
STATUS
approved

Search completed in 0.026 seconds