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.

Previous Showing 11-20 of 60 results. Next

A063496 a(n) = (2*n - 1)*(8*n^2 - 8*n + 3)/3.

Original entry on oeis.org

1, 19, 85, 231, 489, 891, 1469, 2255, 3281, 4579, 6181, 8119, 10425, 13131, 16269, 19871, 23969, 28595, 33781, 39559, 45961, 53019, 60765, 69231, 78449, 88451, 99269, 110935, 123481, 136939, 151341, 166719, 183105, 200531, 219029, 238631, 259369, 281275, 304381
Offset: 1

Views

Author

N. J. A. Sloane, Aug 01 2001

Keywords

Comments

Number of potential flows in a 2 X 2 matrix with integer velocities in -n..n, i.e., number of 2 X 2 matrices with adjacent elements differing by no more than n, counting matrices differing by a constant only once. - R. H. Hardin, Feb 27 2002
Number of ordered quadruples (a,b,c,d), -(n-1) <= a,b,c,d <= n-1, such that a+b+c+d = 0. - Benoit Cloitre, Jun 14 2003
If Y and Z are 2-blocks of a (2n+1)-set X then a(n-1) is the number of 5-subsets of X intersecting both Y and Z. - Milan Janjic, Oct 28 2007
Equals binomial transform of [1, 18, 48, 32, 0, 0, 0, ...]. - Gary W. Adamson, Jul 19 2008

Crossrefs

(1/12)*t*(2*n^3-3*n^2+n)+2*n-1 for t = 2, 4, 6, ... gives A049480, A005894, A063488, A001845, A063489, A005898, A063490, A057813, A063491, A005902, A063492, A005917, A063493, A063494, A063495, A063496.

Programs

  • Magma
    [(2*n-1)*(8*n^2-8*n+3)/3: n in [1..40]]; // Wesley Ivan Hurt, May 09 2014
  • Maple
    A063496:=n->(2*n-1)*(8*n^2-8*n+3)/3; seq(A063496(n), n=1..40); # Wesley Ivan Hurt, May 09 2014
  • Mathematica
    Table[(2*n - 1)*(8*n^2 - 8*n + 3)/3, {n, 40}] (* Wesley Ivan Hurt, May 09 2014 *)
    LinearRecurrence[{4,-6,4,-1}, {1,19,85,231}, 30] (* G. C. Greubel, Dec 01 2017 *)
  • PARI
    a(n) = { (2*n - 1)*(8*n^2 - 8*n + 3)/3 } \\ Harry J. Smith, Aug 23 2009
    
  • PARI
    my(x='x+O('x^30)); Vec(serlaplace((-3+6*x+24*x^2+16*x^3)*exp(x)/3 + 1)) \\ G. C. Greubel, Dec 01 2017
    

Formula

From Peter Bala, Jul 18 2008: (Start)
The following remarks about the C_3 lattice assume the sequence offset is 0.
Partial sums of A010006. So this sequence is the crystal ball sequence for the C_3 lattice - row 3 of A142992. The lattice C_3 consists of all integer lattice points v = (a,b,c) in Z^3 such that a + b + c is even, equipped with the taxicab type norm ||v|| = (1/2) * (|a| + |b| + |c|).
The crystal ball sequence of C_3 gives the number of lattice points v in C_3 with ||v|| <= n for n = 0,1,2,3,... [Bacher et al.].
For example, a(1) = 19 because the origin has norm 0 and the 18 lattice points in Z^3 of norm 1 (as defined above) are +-(2,0,0), +-(0,2,0), +-(0,0,2), +-(1,1,0), +-(1,0,1), +-(0,1,1), +-(1,-1,0), +-(1,0,-1) and +-(0,1,-1). These 18 vectors form a root system of type C_3.
O.g.f.: x*(1 + 15*x + 15*x^2 + x^3)/(1 - x)^4 = x/(1 - x) * T(3, (1 + x)/(1 - x)), where T(n, x) denotes the Chebyshev polynomial of the first kind.
2*log(2) = 4/3 + Sum_{n >= 1} 1/(n*a(n)*a(n+1)). (End)
a(n+1) = (1/Pi) * Integral_{x=0..Pi} (sin((n+1/2)*x)/sin(x/2))^4. - Yalcin Aktar, Nov 02 2011, corrected by R. J. Mathar, Dec 01 2011
From G. C. Greubel, Dec 01 2017: (Start)
G.f.: x*(1 + 15*x + 15*x^2 + x^3)/(1 - x)^4.
E.g.f.: (-3 + 6*x + 24*x^2 + 16*x^3)*exp(x)/3 + 1. (End)
a(n) = A005900(2n-1). - Ivan N. Ianakiev, Mar 27 2022
From Peter Bala, Mar 11 2024: (Start)
Sum_{k = 1..n+1} 1/(k*a(k)*a(k+1)) = 1/(19 - 3/(27 - 60/(43 - 315/(67 - ... -n^2*(4*n^2 - 1)/((2*n + 1)^2 + 2*3^2))))).
E.g.f.: exp(x)*(1 + 18*x + 48*x^2/2! + 32*x^3/3!). Note that -T(6, i*sqrt(x)) = 1 + 18*x + 48*x^2 + 32*x^3, where T(n, x) denotes the n-th Chebyshev polynomial of the first kind. See A008310. (End)

A047969 Square array of nexus numbers a(n,k) = (n+1)^(k+1) - n^(k+1) (n >= 0, k >= 0) read by upwards antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 5, 7, 1, 1, 7, 19, 15, 1, 1, 9, 37, 65, 31, 1, 1, 11, 61, 175, 211, 63, 1, 1, 13, 91, 369, 781, 665, 127, 1, 1, 15, 127, 671, 2101, 3367, 2059, 255, 1, 1, 17, 169, 1105, 4651, 11529, 14197, 6305, 511, 1, 1, 19, 217, 1695, 9031
Offset: 0

Views

Author

Keywords

Comments

If each row started with an initial 0 (i.e., a(n,k) = (n+1)^k - n^k) then each row would be the binomial transform of the preceding row. - Henry Bottomley, May 31 2001
a(n-1, k-1) is the number of ordered k-tuples of positive integers such that the largest of these integers is n. - Alford Arnold, Sep 07 2005
From Alford Arnold, Jul 21 2006: (Start)
The sequences in A047969 can also be calculated using the Eulerian Array (A008292) and Pascal's Triangle (A007318) as illustrated below: (cf. A101095).
1 1 1 1 1 1
1 1 1 1 1 1
-----------------------------------------
1 2 3 4 5 6
1 2 3 4 5
1 3 5 7 9 11
-----------------------------------------
1 3 6 10 15 21
4 12 24 40 60
1 3 6 10
1 7 19 37 61 91
-----------------------------------------
1 4 10 20 35 56
11 44 110 220 385
11 44 110 220
1 4 10
1 15 65 175 369 671
----------------------------------------- (End)
From Peter Bala, Oct 26 2008: (Start)
The above remarks of Alford Arnold may be summarized by saying that (the transpose of) this array is the Hilbert transform of the triangle of Eulerian numbers A008292 (see A145905 for the definition of the Hilbert transform). In this context, A008292 is best viewed as the array of h-vectors of permutohedra of type A. See A108553 for the Hilbert transform of the array of h-vectors of type D permutohedra. Compare this array with A009998.
The polynomials n^k - (n-1)^k, k = 1,2,3,..., which give the nonzero entries in the columns of this array, satisfy a Riemann hypothesis: their zeros lie on the vertical line Re s = 1/2 in the complex plane. See A019538 for the connection between the polynomials n^k - (n-1)^k and the Stirling polynomials of the simplicial complexes dual to the type A permutohedra.
(End)
Empirical: (n+1)^(k+1) - n^(k+1) is the number of first differences of length k+1 arrays of numbers in 0..n, k > 0. - R. H. Hardin, Jun 30 2013
a(n-1, k-1) is the number of bargraphs of width k and height n. Examples: a(1,2) = 7 because we have [1,1,2], [1,2,1], [2,1,1], [1,2,2], [2,1,2], [2,2,1], and [2,2,2]; a(2,1) = 5 because we have [1,3], [2,3], [3,1], [3,2], and [3,3] (bargraphs are given as compositions). This comment is equivalent to A. Arnold's Sep 2005 comment. - Emeric Deutsch, Jan 30 2017

Examples

			Array a begins:
  [n\k][0  1   2    3    4   5  6  ...
  [0]   1  1   1    1    1   1  1  ...
  [1]   1  3   7   15   31  63  ...
  [2]   1  5  19   65  211  ...
  [3]   1  7  37  175  ...
  ...
Triangle T begins:
  n\m   0   1    2     3     4      5      6      7      8     9  10 ...
  0:    1
  1:    1   1
  2:    1   3    1
  3:    1   5    7     1
  4:    1   7   19    15     1
  5:    1   9   37    65    31      1
  6:    1  11   61   175   211     63      1
  7:    1  13   91   369   781    665    127      1
  8:    1  15  127   671  2101   3367   2059    255      1
  9:    1  17  169  1105  4651  11529  14197   6305    511     1
  10:   1  19  217  1695  9031  31031  61741  58975  19171  1023   1
  ...  - _Wolfdieter Lang_, May 07 2021
		

References

  • J. H. Conway and R. K. Guy, The Book of Numbers, Copernicus Press, NY, 1996, p. 54.

Crossrefs

Cf. A047970.
Cf. A009998, A108553 (Hilbert transform of array of h-vectors of type D permutohedra), A145904, A145905.
Row n sequences of array a: A000012, A000225(k+1), A001047(k+1), A005061(k+1), A005060(k+1), A005062(k+1), A016169(k+1), A016177(k+1), A016185(k+1), A016189(k+1), A016195(k+1), A016197(k+1).
Column k sequences of array a: (nexus numbers): A000012, A005408, A003215, A005917(n+1), A022521, A022522, A022523, A022524, A022525, A022526, A022527, A022528.
Cf. A343237 (row reversed triangle).

Programs

  • Mathematica
    Flatten[Table[n = d - e; k = e; (n + 1)^(k + 1) - n^(k + 1), {d, 0, 100}, {e, 0, d}]] (* T. D. Noe, Feb 22 2012 *)
  • Maxima
    T(n,m):=if m=0 then 1 else sum(k!*(-1)^(m+k)*stirling2(m,k)*binomial(n+k-1,n),k,0,m); /* Vladimir Kruchinin, Jan 28 2018 */

Formula

From Vladimir Kruchinin: (Start)
O.g.f. of e.g.f of rows of array: ((1-x)*exp(y))/(1-x*exp(y))^2.
T(n,m) = Sum_{k=0..m} k!*(-1)^(m+k)*Stirling2(m,k)*C(n+k-1,n), T(n,0)=1.(End)
From Wolfdieter Lang, May 07 2021: (Start)
T(n,m) = a(n-m,m) = (n-m+1)^(m+1) - (n-m)^(m+1), n >= 0, m = 0, 1,..., n.
O.g.f. column k of the array: polylog(-(k+1), x)*(1-x)/x. See the Peter Bala comment above, and the Eulerian triangle A008292 formula by Vladeta Jovovic, Sep 02 2002.
E.g.f. of e.g.f. of row of the array: exp(y)*(1 + x*(exp(y) - 1))*exp(x*exp(y)).
O.g.f. of triangle's exponential row polynomials R(n, y) = Sum_{m=0} T(n, m)*(y^m)/m!: G(x, y) = exp(x*y)*(1 - x)/(1 - x*exp(x*y))^2. (End)

A016755 Odd cubes: a(n) = (2*n + 1)^3.

Original entry on oeis.org

1, 27, 125, 343, 729, 1331, 2197, 3375, 4913, 6859, 9261, 12167, 15625, 19683, 24389, 29791, 35937, 42875, 50653, 59319, 68921, 79507, 91125, 103823, 117649, 132651, 148877, 166375, 185193, 205379, 226981, 250047, 274625, 300763, 328509, 357911, 389017, 421875
Offset: 0

Views

Author

Keywords

Comments

Partial sums of A010014. - Jani Melik, May 20 2013
Terms end in the repeating sequence 1, 7, 5, 3, 9, ... - Melvin Peralta, Jul 08 2015

References

  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2003, Section 1.6.3.

Crossrefs

Programs

Formula

Sum_{n >= 0} 1/a(n) = 7 * zeta(3) / 8.
G.f.: (1+23*x+23*x^2+x^3)/(1-4*x+6*x^2-4*x^3+x^4). - Colin Barker, Jan 02 2012
a(n) = A000578(A005408(n)). - Michel Marcus, Jul 09 2015
E.g.f.: exp(x)*(1 + 26*x + 36*x^2 + 8*x^3). See A154537, row n=3. - Wolfdieter Lang, Mar 12 2017
From Bruce J. Nicholson, Dec 08 2019: (Start)
a(n) = 24 * A000330(n) + A005408(n).
a(n) = 2 * A005917(n+1) - A005408(n). (End)
Sum_{n>=0} (-1)^n/a(n) = Pi^3/32 (A153071). - Amiram Eldar, Oct 10 2020
Product_{n>=1} (1 - (-1)^n/a(n)) = (Pi/12)*(1 + sqrt(2)*cosh(sqrt(3)*Pi/4)) (Chamberland and Straub, 2013). - Amiram Eldar, Jan 26 2024

A022521 a(n) = (n+1)^5 - n^5.

Original entry on oeis.org

1, 31, 211, 781, 2101, 4651, 9031, 15961, 26281, 40951, 61051, 87781, 122461, 166531, 221551, 289201, 371281, 469711, 586531, 723901, 884101, 1069531, 1282711, 1526281, 1803001, 2115751, 2467531
Offset: 0

Views

Author

Keywords

Comments

Last digit of a(n) is always 1. Last two digits of a(n) (i.e., a(n) mod 100) are repeated periodically with palindromic part of period 20 {1,31,11,81,1,51,31,61,81,51,51,81,61,31,51,1,81,11,31,1}. Last three digits of a(n) (i.e., a(n) mod 1000) are repeated periodically with palindromic part of period 200. - Alexander Adamchuk, Aug 11 2006
In Conway and Guy, these numbers are called nexus numbers of order 5. - M. F. Hasler, Jan 27 2013
Numbers that can be arranged in a triangular-antitegmatic icosachoron (the 4D version of "rhombic dodecahedal numbers" (A005917)). - Steven Lu, Mar 28 2023

References

  • John H. Conway and Richard K. Guy, The Book of Numbers, Copernicus Press, NY, 1996, p. 54.

Crossrefs

First differences of A000584.
Column k=4 of array A047969.

Programs

Formula

a(n) = A003215(n) + 24 * A006322(n). - Xavier Acloque, Oct 11 2003
G.f.: (-1-x^4-26*x^3-66*x^2-26*x)/(x-1)^5. - Maksym Voznyy (voznyy(AT)mail.ru), Aug 11 2009
G.f.: polylog(-5, x)*(1-x)/x. See the g.f. of the rows of A008292 by Vladeta Jovovic, Sep 02 2002. - Wolfdieter Lang, May 10 2021
Sum_{n>=0} 1/a(n) = c1*tanh(c2/2) - c2*tanh(c1/2), where c1 = tan(3*Pi/10)*Pi and c2 = tan(Pi/10)*Pi. - Amiram Eldar, Jan 27 2022

A057813 a(n) = (2*n+1)*(4*n^2+4*n+3)/3.

Original entry on oeis.org

1, 11, 45, 119, 249, 451, 741, 1135, 1649, 2299, 3101, 4071, 5225, 6579, 8149, 9951, 12001, 14315, 16909, 19799, 23001, 26531, 30405, 34639, 39249, 44251, 49661, 55495, 61769, 68499, 75701, 83391, 91585, 100299, 109549, 119351, 129721, 140675, 152229, 164399
Offset: 0

Views

Author

N. J. A. Sloane, Nov 07 2000

Keywords

Comments

For n>0, 30*a(n) is the sum of the ten distinct products of 2*n-1, 2*n+1, and 2*n+3. For example, when n = 1, we sum the ten distinct products of 1, 3, and 5: 1*1*1 + 1*1*3 + 1*1*5 + 1*3*3 + 1*3*5 + 1*5*5 + 3*3*3 + 3*3*5 + 3*5*5 + 5*5*5 = 330 = 30*11 = 30*a(1). - J. M. Bergot, Apr 06 2014

Crossrefs

1/12*t*(2*n^3-3*n^2+n)+2*n-1 for t = 2, 4, 6, ... gives A049480, A005894, A063488, A001845, A063489, A005898, A063490, A057813, A063491, A005902, A063492, A005917, A063493, A063494, A063495, A063496.

Programs

  • Magma
    [(2*n+1)*(4*n^2+4*n+3)/3 : n in [0..50]] // Wesley Ivan Hurt, Apr 22 2014
  • Maple
    A057813:=n->(2*n + 1)*(4*n^2 + 4*n + 3)/3; seq(A057813(n), n=0..50); # Wesley Ivan Hurt, Apr 06 2014
  • Mathematica
    Table[(2*n + 1)*(4*n^2 + 4*n + 3)/3, {n, 0, 50}] (* David Nacin, Mar 01 2012 *)
  • PARI
    P(x, y, z) = x^3 + x^2*y + x^2*z + x*y^2 + x*y*z + x*z^2 + y^3 + y^2*z + y*z^2 + z^3;
    a(n) = P(2*n-1, 2*n+1, 2*n+3)/30; \\ Michel Marcus, Apr 22 2014
    

Formula

a(n) = 2*A050533(n) + 1. - N. J. A. Sloane, Sep 22 2004
G.f.: (1+7*x+7*x^2+x^3)/(1-x)^4. - Colin Barker, Mar 01 2012
G.f. for sequence with interpolated zeros: 1/(8*x)*sinh(8*arctanh(x)) = 1/(16*x)*( ((1 + x)/(1 - x))^4 - ((1 - x)/(1 + x))^4 ) = 1 + 11*x^2 + 45*x^4 + 119*x^6 + .... Cf. A019560. - Peter Bala, Apr 07 2017
E.g.f.: (3 + 30*x + 36*x^2 + 8*x^3)*exp(x)/3. - G. C. Greubel, Dec 01 2017
From Peter Bala, Mar 26 2024: (Start)
12*a(n) = (2*n + 1)*(a(n + 1) - a(n - 1)).
Sum_{n >= 0} (-1)^n/(a(n)*a(n+1)) = 3*Pi/16 - 1/2. Cf. A016754 and A336266. (End)

A063490 a(n) = (2*n - 1)*(7*n^2 - 7*n + 6)/6.

Original entry on oeis.org

1, 10, 40, 105, 219, 396, 650, 995, 1445, 2014, 2716, 3565, 4575, 5760, 7134, 8711, 10505, 12530, 14800, 17329, 20131, 23220, 26610, 30315, 34349, 38726, 43460, 48565, 54055, 59944, 66246, 72975, 80145, 87770, 95864, 104441, 113515, 123100, 133210, 143859, 155061
Offset: 1

Views

Author

N. J. A. Sloane, Aug 01 2001

Keywords

Comments

From Omar E. Pol, Oct 23 2019: (Start)
a(n) is also the sum of terms that are in the n-th finite row and in the n-th finite column of the square [1,n]x[1,n] of the natural number array A000027; e.g., the [1,3]x[1,3] square is
1..3..6
2..5..9
4..8..13,
so that a(1) = 1, a(2) = 2 + 3 + 5 = 10, a(3) = 4 + 6 + 8 + 9 + 13 = 40.
Hence the partial sums give A185505. (End)

Crossrefs

1/12*t*(2*n^3-3*n^2+n)+2*n-1 for t = 2, 4, 6, ... gives A049480, A005894, A063488, A001845, A063489, A005898, A063490, A057813, A063491, A005902, A063492, A005917, A063493, A063494, A063495, A063496.

Programs

  • Magma
    [(2*n-1)*(7*n^2-7*n+6)/6: n in [1..30]]; // G. C. Greubel, Dec 01 2017
  • Mathematica
    Table[(2*n-1)*(7*n^2-7*n+6)/6, {n,1,50}] (* or *) LinearRecurrence[{4,-6,4,-1}, {1,10,40,105}, 50] (* G. C. Greubel, Dec 01 2017 *)
  • PARI
    a(n) = { (2*n - 1)*(7*n^2 - 7*n + 6)/6 } \\ Harry J. Smith, Aug 23 2009
    
  • PARI
    my(x='x+O('x^30)); Vec(serlaplace((-6 + 12*x + 21*x^2 + 14*x^3 )*exp(x)/6 + 1)) \\ G. C. Greubel, Dec 01 2017
    

Formula

G.f.: x*(1+x)*(1+5*x+x^2)/(1-x)^4. - Colin Barker, Mar 02 2012
a(n) = Sum_{k = n^2-2*n+2..n^2} A064788(k). - Lior Manor, Jan 13 2013
From G. C. Greubel, Dec 01 2017: (Start)
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4).
E.g.f.: (-6 + 12*x + 21*x^2 + 14*x^3)*exp(x)/6 + 1. (End)

A063488 a(n) = (2*n-1)*(n^2 -n +2)/2.

Original entry on oeis.org

1, 6, 20, 49, 99, 176, 286, 435, 629, 874, 1176, 1541, 1975, 2484, 3074, 3751, 4521, 5390, 6364, 7449, 8651, 9976, 11430, 13019, 14749, 16626, 18656, 20845, 23199, 25724, 28426, 31311, 34385, 37654, 41124, 44801, 48691, 52800, 57134, 61699, 66501, 71546, 76840
Offset: 1

Views

Author

N. J. A. Sloane, Aug 01 2001

Keywords

Comments

Sum of two consecutive terms of A006003(n) = n*(n^2+1)/2. a(n) = A006003(n-1) + A006003(n). - Alexander Adamchuk, Jun 03 2006
If a 2-set Y and a 3-set Z are disjoint subsets of an n-set X then a(n-4) is the number of 5-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 08 2007

Crossrefs

1/12*t*n*(2*n^2 - 3*n + 1) + 2*n - 1 for t = 2, 4, 6, ... gives A049480, A005894, A063488, A001845, A063489, A005898, A063490, A057813, A063491, A005902, A063492, A005917, A063493, A063494, A063495, A063496.
Partial sums of A005918.

Programs

  • Magma
    [(2*n-1)*(n^2 -n +2)/2: n in [1..30]]; // G. C. Greubel, Dec 01 2017
  • Mathematica
    Table[(2 n - 1) (n^2 - n + 2)/2, {n, 1, 40}] (* Bruno Berselli, Oct 14 2016 *)
    LinearRecurrence[{4,-6,4,-1}, {1,6,20,49}, 50] (* G. C. Greubel, Dec 01 2017 *)
  • PARI
    a(n) = { (2*n - 1)*(n^2 - n + 2)/2 } \\ Harry J. Smith, Aug 23 2009
    
  • PARI
    my(x='x+O('x^30)); Vec(serlaplace((-2 + 4*x + 3*x^2 + 2*x^3)*exp(x)/2 + 1)) \\ G. C. Greubel, Dec 01 2017
    

Formula

G.f.: (1 + x)*(1 + x + x^2)/(1 - x)^4. - Jaume Oliver Lafont, Aug 30 2009
a(n) = A000217(A000217(n)) - A000217(A000217(n-2)). - Bruno Berselli, Oct 14 2016
E.g.f.: (-2 + 4*x + 3*x^2 + 2*x^3)*exp(x)/2 + 1. - G. C. Greubel, Dec 01 2017

A063494 a(n) = (2*n - 1)*(7*n^2 - 7*n + 3)/3.

Original entry on oeis.org

1, 17, 75, 203, 429, 781, 1287, 1975, 2873, 4009, 5411, 7107, 9125, 11493, 14239, 17391, 20977, 25025, 29563, 34619, 40221, 46397, 53175, 60583, 68649, 77401, 86867, 97075, 108053, 119829, 132431, 145887, 160225, 175473, 191659, 208811, 226957, 246125, 266343, 287639
Offset: 1

Views

Author

N. J. A. Sloane, Aug 01 2001

Keywords

Comments

Interpret A176271 as an infinite square array read by antidiagonals, with rows 1,5,11,19,...; 3,9,17,27,... and so on. The sum of the terms in the n X n upper submatrix are s(n) = 1, 18, 93, 296, ... = n^2*(7*n^2-1)/6, and a(n) = s(n) - s(n-1) are the first differences. - J. M. Bergot, Jun 27 2013

Crossrefs

1/12*t*(2*n^3-3*n^2+n)+2*n-1 for t = 2, 4, 6, ... gives A049480, A005894, A063488, A001845, A063489, A005898, A063490, A057813, A063491, A005902, A063492, A005917, A063493, A063494, A063495, A063496.

Programs

  • Magma
    [(2*n - 1)*(7*n^2 - 7*n + 3)/3: n in [1..30]]; // G. C. Greubel, Dec 01 2017
  • Mathematica
    Table[(2*n - 1)*(7*n^2 - 7*n + 3)/3, {n,1,30}] (* or *) LinearRecurrence[{4,-6,4,-1}, {1,17,75,203}, 30] (* G. C. Greubel, Dec 01 2017 *)
  • PARI
    a(n) = { (2*n - 1)*(7*n^2 - 7*n + 3)/3 } \\ Harry J. Smith, Aug 23 2009
    
  • PARI
    my(x='x+O('x^30)); Vec(serlaplace((-3+6*x+21*x^2+14*x^3)*exp(x)/3 + 1)) \\ G. C. Greubel, Dec 01 2017
    

Formula

G.f.: x*(1+x)*(1+12*x+x^2)/(1-x)^4. - Colin Barker, Mar 02 2012
E.g.f.: (-3 + 6*x + 21*x^2 + 14*x^3)*exp(x)/3 + 1. - G. C. Greubel, Dec 01 2017
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - Wesley Ivan Hurt, May 11 2023

A063491 a(n) = (2*n - 1)*(3*n^2 - 3*n + 2)/2.

Original entry on oeis.org

1, 12, 50, 133, 279, 506, 832, 1275, 1853, 2584, 3486, 4577, 5875, 7398, 9164, 11191, 13497, 16100, 19018, 22269, 25871, 29842, 34200, 38963, 44149, 49776, 55862, 62425, 69483, 77054, 85156, 93807, 103025, 112828, 123234, 134261, 145927, 158250, 171248, 184939
Offset: 1

Views

Author

N. J. A. Sloane, Aug 01 2001

Keywords

Comments

A triangle has sides of lengths 6*n-3, 6*n^2-6*n+4, and 6*n^2-6*n+7; for n>2 its area is 6*sqrt(a(n)^2 - 1). - J. M. Bergot, Aug 30 2013
[The source of this is using (n,n+1), (n+1,n+2), and (n+2,n+3) as (a,b) in the creation of three Pythagorean triangles with sides b^2-a^2, 2*a*b, and a^2+b^2. Combine the three respective sides to create a new larger triangle, then find its area. It is not simply working backwards from the sequence. As well, the sequence has this as its first comment to show that the numbers are actually doing something to find a solution.]

References

  • T. A. Gulliver, Sequences from Arrays of Integers, Int. Math. Journal, Vol. 1, No. 4, pp. 323-332, 2002.

Crossrefs

1/12*t*(2*n^3-3*n^2+n)+2*n-1 for t = 2, 4, 6, ... gives A049480, A005894, A063488, A001845, A063489, A005898, A063490, A057813, A063491, A005902, A063492, A005917, A063493, A063494, A063495, A063496.

Programs

  • Magma
    [(2*n-1)*(3*n^2 -3*n +2)/2: n in [1..30]]; // G. C. Greubel, Dec 01 2017
  • Mathematica
    LinearRecurrence[{4,-6,4,-1},{1,12,50,133},40] (* Harvey P. Dale, Jun 05 2016 *)
    Table[(2*n-1)*(3*n^2 -3*n +2)/2, {n,1,30}] (* G. C. Greubel, Dec 01 2017 *)
  • PARI
    a(n) = { (2*n - 1)*(3*n^2 - 3*n + 2)/2 } \\ Harry J. Smith, Aug 23 2009
    
  • PARI
    my(x='x+O('x^30)); Vec(serlaplace((-2 + 4*x + 9*x^2 + 6*x^3)*exp(x)/2 + 1)) \\ G. C. Greubel, Dec 01 2017
    
  • R
    a <- c(0, 1, 9, 38, 110)
    for(n in (length(a)+1):40)
      a[n] <- +4*a[n-1]-6*a[n-2]+4*a[n-3]-a[n-4]
    a [Yosu Yurramendi, Sep 04 2013]
    

Formula

G.f.: x*(1+x)*(1+7*x+x^2)/(1-x)^4. - Colin Barker, Apr 20 2012
a(n) = +4*a(n-1) -6*a(n-2) +4*a(n-3) -1*a(n-4) n > 3, a(1)=1, a(2)=12, a(3)=50, a(4)=133. - Yosu Yurramendi, Sep 04 2013
E.g.f.: (-2 + 4*x + 9*x^2 + 6*x^3)*exp(x)/2 + 1. - G. C. Greubel, Dec 01 2017
From Bruce J. Nicholson, Jun 17 2020: (Start)
a(n) = A005448(n) * A005408(n-1).
a(n) = A004188(n) + A004188(n-1). (End)

A063492 a(n) = (2*n - 1)*(11*n^2 - 11*n + 6)/6.

Original entry on oeis.org

1, 14, 60, 161, 339, 616, 1014, 1555, 2261, 3154, 4256, 5589, 7175, 9036, 11194, 13671, 16489, 19670, 23236, 27209, 31611, 36464, 41790, 47611, 53949, 60826, 68264, 76285, 84911, 94164, 104066, 114639, 125905, 137886, 150604, 164081, 178339, 193400, 209286, 226019
Offset: 1

Views

Author

N. J. A. Sloane, Aug 01 2001

Keywords

Crossrefs

1/12*t*(2*n^3 - 3*n^2 + n) + 2*n - 1 for t = 2, 4, 6, ... gives A049480, A005894, A063488, A001845, A063489, A005898, A063490, A057813, A063491, A005902, A063492, A005917, A063493, A063494, A063495, A063496.

Programs

Formula

G.f.: x*(1+x)*(1 + 9*x + x^2)/(1-x)^4. - Colin Barker, Apr 24 2012
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) for n>4. - Wesley Ivan Hurt, Dec 16 2015
E.g.f.: (-6 + 12*x + 33*x^2 + 22*x^3)*exp(x)/6 + 1. - G. C. Greubel, Dec 01 2017
Previous Showing 11-20 of 60 results. Next