cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 18 results. Next

A002415 4-dimensional pyramidal numbers: a(n) = n^2*(n^2-1)/12.

Original entry on oeis.org

0, 0, 1, 6, 20, 50, 105, 196, 336, 540, 825, 1210, 1716, 2366, 3185, 4200, 5440, 6936, 8721, 10830, 13300, 16170, 19481, 23276, 27600, 32500, 38025, 44226, 51156, 58870, 67425, 76880, 87296, 98736, 111265, 124950, 139860, 156066, 173641, 192660, 213200, 235340
Offset: 0

Views

Author

Keywords

Comments

Also number of ways to legally insert two pairs of parentheses into a string of m := n-1 letters. (There are initially 2C(m+4,4) (A034827) ways to insert the parentheses, but we must subtract 2(m+1) for illegal clumps of 4 parentheses, 2m(m+1) for clumps of 3 parentheses, C(m+1,2) for 2 clumps of 2 parentheses and (m-1)C(m+1,2) for 1 clump of 2 parentheses, giving m(m+1)^2(m+2)/12 = n^2*(n^2-1)/12.) See also A000217.
E.g., for n=2 there are 6 ways: ((a))b, ((a)b), ((ab)), (a)(b), (a(b)), a((b)).
Let M_n denote the n X n matrix M_n(i,j)=(i+j); then the characteristic polynomial of M_n is x^(n-2) * (x^2-A002378(n)*x - a(n)). - Benoit Cloitre, Nov 09 2002
Let M_n denote the n X n matrix M_n(i,j)=(i-j); then the characteristic polynomial of M_n is x^n + a(n)x^(n-2). - Michael Somos, Nov 14 2002 [See A114327 for the infinite matrix M in triangular form. - Wolfdieter Lang, Feb 05 2018]
Number of permutations of [n] which avoid the pattern 132 and have exactly 2 descents. - Mike Zabrocki, Aug 26 2004
Number of tilings of a <2,n,2> hexagon.
a(n) is the number of squares of side length at least 1 having vertices at the points of an n X n unit grid of points (the vertices of an n-1 X n-1 chessboard). [For a proof, see Comments in A051602. - N. J. A. Sloane, Sep 29 2021] For example, on the 3 X 3 grid (the vertices of a 2 X 2 chessboard) there are four 1 X 1 squares, one (skew) sqrt(2) X sqrt(2) square, and one 3 X 3 square, so a(3)=6. On the 4 X 4 grid (the vertices of a 3 X 3 chessboard) there are 9 1 X 1 squares, 4 2 X 2 squares, 1 3 X 3 square, 4 sqrt(2) X sqrt(2) squares, and 2 sqrt(5) X sqrt(5) squares, so a(4) = 20. See also A024206, A108279. [Comment revised by N. J. A. Sloane, Feb 11 2015]
Kekulé numbers for certain benzenoids. - Emeric Deutsch, Jun 12 2005
Number of distinct components of the Riemann curvature tensor. - Gene Ward Smith, Apr 24 2006
a(n) is the number of 4 X 4 matrices (symmetrical about each diagonal) M = [a,b,c,d;b,e,f,c;c,f,e,b;d,c,b,a] with a+b+c+d=b+e+f+c=n+2; (a,b,c,d,e,f natural numbers). - Philippe Deléham, Apr 11 2007
If a 2-set Y and an (n-2)-set Z are disjoint subsets of an n-set X then a(n-3) is the number of 5-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 19 2007
a(n) is the number of Dyck (n+1)-paths with exactly n-1 peaks. - David Callan, Sep 20 2007
Starting (1,6,20,50,...) = third partial sums of binomial transform of [1,2,0,0,0,...]. a(n) = Sum_{i=0..n} C(n+3,i+3)*b(i), where b(i)=[1,2,0,0,0,...]. - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009
4-dimensional square numbers. - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009
Equals row sums of triangle A177877; a(n), n > 1 = (n-1) terms in (1,2,3,...) dot (...,3,2,1) with additive carryovers. Example: a(4) = 20 = (1,2,3) dot (3,2,1) with carryovers = (1*3) + (2*2 + 3) + (3*1 + 7) = (3 + 7 + 10).
Convolution of the triangular numbers A000217 with the odd numbers A004273.
a(n+2) is the number of 4-tuples (w,x,y,z) with all terms in {0,...,n} and w-x=max{w,x,y,z}-min{w,x,y,z}. - Clark Kimberling, May 28 2012
The second level of finite differences is a(n+2) - 2*a(n+1) + a(n) = (n+1)^2, the squares. - J. M. Bergot, May 29 2012
Because the differences of this sequence give A000330, this is also the number of squares in an n+1 X n+1 grid whose sides are not parallel to the axes.
a(n+2) gives the number of 2*2 arrays that can be populated with 0..n such that rows and columns are nondecreasing. - Jon Perry, Mar 30 2013
For n consecutive numbers 1,2,3,...,n, the sum of all ways of adding the k-tuples of consecutive numbers for n=a(n+1). As an example, let n=4: (1)+(2)+(3)+(4)=10; (1+2)+(2+3)+(3+4)=15; (1+2+3)+(2+3+4)=15; (1+2+3+4)=10 and the sum of these is 50=a(4+1)=a(5). - J. M. Bergot, Apr 19 2013
If P(n,k) = n*(n+1)*(k*n-k+3)/6 is the n-th (k+2)-gonal pyramidal number, then a(n) = P(n,k)*P(n-1,k-1) - P(n-1,k)*P(n,k-1). - Bruno Berselli, Feb 18 2014
For n > 1, a(n) = 1/6 of the area of the trapezoid created by the points (n,n+1), (n+1,n), (1,n^2+n), (n^2+n,1). - J. M. Bergot, May 14 2014
For n > 3, a(n) is twice the area of a triangle with vertices at points (C(n,4),C(n+1,4)), (C(n+1,4),C(n+2,4)), and (C(n+2,4),C(n+3,4)). - J. M. Bergot, Jun 03 2014
a(n) is the dimension of the space of metric curvature tensors (those having the symmetries of the Riemann curvature tensor of a metric) on an n-dimensional real vector space. - Daniel J. F. Fox, Dec 15 2018
Coefficients in the terminating series identity 1 - 6*n/(n + 5) + 20*n*(n - 1)/((n + 5)*(n + 6)) - 50*n*(n - 1)*(n - 2)/((n + 5)*(n + 6)*(n + 7)) + ... = 0 for n = 1,2,3,.... Cf. A000330 and A005585. - Peter Bala, Feb 18 2019

Examples

			a(7) = 6*21 - (6*0 + 4*1 + 2*3 + 0*6 - 2*10 - 4*15) = 196. - _Bruno Berselli_, Jun 22 2013
G.f. = x^2 + 6*x^3 + 20*x^4 + 50*x^5 + 105*x^6 + 196*x^7 + 336*x^8 + ...
		

References

  • O. D. Anderson, Find the next sequence, J. Rec. Math., 8 (No. 4, 1975-1976), 241.
  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 195.
  • S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (p.165).
  • R. Euler and J. Sadek, "The Number of Squares on a Geoboard", Journal of Recreational Mathematics, 251-5 30(4) 1999-2000 Baywood Pub. NY
  • S. Mukai, An Introduction to Invariants and Moduli, Cambridge, 2003; see p. 238.
  • 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).

Crossrefs

a(n) = ((-1)^n)*A053120(2*n, 4)/8 (one-eighth of fifth unsigned column of Chebyshev T-triangle, zeros omitted). Cf. A001296.
Second row of array A103905.
Third column of Narayana numbers A001263.
Partial sums of A000330.
The expression binomial(m+n-1,n)^2-binomial(m+n,n+1)*binomial(m+n-2,n-1) for the values m = 2 through 14 produces sequences A000012, A000217, A002415, A006542, A006857, A108679, A134288, A134289, A134290, A134291, A140925, A140935, A169937.
Cf. A220212 for a list of sequences produced by the convolution of the natural numbers (A000027) with the k-gonal numbers.

Programs

  • GAP
    List([0..45],n->Binomial(n^2,2)/6); # Muniru A Asiru, Dec 15 2018
  • Magma
    [n^2*(n^2-1)/12: n in [0..50]]; // Wesley Ivan Hurt, May 14 2014
    
  • Maple
    A002415 := proc(n) binomial(n^2,2)/6 ; end proc: # Zerinvary Lajos, Jan 07 2008
  • Mathematica
    Table[(n^4 - n^2)/12, {n, 0, 40}] (* Zerinvary Lajos, Mar 21 2007 *)
    LinearRecurrence[{5,-10,10,-5,1},{0,0,1,6,20},40] (* Harvey P. Dale, Nov 29 2011 *)
  • PARI
    a(n) = n^2 * (n^2 - 1) / 12;
    
  • PARI
    x='x+O('x^200); concat([0, 0], Vec(x^2*(1+x)/(1-x)^5)) \\ Altug Alkan, Mar 23 2016
    

Formula

G.f.: x^2*(1+x)/(1-x)^5. - Simon Plouffe in his 1992 dissertation
a(n) = Sum_{i=0..n} (n-i)*i^2 = a(n-1) + A000330(n-1) = A000217(n)*A000292(n-2)/n = A000217(n)*A000217(n-1)/3 = A006011(n-1)/3, convolution of the natural numbers with the squares. - Henry Bottomley, Oct 19 2000
a(n)+1 = A079034(n). - Mario Catalani (mario.catalani(AT)unito.it), Feb 12 2003
a(n) = 2*C(n+2, 4) - C(n+1, 3). - Paul Barry, Mar 04 2003
a(n) = C(n+2, 4) + C(n+1, 4). - Paul Barry, Mar 13 2003
a(n) = Sum_{k=1..n} A000330(n-1). - Benoit Cloitre, Jun 15 2003
a(n) = n*C(n+1,3)/2 = C(n+1,3)*C(n+1,2)/(n+1). - Mitch Harris, Jul 06 2006
a(n) = A006011(n)/3 = A008911(n)/2 = A047928(n-1)/12 = A083374(n)/6. - Zerinvary Lajos, May 09 2007
a(n) = (1/2)*Sum_{1 <= x_1, x_2 <= n} (det V(x_1,x_2))^2 = (1/2)*Sum_{1 <= i,j <= n} (i-j)^2, where V(x_1,x_2) is the Vandermonde matrix of order 2. - Peter Bala, Sep 21 2007
a(n) = C(n+1,3) + 2*C(n+1,4). - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009
a(n) = (1/48)*sinh(2*arccosh(n))^2. - Artur Jasinski, Feb 10 2010
a(n) = n*A000292(n-1)/2. - Tom Copeland, Sep 13 2011
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5), n > 4. - Harvey P. Dale, Nov 29 2011
a(n) = (n-1)*A000217(n-1) - Sum_{i=0..n-2} (n-1-2*i)*A000217(i) for n > 1. - Bruno Berselli, Jun 22 2013
a(n) = C(n,2)*C(n+1,3) - C(n,3)*C(n+1,2). - J. M. Bergot, Sep 17 2013
a(n) = Sum_{k=1..n} ( (2k-n)* k(k+1)/2 ). - Wesley Ivan Hurt, Sep 26 2013
a(n) = floor(n^2/3) + 3*Sum_{k=1..n} k^2*floor((n-k+1)/3). - Mircea Merca, Feb 06 2014
Euler transform of length 2 sequence [6, -1]. - Michael Somos, May 28 2014
G.f. x^2*2F1(3,4;2;x). - R. J. Mathar, Aug 09 2015
Sum_{n>=2} 1/a(n) = 21 - 2*Pi^2 = 1.260791197821282762331... . - Vaclav Kotesovec, Apr 27 2016
a(n) = A080852(2,n-2). - R. J. Mathar, Jul 28 2016
a(n) = A046092(n) * A046092(n-1)/48 = A000217(n) * A000217(n-1)/3. - Bruce J. Nicholson, Jun 06 2017
E.g.f.: (1/12)*exp(x)*x^2*(6 + 6*x + x^2). - Stefano Spezia, Dec 07 2018
Sum_{n>=2} (-1)^n/a(n) = Pi^2 - 9 (See A002388). - Amiram Eldar, Jun 28 2020

Extensions

Typo in link fixed by Matthew Vandermast, Nov 22 2010
Redundant comment deleted and more detail on relationship with A000330 added by Joshua Zucker, Jan 01 2013

A058331 a(n) = 2*n^2 + 1.

Original entry on oeis.org

1, 3, 9, 19, 33, 51, 73, 99, 129, 163, 201, 243, 289, 339, 393, 451, 513, 579, 649, 723, 801, 883, 969, 1059, 1153, 1251, 1353, 1459, 1569, 1683, 1801, 1923, 2049, 2179, 2313, 2451, 2593, 2739, 2889, 3043, 3201, 3363, 3529, 3699, 3873, 4051
Offset: 0

Views

Author

Erich Friedman, Dec 12 2000

Keywords

Comments

Maximal number of regions in the plane that can be formed with n hyperbolas.
Also the number of different 2 X 2 determinants with integer entries from 0 to n.
Number of lattice points in an n-dimensional ball of radius sqrt(2). - David W. Wilson, May 03 2001
Equals A112295(unsigned) * [1, 2, 3, ...]. - Gary W. Adamson, Oct 07 2007
Binomial transform of A166926. - Gary W. Adamson, May 03 2008
a(n) = longest side a of all integer-sided triangles with sides a <= b <= c and inradius n >= 1. Triangle has sides (2n^2 + 1, 2n^2 + 2, 4n^2 + 1).
{a(k): 0 <= k < 3} = divisors of 9. - Reinhard Zumkeller, Jun 17 2009
Number of ways to partition a 3*n X 2 grid into 3 connected equal-area regions. - R. H. Hardin, Oct 31 2009
Let A be the Hessenberg matrix of order n defined by: A[1, j] = 1, A[i, i] := 2, (i > 1), A[i, i - 1] = -1, and A[i, j] = 0 otherwise. Then, for n >= 3, a(n - 1) = coeff(charpoly(A, x), x^(n - 2)). - Milan Janjic, Jan 26 2010
Except for the first term of [A002522] and [A058331] if X = [A058331], Y = [A087113], A = [A002522], we have, for all other terms, Pell's equation: [A058331]^2 - [A002522]*[A087113]^2 = 1; (X^2 - A*Y^2 = 1); e.g., 3^2 -2*2^2 = 1; 9^2 - 5*4^2 = 1; 129^2 - 65*16^2 = 1, and so on. - Vincenzo Librandi, Aug 07 2010
Niven (1961) gives this formula as an example of a formula that does not contain all odd integers, in contrast to 2n + 1 and 2n - 1. - Alonso del Arte, Dec 05 2012
Numbers m such that 2*m-2 is a square. - Vincenzo Librandi, Apr 10 2015
Number of n-tuples from the set {1,0,-1} where at most two elements are nonzero. - Michael Somos, Oct 19 2022
a(n) gives the x-value of the integral solution (x,y) of the Pellian equation x^2 - (n^2 + 1)*y^2 = 1. The y-value is given by 2*n (see Tattersall). - Stefano Spezia, Jul 23 2025

Examples

			a(1) = 3 since (0 0 / 0 0), (1 0 / 0 1) and (0 1 / 1 0) have different determinants.
G.f. = 1 + 3*x + 9*x^2 + 19*x^3 + 33*x^4 + 51*x^5 + 73*x^6 + ... - _Michael Somos_, Oct 19 2022
		

References

  • Ivan Niven, Numbers: Rational and Irrational, New York: Random House for Yale University (1961): 17.
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 256.

Crossrefs

Cf. A000124.
Second row of array A099597.
See A120062 for sequences related to integer-sided triangles with integer inradius n.
Cf. A112295.
Column 2 of array A188645.
Cf. A001105 and A247375. - Bruno Berselli, Sep 16 2014

Programs

  • Haskell
    a058331 = (+ 1) . a001105  -- Reinhard Zumkeller, Dec 13 2014
    
  • Magma
    [2*n^2 + 1 : n in [0..100]]; // Wesley Ivan Hurt, Feb 02 2017
  • Mathematica
    b[g_] := Length[Union[Map[Det, Flatten[ Table[{{i, j}, {k, l}}, {i, 0, g}, {j, 0, g}, {k, 0, g}, {l, 0, g}], 3]]]] Table[b[g], {g, 0, 20}]
    2*Range[0, 49]^2 + 1 (* Alonso del Arte, Dec 05 2012 *)
  • PARI
    a(n)=2*n^2+1 \\ Charles R Greathouse IV, Jun 16 2011
    

Formula

G.f.: (1 + 3x^2)/(1 - x)^3. - Paul Barry, Apr 06 2003
a(n) = M^n * [1 1 1], leftmost term, where M = the 3 X 3 matrix [1 1 1 / 0 1 4 / 0 0 1]. a(0) = 1, a(1) = 3; a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). E.g., a(4) = 33 since M^4 *[1 1 1] = [33 17 1]. - Gary W. Adamson, Nov 11 2004
a(n) = cosh(2*arccosh(n)). - Artur Jasinski, Feb 10 2010
a(n) = 4*n + a(n-1) - 2 for n > 0, a(0) = 1. - Vincenzo Librandi, Aug 07 2010
a(n) = (((n-1)^2 + n^2))/2 + (n^2 + (n+1)^2)/2. - J. M. Bergot, May 31 2012
a(n) = A251599(3*n) for n > 0. - Reinhard Zumkeller, Dec 13 2014
a(n) = sqrt(8*(A000217(n-1)^2 + A000217(n)^2) + 1). - J. M. Bergot, Sep 03 2015
E.g.f.: (2*x^2 + 2*x + 1)*exp(x). - G. C. Greubel, Jul 14 2017
a(n) = A002378(n) + A002061(n). - Bruce J. Nicholson, Aug 06 2017
From Amiram Eldar, Jul 15 2020: (Start)
Sum_{n>=0} 1/a(n) = (1 + (Pi/sqrt(2))*coth(Pi/sqrt(2)))/2.
Sum_{n>=0} (-1)^n/a(n) = (1 + (Pi/sqrt(2))*csch(Pi/sqrt(2)))/2. (End)
From Amiram Eldar, Feb 05 2021: (Start)
Product_{n>=0} (1 + 1/a(n)) = sqrt(2)*csch(Pi/sqrt(2))*sinh(Pi).
Product_{n>=1} (1 - 1/a(n)) = (Pi/sqrt(2))*csch(Pi/sqrt(2)). (End)
From Leo Tavares, May 23 2022: (Start)
a(n) = A000384(n+1) - 3*n.
a(n) = 3*A000217(n) + A000217(n-2). (End)
a(n) = a(-n) for all n in Z and A037235(n) = Sum_{k=0..n-1} a(k). - Michael Somos, Oct 19 2022

Extensions

Revised description from Noam Katz (noamkj(AT)hotmail.com), Jan 28 2001

A132592 X-values of solutions to the equation X*(X + 1) - 8*Y^2 = 0.

Original entry on oeis.org

0, 8, 288, 9800, 332928, 11309768, 384199200, 13051463048, 443365544448, 15061377048200, 511643454094368, 17380816062160328, 590436102659356800, 20057446674355970888, 681362750825443653408, 23146276081390728245000, 786292024016459316676608, 26710782540478226038759688
Offset: 0

Views

Author

Mohamed Bouhamida, Nov 14 2007

Keywords

Comments

Equivalently, numbers k such that both k/2 and k+1 are squares. - Karl-Heinz Hofmann, Sep 20 2022
Equivalently, numbers k such that the k-dimensional volume and total (k-1)-dimensional volume are equal, with side length being a positive integer, for all regular polyhedra constructible in k dimensions. - Matt Moir, Jul 09 2024

Crossrefs

Programs

Formula

a(0)=0, a(1)=8 and a(n) = 34*a(n-1) - a(n-2) + 16.
a(n) = (A056771(n) - 1)/2. - Max Alekseyev, Nov 13 2009
a(n) = sinh(2*n*arccosh(sqrt(2))^2) (n=0,1,2,3,...). - Artur Jasinski, Feb 10 2010
G.f.: -8*x*(x+1)/((x-1)*(x^2-34*x+1)). - Colin Barker, Oct 24 2012
a(n) = A055792(n+1)-1 = A001541(n)^2 - 1. - Antti Karttunen, Oct 03 2016

Extensions

More terms from Max Alekseyev, Nov 13 2009

A173129 a(n) = cosh(2 * n * arccosh(n)).

Original entry on oeis.org

1, 1, 97, 19601, 7380481, 4517251249, 4097989415521, 5170128475599457, 8661355881006882817, 18605234632923999244961, 49862414878754347585980001, 163104845048002042971670685041, 639582975902942936737758325440001
Offset: 0

Views

Author

Artur Jasinski, Feb 10 2010

Keywords

Crossrefs

Programs

  • Maple
    seq(orthopoly[T](2*n,n), n=0..50); # Robert Israel, Dec 27 2018
  • Mathematica
    Table[Round[Cosh[2 n ArcCosh[n]]], {n, 0, 20}] (* Artur Jasinski, Feb 10 2010 *)
    Round[Table[1/2 (x - Sqrt[ -1 + x^2])^(2 x) + 1/2 (x + Sqrt[ -1 + x^2])^(2 x), {x, 0, 10}]] (* Artur Jasinski, Feb 14 2010 *)
    Table[ChebyshevT[2*n, n], {n, 0, 15}] (* Vaclav Kotesovec, Nov 07 2021 *)
  • PARI
    {a(n) = sum(k=0, n, binomial(2*n, 2*k)*(n^2-1)^(n-k)*n^(2*k))} \\ Seiichi Manyama, Dec 27 2018
    
  • PARI
    {a(n) = polchebyshev(2*n, 1, n)} \\ Seiichi Manyama, Dec 28 2018
    
  • PARI
    {a(n) = polchebyshev(n, 1, 2*n^2-1)} \\ Seiichi Manyama, Dec 29 2018

Formula

a(n) = (1/2)*((n+sqrt(n^2-1))^(2*n) + (n-sqrt(n^2-1))^(2*n)). - Artur Jasinski, Feb 14 2010, corrected by Vaclav Kotesovec, Apr 05 2016
a(n) = Sum_{k=0..n} binomial(2*n,2*k)*(n^2-1)^(n-k)*n^(2*k). - Seiichi Manyama, Dec 27 2018
a(n) = T_{2n}(n) where T_{2n} is a Chebyshev polynomial of the first kind. - Robert Israel, Dec 27 2018
a(n) = T_{n}(2*n^2-1) where T_{n}(x) is a Chebyshev polynomial of the first kind. - Seiichi Manyama, Dec 29 2018

A173127 a(n) = sinh((2n-1)*arcsinh(3)).

Original entry on oeis.org

-3, 3, 117, 4443, 168717, 6406803, 243289797, 9238605483, 350823718557, 13322062699683, 505887558869397, 19210405174337403, 729489509065951917, 27701390939331835443, 1051923366185543794917, 39945386524111332371403
Offset: 0

Views

Author

Artur Jasinski, Feb 10 2010

Keywords

Comments

Numbers n such that ((n^2 + 1)/10) is a square. - Vincenzo Librandi, Jan 02 2012

Crossrefs

Programs

  • Magma
    [-3] cat [n: n in [0..10^7]|IsSquare((n^2+1)/10)]; // Vincenzo Librandi, Jan 02 2012
    
  • Mathematica
    LinearRecurrence[{38,-1},{-3,3},30] (* Harvey P. Dale, Jan 14 2015 *)
  • Python
    from itertools import islice
    def A173127_gen(): # generator of terms
        x, y = -30, 10
        while True:
            yield x//10
            x, y = x*19+y*60, x*6+y*19
    A173127_list = list(islice(A173127_gen(),20)) # Chai Wah Wu, Apr 24 2025

Formula

a(n) = (1/2)*((-3+sqrt(10))*(19+6*sqrt(10))^n + (-3-sqrt(10))*(19-6*sqrt(10))^n).
a(n) = -a(-n+1).
G.f.: -3*(1-39*x)/(1-38*x+x^2). - Bruno Berselli, Jan 03 2011
E.g.f.: exp(19*x)*(-3*cosh(6*sqrt(10)*x) + sqrt(10)*sinh(6*sqrt(10)*x)). - Stefano Spezia, Apr 24 2025

A173128 a(n) = cosh(2*n*arcsinh(n)).

Original entry on oeis.org

1, 3, 161, 27379, 9478657, 5517751251, 4841332221601, 5964153172084899, 9814664424981012481, 20791777842234580902499, 55106605639755476546020001, 178627672869645203363556318483, 695165908550906808156689590141441
Offset: 0

Views

Author

Artur Jasinski, Feb 10 2010

Keywords

Crossrefs

Programs

  • Maple
    seq(expand( (1/2)*((n + sqrt(n^2 + 1))^(2*n) + (n - sqrt(n^2 + 1))^(2*n))), n=0..30); # Robert Israel, Apr 05 2016
  • Mathematica
    Round[Table[Cosh[2 n ArcSinh[n]], {n, 0, 20}]] (* Artur Jasinski *)
    Round[Table[1/2 (x - Sqrt[1 + x^2])^(2 x) + 1/2 (x + Sqrt[1 + x^2])^(2 x), {x, 0, 20}]] (* Artur Jasinski, Feb 14 2010 *)
  • PARI
    {a(n) = sum(k=0, n, binomial(2*n, 2*k)*(n^2+1)^(n-k)*n^(2*k))} \\ Seiichi Manyama, Dec 27 2018
    
  • PARI
    {a(n) = polchebyshev(n, 1, 2*n^2+1)} \\ Seiichi Manyama, Dec 29 2018

Formula

a(n) = (1/2)*((n + sqrt(n^2 + 1))^(2*n) + (n - sqrt(n^2 + 1))^(2*n)). - Artur Jasinski, Feb 14 2010, corrected by Vaclav Kotesovec, Apr 05 2016
a(n) = Sum_{k=0..n} binomial(2*n,2*k)*(n^2+1)^(n-k)*n^(2*k). - Seiichi Manyama, Dec 27 2018
a(n) = T_{n}(2*n^2+1) where T_{n}(x) is a Chebyshev polynomial of the first kind. - Seiichi Manyama, Dec 29 2018

A173130 a(n) = Cosh[(2 n - 1) ArcCosh[n]].

Original entry on oeis.org

0, 1, 26, 3363, 937444, 456335045, 343904160606, 371198523608647, 543466014742175624, 1036834190110356583689, 2499384905955651114739810, 7429238104512325157021090411, 26695718139185294187938997247212
Offset: 0

Views

Author

Artur Jasinski, Feb 10 2010

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Round[Cosh[(2 n - 1) ArcCosh[n]]], {n, 0, 20}] (* Artur Jasinski *)

Formula

a(n) ~ 2^(2*n-2) * n^(2*n-1). - Vaclav Kotesovec, Apr 05 2016

A173131 a(n) = (Cosh[(2n-1)ArcSinh[n]])^2.

Original entry on oeis.org

1, 2, 1445, 19740250, 1361599599377, 298514762397852026, 160545187370375075046277, 179656719395983409634002348450, 373368546362937441101158606899394625
Offset: 0

Views

Author

Artur Jasinski, Feb 10 2010

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Round[Cosh[(2 n - 1) ArcSinh[n]]^2], {n, 0, 10}] (* Artur Jasinski *)

Formula

a(n) ~ 2^(4*n-4) * n^(4*n-2). - Vaclav Kotesovec, Apr 05 2016

A173148 a(n) = cos(2*n*arccos(sqrt(n))).

Original entry on oeis.org

1, 1, 17, 485, 18817, 930249, 55989361, 3974443213, 325142092801, 30122754096401, 3117419602578001, 356452534779818421, 44627167107085622401, 6071840759403431812825, 892064955046043465408177, 140751338790698080509966749, 23737154316161495960243527681
Offset: 0

Views

Author

Artur Jasinski, Feb 11 2010

Keywords

Comments

The Chebyshev polynomial T_n is defined by cos(nx) = T_n(cos(x)). So T_2n(cos(x)) = cos(2nx) = cos^2(nx) - 1 = (T_n(x))^2 - 1 consists of only even powers of x. As a result, a(n) = T_2n(sqrt(n)) is an integer. - Michael B. Porter, Jan 01 2019

Crossrefs

Programs

  • GAP
    a:=List([0..20],n->Sum([0..n],k->Binomial(2*n,2*k)*(n-1)^(n-k)*n^k));; Print(a); # Muniru A Asiru, Jan 03 2019
    
  • Magma
    [&+[Binomial(2*n,2*k)*(n-1)^(n-k)*n^k: k in [0..n]]: n in [0..20]]; // Vincenzo Librandi, Jan 03 2019
  • Mathematica
    Table[Round[Cos[2 n ArcCos[Sqrt[n]]]], {n, 0, 30}] (* Artur Jasinski, Feb 11 2010 *)
  • PARI
    {a(n) = sum(k=0, n, binomial(2*n, 2*k)*(n-1)^(n-k)*n^k)} \\ Seiichi Manyama, Dec 27 2018
    
  • PARI
    {a(n) = round(cosh(2*n*acosh(sqrt(n))))} \\ Seiichi Manyama, Dec 27 2018
    
  • PARI
    {a(n) = polchebyshev(n, 1, 2*n-1)} \\ Seiichi Manyama, Dec 29 2018
    

Formula

a(n) ~ exp(-1/2) * 2^(2*n-1) * n^n. - Vaclav Kotesovec, Apr 05 2016
a(n) = Sum_{k=0..n} binomial(2*n,2*k)*(n-1)^(n-k)*n^k. - Seiichi Manyama, Dec 27 2018
a(n) = cosh(2*n*arccosh(sqrt(n))). - Seiichi Manyama, Dec 27 2018
a(n) = T_{2*n}(sqrt(n)) = T_{n}(2*n-1) where T_{n}(x) is a Chebyshev polynomial of the first kind. - Seiichi Manyama, Dec 29 2018
a(n) = A322790(n-1, n) for n > 0. - Seiichi Manyama, Dec 29 2018

Extensions

Minor edits by Vaclav Kotesovec, Apr 05 2016

A173133 a(n) = Sinh[(2n-1) ArcSinh[n]].

Original entry on oeis.org

0, 1, 38, 4443, 1166876, 546365045, 400680904674, 423859315570607, 611038907405197432, 1151555487914640463209, 2748476184146759127540190, 8102732939160371170806346243, 28915133156938367486730067779348
Offset: 0

Views

Author

Artur Jasinski, Feb 10 2010

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Round[Sinh[(2 n - 1) ArcSinh[n]]], {n, 0, 20}] (* Artur Jasinski *)
    Round[Table[1/2 (n - Sqrt[1 + n^2])^(2 n - 1) + 1/2 (n + Sqrt[1 + n^2])^(2 n - 1), {n, 0, 10}]] (* Artur Jasinski, Feb 14 2010 *)

Formula

a(n) = 1/2 (n - sqrt(1 + n^2))^(2 n - 1) + 1/2 (n + sqrt(1 + n^2))^(2 n - 1). - Artur Jasinski, Feb 14 2010

Extensions

Minor edits by Vaclav Kotesovec, Apr 05 2016
Showing 1-10 of 18 results. Next