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 69 results. Next

A186339 a(n)=A006125(n+1)*2^A001840(n).

Original entry on oeis.org

1, 4, 32, 512, 32768, 4194304, 1073741824, 1099511627776, 2251799813685248, 9223372036854775808, 151115727451828646838272, 4951760157141521099596496896, 324518553658426726783156020576256, 85070591730234615865843651857942052864
Offset: 0

Views

Author

Paul Barry, Feb 18 2011

Keywords

Comments

Hankel transform of A186338.

Formula

a(n)=2^binomial(n+1,2)*2^floor((n+1)(n+2)/6).

A091361 Numbers n such that A001840(n) == 0 (mod n).

Original entry on oeis.org

1, 2, 3, 9, 15, 21, 27, 33, 39, 45, 51, 57, 63, 69, 75, 81, 87, 93, 99, 105, 111, 117, 123, 129, 135, 141, 147, 153, 159, 165, 171, 177, 183, 189, 195, 201, 207, 213, 219, 225, 231, 237, 243, 249, 255, 261, 267, 273, 279, 285, 291, 297, 303, 309, 315, 321, 327
Offset: 1

Views

Author

Jon Perry, Mar 01 2004

Keywords

Comments

Apart from 1 and 2 it is conjectured that the only values present are congruent to 3 mod 6 (all these values are present).
From R. J. Mathar, Feb 25 2008: (Start)
Proof of the conjecture that this is 1 and 2 followed by A016945 follows by considering the 6 cases n=6k-1, 6k, 6k+1, 6k+2, 6k+3 or 6k+4, individual evaluation of A001840(n) with their corresponding 3 formulas quoted in A001840 in each case and searching for solutions of the form A001840(n) = t*n for integer t.
Example: A001840[6k+4]=A001840[3(2k+1)+1]=(2k+2)(6k+5)/2=t*(6k+4) implies t=k+7/6+1/[6(3k+2)] which cannot be solved in integers t and k. So numbers of the form 6k+4 are not members here. (End)

Examples

			A001840(9)=18, so 9 is in the sequence.
		

Crossrefs

Cf. A001840.

Programs

  • Mathematica
    (* b = A001840 *) b[0] = 0; b[1] = 1; b[n_] := b[n] = n (n + 1)/2 - b[n - 1] - b[n - 2]; Reap[For[n = 1, n <= 400, n++, If[Mod[b[n], n] == 0, Sow[n]]]][[2, 1]] (* Jean-François Alcover, Feb 09 2019 *)

Formula

G.f.: conjecture: 2*(1+x)/(1-x)/G(0) +x, where G(k)= 1 + 1/(1 - x*(3*k+1)/(x*(3*k+4) + 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 06 2013

A165189 Partial sums of partial sums of (A001840 interleaved with zeros).

Original entry on oeis.org

1, 2, 5, 8, 14, 20, 31, 42, 60, 78, 105, 132, 171, 210, 264, 318, 390, 462, 556, 650, 770, 890, 1040, 1190, 1375, 1560, 1785, 2010, 2280, 2550, 2871, 3192, 3570, 3948, 4389, 4830, 5341, 5852, 6440, 7028, 7700, 8372, 9136, 9900, 10764, 11628, 12600, 13572
Offset: 1

Views

Author

Alford Arnold, Sep 16 2009

Keywords

Comments

Also convolution of period six sequence 1,0,0,0,0,0,1,... (A079979) with sequence 1,2,5,8,14,20,30,40,... (A006918 without initial zero).

Examples

			A001840 interleaved with zeros is
1, 0, 2, 0, 3, 0, 5, 0, 7, 0, 9, 0, 12, 0, 15, 0, ...
Partial sums thereof are
1, 1, 3, 3, 6, 6, 11, 11, 18, 18, 27, 27, 39, 39, 54, 54, ...
This equals A014125 interleaved with itself.
Partial sums thereof are
1, 2, 5, 8, 14, 20, 31, 42, 60, 78, 105, 132, 171, 210, 264, 318, ...
		

Crossrefs

Cf. A001840 (expansion of x/((1-x)^3*(1+x+x^2))), A001840 (expansion of x/((1-x)^2*(1-x^3))), A079979, A006918, A014125.

Programs

  • Mathematica
    Drop[Accumulate[Accumulate[Riffle[LinearRecurrence[{2,-1,1,-2,1},{0,1,2,3,5},30],0]]],2] (* or *) LinearRecurrence[{2,1,-4,1,2,0,-2,-1,4,-1,-2,1},{1,2,5,8,14,20,31,42,60,78,105,132},50] (* Harvey P. Dale, Jun 08 2018 *)
  • PARI
    /* first computes u = A001840 interleaved with zeros, then v = partial sums, then w = second partial sums */ {m=50; u=vector(m, n, polcoeff(x/((1-x^2)^3*(1+x^2+x^4))+x*O(x^(n)),n)); v=vector(m); a=u[1]; v[1]=a; for(n=2, m, a+=u[n]; v[n]=a); w=vector(m-1); a=v[1]; w[1]=a; for(n=2, m-1, a+=v[n]; w[n]=a); w} \\ Klaus Brockhaus, Sep 21 2009

Formula

G.f.: x/((1-x)^5*(1+x)^3*(1-x+x^2)*(1+x+x^2)).
54*a(n) = 631/64 +405/16*n +3/32*n^4 +15/8*n^3 +381/32*n^2 -(-1)^n*( 9/32*n^2 +45/16*n +375/64) -A131713(n) -3*A057079(n). - R. J. Mathar, Jun 16 2018

Extensions

Edited and corrected by R. J. Mathar, Klaus Brockhaus and N. J. A. Sloane, Sep 21 2009 - Sep 25 2009

A186340 a(n) = 2^A001840(n).

Original entry on oeis.org

1, 2, 4, 8, 32, 128, 512, 4096, 32768, 262144, 4194304, 67108864, 1073741824, 34359738368, 1099511627776, 35184372088832, 2251799813685248, 144115188075855872, 9223372036854775808, 1180591620717411303424, 151115727451828646838272
Offset: 0

Views

Author

Paul Barry, Feb 18 2011

Keywords

Comments

A factor in the Hankel transform A186339 of A186338.
a(n)*a(n-4) = 2*a(n-1)*a(n-3) = a(n-1)*a(n-3) + c(n)*a(n-2)^2, where c(3*n+2) = 2, c(3*n) = c(3*n+1) = 1 for all n in Z. - Michael Somos, Oct 19 2018

Examples

			G.f. = 1 + 2*x + 4*x^2 + 8*x^3 + 32*x^4 + 128*x^5 + 512*x^6 + ... - _Michael Somos_, Oct 19 2018
		

Crossrefs

Cf. A058937.

Programs

  • Mathematica
    a[ n_] := 2^Quotient[ Binomial[n + 2, 2], 3]; (* Michael Somos, Oct 19 2018 *)
  • PARI
    {a(n) = 2^(binomial(n+2, 2)\3)}; /* Michael Somos, Oct 19 2018 */

Formula

a(n)=2^floor((n+1)(n+2)/6).

A000217 Triangular numbers: a(n) = binomial(n+1,2) = n*(n+1)/2 = 0 + 1 + 2 + ... + n.

Original entry on oeis.org

0, 1, 3, 6, 10, 15, 21, 28, 36, 45, 55, 66, 78, 91, 105, 120, 136, 153, 171, 190, 210, 231, 253, 276, 300, 325, 351, 378, 406, 435, 465, 496, 528, 561, 595, 630, 666, 703, 741, 780, 820, 861, 903, 946, 990, 1035, 1081, 1128, 1176, 1225, 1275, 1326, 1378, 1431
Offset: 0

Views

Author

Keywords

Comments

Also referred to as T(n) or C(n+1, 2) or binomial(n+1, 2) (preferred).
Also generalized hexagonal numbers: n*(2*n-1), n=0, +-1, +-2, +-3, ... Generalized k-gonal numbers are second k-gonal numbers and positive terms of k-gonal numbers interleaved, k >= 5. In this case k = 6. - Omar E. Pol, Sep 13 2011 and Aug 04 2012
Number of edges in complete graph of order n+1, K_{n+1}.
Number of legal ways to insert a pair of parentheses in a string of n letters. E.g., there are 6 ways for three letters: (a)bc, (ab)c, (abc), a(b)c, a(bc), ab(c). Proof: there are C(n+2,2) ways to choose where the parentheses might go, but n + 1 of them are illegal because the parentheses are adjacent. Cf. A002415.
For n >= 1, a(n) is also the genus of a nonsingular curve of degree n+2, such as the Fermat curve x^(n+2) + y^(n+2) = 1. - Ahmed Fares (ahmedfares(AT)my_deja.com), Feb 21 2001
From Harnack's theorem (1876), the number of branches of a nonsingular curve of order n is bounded by a(n-1)+1, and the bound can be achieved. See also A152947. - Benoit Cloitre, Aug 29 2002. Corrected by Robert McLachlan, Aug 19 2024
Number of tiles in the set of double-n dominoes. - Scott A. Brown, Sep 24 2002
Number of ways a chain of n non-identical links can be broken up. This is based on a similar problem in the field of proteomics: the number of ways a peptide of n amino acid residues can be broken up in a mass spectrometer. In general, each amino acid has a different mass, so AB and BC would have different masses. - James A. Raymond, Apr 08 2003
Triangular numbers - odd numbers = shifted triangular numbers; 1, 3, 6, 10, 15, 21, ... - 1, 3, 5, 7, 9, 11, ... = 0, 0, 1, 3, 6, 10, ... - Xavier Acloque, Oct 31 2003 [Corrected by Derek Orr, May 05 2015]
Centered polygonal numbers are the result of [number of sides * A000217 + 1]. E.g., centered pentagonal numbers (1,6,16,31,...) = 5 * (0,1,3,6,...) + 1. Centered heptagonal numbers (1,8,22,43,...) = 7 * (0,1,3,6,...) + 1. - Xavier Acloque, Oct 31 2003
Maximum number of lines formed by the intersection of n+1 planes. - Ron R. King, Mar 29 2004
Number of permutations of [n] which avoid the pattern 132 and have exactly 1 descent. - Mike Zabrocki, Aug 26 2004
Number of ternary words of length n-1 with subwords (0,1), (0,2) and (1,2) not allowed. - Olivier Gérard, Aug 28 2012
Number of ways two different numbers can be selected from the set {0,1,2,...,n} without repetition, or, number of ways two different numbers can be selected from the set {1,2,...,n} with repetition.
Conjecturally, 1, 6, 120 are the only numbers that are both triangular and factorial. - Christopher M. Tomaszewski (cmt1288(AT)comcast.net), Mar 30 2005
Binomial transform is {0, 1, 5, 18, 56, 160, 432, ...}, A001793 with one leading zero. - Philippe Deléham, Aug 02 2005
Each pair of neighboring terms adds to a perfect square. - Zak Seidov, Mar 21 2006
Number of transpositions in the symmetric group of n+1 letters, i.e., the number of permutations that leave all but two elements fixed. - Geoffrey Critzer, Jun 23 2006
With rho(n):=exp(i*2*Pi/n) (an n-th root of 1) one has, for n >= 1, rho(n)^a(n) = (-1)^(n+1). Just use the triviality a(2*k+1) == 0 (mod (2*k+1)) and a(2*k) == k (mod (2*k)).
a(n) is the number of terms in the expansion of (a_1 + a_2 + a_3)^(n-1). - Sergio Falcon, Feb 12 2007
a(n+1) is the number of terms in the complete homogeneous symmetric polynomial of degree n in 2 variables. - Richard Barnes, Sep 06 2017
The number of distinct handshakes in a room with n+1 people. - Mohammad K. Azarian, Apr 12 2007 [corrected, Joerg Arndt, Jan 18 2016]
Equal to the rank (minimal cardinality of a generating set) of the semigroup PT_n\S_n, where PT_n and S_n denote the partial transformation semigroup and symmetric group on [n]. - James East, May 03 2007
a(n) gives the total number of triangles found when cevians are drawn from a single vertex on a triangle to the side opposite that vertex, where n = the number of cevians drawn+1. For instance, with 1 cevian drawn, n = 1+1 = 2 and a(n)= 2*(2+1)/2 = 3 so there is a total of 3 triangles in the figure. If 2 cevians are drawn from one point to the opposite side, then n = 1+2 = 3 and a(n) = 3*(3+1)/2 = 6 so there is a total of 6 triangles in the figure. - Noah Priluck (npriluck(AT)gmail.com), Apr 30 2007
For n >= 1, a(n) is the number of ways in which n-1 can be written as a sum of three nonnegative integers if representations differing in the order of the terms are considered to be different. In other words, for n >= 1, a(n) is the number of nonnegative integral solutions of the equation x + y + z = n-1. - Amarnath Murthy, Apr 22 2001 (edited by Robert A. Beeler)
a(n) is the number of levels with energy n + 3/2 (in units of h*f0, with Planck's constant h and the oscillator frequency f0) of the three-dimensional isotropic harmonic quantum oscillator. See the comment by A. Murthy above: n = n1 + n2 + n3 with positive integers and ordered. Proof from the o.g.f. See the A. Messiah reference. - Wolfdieter Lang, Jun 29 2007
From Hieronymus Fischer, Aug 06 2007: (Start)
Numbers m >= 0 such that round(sqrt(2m+1)) - round(sqrt(2m)) = 1.
Numbers m >= 0 such that ceiling(2*sqrt(2m+1)) - 1 = 1 + floor(2*sqrt(2m)).
Numbers m >= 0 such that fract(sqrt(2m+1)) > 1/2 and fract(sqrt(2m)) < 1/2, where fract(x) is the fractional part of x (i.e., x - floor(x), x >= 0). (End)
If Y and Z are 3-blocks of an n-set X, then, for n >= 6, a(n-1) is the number of (n-2)-subsets of X intersecting both Y and Z. - Milan Janjic, Nov 09 2007
Equals row sums of triangle A143320, n > 0. - Gary W. Adamson, Aug 07 2008
a(n) is also an even perfect number in A000396 iff n is a Mersenne prime A000668. - Omar E. Pol, Sep 05 2008. Unnecessary assumption removed and clarified by Rick L. Shepherd, Apr 14 2025
Equals row sums of triangle A152204. - Gary W. Adamson, Nov 29 2008
The number of matches played in a round robin tournament: n*(n-1)/2 gives the number of matches needed for n players. Everyone plays against everyone else exactly once. - Georg Wrede (georg(AT)iki.fi), Dec 18 2008
-a(n+1) = E(2)*binomial(n+2,2) (n >= 0) where E(n) are the Euler numbers in the enumeration A122045. Viewed this way, a(n) is the special case k=2 in the sequence of diagonals in the triangle A153641. - Peter Luschny, Jan 06 2009
Equivalent to the first differences of successive tetrahedral numbers. See A000292. - Jeremy Cahill (jcahill(AT)inbox.com), Apr 15 2009
The general formula for alternating sums of powers is in terms of the Swiss-Knife polynomials P(n,x) A153641 2^(-n-1)(P(n,1)-(-1)^k P(n,2k+1)). Thus a(k) = |2^(-3)(P(2,1)-(-1)^k P(2,2k+1))|. - Peter Luschny, Jul 12 2009
a(n) is the smallest number > a(n-1) such that gcd(n,a(n)) = gcd(n,a(n-1)). If n is odd this gcd is n; if n is even it is n/2. - Franklin T. Adams-Watters, Aug 06 2009
Partial sums of A001477. - Juri-Stepan Gerasimov, Jan 25 2010. [A-number corrected by Omar E. Pol, Jun 05 2012]
The numbers along the right edge of Floyd's triangle are 1, 3, 6, 10, 15, .... - Paul Muljadi, Jan 25 2010
From Charlie Marion, Dec 03 2010: (Start)
More generally, a(2k+1) == j*(2j-1) (mod 2k+2j+1) and
a(2k) == [-k + 2j*(j-1)] (mod 2k+2j).
Column sums of:
1 3 5 7 9 ...
1 3 5 ...
1 ...
...............
---------------
1 3 6 10 15 ...
Sum_{n>=1} 1/a(n)^2 = 4*Pi^2/3-12 = 12 less than the volume of a sphere with radius Pi^(1/3).
(End)
A004201(a(n)) = A000290(n); A004202(a(n)) = A002378(n). - Reinhard Zumkeller, Feb 12 2011
1/a(n+1), n >= 0, has e.g.f. -2*(1+x-exp(x))/x^2, and o.g.f. 2*(x+(1-x)*log(1-x))/x^2 (see the Stephen Crowley formula line). -1/(2*a(n+1)) is the z-sequence for the Sheffer triangle of the coefficients of the Bernoulli polynomials A196838/A196839. - Wolfdieter Lang, Oct 26 2011
From Charlie Marion, Feb 23 2012: (Start)
a(n) + a(A002315(k)*n + A001108(k+1)) = (A001653(k+1)*n + A001109(k+1))^2. For k=0 we obtain a(n) + a(n+1) = (n+1)^2 (identity added by N. J. A. Sloane on Feb 19 2004).
a(n) + a(A002315(k)*n - A055997(k+1)) = (A001653(k+1)*n - A001109(k))^2.
(End)
Plot the three points (0,0), (a(n), a(n+1)), (a(n+1), a(n+2)) to form a triangle. The area will be a(n+1)/2. - J. M. Bergot, May 04 2012
The sum of four consecutive triangular numbers, beginning with a(n)=n*(n+1)/2, minus 2 is 2*(n+2)^2. a(n)*a(n+2)/2 = a(a(n+1)-1). - J. M. Bergot, May 17 2012
(a(n)*a(n+3) - a(n+1)*a(n+2))*(a(n+1)*a(n+4) - a(n+2)*a(n+3))/8 = a((n^2+5*n+4)/2). - J. M. Bergot, May 18 2012
a(n)*a(n+1) + a(n+2)*a(n+3) + 3 = a(n^2 + 4*n + 6). - J. M. Bergot, May 22 2012
In general, a(n)*a(n+1) + a(n+k)*a(n+k+1) + a(k-1)*a(k) = a(n^2 + (k+2)*n + k*(k+1)). - Charlie Marion, Sep 11 2012
a(n)*a(n+3) + a(n+1)*a(n+2) = a(n^2 + 4*n + 2). - J. M. Bergot, May 22 2012
In general, a(n)*a(n+k) + a(n+1)*a(n+k-1) = a(n^2 + (k+1)*n + k-1). - Charlie Marion, Sep 11 2012
a(n)*a(n+2) + a(n+1)*a(n+3) = a(n^2 + 4*n + 3). - J. M. Bergot, May 22 2012
Three points (a(n),a(n+1)), (a(n+1),a(n)) and (a(n+2),a(n+3)) form a triangle with area 4*a(n+1). - J. M. Bergot, May 23 2012
a(n) + a(n+k) = (n+k)^2 - (k^2 + (2n-1)*k -2n)/2. For k=1 we obtain a(n) + a(n+1) = (n+1)^2 (see below). - Charlie Marion, Oct 02 2012
In n-space we can define a(n-1) nontrivial orthogonal projections. For example, in 3-space there are a(2)=3 (namely point onto line, point onto plane, line onto plane). - Douglas Latimer, Dec 17 2012
From James East, Jan 08 2013: (Start)
For n >= 1, a(n) is equal to the rank (minimal cardinality of a generating set) and idempotent rank (minimal cardinality of an idempotent generating set) of the semigroup P_n\S_n, where P_n and S_n denote the partition monoid and symmetric group on [n].
For n >= 3, a(n-1) is equal to the rank and idempotent rank of the semigroup T_n\S_n, where T_n and S_n denote the full transformation semigroup and symmetric group on [n].
(End)
For n >= 3, a(n) is equal to the rank and idempotent rank of the semigroup PT_n\S_n, where PT_n and S_n denote the partial transformation semigroup and symmetric group on [n]. - James East, Jan 15 2013
Conjecture: For n > 0, there is always a prime between A000217(n) and A000217(n+1). Sequence A065383 has the first 1000 of these primes. - Ivan N. Ianakiev, Mar 11 2013
The formula, a(n)*a(n+4k+2)/2 + a(k) = a(a(n+2k+1) - (k^2+(k+1)^2)), is a generalization of the formula a(n)*a(n+2)/2 = a(a(n+1)-1) in Bergot's comment dated May 17 2012. - Charlie Marion, Mar 28 2013
The series Sum_{k>=1} 1/a(k) = 2, given in a formula below by Jon Perry, Jul 13 2003, has partial sums 2*n/(n+1) (telescopic sum) = A022998(n)/A026741(n+1). - Wolfdieter Lang, Apr 09 2013
For odd m = 2k+1, we have the recurrence a(m*n + k) = m^2*a(n) + a(k). Corollary: If number T is in the sequence then so is 9*T+1. - Lekraj Beedassy, May 29 2013
Euler, in Section 87 of the Opera Postuma, shows that whenever T is a triangular number then 9*T + 1, 25*T + 3, 49*T + 6 and 81*T + 10 are also triangular numbers. In general, if T is a triangular number then (2*k + 1)^2*T + k*(k + 1)/2 is also a triangular number. - Peter Bala, Jan 05 2015
Using 1/b and 1/(b+2) will give a Pythagorean triangle with sides 2*b + 2, b^2 + 2*b, and b^2 + 2*b + 2. Set b=n-1 to give a triangle with sides of lengths 2*n,n^2-1, and n^2 + 1. One-fourth the perimeter = a(n) for n > 1. - J. M. Bergot, Jul 24 2013
a(n) = A028896(n)/6, where A028896(n) = s(n) - s(n-1) are the first differences of s(n) = n^3 + 3*n^2 + 2*n - 8. s(n) can be interpreted as the sum of the 12 edge lengths plus the sum of the 6 face areas plus the volume of an n X (n-1) X (n-2) rectangular prism. - J. M. Bergot, Aug 13 2013
Dimension of orthogonal group O(n+1). - Eric M. Schmidt, Sep 08 2013
Number of positive roots in the root system of type A_n (for n > 0). - Tom Edgar, Nov 05 2013
A formula for the r-th successive summation of k, for k = 1 to n, is binomial(n+r,r+1) [H. W. Gould]. - Gary Detlefs, Jan 02 2014
Also the alternating row sums of A095831. Also the alternating row sums of A055461, for n >= 1. - Omar E. Pol, Jan 26 2014
For n >= 3, a(n-2) is the number of permutations of 1,2,...,n with the distribution of up (1) - down (0) elements 0...011 (n-3 zeros), or, the same, a(n-2) is up-down coefficient {n,3} (see comment in A060351). - Vladimir Shevelev, Feb 14 2014
a(n) is the dimension of the vector space of symmetric n X n matrices. - Derek Orr, Mar 29 2014
Non-vanishing subdiagonal of A132440^2/2, aside from the initial zero. First subdiagonal of unsigned A238363. Cf. A130534 for relations to colored forests, disposition of flags on flagpoles, and colorings of the vertices of complete graphs. - Tom Copeland, Apr 05 2014
The number of Sidon subsets of {1,...,n+1} of size 2. - Carl Najafi, Apr 27 2014
Number of factors in the definition of the Vandermonde determinant V(x_1,x_2,...,x_n) = Product_{1 <= i < k <= n} x_i - x_k. - Tom Copeland, Apr 27 2014
Number of weak compositions of n into three parts. - Robert A. Beeler, May 20 2014
Suppose a bag contains a(n) red marbles and a(n+1) blue marbles, where a(n), a(n+1) are consecutive triangular numbers. Then, for n > 0, the probability of choosing two marbles at random and getting two red or two blue is 1/2. In general, for k > 2, let b(0) = 0, b(1) = 1 and, for n > 1, b(n) = (k-1)*b(n-1) - b(n-2) + 1. Suppose, for n > 0, a bag contains b(n) red marbles and b(n+1) blue marbles. Then the probability of choosing two marbles at random and getting two red or two blue is (k-1)/(k+1). See also A027941, A061278, A089817, A053142, A092521. - Charlie Marion, Nov 03 2014
Let O(n) be the oblong number n(n+1) = A002378 and S(n) the square number n^2 = A000290(n). Then a(4n) = O(3n) - O(n), a(4n+1) = S(3n+1) - S(n), a(4n+2) = S(3n+2) - S(n+1) and a(4n+3) = O(3n+2) - O(n). - Charlie Marion, Feb 21 2015
Consider the partition of the natural numbers into parts from the set S=(1,2,3,...,n). The length (order) of the signature of the resulting sequence is given by the triangular numbers. E.g., for n=10, the signature length is 55. - David Neil McGrath, May 05 2015
a(n) counts the partitions of (n-1) unlabeled objects into three (3) parts (labeled a,b,c), e.g., a(5)=15 for (n-1)=4. These are (aaaa),(bbbb),(cccc),(aaab),(aaac),(aabb),(aacc),(aabc),(abbc),(abcc),(abbb),(accc),(bbcc),(bccc),(bbbc). - David Neil McGrath, May 21 2015
Conjecture: the sequence is the genus/deficiency of the sinusoidal spirals of index n which are algebraic curves. The value 0 corresponds to the case of the Bernoulli Lemniscate n=2. So the formula conjectured is (n-1)(n-2)/2. - Wolfgang Tintemann, Aug 02 2015
Conjecture: Let m be any positive integer. Then, for each n = 1,2,3,... the set {Sum_{k=s..t} 1/k^m: 1 <= s <= t <= n} has cardinality a(n) = n*(n+1)/2; in other words, all the sums Sum_{k=s..t} 1/k^m with 1 <= s <= t are pairwise distinct. (I have checked this conjecture via a computer and found no counterexample.) - Zhi-Wei Sun, Sep 09 2015
The Pisano period lengths of reading the sequence modulo m seem to be A022998(m). - R. J. Mathar, Nov 29 2015
For n >= 1, a(n) is the number of compositions of n+4 into n parts avoiding the part 2. - Milan Janjic, Jan 07 2016
In this sequence only 3 is prime. - Fabian Kopp, Jan 09 2016
Suppose you are playing Bulgarian Solitaire (see A242424 and Chamberland's and Gardner's books) and, for n > 0, you are starting with a single pile of a(n) cards. Then the number of operations needed to reach the fixed state {n, n-1,...,1} is a(n-1). For example, {6}->{5,1}->{4,2}->{3,2,1}. - Charlie Marion, Jan 14 2016
Numbers k such that 8k + 1 is a square. - Juri-Stepan Gerasimov, Apr 09 2016
Every perfect cube is the difference of the squares of two consecutive triangular numbers. 1^2-0^2 = 1^3, 3^2-1^2 = 2^3, 6^2-3^2 = 3^3. - Miquel Cerda, Jun 26 2016
For n > 1, a(n) = tau_n(k*) where tau_n(k) is the number of ordered n-factorizations of k and k* is the square of a prime. For example, tau_3(4) = tau_3(9) = tau_3(25) = tau_3(49) = 6 (see A007425) since the number of divisors of 4, 9, 25, and 49's divisors is 6, and a(3) = 6. - Melvin Peralta, Aug 29 2016
In an (n+1)-dimensional hypercube, number of two-dimensional faces congruent with a vertex (see also A001788). - Stanislav Sykora, Oct 23 2016
Generalizations of the familiar formulas, a(n) + a(n+1) = (n+1)^2 (Feb 19 2004) and a(n)^2 + a(n+1)^2 = a((n+1)^2) (Nov 22 2006), follow: a(n) + a(n+2k-1) + 4a(k-1) = (n+k)^2 + 6a(k-1) and a(n)^2 + a(n+2k-1)^2 + (4a(k-1))^2 + 3a(k-1) = a((n+k)^2 + 6a(k-1)). - Charlie Marion, Nov 27 2016
a(n) is also the greatest possible number of diagonals in a polyhedron with n+4 vertices. - Vladimir Letsko, Dec 19 2016
For n > 0, 2^5 * (binomial(n+1,2))^2 represents the first integer in a sum of 2*(2*n + 1)^2 consecutive integers that equals (2*n + 1)^6. - Patrick J. McNab, Dec 25 2016
Does not satisfy Benford's law (cf. Ross, 2012). - N. J. A. Sloane, Feb 12 2017
Number of ordered triples (a,b,c) of positive integers not larger than n such that a+b+c = 2n+1. - Aviel Livay, Feb 13 2017
Number of inequivalent tetrahedral face colorings using at most n colors so that no color appears only once. - David Nacin, Feb 22 2017
Also the Wiener index of the complete graph K_{n+1}. - Eric W. Weisstein, Sep 07 2017
Number of intersections between the Bernstein polynomials of degree n. - Eric Desbiaux, Apr 01 2018
a(n) is the area of a triangle with vertices at (1,1), (n+1,n+2), and ((n+1)^2, (n+2)^2). - Art Baker, Dec 06 2018
For n > 0, a(n) is the smallest k > 0 such that n divides numerator of (1/a(1) + 1/a(2) + ... + 1/a(n-1) + 1/k). It should be noted that 1/1 + 1/3 + 1/6 + ... + 2/(n(n+1)) = 2n/(n+1). - Thomas Ordowski, Aug 04 2019
Upper bound of the number of lines in an n-homogeneous supersolvable line arrangement (see Theorem 1.1 in Dimca). - Stefano Spezia, Oct 04 2019
For n > 0, a(n+1) is the number of lattice points on a triangular grid with side length n. - Wesley Ivan Hurt, Aug 12 2020
From Michael Chu, May 04 2022: (Start)
Maximum number of distinct nonempty substrings of a string of length n.
Maximum cardinality of the sumset A+A, where A is a set of n numbers. (End)
a(n) is the number of parking functions of size n avoiding the patterns 123, 132, and 312. - Lara Pudwell, Apr 10 2023
Suppose two rows, each consisting of n evenly spaced dots, are drawn in parallel. Suppose we bijectively draw lines between the dots of the two rows. For n >= 1, a(n - 1) is the maximal possible number of intersections between the lines. Equivalently, the maximal number of inversions in a permutation of [n]. - Sela Fried, Apr 18 2023
The following equation complements the generalization in Bala's Comment (Jan 05 2015). (2k + 1)^2*a(n) + a(k) = a((2k + 1)*n + k). - Charlie Marion, Aug 28 2023
a(n) + a(n+k) + a(k-1) + (k-1)*n = (n+k)^2. For k = 1, we have a(n) + a(n+1) = (n+1)^2. - Charlie Marion, Nov 17 2023
a(n+1)/3 is the expected number of steps to escape from a linear row of n positions starting at a random location and randomly performing steps -1 or +1 with equal probability. - Hugo Pfoertner, Jul 22 2025
a(n+1) is the number of nonnegative integer solutions to p + q + r = n. By Sylvester's law of inertia, it is also the number of congruence classes of real symmetric n-by-n matrices or equivalently, the number of symmetric bilinear forms on a real n-dimensional vector space. - Paawan Jethva, Jul 24 2025

Examples

			G.f.: x + 3*x^2 + 6*x^3 + 10*x^4 + 15*x^5 + 21*x^6 + 28*x^7 + 36*x^8 + 45*x^9 + ...
When n=3, a(3) = 4*3/2 = 6.
Example(a(4)=10): ABCD where A, B, C and D are different links in a chain or different amino acids in a peptide possible fragments: A, B, C, D, AB, ABC, ABCD, BC, BCD, CD = 10.
a(2): hollyhock leaves on the Tokugawa Mon, a(4): points in Pythagorean tetractys, a(5): object balls in eight-ball billiards. - _Bradley Klee_, Aug 24 2015
From _Gus Wiseman_, Oct 28 2020: (Start)
The a(1) = 1 through a(5) = 15 ordered triples of positive integers summing to n + 2 [Beeler, McGrath above] are the following. These compositions are ranked by A014311.
  (111)  (112)  (113)  (114)  (115)
         (121)  (122)  (123)  (124)
         (211)  (131)  (132)  (133)
                (212)  (141)  (142)
                (221)  (213)  (151)
                (311)  (222)  (214)
                       (231)  (223)
                       (312)  (232)
                       (321)  (241)
                       (411)  (313)
                              (322)
                              (331)
                              (412)
                              (421)
                              (511)
The unordered version is A001399(n-3) = A069905(n), with Heinz numbers A014612.
The strict case is A001399(n-6)*6, ranked by A337453.
The unordered strict case is A001399(n-6), with Heinz numbers A007304.
(End)
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 828.
  • C. Alsina and R. B. Nelson, Charming Proofs: A Journey into Elegant Mathematics, MAA, 2010. See Chapter 1.
  • T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 2.
  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 189.
  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, p. 109ff.
  • Marc Chamberland, Single Digits: In Praise of Small Numbers, Chapter 3, The Number Three, p. 72, Princeton University Press, 2015.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 155.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 33, 38, 40, 70.
  • J. M. De Koninck and A. Mercier, 1001 Problèmes en Théorie Classique des Nombres, Problème 309 pp 46-196, Ellipses, Paris, 2004
  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 6.
  • 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. 1.
  • Martin Gardner, Colossal Book of Mathematics, Chapter 34, Bulgarian Solitaire and Other Seemingly Endless Tasks, pp. 455-467, W. W. Norton & Company, 2001.
  • James Gleick, The Information: A History, A Theory, A Flood, Pantheon, 2011. [On page 82 mentions a table of the first 19999 triangular numbers published by E. de Joncort in 1762.]
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §4.6 Mathematical Proof and §8.6 Figurate Numbers, pp. 158-159, 289-290.
  • Cay S. Horstmann, Scala for the Impatient. Upper Saddle River, New Jersey: Addison-Wesley (2012): 171.
  • Elemer Labos, On the number of RGB-colors we can distinguish. Partition Spectra. Lecture at 7th Hungarian Conference on Biometry and Biomathematics. Budapest. Jul 06 2005.
  • A. Messiah, Quantum Mechanics, Vol.1, North Holland, Amsterdam, 1965, p. 457.
  • J. C. P. Miller, editor, Table of Binomial Coefficients. Royal Society Mathematical Tables, Vol. 3, Cambridge Univ. Press, 1954.
  • Alfred S. Posamentier, Math Charmers, Tantalizing Tidbits for the Mind, Prometheus Books, NY, 2003, pages 52-53, 129-132, 274.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 2-6, 13.
  • T. Trotter, Some Identities for the Triangular Numbers, Journal of Recreational Mathematics, Spring 1973, 6(2).
  • D. Wells, The Penguin Dictionary of Curious and Interesting Numbers, pp. 91-93 Penguin Books 1987.

Crossrefs

The figurate numbers, with parameter k as in the second Python program: A001477 (k=0), this sequence (k=1), A000290 (k=2), A000326 (k=3), A000384 (k=4), A000566 (k=5), A000567 (k=6), A001106 (k=7), A001107 (k=8).
a(n) = A110449(n, 0).
a(n) = A110555(n+2, 2).
A diagonal of A008291.
Column 2 of A195152.
Numbers of the form n*t(n+k,h)-(n+k)*t(n,h), where t(i,h) = i*(i+2*h+1)/2 for any h (for A000217 is k=1): A005563, A067728, A140091, A140681, A212331.
Boustrophedon transforms: A000718, A000746.
Iterations: A007501 (start=2), A013589 (start=4), A050542 (start=5), A050548 (start=7), A050536 (start=8), A050909 (start=9).
Cf. A002817 (doubly triangular numbers), A075528 (solutions of a(n)=a(m)/2).
Cf. A104712 (first column, starting with a(1)).
Some generalized k-gonal numbers are A001318 (k=5), this sequence (k=6), A085787 (k=7), etc.
A001399(n-3) = A069905(n) = A211540(n+2) counts 3-part partitions.
A001399(n-6) = A069905(n-3) = A211540(n-1) counts 3-part strict partitions.
A011782 counts compositions of any length.
A337461 counts pairwise coprime triples, with unordered version A307719.

Programs

  • Haskell
    a000217 n = a000217_list !! n
    a000217_list = scanl1 (+) [0..] -- Reinhard Zumkeller, Sep 23 2011
    
  • J
    a000217=: *-:@>: NB. Stephen Makdisi, May 02 2018
    
  • Magma
    [n*(n+1)/2: n in [0..60]]; // Bruno Berselli, Jul 11 2014
    
  • Magma
    [n: n in [0..1500] | IsSquare(8*n+1)]; // Juri-Stepan Gerasimov, Apr 09 2016
    
  • Maple
    A000217 := proc(n) n*(n+1)/2; end;
    istriangular:=proc(n) local t1; t1:=floor(sqrt(2*n)); if n = t1*(t1+1)/2 then return true else return false; end if; end proc; # N. J. A. Sloane, May 25 2008
    ZL := [S, {S=Prod(B, B, B), B=Set(Z, 1 <= card)}, unlabeled]:
    seq(combstruct[count](ZL, size=n), n=2..55); # Zerinvary Lajos, Mar 24 2007
    isA000217 := proc(n)
        issqr(1+8*n) ;
    end proc: # R. J. Mathar, Nov 29 2015 [This is the recipe Leonhard Euler proposes in chapter VII of his "Vollständige Anleitung zur Algebra", 1765. Peter Luschny, Sep 02 2022]
  • Mathematica
    Array[ #*(# - 1)/2 &, 54] (* Zerinvary Lajos, Jul 10 2009 *)
    FoldList[#1 + #2 &, 0, Range@ 50] (* Robert G. Wilson v, Feb 02 2011 *)
    Accumulate[Range[0,70]] (* Harvey P. Dale, Sep 09 2012 *)
    CoefficientList[Series[x / (1 - x)^3, {x, 0, 50}], x] (* Vincenzo Librandi, Jul 30 2014 *)
    (* For Mathematica 10.4+ *) Table[PolygonalNumber[n], {n, 0, 53}] (* Arkadiusz Wesolowski, Aug 27 2016 *)
    LinearRecurrence[{3, -3, 1}, {0, 1, 3}, 54] (* Robert G. Wilson v, Dec 04 2016 *)
    (* The following Mathematica program, courtesy of Steven J. Miller, is useful for testing if a sequence is Benford. To test a different sequence only one line needs to be changed. This strongly suggests that the triangular numbers are not Benford, since the second and third columns of the output disagree. - N. J. A. Sloane, Feb 12 2017 *)
    fd[x_] := Floor[10^Mod[Log[10, x], 1]]
    benfordtest[num_] := Module[{},
       For[d = 1, d <= 9, d++, digit[d] = 0];
       For[n = 1, n <= num, n++,
        {
         d = fd[n(n+1)/2];
         If[d != 0, digit[d] = digit[d] + 1];
         }];
       For[d = 1, d <= 9, d++, digit[d] = 1.0 digit[d]/num];
       For[d = 1, d <= 9, d++,
        Print[d, " ", 100.0 digit[d], " ", 100.0 Log[10, (d + 1)/d]]];
       ];
    benfordtest[20000]
    Table[Length[Join@@Permutations/@IntegerPartitions[n,{3}]],{n,0,15}] (* Gus Wiseman, Oct 28 2020 *)
  • PARI
    A000217(n) = n * (n + 1) / 2;
    
  • PARI
    is_A000217(n)=n*2==(1+n=sqrtint(2*n))*n \\ M. F. Hasler, May 24 2012
    
  • PARI
    is(n)=ispolygonal(n,3) \\ Charles R Greathouse IV, Feb 28 2014
    
  • PARI
    list(lim)=my(v=List(),n,t); while((t=n*n++/2)<=lim,listput(v,t)); Vec(v) \\ Charles R Greathouse IV, Jun 18 2021
    
  • Python
    for n in range(0,60): print(n*(n+1)//2, end=', ') # Stefano Spezia, Dec 06 2018
    
  • Python
    # Intended to compute the initial segment of the sequence, not
    # isolated terms. If in the iteration the line "x, y = x + y + 1, y + 1"
    # is replaced by "x, y = x + y + k, y + k" then the figurate numbers are obtained,
    # for k = 0 (natural A001477), k = 1 (triangular), k = 2 (squares), k = 3 (pentagonal), k = 4 (hexagonal), k = 5 (heptagonal), k = 6 (octagonal), etc.
    def aList():
        x, y = 1, 1
        yield 0
        while True:
            yield x
            x, y = x + y + 1, y + 1
    A000217 = aList()
    print([next(A000217) for i in range(54)]) # Peter Luschny, Aug 03 2019
  • SageMath
    [n*(n+1)/2 for n in (0..60)] # Bruno Berselli, Jul 11 2014
    
  • Scala
    (1 to 53).scanLeft(0)( + ) // Horstmann (2012), p. 171
    
  • Scheme
    (define (A000217 n) (/ (* n (+ n 1)) 2)) ;; Antti Karttunen, Jul 08 2017
    

Formula

G.f.: x/(1-x)^3. - Simon Plouffe in his 1992 dissertation
E.g.f.: exp(x)*(x+x^2/2).
a(n) = a(-1-n).
a(n) + a(n-1)*a(n+1) = a(n)^2. - Terrel Trotter, Jr., Apr 08 2002
a(n) = (-1)^n*Sum_{k=1..n} (-1)^k*k^2. - Benoit Cloitre, Aug 29 2002
a(n+1) = ((n+2)/n)*a(n), Sum_{n>=1} 1/a(n) = 2. - Jon Perry, Jul 13 2003
For n > 0, a(n) = A001109(n) - Sum_{k=0..n-1} (2*k+1)*A001652(n-1-k); e.g., 10 = 204 - (1*119 + 3*20 + 5*3 + 7*0). - Charlie Marion, Jul 18 2003
With interpolated zeros, this is n*(n+2)*(1+(-1)^n)/16. - Benoit Cloitre, Aug 19 2003
a(n+1) is the determinant of the n X n symmetric Pascal matrix M_(i, j) = binomial(i+j+1, i). - Benoit Cloitre, Aug 19 2003
a(n) = ((n+1)^3 - n^3 - 1)/6. - Xavier Acloque, Oct 24 2003
a(n) = a(n-1) + (1 + sqrt(1 + 8*a(n-1)))/2. This recursive relation is inverted when taking the negative branch of the square root, i.e., a(n) is transformed into a(n-1) rather than a(n+1). - Carl R. White, Nov 04 2003
a(n) = Sum_{k=1..n} phi(k)*floor(n/k) = Sum_{k=1..n} A000010(k)*A010766(n, k) (R. Dedekind). - Vladeta Jovovic, Feb 05 2004
a(n) + a(n+1) = (n+1)^2. - N. J. A. Sloane, Feb 19 2004
a(n) = a(n-2) + 2*n - 1. - Paul Barry, Jul 17 2004
a(n) = sqrt(Sum_{i=1..n} Sum_{j=1..n} (i*j)) = sqrt(A000537(n)). - Alexander Adamchuk, Oct 24 2004
a(n) = sqrt(sqrt(Sum_{i=1..n} Sum_{j=1..n} (i*j)^3)) = (Sum_{i=1..n} Sum_{j=1..n} Sum_{k=1..n} (i*j*k)^3)^(1/6). - Alexander Adamchuk, Oct 26 2004
a(n) == 1 (mod n+2) if n is odd and a(n) == n/2+2 (mod n+2) if n is even. - Jon Perry, Dec 16 2004
a(0) = 0, a(1) = 1, a(n) = 2*a(n-1) - a(n-2) + 1. - Miklos Kristof, Mar 09 2005
a(n) = a(n-1) + n. - Zak Seidov, Mar 06 2005
a(n) = A108299(n+3,4) = -A108299(n+4,5). - Reinhard Zumkeller, Jun 01 2005
a(n) = A111808(n,2) for n > 1. - Reinhard Zumkeller, Aug 17 2005
a(n)*a(n+1) = A006011(n+1) = (n+1)^2*(n^2+2)/4 = 3*A002415(n+1) = 1/2*a(n^2+2*n). a(n-1)*a(n) = (1/2)*a(n^2-1). - Alexander Adamchuk, Apr 13 2006 [Corrected and edited by Charlie Marion, Nov 26 2010]
a(n) = floor((2*n+1)^2/8). - Paul Barry, May 29 2006
For positive n, we have a(8*a(n))/a(n) = 4*(2*n+1)^2 = (4*n+2)^2, i.e., a(A033996(n))/a(n) = 4*A016754(n) = (A016825(n))^2 = A016826(n). - Lekraj Beedassy, Jul 29 2006
a(n)^2 + a(n+1)^2 = a((n+1)^2) [R B Nelsen, Math Mag 70 (2) (1997), p. 130]. - R. J. Mathar, Nov 22 2006
a(n) = A126890(n,0). - Reinhard Zumkeller, Dec 30 2006
a(n)*a(n+k)+a(n+1)*a(n+1+k) = a((n+1)*(n+1+k)). Generalizes previous formula dated Nov 22 2006 [and comments by J. M. Bergot dated May 22 2012]. - Charlie Marion, Feb 04 2011
(sqrt(8*a(n)+1)-1)/2 = n. - David W. Cantrell (DWCantrell(AT)sigmaxi.net), Feb 26 2007
a(n) = A023896(n) + A067392(n). - Lekraj Beedassy, Mar 02 2007
Sum_{k=0..n} a(k)*A039599(n,k) = A002457(n-1), for n >= 1. - Philippe Deléham, Jun 10 2007
8*a(n)^3 + a(n)^2 = Y(n)^2, where Y(n) = n*(n+1)*(2*n+1)/2 = 3*A000330(n). - Mohamed Bouhamida, Nov 06 2007 [Edited by Derek Orr, May 05 2015]
A general formula for polygonal numbers is P(k,n) = (k-2)*(n-1)n/2 + n = n + (k-2)*A000217(n-1), for n >= 1, k >= 3. - Omar E. Pol, Apr 28 2008 and Mar 31 2013
a(3*n) = A081266(n), a(4*n) = A033585(n), a(5*n) = A144312(n), a(6*n) = A144314(n). - Reinhard Zumkeller, Sep 17 2008
a(n) = A022264(n) - A049450(n). - Reinhard Zumkeller, Oct 09 2008
If we define f(n,i,a) = Sum_{j=0..k-1} (binomial(n,k)*Stirling1(n-k,i)*Product_{j=0..k-1} (-a-j)), then a(n) = -f(n,n-1,1), for n >= 1. - Milan Janjic, Dec 20 2008
4*a(x) + 4*a(y) + 1 = (x+y+1)^2 + (x-y)^2. - Vladimir Shevelev, Jan 21 2009
a(n) = A000124(n-1) + n-1 for n >= 2. a(n) = A000124(n) - 1. - Jaroslav Krizek, Jun 16 2009
An exponential generating function for the inverse of this sequence is given by Sum_{m>=0} ((Pochhammer(1, m)*Pochhammer(1, m))*x^m/(Pochhammer(3, m)*factorial(m))) = ((2-2*x)*log(1-x)+2*x)/x^2, the n-th derivative of which has a closed form which must be evaluated by taking the limit as x->0. A000217(n+1) = (lim_{x->0} d^n/dx^n (((2-2*x)*log(1-x)+2*x)/x^2))^-1 = (lim_{x->0} (2*Gamma(n)*(-1/x)^n*(n*(x/(-1+x))^n*(-x+1+n)*LerchPhi(x/(-1+x), 1, n) + (-1+x)*(n+1)*(x/(-1+x))^n + n*(log(1-x)+log(-1/(-1+x)))*(-x+1+n))/x^2))^-1. - Stephen Crowley, Jun 28 2009
a(n) = A034856(n+1) - A005408(n) = A005843(n) + A000124(n) - A005408(n). - Jaroslav Krizek, Sep 05 2009
a(A006894(n)) = a(A072638(n-1)+1) = A072638(n) = A006894(n+1)-1 for n >= 1. For n=4, a(11) = 66. - Jaroslav Krizek, Sep 12 2009
With offset 1, a(n) = floor(n^3/(n+1))/2. - Gary Detlefs, Feb 14 2010
a(n) = 4*a(floor(n/2)) + (-1)^(n+1)*floor((n+1)/2). - Bruno Berselli, May 23 2010
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3); a(0)=0, a(1)=1. - Mark Dols, Aug 20 2010
From Charlie Marion, Oct 15 2010: (Start)
a(n) + 2*a(n-1) + a(n-2) = n^2 + (n-1)^2; and
a(n) + 3*a(n-1) + 3*a(n-2) + a(n-3) = n^2 + 2*(n-1)^2 + (n-2)^2.
In general, for n >= m > 2, Sum_{k=0..m} binomial(m,m-k)*a(n-k) = Sum_{k=0..m-1} binomial(m-1,m-1-k)*(n-k)^2.
a(n) - 2*a(n-1) + a(n-2) = 1, a(n) - 3*a(n-1) + 3*a(n-2) - a(n-3) = 0 and a(n) - 4*a(n-1) + 6*a(n-2) - 4*(a-3) + a(n-4) = 0.
In general, for n >= m > 2, Sum_{k=0..m} (-1)^k*binomial(m,m-k)*a(n-k) = 0.
(End)
a(n) = sqrt(A000537(n)). - Zak Seidov, Dec 07 2010
For n > 0, a(n) = 1/(Integral_{x=0..Pi/2} 4*(sin(x))^(2*n-1)*(cos(x))^3). - Francesco Daddi, Aug 02 2011
a(n) = A110654(n)*A008619(n). - Reinhard Zumkeller, Aug 24 2011
a(2*k-1) = A000384(k), a(2*k) = A014105(k), k > 0. - Omar E. Pol, Sep 13 2011
a(n) = A026741(n)*A026741(n+1). - Charles R Greathouse IV, Apr 01 2012
a(n) + a(a(n)) + 1 = a(a(n)+1). - J. M. Bergot, Apr 27 2012
a(n) = -s(n+1,n), where s(n,k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
a(n)*a(n+1) = a(Sum_{m=1..n} A005408(m))/2, for n >= 1. For example, if n=8, then a(8)*a(9) = a(80)/2 = 1620. - Ivan N. Ianakiev, May 27 2012
a(n) = A002378(n)/2 = (A001318(n) + A085787(n))/2. - Omar E. Pol, Jan 11 2013
G.f.: x * (1 + 3x + 6x^2 + ...) = x * Product_{j>=0} (1+x^(2^j))^3 = x * A(x) * A(x^2) * A(x^4) * ..., where A(x) = (1 + 3x + 3x^2 + x^3). - Gary W. Adamson, Jun 26 2012
G.f.: G(0) where G(k) = 1 + (2*k+3)*x/(2*k+1 - x*(k+2)*(2*k+1)/(x*(k+2) + (k+1)/G(k+1))); (continued fraction, 3rd kind, 3-step). - Sergei N. Gladkovskii, Nov 23 2012
a(n) = A002088(n) + A063985(n). - Reinhard Zumkeller, Jan 21 2013
G.f.: x + 3*x^2/(Q(0)-3*x) where Q(k) = 1 + k*(x+1) + 3*x - x*(k+1)*(k+4)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Mar 14 2013
a(n) + a(n+1) + a(n+2) + a(n+3) + n = a(2*n+4). - Ivan N. Ianakiev, Mar 16 2013
a(n) + a(n+1) + ... + a(n+8) + 6*n = a(3*n+15). - Charlie Marion, Mar 18 2013
a(n) + a(n+1) + ... + a(n+20) + 2*n^2 + 57*n = a(5*n+55). - Charlie Marion, Mar 18 2013
3*a(n) + a(n-1) = a(2*n), for n > 0. - Ivan N. Ianakiev, Apr 05 2013
In general, a(k*n) = (2*k-1)*a(n) + a((k-1)*n-1). - Charlie Marion, Apr 20 2015
Also, a(k*n) = a(k)*a(n) + a(k-1)*a(n-1). - Robert Israel, Apr 20 2015
a(n+1) = det(binomial(i+2,j+1), 1 <= i,j <= n). - Mircea Merca, Apr 06 2013
a(n) = floor(n/2) + ceiling(n^2/2) = n - floor(n/2) + floor(n^2/2). - Wesley Ivan Hurt, Jun 15 2013
a(n) = floor((n+1)/(exp(2/(n+1))-1)). - Richard R. Forberg, Jun 22 2013
Sum_{n>=1} a(n)/n! = 3*exp(1)/2 by the e.g.f. Also see A067764 regarding ratios calculated this way for binomial coefficients in general. - Richard R. Forberg, Jul 15 2013
Sum_{n>=1} (-1)^(n+1)/a(n) = 4*log(2) - 2 = 0.7725887... . - Richard R. Forberg, Aug 11 2014
2/(Sum_{n>=m} 1/a(n)) = m, for m > 0. - Richard R. Forberg, Aug 12 2014
A228474(a(n))=n; A248952(a(n))=0; A248953(a(n))=a(n); A248961(a(n))=A000330(n). - Reinhard Zumkeller, Oct 20 2014
a(a(n)-1) + a(a(n+2)-1) + 1 = A000124(n+1)^2. - Charlie Marion, Nov 04 2014
a(n) = 2*A000292(n) - A000330(n). - Luciano Ancora, Mar 14 2015
a(n) = A007494(n-1) + A099392(n) for n > 0. - Bui Quang Tuan, Mar 27 2015
Sum_{k=0..n} k*a(k+1) = a(A000096(n+1)). - Charlie Marion, Jul 15 2015
Let O(n) be the oblong number n(n+1) = A002378(n) and S(n) the square number n^2 = A000290(n). Then a(n) + a(n+2k) = O(n+k) + S(k) and a(n) + a(n+2k+1) = S(n+k+1) + O(k). - Charlie Marion, Jul 16 2015
A generalization of the Nov 22 2006 formula, a(n)^2 + a(n+1)^2 = a((n+1)^2), follows. Let T(k,n) = a(n) + k. Then for all k, T(k,n)^2 + T(k,n+1)^2 = T(k,(n+1)^2 + 2*k) - 2*k. - Charlie Marion, Dec 10 2015
a(n)^2 + a(n+1)^2 = a(a(n) + a(n+1)). Deducible from N. J. A. Sloane's a(n) + a(n+1) = (n+1)^2 and R. B. Nelson's a(n)^2 + a(n+1)^2 = a((n+1)^2). - Ben Paul Thurston, Dec 28 2015
Dirichlet g.f.: (zeta(s-2) + zeta(s-1))/2. - Ilya Gutkovskiy, Jun 26 2016
a(n)^2 - a(n-1)^2 = n^3. - Miquel Cerda, Jun 29 2016
a(n) = A080851(0,n-1). - R. J. Mathar, Jul 28 2016
a(n) = A000290(n-1) - A034856(n-4). - Peter M. Chema, Sep 25 2016
a(n)^2 + a(n+3)^2 + 19 = a(n^2 + 4*n + 10). - Charlie Marion, Nov 23 2016
2*a(n)^2 + a(n) = a(n^2+n). - Charlie Marion, Nov 29 2016
G.f.: x/(1-x)^3 = (x * r(x) * r(x^3) * r(x^9) * r(x^27) * ...), where r(x) = (1 + x + x^2)^3 = (1 + 3*x + 6*x^2 + 7*x^3 + 6*x^4 + 3*x^5 + x^6). - Gary W. Adamson, Dec 03 2016
a(n) = sum of the elements of inverse of matrix Q(n), where Q(n) has elements q_i,j = 1/(1-4*(i-j)^2). So if e = appropriately sized vector consisting of 1's, then a(n) = e'.Q(n)^-1.e. - Michael Yukish, Mar 20 2017
a(n) = Sum_{k=1..n} ((2*k-1)!!*(2*n-2*k-1)!!)/((2*k-2)!!*(2*n-2*k)!!). - Michael Yukish, Mar 20 2017
Sum_{i=0..k-1} a(n+i) = (3*k*n^2 + 3*n*k^2 + k^3 - k)/6. - Christopher Hohl, Feb 23 2019
a(n) = A060544(n + 1) - A016754(n). - Ralf Steiner, Nov 09 2019
a(n) == 0 (mod n) iff n is odd (see De Koninck reference). - Bernard Schott, Jan 10 2020
8*a(k)*a(n) + ((a(k)-1)*n + a(k))^2 = ((a(k)+1)*n + a(k))^2. This formula reduces to the well-known formula, 8*a(n) + 1 = (2*n+1)^2, when k = 1. - Charlie Marion, Jul 23 2020
a(k)*a(n) = Sum_{i = 0..k-1} (-1)^i*a((k-i)*(n-i)). - Charlie Marion, Dec 04 2020
From Amiram Eldar, Jan 20 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = cosh(sqrt(7)*Pi/2)/(2*Pi).
Product_{n>=2} (1 - 1/a(n)) = 1/3. (End)
a(n) = Sum_{k=1..2*n-1} (-1)^(k+1)*a(k)*a(2*n-k). For example, for n = 4, 1*28 - 3*21 + 6*15 - 10*10 + 15*6 - 21*3 + 28*1 = 10. - Charlie Marion, Mar 23 2022
2*a(n) = A000384(n) - n^2 + 2*n. In general, if P(k,n) = the n-th k-gonal number, then (j+1)*a(n) = P(5 + j, n) - n^2 + (j+1)*n. More generally, (j+1)*P(k,n) = P(2*k + (k-2)*(j-1),n) - n^2 + (j+1)*n. - Charlie Marion, Mar 14 2023
a(n) = A109613(n) * A004526(n+1). - Torlach Rush, Nov 10 2023
a(n) = (1/6)* Sum_{k = 0..3*n} (-1)^(n+k+1) * k*(k + 1) * binomial(3*n+k, 2*k). - Peter Bala, Nov 03 2024
From Peter Bala, Jul 05 2025: (Start)
The following series telescope: for k >= 0,
Sum_{n >= 1} a(n)*a(n+2)*...*a(n+2*k)/(a(n+1)*a(n+3)*...*a(n+2*k+3)) = 1/(2*k + 3);
Sum_{n >= 1} a(n+1)*a(n+3)*...*a(n+2*k+1)/(a(n)*a(n+2)*...*a(n+2*k+2)) = 2/(2*k + 3) * Sum_{i = 1..2*k+3} 1/i. (End)

Extensions

Edited by Derek Orr, May 05 2015

A002260 Triangle read by rows: T(n,k) = k for n >= 1, k = 1..n.

Original entry on oeis.org

1, 1, 2, 1, 2, 3, 1, 2, 3, 4, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 6, 1, 2, 3, 4, 5, 6, 7, 1, 2, 3, 4, 5, 6, 7, 8, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14
Offset: 1

Views

Author

Angele Hamel (amh(AT)maths.soton.ac.uk)

Keywords

Comments

Old name: integers 1 to k followed by integers 1 to k+1 etc. (a fractal sequence).
Start counting again and again.
This is a "doubly fractal sequence" - see the Franklin T. Adams-Watters link.
The PARI functions t1, t2 can be used to read a square array T(n,k) (n >= 1, k >= 1) by antidiagonals downwards: n -> T(t1(n), t2(n)). - Michael Somos, Aug 23 2002
Reading this sequence as the antidiagonals of a rectangular array, row n is (n,n,n,...); this is the weight array (Cf. A144112) of the array A127779 (rectangular). - Clark Kimberling, Sep 16 2008
The upper trim of an arbitrary fractal sequence s is s, but the lower trim of s, although a fractal sequence, need not be s itself. However, the lower trim of A002260 is A002260. (The upper trim of s is what remains after the first occurrence of each term is deleted; the lower trim of s is what remains after all 0's are deleted from the sequence s-1.) - Clark Kimberling, Nov 02 2009
Eigensequence of the triangle = A001710 starting (1, 3, 12, 60, 360, ...). - Gary W. Adamson, Aug 02 2010
The triangle sums, see A180662 for their definitions, link this triangle of natural numbers with twenty-three different sequences, see the crossrefs. The mirror image of this triangle is A004736. - Johannes W. Meijer, Sep 22 2010
A002260 is the self-fission of the polynomial sequence (q(n,x)), where q(n,x) = x^n + x^(n-1) + ... + x + 1. See A193842 for the definition of fission. - Clark Kimberling, Aug 07 2011
Sequence B is called a reluctant sequence of sequence A, if B is triangle array read by rows: row number k coincides with first k elements of the sequence A. Sequence A002260 is reluctant sequence of sequence 1,2,3,... (A000027). - Boris Putievskiy, Dec 12 2012
This is the maximal sequence of positive integers, such that once an integer k has occurred, the number of k's always exceeds the number of (k+1)'s for the remainder of the sequence, with the first occurrence of the integers being in order. - Franklin T. Adams-Watters, Oct 23 2013
A002260 are the k antidiagonal numerators of rationals in Cantor's proof of 1-to-1 correspondence between rationals and naturals; the denominators are k-numerator+1. - Adriano Caroli, Mar 24 2015
T(n,k) gives the distance to the largest triangular number < n. - Ctibor O. Zizka, Apr 09 2020

Examples

			First six rows:
  1
  1   2
  1   2   3
  1   2   3   4
  1   2   3   4   5
  1   2   3   4   5   6
		

References

  • Clark Kimberling, "Fractal sequences and interspersions," Ars Combinatoria 45 (1997) 157-168. (Introduces upper trimming, lower trimming, and signature sequences.)
  • M. Myers, Smarandache Crescendo Subsequences, R. H. Wilde, An Anthology in Memoriam, Bristol Banner Books, Bristol, 1998, p. 19.
  • F. Smarandache, Sequences of Numbers Involved in Unsolved Problems, Hexis, Phoenix, 2006.

Crossrefs

Cf. A140756 (alternating signs).
Triangle sums (see the comments): A000217 (Row1, Kn11); A004526 (Row2); A000096 (Kn12); A055998 (Kn13); A055999 (Kn14); A056000 (Kn15); A056115 (Kn16); A056119 (Kn17); A056121 (Kn18); A056126 (Kn19); A051942 (Kn110); A101859 (Kn111); A132754 (Kn112); A132755 (Kn113); A132756 (Kn114); A132757 (Kn115); A132758 (Kn116); A002620 (Kn21); A000290 (Kn3); A001840 (Ca2); A000326 (Ca3); A001972 (Gi2); A000384 (Gi3).
Cf. A108872.

Programs

  • Haskell
    a002260 n k = k
    a002260_row n = [1..n]
    a002260_tabl = iterate (\row -> map (+ 1) (0 : row)) [1]
    -- Reinhard Zumkeller, Aug 04 2014, Jul 03 2012
    
  • Maple
    at:=0; for n from 1 to 150 do for i from 1 to n do at:=at+1; lprint(at,i); od: od: # N. J. A. Sloane, Nov 01 2006
    seq(seq(i,i=1..k),k=1..13); # Peter Luschny, Jul 06 2009
  • Mathematica
    FoldList[{#1, #2} &, 1, Range[2, 13]] // Flatten (* Robert G. Wilson v, May 10 2011 *)
    Flatten[Table[Range[n],{n,20}]] (* Harvey P. Dale, Jun 20 2013 *)
  • Maxima
    T(n,k):=sum((i+k)*binomial(i+k-1,i)*binomial(k,n-i-k+1)*(-1)^(n-i-k+1),i,max(0,n+1-2*k),n-k+1); /* Vladimir Kruchinin, Oct 18 2013 */
    
  • PARI
    t1(n)=n-binomial(floor(1/2+sqrt(2*n)),2) /* this sequence */
    
  • PARI
    A002260(n)=n-binomial((sqrtint(8*n)+1)\2,2) \\ M. F. Hasler, Mar 10 2014
    
  • Python
    from math import isqrt, comb
    def A002260(n): return n-comb((m:=isqrt(k:=n<<1))+(k>m*(m+1)),2) # Chai Wah Wu, Nov 08 2024

Formula

a(n) = 1 + A002262(n).
n-th term is n - m*(m+1)/2 + 1, where m = floor((sqrt(8*n+1) - 1) / 2).
The above formula is for offset 0; for offset 1, use a(n) = n-m*(m+1)/2 where m = floor((-1+sqrt(8*n-7))/2). - Clark Kimberling, Jun 14 2011
a(k * (k + 1) / 2 + i) = i for k >= 0 and 0 < i <= k + 1. - Reinhard Zumkeller, Aug 14 2001
a(n) = (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2. - Brian Tenneson, Oct 11 2003
a(n) = n - binomial(floor((1+sqrt(8*n))/2), 2). - Paul Barry, May 25 2004
T(n,k) = A001511(A118413(n,k)); T(n,k) = A003602(A118416(n,k)). - Reinhard Zumkeller, Apr 27 2006
a(A000217(n)) = A000217(n) - A000217(n-1), a(A000217(n-1) + 1) = 1, a(A000217(n) - 1) = A000217(n) - A000217(n-1) - 1. - Alexander R. Povolotsky, May 28 2008
a(A169581(n)) = A038566(n). - Reinhard Zumkeller, Dec 02 2009
T(n,k) = Sum_{i=1..k} i*binomial(k,i)*binomial(n-k,n-i) (regarded as triangle, see the example). - Mircea Merca, Apr 11 2012
T(n,k) = Sum_{i=max(0,n+1-2*k)..n-k+1} (i+k)*binomial(i+k-1,i)*binomial(k,n-i-k+1)*(-1)^(n-i-k+1). - Vladimir Kruchinin, Oct 18 2013
G.f.: x*y / ((1 - x) * (1 - x*y)^2) = Sum_{n,k>0} T(n,k) * x^n * y^k. - Michael Somos, Sep 17 2014
a(n) = n - S(n) where S(n) = sum of distinct terms in {a(1), a(2), ..., a(n-1)}. - David James Sycamore, Mar 10 2025

Extensions

More terms from Reinhard Zumkeller, Apr 27 2006
Incorrect program removed by Franklin T. Adams-Watters, Mar 19 2010
New name from Omar E. Pol, Jul 15 2012

A007304 Sphenic numbers: products of 3 distinct primes.

Original entry on oeis.org

30, 42, 66, 70, 78, 102, 105, 110, 114, 130, 138, 154, 165, 170, 174, 182, 186, 190, 195, 222, 230, 231, 238, 246, 255, 258, 266, 273, 282, 285, 286, 290, 310, 318, 322, 345, 354, 357, 366, 370, 374, 385, 399, 402, 406, 410, 418, 426, 429, 430, 434, 435, 438
Offset: 1

Views

Author

Keywords

Comments

Note the distinctions between this and "n has exactly three prime factors" (A014612) or "n has exactly three distinct prime factors." (A033992). The word "sphenic" also means "shaped like a wedge" [American Heritage Dictionary] as in dentation with "sphenic molars." - Jonathan Vos Post, Sep 11 2005
Also the volume of a sphenic brick. A sphenic brick is a rectangular parallelepiped whose sides are components of a sphenic number, namely whose sides are three distinct primes. Example: The distinct prime triple (3,5,7) produces a 3x5x7 unit brick which has volume 105 cubic units. 3-D analog of 2-D A037074 Product of twin primes, per Cino Hilliard's comment. Compare with 3-D A107768 Golden 3-almost primes = Volumes of bricks (rectangular parallelepipeds) each of whose faces has golden semiprime area. - Jonathan Vos Post, Jan 08 2007
Sum(n>=1, 1/a(n)^s) = (1/6)*(P(s)^3 - P(3*s) - 3*(P(s)*P(2*s)-P(3*s))), where P is prime zeta function. - Enrique Pérez Herrero, Jun 28 2012
Also numbers n with A001222(n)=3 and A001221(n)=3. - Enrique Pérez Herrero, Jun 28 2012
n = 265550 is the smallest n with a(n) (=1279789) < A006881(n) (=1279793). - Peter Dolland, Apr 11 2020

Examples

			From _Gus Wiseman_, Nov 05 2020: (Start)
Also Heinz numbers of strict integer partitions into three parts, where the Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). These partitions are counted by A001399(n-6) = A069905(n-3), with ordered version A001399(n-6)*6. The sequence of terms together with their prime indices begins:
     30: {1,2,3}     182: {1,4,6}     286: {1,5,6}
     42: {1,2,4}     186: {1,2,11}    290: {1,3,10}
     66: {1,2,5}     190: {1,3,8}     310: {1,3,11}
     70: {1,3,4}     195: {2,3,6}     318: {1,2,16}
     78: {1,2,6}     222: {1,2,12}    322: {1,4,9}
    102: {1,2,7}     230: {1,3,9}     345: {2,3,9}
    105: {2,3,4}     231: {2,4,5}     354: {1,2,17}
    110: {1,3,5}     238: {1,4,7}     357: {2,4,7}
    114: {1,2,8}     246: {1,2,13}    366: {1,2,18}
    130: {1,3,6}     255: {2,3,7}     370: {1,3,12}
    138: {1,2,9}     258: {1,2,14}    374: {1,5,7}
    154: {1,4,5}     266: {1,4,8}     385: {3,4,5}
    165: {2,3,5}     273: {2,4,6}     399: {2,4,8}
    170: {1,3,7}     282: {1,2,15}    402: {1,2,19}
    174: {1,2,10}    285: {2,3,8}     406: {1,4,10}
(End)
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • "Sphenic", The American Heritage Dictionary of the English Language, Fourth Edition, Houghton Mifflin Company, 2000.

Crossrefs

Products of exactly k distinct primes, for k = 1 to 6: A000040, A006881. A007304, A046386, A046387, A067885.
Cf. A162143 (a(n)^2).
For the following, NNS means "not necessarily strict".
A014612 is the NNS version.
A046389 is the restriction to odds (NNS: A046316).
A075819 is the restriction to evens (NNS: A075818).
A239656 gives first differences.
A285508 lists terms of A014612 that are not squarefree.
A307534 is the case where all prime indices are odd (NNS: A338471).
A337453 is a different ranking of ordered triples (NNS: A014311).
A338557 is the case where all prime indices are even (NNS: A338556).
A001399(n-6) counts strict 3-part partitions (NNS: A001399(n-3)).
A005117 lists squarefree numbers.
A008289 counts strict partitions by sum and length.
A220377 counts 3-part pairwise coprime strict partitions (NNS: A307719).

Programs

  • Haskell
    a007304 n = a007304_list !! (n-1)
    a007304_list = filter f [1..] where
    f u = p < q && q < w && a010051 w == 1 where
    p = a020639 u; v = div u p; q = a020639 v; w = div v q
    -- Reinhard Zumkeller, Mar 23 2014
    
  • Maple
    with(numtheory): a:=proc(n) if bigomega(n)=3 and nops(factorset(n))=3 then n else fi end: seq(a(n),n=1..450); # Emeric Deutsch
    A007304 := proc(n)
        option remember;
        local a;
        if n =1 then
            30;
        else
            for a from procname(n-1)+1 do
                if bigomega(a)=3 and nops(factorset(a))=3 then
                    return a;
                end if;
            end do:
        end if;
    end proc: # R. J. Mathar, Dec 06 2016
    is_a := proc(n) local P; P := NumberTheory:-PrimeFactors(n); nops(P) = 3 and n = mul(P) end:
    A007304List := upto -> select(is_a, [seq(1..upto)]):  # Peter Luschny, Apr 14 2025
  • Mathematica
    Union[Flatten[Table[Prime[n]*Prime[m]*Prime[k], {k, 20}, {n, k+1, 20}, {m, n+1, 20}]]]
    Take[ Sort@ Flatten@ Table[ Prime@i Prime@j Prime@k, {i, 3, 21}, {j, 2, i - 1}, {k, j - 1}], 53] (* Robert G. Wilson v *)
    With[{upto=500},Sort[Select[Times@@@Subsets[Prime[Range[Ceiling[upto/6]]],{3}],#<=upto&]]] (* Harvey P. Dale, Jan 08 2015 *)
    Select[Range[100],SquareFreeQ[#]&&PrimeOmega[#]==3&] (* Gus Wiseman, Nov 05 2020 *)
  • PARI
    for(n=1,1e4,if(bigomega(n)==3 && omega(n)==3,print1(n", "))) \\ Charles R Greathouse IV, Jun 10 2011
    
  • PARI
    list(lim)=my(v=List(),t);forprime(p=2,(lim)^(1/3),forprime(q=p+1,sqrt(lim\p),t=p*q;forprime(r=q+1,lim\t,listput(v,t*r))));vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 20 2011
    
  • PARI
    list(lim)=my(v=List(), t); forprime(p=2, sqrtnint(lim\=1,3), forprime(q=p+1, sqrtint(lim\p), t=p*q; forprime(r=q+1, lim\t, listput(v, t*r)))); Set(v) \\ Charles R Greathouse IV, Jan 21 2025
    
  • Python
    from math import isqrt
    from sympy import primepi, primerange, integer_nthroot
    def A007304(n):
        def f(x): return int(n+x-sum(primepi(x//(k*m))-b for a,k in enumerate(primerange(integer_nthroot(x,3)[0]+1),1) for b,m in enumerate(primerange(k+1,isqrt(x//k)+1),a+1)))
        kmin, kmax = 0,1
        while f(kmax) > kmax:
            kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax # Chai Wah Wu, Aug 29 2024
    
  • SageMath
    def is_a(n):
        P = prime_divisors(n)
        return len(P) == 3 and prod(P) == n
    print([n for n in range(1, 439) if is_a(n)]) # Peter Luschny, Apr 14 2025

Formula

A008683(a(n)) = -1.
A000005(a(n)) = 8. - R. J. Mathar, Aug 14 2009
A002033(a(n)-1) = 13. - Juri-Stepan Gerasimov, Oct 07 2009, R. J. Mathar, Oct 14 2009
A178254(a(n)) = 36. - Reinhard Zumkeller, May 24 2010
A050326(a(n)) = 5, subsequence of A225228. - Reinhard Zumkeller, May 03 2013
a(n) ~ 2n log n/(log log n)^2. - Charles R Greathouse IV, Sep 14 2015

Extensions

More terms from Robert G. Wilson v, Jan 04 2006
Comment concerning number of divisors corrected by R. J. Mathar, Aug 14 2009

A001399 a(n) is the number of partitions of n into at most 3 parts; also partitions of n+3 in which the greatest part is 3; also number of unlabeled multigraphs with 3 nodes and n edges.

Original entry on oeis.org

1, 1, 2, 3, 4, 5, 7, 8, 10, 12, 14, 16, 19, 21, 24, 27, 30, 33, 37, 40, 44, 48, 52, 56, 61, 65, 70, 75, 80, 85, 91, 96, 102, 108, 114, 120, 127, 133, 140, 147, 154, 161, 169, 176, 184, 192, 200, 208, 217, 225, 234, 243, 252, 261, 271, 280, 290, 300, 310, 320, 331, 341
Offset: 0

Views

Author

Keywords

Comments

Also number of tripods (trees with exactly 3 leaves) on n vertices. - Eric W. Weisstein, Mar 05 2011
Also number of partitions of n+3 into exactly 3 parts; number of partitions of n in which the greatest part is less than or equal to 3; and the number of nonnegative solutions to b + 2c + 3d = n.
Also a(n) gives number of partitions of n+6 into 3 distinct parts and number of partitions of 2n+9 into 3 distinct and odd parts, e.g., 15 = 11 + 3 + 1 = 9 + 5 + 1 = 7 + 5 + 3. - Jon Perry, Jan 07 2004
Also bracelets with n+3 beads 3 of which are red (so there are 2 possibilities with 5 beads).
More generally, the number of partitions of n into at most k parts is also the number of partitions of n+k into k positive parts, the number of partitions of n+k in which the greatest part is k, the number of partitions of n in which the greatest part is less than or equal to k, the number of partitions of n+k(k+1)/2 into exactly k distinct positive parts, the number of nonnegative solutions to b + 2c + 3d + ... + kz = n and the number of nonnegative solutions to 2c + 3d + ... + kz <= n. - Henry Bottomley, Apr 17 2001
Also coefficient of q^n in the expansion of (m choose 3)_q as m goes to infinity. - Y. Kelly Itakura (yitkr(AT)mta.ca), Aug 21 2002
From Winston C. Yang (winston(AT)cs.wisc.edu), Apr 30 2002: (Start)
Write 1,2,3,4,... in a hexagonal spiral around 0, then a(n) for n > 0 is formed by the folding points (including the initial 1). The spiral begins:
.
85--84--83--82--81--80
/ \
86 56--55--54--53--52 79
/ / \ \
87 57 33--32--31--30 51 78
/ / / \ \ \
88 58 34 16--15--14 29 50 77
/ / / / \ \ \ \
89 59 35 17 5---4 13 28 49 76
/ / / / / \ \ \ \ \
90 60 36 18 6 0 3 12 27 48 75
/ / / / / / / / / / /
91 61 37 19 7 1---2 11 26 47 74
\ \ \ \ / / / /
62 38 20 8---9--10 25 46 73
\ \ \ / / /
63 39 21--22--23--24 45 72
\ \ / /
64 40--41--42--43--44 71
\ /
65--66--67--68--69--70
.
a(p) is maximal number of hexagons in a polyhex with perimeter at most 2p + 6. (End)
a(n-3) is the number of partitions of n into 3 distinct parts, where 0 is allowed as a part. E.g., at n=9, we can write 8+1+0, 7+2+0, 6+3+0, 4+5+0, 1+2+6, 1+3+5 and 2+3+4, which is a(6)=7. - Jon Perry, Jul 08 2003
a(n) gives number of partitions of n+6 into parts <=3 where each part is used at least once (subtract 6=1+2+3 from n). - Jon Perry, Jul 03 2004
This is also the number of partitions of n+3 into exactly 3 parts (there is a 1-to-1 correspondence between the number of partitions of n+3 in which the greatest part is 3 and the number of partitions of n+3 into exactly three parts). - Graeme McRae, Feb 07 2005
Apply the Riordan array (1/(1-x^3),x) to floor((n+2)/2). - Paul Barry, Apr 16 2005
Also, number of triangles that can be created with odd perimeter 3,5,7,9,11,... with all sides whole numbers. Note that triangles with even perimeter can be generated from the odd ones by increasing each side by 1. E.g., a(1) = 1 because perimeter 3 can make {1,1,1} 1 triangle. a(4) = 3 because perimeter 9 can make {1,4,4} {2,3,4} {3,3,3} 3 possible triangles. - Bruce Love (bruce_love(AT)ofs.edu.sg), Nov 20 2006
Also number of nonnegative solutions of the Diophantine equation x+2*y+3*z=n, cf. Pólya/Szegő reference.
From Vladimir Shevelev, Apr 23 2011: (Start)
Also a(n-3), n >= 3, is the number of non-equivalent necklaces of 3 beads each of them painted by one of n colors.
The sequence {a(n-3), n >= 3} solves the so-called Reis problem about convex k-gons in case k=3 (see our comment to A032279).
a(n-3) (n >= 3) is an essentially unimprovable upper estimate for the number of distinct values of the permanent in (0,1)-circulants of order n with three 1's in every row. (End)
A001399(n) is the number of 3-tuples (w,x,y) having all terms in {0,...,n} and w = 2*x+3*y. - Clark Kimberling, Jun 04 2012
Also, for n >= 3, a(n-3) is the number of the distinct triangles in an n-gon, see the Ngaokrajang links. - Kival Ngaokrajang, Mar 16 2013
Also, a(n) is the total number of 5-curve coin patterns (5C4S type: 5 curves covering full 4 coins and symmetry) packing into fountain of coins base (n+3). See illustration in links. - Kival Ngaokrajang, Oct 16 2013
Also a(n) = half the number of minimal zero sequences for Z_n of length 3 [Ponomarenko]. - N. J. A. Sloane, Feb 25 2014
Also, a(n) equals the number of linearly-independent terms at 2n-th order in the power series expansion of an Octahedral Rotational Energy Surface (cf. Harter & Patterson). - Bradley Klee, Jul 31 2015
Also Molien series for invariants of finite Coxeter groups D_3 and A_3. - N. J. A. Sloane, Jan 10 2016
Number of different distributions of n+6 identical balls in 3 boxes as x,y,z where 0 < x < y < z. - Ece Uslu and Esin Becenen, Jan 11 2016
a(n) is also the number of partitions of 2*n with <= n parts and no part >= 4. The bijection to partitions of n with no part >= 4 is: 1 <-> 2, 2 <-> 1 + 3, 3 <-> 3 + 3 (observing the order of these rules). The <- direction uses the following fact for partitions of 2*n with <= n parts and no part >=4: for each part 1 there is a part 3, and an even number (including 0) of remaining parts 3. - Wolfdieter Lang, May 21 2019
List of the terms in A000567(n>=1), A049450(n>=1), A033428(n>=1), A049451(n>=1), A045944(n>=1), and A003215(n) in nondecreasing order. List of the numbers A056105(n)-1, A056106(n)-1, A056107(n)-1, A056108(n)-1, A056109(n)-1, and A003215(m) with n >= 1 and m >= 0 in nondecreasing order. Numbers of the forms 3n*(n-1)+1, n*(3n-2), n*(3n-1), 3n^2, n*(3n+1), n*(3n+2) with n >= 1 listed in nondecreasing order. Integers m such that lattice points from 1 through m on a hexagonal spiral starting at 1 forms a convex polygon. - Ya-Ping Lu, Jan 24 2024

Examples

			G.f. = 1 + x + 2*x^2 + 3*x^3 + 4*x^4 + 5*x^5 + 7*x^6 + 8*x^7 + 10*x^8 + 12*x^9 + ...
Recall that in a necklace the adjacent beads have distinct colors. Suppose we have n colors with labels 1,...,n. Two colorings of the beads are equivalent if the cyclic sequences of the distances modulo n between labels of adjacent colors have the same period. If n=4, all colorings are equivalent. E.g., for the colorings {1,2,3} and {1,2,4} we have the same period {1,1,2} of distances modulo 4. So, a(n-3)=a(1)=1. If n=5, then we have two such periods {1,1,3} and {1,2,2} modulo 5. Thus a(2)=2. - _Vladimir Shevelev_, Apr 23 2011
a(0) = 1, i.e., {1,2,3} Number of different distributions of 6 identical balls to 3 boxes as x,y and z where 0 < x < y < z. - _Ece Uslu_, Esin Becenen, Jan 11 2016
a(3) = 3, i.e., {1,2,6}, {1,3,5}, {2,3,4} Number of different distributions of 9 identical balls in 3 boxes as x,y and z where 0 < x < y < z. - _Ece Uslu_, Esin Becenen, Jan 11 2016
From _Gus Wiseman_, Apr 15 2019: (Start)
The a(0) = 1 through a(8) = 10 integer partitions of n with at most three parts are the following. The Heinz numbers of these partitions are given by A037144.
  ()  (1)  (2)   (3)    (4)    (5)    (6)    (7)    (8)
           (11)  (21)   (22)   (32)   (33)   (43)   (44)
                 (111)  (31)   (41)   (42)   (52)   (53)
                        (211)  (221)  (51)   (61)   (62)
                               (311)  (222)  (322)  (71)
                                      (321)  (331)  (332)
                                      (411)  (421)  (422)
                                             (511)  (431)
                                                    (521)
                                                    (611)
The a(0) = 1 through a(7) = 8 integer partitions of n + 3 whose greatest part is 3 are the following. The Heinz numbers of these partitions are given by A080193.
  (3)  (31)  (32)   (33)    (322)    (332)     (333)      (3322)
             (311)  (321)   (331)    (3221)    (3222)     (3331)
                    (3111)  (3211)   (3311)    (3321)     (32221)
                            (31111)  (32111)   (32211)    (33211)
                                     (311111)  (33111)    (322111)
                                               (321111)   (331111)
                                               (3111111)  (3211111)
                                                          (31111111)
Non-isomorphic representatives of the a(0) = 1 through a(5) = 5 unlabeled multigraphs with 3 vertices and n edges are the following.
  {}  {12}  {12,12}  {12,12,12}  {12,12,12,12}  {12,12,12,12,12}
            {13,23}  {12,13,23}  {12,13,23,23}  {12,13,13,23,23}
                     {13,23,23}  {13,13,23,23}  {12,13,23,23,23}
                                 {13,23,23,23}  {13,13,23,23,23}
                                                {13,23,23,23,23}
The a(0) = 1 through a(8) = 10 strict integer partitions of n - 6 with three parts are the following (A = 10, B = 11). The Heinz numbers of these partitions are given by A007304.
  (321)  (421)  (431)  (432)  (532)  (542)  (543)  (643)   (653)
                (521)  (531)  (541)  (632)  (642)  (652)   (743)
                       (621)  (631)  (641)  (651)  (742)   (752)
                              (721)  (731)  (732)  (751)   (761)
                                     (821)  (741)  (832)   (842)
                                            (831)  (841)   (851)
                                            (921)  (931)   (932)
                                                   (A21)   (941)
                                                           (A31)
                                                           (B21)
The a(0) = 1 through a(8) = 10 integer partitions of n + 3 with three parts are the following. The Heinz numbers of these partitions are given by A014612.
  (111)  (211)  (221)  (222)  (322)  (332)  (333)  (433)  (443)
                (311)  (321)  (331)  (422)  (432)  (442)  (533)
                       (411)  (421)  (431)  (441)  (532)  (542)
                              (511)  (521)  (522)  (541)  (551)
                                     (611)  (531)  (622)  (632)
                                            (621)  (631)  (641)
                                            (711)  (721)  (722)
                                                   (811)  (731)
                                                          (821)
                                                          (911)
The a(0) = 1 through a(8) = 10 integer partitions of n whose greatest part is <= 3 are the following. The Heinz numbers of these partitions are given by A051037.
  ()  (1)  (2)   (3)    (22)    (32)     (33)      (322)      (332)
           (11)  (21)   (31)    (221)    (222)     (331)      (2222)
                 (111)  (211)   (311)    (321)     (2221)     (3221)
                        (1111)  (2111)   (2211)    (3211)     (3311)
                                (11111)  (3111)    (22111)    (22211)
                                         (21111)   (31111)    (32111)
                                         (111111)  (211111)   (221111)
                                                   (1111111)  (311111)
                                                              (2111111)
                                                              (11111111)
The a(0) = 1 through a(6) = 7 strict integer partitions of 2n+9 with 3 parts, all of which are odd, are the following. The Heinz numbers of these partitions are given by A307534.
  (5,3,1)  (7,3,1)  (7,5,1)  (7,5,3)   (9,5,3)   (9,7,3)   (9,7,5)
                    (9,3,1)  (9,5,1)   (9,7,1)   (11,5,3)  (11,7,3)
                             (11,3,1)  (11,5,1)  (11,7,1)  (11,9,1)
                                       (13,3,1)  (13,5,1)  (13,5,3)
                                                 (15,3,1)  (13,7,1)
                                                           (15,5,1)
                                                           (17,3,1)
The a(0) = 1 through a(8) = 10 strict integer partitions of n + 3 with 3 parts where 0 is allowed as a part (A = 10):
  (210)  (310)  (320)  (420)  (430)  (530)  (540)  (640)  (650)
                (410)  (510)  (520)  (620)  (630)  (730)  (740)
                       (321)  (610)  (710)  (720)  (820)  (830)
                              (421)  (431)  (810)  (910)  (920)
                                     (521)  (432)  (532)  (A10)
                                            (531)  (541)  (542)
                                            (621)  (631)  (632)
                                                   (721)  (641)
                                                          (731)
                                                          (821)
The a(0) = 1 through a(7) = 7 integer partitions of n + 6 whose distinct parts are 1, 2, and 3 are the following. The Heinz numbers of these partitions are given by A143207.
  (321)  (3211)  (3221)   (3321)    (32221)    (33221)     (33321)
                 (32111)  (32211)   (33211)    (322211)    (322221)
                          (321111)  (322111)   (332111)    (332211)
                                    (3211111)  (3221111)   (3222111)
                                               (32111111)  (3321111)
                                                           (32211111)
                                                           (321111111)
(End)
Partitions of 2*n with <= n parts and no part >= 4: a(3) = 3 from (2^3), (1,2,3), (3^2) mapping to (1^3), (1,2), (3), the partitions of 3 with no part >= 4, respectively. - _Wolfdieter Lang_, May 21 2019
		

References

  • R. Ayoub, An Introduction to the Analytic Theory of Numbers, Amer. Math. Soc., 1963; Chapter III, Problem 33.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 110, D(n); page 263, #18, P_n^{3}.
  • J. L. Gross and J. Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 517.
  • H. Gupta et al., Tables of Partitions. Royal Society Mathematical Tables, Vol. 4, Cambridge Univ. Press, 1958, p. 2.
  • F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, p. 88, (4.1.18).
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 275.
  • R. Honsberger, Mathematical Gems III, Math. Assoc. Amer., 1985, p. 39.
  • J. H. van Lint, Combinatorial Seminar Eindhoven, Lecture Notes Math., 382 (1974), see pp. 33-34.
  • G. Pólya and G. Szegő, Problems and Theorems in Analysis I (Springer 1924, reprinted 1972), Part One, Chap. 1, Sect. 1, Problem 25.
  • 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

Programs

  • Haskell
    a001399 = p [1,2,3] where
       p _      0 = 1
       p []     _ = 0
       p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Feb 28 2013
    
  • Magma
    I:=[1,1,2,3,4,5]; [n le 6 select I[n] else Self(n-1)+Self(n-2)-Self(n-4)-Self(n-5)+Self(n-6): n in [1..80]]; // Vincenzo Librandi, Feb 14 2015
    
  • Magma
    [#RestrictedPartitions(n,{1,2,3}): n in [0..62]]; // Marius A. Burtea, Jan 06 2019
    
  • Magma
    [Round((n+3)^2/12): n in [0..70]]; // Marius A. Burtea, Jan 06 2019
    
  • Maple
    A001399 := proc(n)
        round( (n+3)^2/12) ;
    end proc:
    seq(A001399(n),n=0..40) ;
    with(combstruct):ZL4:=[S,{S=Set(Cycle(Z,card<4))}, unlabeled]:seq(count(ZL4,size=n),n=0..61); # Zerinvary Lajos, Sep 24 2007
    B:=[S,{S = Set(Sequence(Z,1 <= card),card <=3)},unlabelled]: seq(combstruct[count](B, size=n), n=0..61); # Zerinvary Lajos, Mar 21 2009
  • Mathematica
    CoefficientList[ Series[ 1/((1 - x)*(1 - x^2)*(1 - x^3)), {x, 0, 65} ], x ]
    Table[ Length[ IntegerPartitions[n, 3]], {n, 0, 61} ] (* corrected by Jean-François Alcover, Aug 08 2012 *)
    k = 3; Table[(Apply[Plus, Map[EulerPhi[ # ]Binomial[n/#, k/# ] &, Divisors[GCD[n, k]]]]/n + Binomial[If[OddQ[n], n - 1, n - If[OddQ[k], 2, 0]]/2, If[OddQ[k], k - 1, k]/2])/2, {n, k, 50}] (* Robert A. Russell, Sep 27 2004 *)
    LinearRecurrence[{1,1,0,-1,-1,1},{1,1,2,3,4,5},70] (* Harvey P. Dale, Jun 21 2012 *)
    a[ n_] := With[{m = Abs[n + 3] - 3}, Length[ IntegerPartitions[ m, 3]]]; (* Michael Somos, Dec 25 2014 *)
    k=3 (* Number of red beads in bracelet problem *);CoefficientList[Series[(1/k Plus@@(EulerPhi[#] (1-x^#)^(-(k/#))&/@Divisors[k])+(1+x)/(1-x^2)^Floor[(k+2)/2])/2,{x,0,50}],x] (* Herbert Kociemba, Nov 04 2016 *)
    Table[Length[Select[IntegerPartitions[n,{3}],UnsameQ@@#&]],{n,0,30}] (* Gus Wiseman, Apr 15 2019 *)
  • PARI
    {a(n) = round((n + 3)^2 / 12)}; /* Michael Somos, Sep 04 2006 */
    
  • Python
    [round((n+3)**2 / 12) for n in range(0,62)] # Ya-Ping Lu, Jan 24 2024

Formula

G.f.: 1/((1 - x) * (1 - x^2) * (1 - x^3)) = -1/((x+1)*(x^2+x+1)*(x-1)^3); Simon Plouffe in his 1992 dissertation
a(n) = round((n + 3)^2/12). Note that this cannot be of the form (2*i + 1)/2, so ties never arise.
a(n) = A008284(n+3, 3), n >= 0.
a(n) = 1 + a(n-2) + a(n-3) - a(n-5) for all n in Z. - Michael Somos, Sep 04 2006
a(n) = a(-6 - n) for all n in Z. - Michael Somos, Sep 04 2006
a(6*n) = A003215(n), a(6*n + 1) = A000567(n + 1), a(6*n + 2) = A049450(n + 1), a(6*n + 3) = A033428(n + 1), a(6*n + 4) = A049451(n + 1), a(6*n + 5) = A045944(n + 1).
a(n) = a(n-1) + A008615(n+2) = a(n-2) + A008620(n) = a(n-3) + A008619(n) = A001840(n+1) - a(n-1) = A002620(n+2) - A001840(n) = A000601(n) - A000601(n-1). - Henry Bottomley, Apr 17 2001
P(n, 3) = (1/72) * (6*n^2 - 7 - 9*pcr{1, -1}(2, n) + 8*pcr{2, -1, -1}(3, n)) (see Comtet). [Here "pcr" stands for "prime circulator" and it is defined on p. 109 of Comtet, while the formula appears on p. 110. - Petros Hadjicostas, Oct 03 2019]
Let m > 0 and -3 <= p <= 2 be defined by n = 6*m+p-3; then for n > -3, a(n) = 3*m^2 + p*m, and for n = -3, a(n) = 3*m^2 + p*m + 1. - Floor van Lamoen, Jul 23 2001
72*a(n) = 17 + 6*(n+1)*(n+5) + 9*(-1)^n - 8*A061347(n). - Benoit Cloitre, Feb 09 2003
From Jon Perry, Jun 17 2003: (Start)
a(n) = 6*t(floor(n/6)) + (n%6) * (floor(n/6) + 1) + (n mod 6 == 0?1:0), where t(n) = n*(n+1)/2.
a(n) = ceiling(1/12*n^2 + 1/2*n) + (n mod 6 == 0?1:0).
[Here "n%6" means "n mod 6" while "(n mod 6 == 0?1:0)" means "if n mod 6 == 0 then 1, else 0" (as in C).]
(End)
a(n) = Sum_{i=0..floor(n/3)} 1 + floor((n - 3*i)/2). - Jon Perry, Jun 27 2003
a(n) = Sum_{k=0..n} floor((k + 2)/2) * (cos(2*Pi*(n - k)/3 + Pi/3)/3 + sqrt(3) * sin(2*Pi*(n-k)/3 + Pi/3)/3 + 1/3). - Paul Barry, Apr 16 2005
(m choose 3)_q = (q^m-1) * (q^(m-1) - 1) * (q^(m-2) - 1)/((q^3 - 1) * (q^2 - 1) * (q - 1)).
a(n) = Sum_{k=0..floor(n/2)} floor((3 + n - 2*k)/3). - Paul Barry, Nov 11 2003
A117220(n) = a(A003586(n)). - Reinhard Zumkeller, Mar 04 2006
a(n) = 3 * Sum_{i=2..n+1} floor(i/2) - floor(i/3). - Thomas Wieder, Feb 11 2007
Identical to the number of points inside or on the boundary of the integer grid of {I, J}, bounded by the three straight lines I = 0, I - J = 0 and I + 2J = n. - Jonathan Vos Post, Jul 03 2007
a(n) = A026820(n,3) for n > 2. - Reinhard Zumkeller, Jan 21 2010
Euler transform of length 3 sequence [ 1, 1, 1]. - Michael Somos, Feb 25 2012
a(n) = A005044(2*n + 3) = A005044(2*n + 6). - Michael Somos, Feb 25 2012
a(n) = A000212(n+3) - A002620(n+3). - Richard R. Forberg, Dec 08 2013
a(n) = a(n-1) + a(n-2) - a(n-4) - a(n-5) + a(n-6). - David Neil McGrath, Feb 14 2015
a(n) = floor((n^2+3)/12) + floor((n+2)/2). - Giacomo Guglieri, Apr 02 2019
From Devansh Singh, May 28 2020: (Start)
Let p(n, 3) be the number of 3-part integer partitions in which every part is > 0.
Then for n >= 3, p(n, 3) is equal to:
(n^2 - 1)/12 when n is odd and 3 does not divide n.
(n^2 + 3)/12 when n is odd and 3 divides n.
(n^2 - 4)/12 when n is even and 3 does not divide n.
(n^2)/12 when n is even and 3 divides n.
For n >= 3, p(n, 3) = a(n-3). (End)
a(n) = floor(((n+3)^2 + 4)/12). - Vladimír Modrák, Zuzana Soltysova, Dec 08 2020
Sum_{n>=0} 1/a(n) = 15/4 - Pi/(2*sqrt(3)) + Pi^2/18 + tanh(Pi/(2*sqrt(3)))*Pi/sqrt(3). - Amiram Eldar, Sep 29 2022
E.g.f.: exp(-x)*(9 + exp(2*x)*(47 + 42*x + 6*x^2) + 16*exp(x/2)*cos(sqrt(3)*x/2))/72. - Stefano Spezia, Mar 05 2023
a(6n) = 1+6*A000217(n); Sum_{i=1..n} a(6*i) = A000578(n+1). - David García Herrero, May 05 2024

Extensions

Name edited by Gus Wiseman, Apr 15 2019

A000389 Binomial coefficients C(n,5).

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 6, 21, 56, 126, 252, 462, 792, 1287, 2002, 3003, 4368, 6188, 8568, 11628, 15504, 20349, 26334, 33649, 42504, 53130, 65780, 80730, 98280, 118755, 142506, 169911, 201376, 237336, 278256, 324632, 376992, 435897, 501942, 575757, 658008, 749398
Offset: 0

Views

Author

Keywords

Comments

a(n+4) is the number of inequivalent ways of coloring the vertices of a regular 4-dimensional simplex with n colors, under the full symmetric group S_5 of order 120, with cycle index (x1^5 + 10*x1^3*x2 + 20*x1^2*x3 + 15*x1*x2^2 + 30*x1*x4 + 20*x2*x3 + 24*x5)/120.
Figurate numbers based on 5-dimensional regular simplex. According to Hyun Kwang Kim, it appears that every nonnegative integer can be represented as the sum of g = 10 of these 5-simplex(n) numbers (compared with g=3 for triangular numbers, g=5 for tetrahedral numbers and g=8 for pentatope numbers). - Jonathan Vos Post, Nov 28 2004
The convolution of the nonnegative integers (A001477) with the tetrahedral numbers (A000292), which are the convolution of the nonnegative integers with themselves (making appropriate allowances for offsets of all sequences). - Graeme McRae, Jun 07 2006
a(n) is the number of terms in the expansion of (a_1 + a_2 + a_3 + a_4 + a_5 + a_6)^n. - Sergio Falcon, Feb 12 2007
Product of five consecutive numbers divided by 120. - Artur Jasinski, Dec 02 2007
Equals binomial transform of [1, 5, 10, 10, 5, 1, 0, 0, 0, ...]. - Gary W. Adamson, Feb 02 2009
Equals INVERTi transform of A099242 (1, 7, 34, 153, 686, 3088, ...). - Gary W. Adamson, Feb 02 2009
For a team with n basketball players (n>=5), this sequence is the number of possible starting lineups of 5 players, without regard to the positions (center, forward, guard) of the players. - Mohammad K. Azarian, Sep 10 2009
a(n) is the number of different patterns, regardless of order, when throwing (n-5) 6-sided dice. For example, one die can display the 6 numbers 1, 2, ..., 6; two dice can display the 21 digit-pairs 11, 12, ..., 56, 66. - Ian Duff, Nov 16 2009
Sum of the first n pentatope numbers (1, 5, 15, 35, 70, 126, 210, ...), see A000332. - Paul Muljadi, Dec 16 2009
Sum_{n>=0} a(n)/n! = e/120. Sum_{n>=4} a(n)/(n-4)! = 501*e/120. See A067764 regarding the second ratio. - Richard R. Forberg, Dec 26 2013
For a set of integers {1,2,...,n}, a(n) is the sum of the 2 smallest elements of each subset with 4 elements, which is 3*C(n+1,5) (for n>=4), hence a(n) = 3*C(n+1,5) = 3*A000389(n+1). - Serhat Bulut, Mar 11 2015
a(n) = fallfac(n,5)/5! is also the number of independent components of an antisymmetric tensor of rank 5 and dimension n >= 1. Here fallfac is the falling factorial. - Wolfdieter Lang, Dec 10 2015
Number of compositions (ordered partitions) of n+1 into exactly 6 parts. - Juergen Will, Jan 02 2016
Number of weak compositions (ordered weak partitions) of n-5 into exactly 6 parts. - Juergen Will, Jan 02 2016
a(n+3) could be the general number of all geodetic graphs of diameter n>=2 homeomorphic to the Petersen Graph. - Carlos Enrique Frasser, May 24 2018
From Robert A. Russell, Dec 24 2020: (Start)
a(n) is the number of chiral pairs of colorings of the 5 tetrahedral facets (or vertices) of the regular 4-D simplex (5-cell, pentachoron, Schläfli symbol {3,3,3}) using subsets of a set of n colors. Each member of a chiral pair is a reflection but not a rotation of the other.
a(n+4) is the number of unoriented colorings of the 5 tetrahedral facets of the regular 4-D simplex (5-cell, pentachoron) using subsets of a set of n colors. Each chiral pair is counted as one when enumerating unoriented arrangements. (End)
For integer m and positive integer r >= 4, the polynomial a(n) + a(n + m) + a(n + 2*m) + ... + a(n + r*m) in n has its zeros on the vertical line Re(n) = (4 - r*m)/2 in the complex plane. - Peter Bala, Jun 02 2024

Examples

			G.f. = x^5 + 6*x^6 + 21*x^7 + 56*x^8 + 126*x^9 + 252*x^10 + 462*x^11 + ...
For A={1,2,3,4}, the only subset with 4 elements is {1,2,3,4}; sum of 2 minimum elements of this subset: a(4) = 1+2 = 3 = 3*C(4+1,5).
For A={1,2,3,4,5}, the subsets with 4 elements are {1,2,3,4}, {1,2,3,5}, {1,2,4,5}, {1,3,4,5}, {2,3,4,5}; sum of 2 smallest elements of each subset: a(5) = (1+2)+(1+2)+(1+2)+(1+3)+(2+3) = 18 = 3*C(5+1,5). - _Serhat Bulut_, Mar 11 2015
a(6) = 6 from the six independent components of an antisymmetric tensor A of rank 5 and dimension 6: A(1,2,3,4,5), A(1,2,3,4,6), A(1,2,3,5,6), A(1,2,4,5,6), A(1,3,4,5,6), A(2,3,4,5,6). See the Dec 10 2015 comment. - _Wolfdieter Lang_, Dec 10 2015
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 828.
  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 196.
  • 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. 7.
  • Gupta, Hansraj; Partitions of j-partite numbers into twelve or a smaller number of parts. Collection of articles dedicated to Professor P. L. Bhatnagar on his sixtieth birthday. Math. Student 40 (1972), 401-441 (1974).
  • J. C. P. Miller, editor, Table of Binomial Coefficients. Royal Society Mathematical Tables, Vol. 3, Cambridge Univ. Press, 1954.
  • 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

Cf. A099242. - Gary W. Adamson, Feb 02 2009
Cf. A242023. A104712 (fourth column, k=5).
5-cell colorings: A337895 (oriented), A132366(n-1) (achiral).
Unoriented colorings: A063843 (5-cell edges, faces), A128767 (8-cell vertices, 16-cell facets), A337957 (16-cell vertices, 8-cell facets), A338949 (24-cell), A338965 (600-cell vertices, 120-cell facets).
Chiral colorings: A331352 (5-cell edges, faces), A337954 (8-cell vertices, 16-cell facets), A234249 (16-cell vertices, 8-cell facets), A338950 (24-cell), A338966 (600-cell vertices, 120-cell facets).

Programs

  • Haskell
    a000389 n = a000389_list !! n
    a000389_list = 0 : 0 : f [] a000217_list where
       f xs (t:ts) = (sum $ zipWith (*) xs a000217_list) : f (t:xs) ts
    -- Reinhard Zumkeller, Mar 03 2015, Apr 13 2012
    
  • Magma
    [Binomial(n, 5): n in [0..40]]; // Vincenzo Librandi, Mar 12 2015
  • Maple
    f:=n->(1/120)*(n^5-10*n^4+35*n^3-50*n^2+24*n): seq(f(n), n=0..60);
    ZL := [S, {S=Prod(B,B,B,B,B,B), B=Set(Z, 1 <= card)}, unlabeled]: seq(combstruct[count](ZL, size=n+1), n=0..42); # Zerinvary Lajos, Mar 13 2007
    A000389:=1/(z-1)**6; # Simon Plouffe, 1992 dissertation
  • Mathematica
    Table[Binomial[n, 5], {n, 5, 50}] (* Stefan Steinerberger, Apr 02 2006 *)
    CoefficientList[Series[x^5 / (1 - x)^6, {x, 0, 40}], x] (* Vincenzo Librandi, Mar 12 2015 *)
    LinearRecurrence[{6,-15,20,-15,6,-1},{0,0,0,0,0,1},50] (* Harvey P. Dale, Jul 17 2016 *)
  • PARI
    (conv(u,v)=local(w); w=vector(length(u),i,sum(j=1,i,u[j]*v[i+1-j])); w);
    (t(n)=n*(n+1)/2); u=vector(10,i,t(i)); conv(u,u)
    

Formula

G.f.: x^5/(1-x)^6.
a(n) = n*(n-1)*(n-2)*(n-3)*(n-4)/120.
a(n) = (n^5-10*n^4+35*n^3-50*n^2+24*n)/120. (Replace all x_i's in the cycle index with n.)
a(n+2) = Sum_{i+j+k=n} i*j*k. - Benoit Cloitre, Nov 01 2002
Convolution of triangular numbers (A000217) with themselves.
Partial sums of A000332. - Alexander Adamchuk, Dec 19 2004
a(n) = -A110555(n+1,5). - Reinhard Zumkeller, Jul 27 2005
a(n+3) = (1/2!)*(d^2/dx^2)S(n,x)|A049310.%20-%20_Wolfdieter%20Lang">{x=2}, n>=2, one half of second derivative of Chebyshev S-polynomials evaluated at x=2. See A049310. - _Wolfdieter Lang, Apr 04 2007
a(n) = A052787(n+5)/120. - Zerinvary Lajos, Apr 26 2007
Sum_{n>=5} 1/a(n) = 5/4. - R. J. Mathar, Jan 27 2009
For n>4, a(n) = 1/(Integral_{x=0..Pi/2} 10*(sin(x))^(2*n-9)*(cos(x))^9). - Francesco Daddi, Aug 02 2011
Sum_{n>=5} (-1)^(n + 1)/a(n) = 80*log(2) - 655/12 = 0.8684411114... - Richard R. Forberg, Aug 11 2014
a(n) = -a(4-n) for all n in Z. - Michael Somos, Oct 07 2014
0 = a(n)*(+a(n+1) + 4*a(n+2)) + a(n+1)*(-6*a(n+1) + a(n+2)) for all n in Z. - Michael Somos, Oct 07 2014
a(n) = 3*C(n+1, 5) = 3*A000389(n+1). - Serhat Bulut, Mar 11 2015
From Ilya Gutkovskiy, Jul 23 2016: (Start)
E.g.f.: x^5*exp(x)/120.
Inverse binomial transform of A054849. (End)
From Robert A. Russell, Dec 24 2020: (Start)
a(n) = A337895(n) - a(n+4) = (A337895(n) - A132366(n-1)) / 2 = a(n+4) - A132366(n-1).
a(n+4) = A337895(n) - a(n) = (A337895(n) + A132366(n-1)) / 2 = a(n) + A132366(n-1).
a(n+4) = 1*C(n,1) + 4*C(n,2) + 6*C(n,3) + 4*C(n,4) + 1*C(n,5), where the coefficient of C(n,k) is the number of unoriented pentachoron colorings using exactly k colors. (End)

Extensions

Corrected formulas that had been based on other offsets. - R. J. Mathar, Jun 16 2009
I changed the offset to 0. This will require some further adjustments to the formulas. - N. J. A. Sloane, Aug 01 2010

A002623 Expansion of 1/((1-x)^4*(1+x)).

Original entry on oeis.org

1, 3, 7, 13, 22, 34, 50, 70, 95, 125, 161, 203, 252, 308, 372, 444, 525, 615, 715, 825, 946, 1078, 1222, 1378, 1547, 1729, 1925, 2135, 2360, 2600, 2856, 3128, 3417, 3723, 4047, 4389, 4750, 5130, 5530, 5950, 6391, 6853, 7337, 7843, 8372, 8924, 9500
Offset: 0

Views

Author

Keywords

Comments

Also a(n) is the number of nondegenerate triangles that can be made from rods of lengths 1 to n+1. - Alfred Bruckstein; corrected by Hans Rudolf Widmer, Nov 02 2023
Also number of circumscribable (or escrible) quadrilaterals that can be made from rods of length 1,2,3,4,...,n. - Antreas P. Hatzipolakis (xpolakis(AT)otenet.gr)
Also number of 2 X n binary matrices up to row and column permutation (see the link: Binary matrices up to row and column permutations). - Vladeta Jovovic
Also partial sum of alternate triangular numbers (1, 3, 1+6, 3+10, 1+6+15, 3+10+21, etc.); and also number of triangles pointing in opposite direction to largest triangle in triangular matchstick arrangement of side n+2 (cf. A002717, also the Larsen article). - Henry Bottomley, Aug 08 2000
Ordered union of A002412(n+1) and A016061(n+1). - Lekraj Beedassy, Oct 13 2003
Also Molien series for certain 4-D representation of cyclic group of order 2. - N. J. A. Sloane, Jun 12 2004
From Radu Grigore (radugrigore(AT)gmail.com), Jun 19 2004: (Start)
a(n) = floor( (n+2)*(n+4)*(2n+3) / 24 ). E.g., a(2) = floor(4*6*7/24) = 7 because there are 7 upside down triangles (6 of size one and 1 of size two) in the matchstick figure:
/\
/\/\
/\/\/\
/\/\/\/\
(End)
Number of non-congruent non-parallelogram trapezoids with positive integer sides (trapezints) and perimeter 2n+5. Also with perimeter 2n+8. - Michael Somos, May 12 2005
a(n) = A108561(n+4,n) for n > 0. - Reinhard Zumkeller, Jun 10 2005
Also number of nonisomorphic planes with n points and 2 lines. E.g., a(0)=1 because with no points, we just have two empty lines. a(1)=3 because the one point may belong to 0, 1 or 2 lines. a(2)=7 because there are 7 ways to determine which of 2 points belong to which of 2 lines, up to isomorphism, i.e., up to a bijection f on the sets of points and a bijection g on the sets of lines, such that A belongs to a iff f(A) belongs to g(a). - Bjorn Kjos-Hanssen (bjorn(AT)math.uconn.edu), Nov 10 2005
a(n-2) is the number of ways to pick two non-overlapping subwords of equal nonzero length from a word of length n. E.g., a(5-2)=a(3)=13 since the word 12345 of length 5 has the following subword pairs: 1,2; 1,3; 1,4; 1,5; 2,3; 2,4; 2,5; 3,4; 3,5; 4,5; 12,34; 12,45; 23,45. - Michael Somos, Oct 22 2006
Partial sums of A002620. - G.H.J. van Rees (vanrees(AT)cs.umanitoba.ca), Feb 16 2007
From Philippe LALLOUET (philip.lallouet(AT)orange.fr), Oct 19 2007: (Start)
Also number of squares of any size in a staircase of n steps built with unit squares:
||__
||__|
||__||
For a staircase of 3 steps 6 squares of size 1 and 1 square of size 2, hence c(3)=7.
Columns sums of:
1 3 6 10 15 21 28 ...
1 3 6 10 15 ...
1 3 6 ...
1 ...
---------------------
1 3 7 13 22 34 50 ...
(End)
a(n) = sum of row n+1 of triangle A134446. Also, binomial transform of [1, 2, 2, 0, 1, -2, 4, -8, 16, -32, ...]. - Gary W. Adamson, Oct 25 2007
Let b(n) be the number of 4-tuples (w,x,y,z) having all terms in {1,...,n} and 2w=x+y+z+n; then b(n+3) = a(n) for n >= 0. - Clark Kimberling, May 08 2012
a(n) is the number of 3-tuples (w,x,y) having all terms in {0,...,n} and w >= x+y and x <= y. - Clark Kimberling, Jun 04 2012
Also, number of unlabeled bipartite graphs with two left vertices and n right vertices. - Yavuz Oruc, Jan 14 2018
Also number of triples (x,y,z) with 0 < x <= y <= z <= n + 1, x + y > z. - Ralf Steiner, Feb 06 2020
Bisections A002412 and A016061: a(2*k) = k*(k+1)*(4*k-1)/3! and a(2*k+1) = (k+1)*(k+2)*(4*k+9)/3!, for k >= 0. See the Woolhouse link, II. Solution by Stephen Watson, p. 65, with index shifts. - Mo Li, Apr 02 2020
Also, Wiener index of the square of the path graph P_(n+2). - Allan Bickle, Aug 01 2020
Maximum Wiener index of all maximal 2-degenerate graphs with n+2 vertices. (A maximal 2-degenerate graph can be constructed from a 2-clique by iteratively adding a new 2-leaf (vertex of degree 2) adjacent to two existing vertices.) The extremal graphs are squares of paths, so the bound also applies to 2-trees and maximal outerplanar graphs. - Allan Bickle, Sep 15 2022

Examples

			G.f. = 1 + 3*x + 7*x^2 + 13*x^3 + 22*x^4 + 34*x^5 + 50*x^6 + 70*x^7 + 95*x^8 + ...
		

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 74, Problem 7.
  • P. Diaconis, R. L. Graham and B. Sturmfels, Primitive partition identities, in Combinatorics: Paul Erdős is Eighty, Vol. 2, Bolyai Soc. Math. Stud., 2, 1996, pp. 173-192.
  • H. Gupta, Partitions of j-partite numbers into twelve or a smaller number of parts. Collection of articles dedicated to Professor P. L. Bhatnagar on his sixtieth birthday. Math. Student 40 (1972), 401-441 (1974).
  • I. Siap, Linear codes over F_2 + u*F_2 and their complete weight enumerators, in Codes and Designs (Ohio State, May 18, 2000), pp. 259-271. De Gruyter, 2002.
  • 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

Cf. A002620 (first differences), A000292, A001752 (partial sums), A062109 (binomial transf.).
Bisections A002412, A016061.
Cf. also A002717 (a companion sequence), A002727, A006148, A057524, A134446, A014125, A122046, A122047.
The maximum Wiener index of all maximal k-degenerate graphs for k=1..6 are given in A000292, A002623 (this sequence), A014125, A122046, A122047, A175724, respectively.

Programs

  • Maple
    A002623 := n->(1/16)*(1+(-1)^n)+(n+1)/8+binomial(n+2,2)/4+binomial(n+3,3)/2;
    seq( ((2*n+3)*(n+2)*(n+1)/6-floor((n+2)/2))/4,n=1..47); # Lewis
    a := n -> ((-1)^n*3 + 45 + 68*n + 30*n^2 + 4*n^3) / 48:
    seq(a(n), n=0..46); # Peter Luschny, Jan 22 2018
  • Mathematica
    CoefficientList[Series[1/((1-x)^3(1-x^2)),{x,0,50}],x] (* or *) LinearRecurrence[{3,-2,-2,3,-1},{1,3,7,13,22},50] (* Harvey P. Dale, Jul 19 2011 *)
    Table[((2 n^3 + 15 n^2 + 34 n + 45 / 2 + (3/2) (-1)^n) / 24), {n, 0, 100}] (* Vincenzo Librandi, Jan 15 2018 *)
    a[ n_] := Floor[(n + 2)*(n + 4)*(2*n + 3)/24]; (* Michael Somos, Feb 19 2024 *)
  • PARI
    {a(n) = (8 + 34/3*n + 5*n^2 + 2/3*n^3) \ 8}; /* Michael Somos, Sep 04 1999 */
    
  • PARI
    x='x+O('x^50); Vec(1/((1 - x)^3 * (1 - x^2))) \\ Indranil Ghosh, Apr 04 2017
    
  • Python
    def A002623(n): return ((n+2)*(n+4)*((n<<1)+3)>>3)//3 # Chai Wah Wu, Mar 25 2024

Formula

a(n+1) = a(n) + {(k-1)*k if n=2*k} or {k*k if n=2*k+1}.
a(n)+a(n+1) = A000292(n+1).
a(n) = a(n-2) + A000217(n+1) = A002717(n+2) - A000292(n+1).
Also: a(n) = C(n+3, 3) - a(n-1) with a(0)=1. - Labos Elemer, Apr 26 2003
From Paul Barry, Jul 01 2003: (Start)
a(n) = Sum_{k=0..n} (-1)^(n-k)*C(k+3,3).
The signed version 1, -3, 7, ... has the formula:
a(n) = (4*n^3 + 30*n^2 + 68*n + 45)*(-1)^n/48 + 1/16.
This is the partial sums of the signed version of A000292. (End)
From Paul Barry, Jul 21 2003: (Start)
a(n) = Sum_{k=0..n} floor((k+2)^2/4).
a(n) = Sum_{k=0..n} Sum_{j=0..k} Sum_{i=0..j} (1+(-1)^i)/2. (End)
a(n) = a(n - 2) + (n*(n - 1))/2, with n>2, a(1)=0, a(2)=1; a(n) = (4*n^3+6*n^2-4*n+3*(-1)^n-3)/48, with offset 2. - Cecilia Rossiter (cecilia(AT)noticingnumbers.net), Dec 14 2004 (formula simplified by Bruno Berselli, Aug 29 2013)
a(n) = ((2*n+3)*(n+2)*(n+1)/6-floor((n+2)/2))/4, with offset 1. - Jerry W. Lewis (JLewis(AT)wyeth.com), Mar 23 2005
a(n) = 2*a(n-1) - a(n-2) + 1 + floor(n/2). - Bjorn Kjos-Hanssen (bjorn(AT)math.uconn.edu), Nov 10 2005
A002620(n+3) = a(n+1) - a(n). - Michael Somos, Sep 04 1999
Euler transform of length 2 sequence [ 3, 1]. - Michael Somos, Sep 04 2006
a(n) = -a(-5-n) for all n in Z. - Michael Somos, Sep 04 2006
Let P(i,k) be the number of integer partitions of n into k parts, then with k=2 we have a(n) = sum_{m=1}^{n} sum_{i=k}^{m} P(i,k). For k=1 we get A000217 = triangular numbers. - Thomas Wieder, Feb 18 2007
a(n) = (n+(3+(-1)^n)/2)*(n+(7+(-1)^n)/2)*(2*n+5-2*(-1)^n)/24. - Philippe LALLOUET (philip.lallouet(AT)orange.fr), Oct 19 2007 (corrected by Bruno Berselli, Aug 30 2013)
From Johannes W. Meijer, May 20 2011: (Start)
a(n) = A006918(n+1) + A006918(n).
a(n) = A058187(n-2) + 2*A058187(n-1) + A058187(n). (End)
a(0)=1, a(1)=3, a(2)=7, a(3)=13, a(4)=22; for n > 4, a(n) = 3*a(n-1) - 2*a(n-2) - 2*a(n-3) + 3*a(n-4) - a(n-5). - Harvey P. Dale, Jul 19 2011
a(n) = Sum_{i=0..n+2} floor(i/2)*ceiling(i/2). - Bruno Berselli, Aug 30 2013
a(n) = 15/16 + (1/16)*(-1)^n + (17/12)*n + (5/8)*n^2 + (1/12)*n^3. - Robert Israel, Jul 07 2014
a(n) = Sum_{i=0..n+2} (n+1-i)*floor(i/2+1). - Bruno Berselli, Apr 04 2017
a(n) = 1 + floor((2*n^3 + 15*n^2 + 34*n) / 24). - Allan Bickle, Aug 01 2020
E.g.f.: ((24 + 51*x + 21*x^2 + 2*x^3)*cosh(x) + (21 + 51*x + 21*x^2 + 2*x^3)*sinh(x))/24. - Stefano Spezia, Jun 02 2021
Showing 1-10 of 69 results. Next