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

A014125 Bisection of A001400.

Original entry on oeis.org

1, 3, 6, 11, 18, 27, 39, 54, 72, 94, 120, 150, 185, 225, 270, 321, 378, 441, 511, 588, 672, 764, 864, 972, 1089, 1215, 1350, 1495, 1650, 1815, 1991, 2178, 2376, 2586, 2808, 3042, 3289, 3549, 3822, 4109, 4410, 4725, 5055, 5400, 5760, 6136, 6528, 6936, 7361, 7803
Offset: 0

Views

Author

Keywords

Comments

Also Schoenheim bound L_1(n,5,4).
Degrees of polynomials defined by p(n) = (x^(n+1)*p(n-1)p(n-3) + p(n-2)^2)/p(n-4), p(-4)=p(-3)=p(-2)=p(-1)=1. - Michael Somos, Jul 21 2004
Degrees of polynomial tau-functions of q-discrete Painlevé I, which generate sequence A095708 when q=2 (up to an offset of 3). - Andrew Hone, Jul 29 2004
Because of the Laurent phenomenon for the general q-discrete Painlevé I tau-function recurrence p(n) = (a*x^(n+1)*p(n-1)*p(n-3) + b*p(n-2)^2)/p(n-4), p(n) for n > -1 will always be a polynomial in x and a Laurent polynomial in a,b and the initial data p(-4),p(-3),p(-2),p(-1). - Andrew Hone, Jul 29 2004
Create the sequence 0,0,0,0,0,6,18,36,66,108,... so that the sum of three consecutive terms b(n) + b(n+1) + b(n+2) = A007531(n), with b(0)=0; then a(n) = b(n+5)/6. - J. M. Bergot, Jul 30 2013
Number of partitions of n into three kinds of part 1 and one kind of part 3. - Joerg Arndt, Sep 28 2015
First differences are A001840(k) starting with k=2; second differences are A086161(k) starting with k=1. - Bob Selcoe, Sep 28 2015
Maximum Wiener index of all maximal planar graphs with n+2 vertices. The extremal graphs are cubes of paths. - Allan Bickle, Jul 09 2022
Maximum Wiener index of all maximal 3-degenerate graphs with n+2 vertices. (A maximal 3-degenerate graph can be constructed from a 3-clique by iteratively adding a new 3-leaf (vertex of degree 3) adjacent to three existing vertices.) The extremal graphs are cubes of paths, so the bound also applies to 3-trees. - Allan Bickle, Sep 18 2022

Examples

			Polynomials: p(0)=x+1, p(1)=x^3+x^2+1, p(2)=x^6+x^5+x^3+x^2+2x+1, ...
a(12)=185:  A000217(13)=91 + a(9)=94 == 91+55+28+10+1 = 185. - _Bob Selcoe_, Sep 27 2015
a(3)=11: the 11 partitions of 3 are {1a,1a,1a}, {1a,1a,1b}, {1a,1a,1c}, {1a,1b,1b}, {1a,1b,1c}, {1a,1c,1c}, {1b,1b,1b}, {1b,1b,1c}, {1b,1c,1c}, {1c,1c,1c}, {3}. - _Bob Selcoe_, Oct 04 2015
		

References

  • W. H. Mills and R. C. Mullin, Coverings and packings, pp. 371-399 of Jeffrey H. Dinitz and D. R. Stinson, editors, Contemporary Design Theory, Wiley, 1992. See Eq. 1.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • L. Smiley, Hidden Hexagons (preprint).

Crossrefs

A column of A036838.
Maximum Wiener index of all maximal k-degenerate graphs for k=1..6: A000292, A002623, A014125, A122046, A122047, A175724.

Programs

  • Magma
    [n^3/18+n^2/2+4*n/3+1+(((n+1) mod 3)-1)/9 : n in [0..50]]; // Wesley Ivan Hurt, Apr 14 2015
    
  • Magma
    I:=[1,3,6,11,18,27]; [n le 6 select I[n] else 3*Self(n-1) -3*Self(n-2) +2*Self(n-3)-3*Self(n-4)+3*Self(n-5)-Self(n-6): n in [1..50]]; // Vincenzo Librandi, Apr 15 2015
    
  • Maple
    L := proc(v,k,t,l) local i,t1; t1 := l; for i from v-t+1 to v do t1 := ceil(t1*i/(i-(v-k))); od: t1; end; # gives Schoenheim bound L_l(v,k,t). Current sequence is L_1(n,n-3,n-4,1).
  • Mathematica
    CoefficientList[Series[1/((1 - x)^3*(1 - x^3)), {x, 0, 50}], x] (* Wesley Ivan Hurt, Apr 14 2015 *)
  • PARI
    a(n)=if(n<-5,-a(-6-n),polcoeff(1/(1-x)^3/(1-x^3)+x^n*O(x),n)) /* Michael Somos, Jul 21 2004 */
    
  • PARI
    my(x='x+O('x^50)); Vec(1/((1-x)^3*(1-x^3))) \\ Altug Alkan, Oct 16 2015
    
  • PARI
    a(n)=(n^3 + 9*n^2 + 24*n + 19)\/18 \\ Charles R Greathouse IV, Jun 29 2020
    
  • Sage
    [(binomial(n+4,3) - ((n+4)//3))/3 for n in (0..50)] # G. C. Greubel, Apr 28 2019

Formula

G.f.: 1/((1-x)^3*(1-x^3)).
a(n) = -a(-6-n) = 3*a(n-1) -3*a(n-2) +2*a(n-3) -3*a(n-4) +3*a(n-5) -a(n-6).
The simplest recurrence is fourth order: a(n) = a(n-1) + a(n-3) - a(n-4) + n + 1, which gives the g.f.: 1/((1-x)^3*(1-x^3)), with a(-4) = a(-3) = a(-2) = a(-1) = 0.
a(n) = n^3/18 + n^2/2 + 4*n/3 + 1 + 2/(9*sqrt(3))*sin(2*Pi*n/3). - Andrew Hone, Jul 29 2004
a(n) = n^3/18 + n^2/2 + 4*n/3 + 1 + (((n+1) mod 3) - 1)/9. - same formula, simplified by Gerald Hillier, Apr 14 2015
a(n) = (2*A000027(n+1) + 3*A000292(n+1) + A049347(n-1) + 1 + 3*A000217(n+1))/9. - R. J. Mathar, Nov 16 2007
From Johannes W. Meijer, May 20 2011: (Start)
a(n) = A144677(n) + A144677(n-1) + A144677(n-2).
a(n) = A190717(n-4) + 2*A190717(n-3) + 3*A190717(n-2) + 2*A190717(n-1) + A190717(n). (End)
3*a(n) = binomial(n+4,3) - floor((n+4)/3). - Bruno Berselli, Nov 08 2013
a(n) = A000217(n+1) + a(n-3) = Sum_{j>=0, n>=3*j} (n-3*j+1)*(n-3*j+2)/2. - Bob Selcoe, Sep 27 2015
a(n) = round(((2*n+5)^3 + 3*(2*n+5)^2 - 9*(2*n+5))/144). - Giacomo Guglieri, Jun 28 2020
a(n) = floor(((n+2)^3 + 3*(n+2)^2)/18). - Allan Bickle, Aug 01 2020
a(n) = Sum_{j=0..n} (n-j+1)*floor((j+3)/3). - G. C. Greubel, Oct 18 2021
E.g.f.: exp(x) + exp(x)*x*(34 + 12*x + x^2)/18 + 2*exp(-x/2)*sin(sqrt(3)*x/2)/(9*sqrt(3)). - Stefano Spezia, Apr 05 2023

Extensions

More terms from James Sellers, Dec 24 1999

A139672 Convolution of A008619 and A001400.

Original entry on oeis.org

1, 2, 5, 9, 17, 27, 44, 65, 97, 136, 191, 257, 346, 451, 587, 746, 946, 1177, 1461, 1786, 2178, 2623, 3151, 3746, 4443, 5223, 6126, 7131, 8283, 9558, 11007, 12603, 14403, 16377, 18588, 21003, 23692, 26618, 29858, 33372, 37244, 41430, 46022, 50972
Offset: 1

Views

Author

Alford Arnold, Apr 29 2008, May 01 2008

Keywords

Comments

This is row 21 of a table of values related to Molien series. It is the product of the sequence on row 3 (A008619) with the sequence on row 7 (A001400).
This table may be constructed by moving the rows of table A008284 to prime locations and generating the composite locations by multiplication in a manner similar to the calculation illustrated in the present sequence.
Rows 1 thru 20 and 22 thru 25 are as follows:

Programs

  • Maple
    a:= proc(n) local m, r; m:= iquo (n, 12, 'r'); r:= r+1; (19+ (145+ (260+ 15* (r+9)*r+ (405+ 90*r+ 216*m) *m) *m) *m) *m/5+ [0, 1, 2, 5, 9, 17, 27, 44, 65, 97, 136, 191][r]+ [0, 16, 37, 77, 128, 208, 307, 447, 616, 840, 1105, 1441][r]*m/2+ [0, 52, 119, 213, 328, 476, 651, 865, 1112, 1404, 1735, 2117][r]*m^2/2 end: seq (a(n), n=1..50); # Alois P. Heinz, Nov 10 2008
  • Mathematica
    CoefficientList[Series[x/((x^2+x+1)(x^2+1)(x+1)^3 (x-1)^6),{x,0,50}],x] (* or *) LinearRecurrence[{2,1,-3,0,-1,2,2,-1,0,-3,1,2,-1},{0,1,2,5,9,17,27,44,65,97,136,191,257},50] (* Harvey P. Dale, Feb 17 2016 *)

Formula

G.f.: x/((x^2+x+1)*(x^2+1)*(x+1)^3*(x-1)^6). - Alois P. Heinz, Nov 10 2008
a(n)= -A049347(n)/27 +(2*n+11)*(6*n^4+132*n^3+914*n^2+2068*n+1055)/69120 -(-1)^n*(51/512+n^2/256+11*n/256+A057077(n)/32 ). - R. J. Mathar, Nov 21 2008

Extensions

More terms from Alois P. Heinz, Nov 10 2008
Corrected A-number in definition. Added formula. - R. J. Mathar, Nov 21 2008

A008619 Positive integers repeated.

Original entry on oeis.org

1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8, 8, 9, 9, 10, 10, 11, 11, 12, 12, 13, 13, 14, 14, 15, 15, 16, 16, 17, 17, 18, 18, 19, 19, 20, 20, 21, 21, 22, 22, 23, 23, 24, 24, 25, 25, 26, 26, 27, 27, 28, 28, 29, 29, 30, 30, 31, 31, 32, 32, 33, 33, 34, 34, 35, 35, 36, 36, 37, 37, 38
Offset: 0

Views

Author

Keywords

Comments

The floor of the arithmetic mean of the first n+1 positive integers. - Cino Hilliard, Sep 06 2003
Number of partitions of n into powers of 2 where no power is used more than three times, or 4th binary partition function (see A072170).
Number of partitions of n in which the greatest part is at most 2. - Robert G. Wilson v, Jan 11 2002
Number of partitions of n into at most 2 parts. - Jon Perry, Jun 16 2003
a(n) = #{k=0..n: k+n is even}. - Paul Barry, Sep 13 2003
Number of symmetric Dyck paths of semilength n+2 and having two peaks. E.g., a(6)=4 because we have UUUUUUU*DU*DDDDDDD, UUUUUU*DDUU*DDDDDD, UUUUU*DDDUUU*DDDDD and UUUU*DDDDUUUU*DDDD, where U=(1,1), D=(1,-1) and * indicates a peak. - Emeric Deutsch, Jan 12 2004
Smallest positive integer whose harmonic mean with another positive integer is n (for n > 0). For example, a(6)=4 is already given (as 4 is the smallest positive integer such that the harmonic mean of 4 (with 12) is 6) - but the harmonic mean of 2 (with -6) is also 6 and 2 < 4, so the two positive integer restrictions need to be imposed to rule out both 2 and -6.
Second outermost diagonal of Losanitsch's triangle (A034851). - Alonso del Arte, Mar 12 2006
Arithmetic mean of n-th row of A080511. - Amarnath Murthy, Mar 20 2003
a(n) is the number of ways to pay n euros (or dollars) with coins of one and two euros (respectively dollars). - Richard Choulet and Robert G. Wilson v, Dec 31 2007
Inverse binomial transform of A045623. - Philippe Deléham, Dec 30 2008
Coefficient of q^n in the expansion of (m choose 2)_q as m goes to infinity. - Y. Kelly Itakura (yitkr(AT)mta.ca), Aug 21 2002
Binomial transform of (-1)^n*A034008(n) = [1,0,1,-2,4,-8,16,-32,...]. - Philippe Deléham, Nov 15 2009
From Jon Perry_, Nov 16 2010: (Start)
Column sums of:
1 1 1 1 1 1...
1 1 1 1...
1 1...
..............
--------------
1 1 2 2 3 3... (End)
This sequence is also the half-convolution of the powers of 1 sequence A000012 with itself. For the definition of half-convolution see a comment on A201204, where also the rule for the o.g.f. is given. - Wolfdieter Lang, Jan 09 2012
a(n) is also the number of roots of the n-th Bernoulli polynomial in the right half-plane for n>0. - Michel Lagneau, Nov 08 2012
a(n) is the number of symmetry-allowed, linearly-independent terms at n-th order in the series expansion of the Exe vibronic perturbation matrix, H(Q) (cf. Viel & Eisfeld). - Bradley Klee, Jul 21 2015
a(n) is the number of distinct integers in the n-th row of Pascal's triangle. - Melvin Peralta, Feb 03 2016
a(n+1) for n >= 3 is the diameter of the Generalized Petersen Graph G(n, 1). - Nick Mayers, Jun 06 2016
The arithmetic function v_1(n,2) as defined in A289198. - Robert Price, Aug 22 2017
Also, this sequence is the second column in the triangle of the coefficients of the sum of two consecutive Fibonacci polynomials F(n+1, x) and F(n, x) (n>=0) in ascending powers of x. - Mohammad K. Azarian, Jul 18 2018
a(n+2) is the least k such that given any k integers, there exist two of them whose sum or difference is divisible by n. - Pablo Hueso Merino, May 09 2020
Column k = 2 of A051159. - John Keith, Jun 28 2021

References

  • D. J. Benson, Polynomial Invariants of Finite Groups, Cambridge, 1993, p. 100.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 109, Eq. [6c]; p. 116, P(n,2).
  • D. Parisse, 'The tower of Hanoi and the Stern-Brocot Array', Thesis, Munich 1997

Crossrefs

Essentially same as A004526.
Harmonic mean of a(n) and A056136 is n.
a(n)=A010766(n+2, 2).
Cf. A010551 (partial products).
Cf. A263997 (a block spiral).
Cf. A289187.
Column 2 of A235791.

Programs

  • Haskell
    a008619 = (+ 1) . (`div` 2)
    a008619_list = concatMap (\x -> [x,x]) [1..]
    -- Reinhard Zumkeller, Apr 02 2012
    
  • Magma
    I:=[1,1,2]; [n le 3 select I[n] else Self(n-1)+Self(n-2)-Self(n-3): n in [1..100]]; // Vincenzo Librandi, Feb 04 2015
    
  • Maple
    a:= n-> iquo(n+2, 2): seq(a(n), n=0..75);
  • Mathematica
    Flatten[Table[{n,n},{n,35}]] (* Harvey P. Dale, Sep 20 2011 *)
    With[{c=Range[40]},Riffle[c,c]] (* Harvey P. Dale, Feb 23 2013 *)
    CoefficientList[Series[1/(1 - x - x^2 + x^3), {x, 0, 75}], x] (* Robert G. Wilson v, Feb 05 2015 *)
    LinearRecurrence[{1, 1, -1}, {1, 1, 2}, 75] (* Robert G. Wilson v, Feb 05 2015 *)
    Table[QBinomial[n, 2, -1], {n, 2, 75}] (* John Keith, Jun 28 2021 *)
  • PARI
    a(n)=n\2+1
    
  • Python
    def A008619(n): return (n>>1)+1 # Chai Wah Wu, Jul 07 2022
  • Sage
    a = lambda n: 1 if n==0 else a(n-1)+1 if 2.divides(n) else a(n-1) # Peter Luschny, Feb 05 2015
    
  • Scala
    (2 to 99).map( / 2) // _Alonso del Arte, May 09 2020
    

Formula

Euler transform of [1, 1].
a(n) = 1 + floor(n/2).
G.f.: 1/((1-x)(1-x^2)).
E.g.f.: ((3+2*x)*exp(x) + exp(-x))/4.
a(n) = a(n-1) + a(n-2) - a(n-3) = -a(-3-n).
a(0) = a(1) = 1 and a(n) = floor( (a(n-1) + a(n-2))/2 + 1 ).
a(n) = (2*n + 3 + (-1)^n)/4. - Paul Barry, May 27 2003
a(n) = Sum_{k=0..n} Sum_{j=0..k} Sum_{i=0..j} binomial(j, i)*(-2)^i. - Paul Barry, Aug 26 2003
E.g.f.: ((1+x)*exp(x) + cosh(x))/2. - Paul Barry, Sep 13 2003
a(n) = A108299(n-1,n)*(-1)^floor(n/2) for n > 0. - Reinhard Zumkeller, Jun 01 2005
a(n) = A108561(n+2,n) for n > 0. - Reinhard Zumkeller, Jun 10 2005
a(n) = A125291(A125293(n)) for n>0. - Reinhard Zumkeller, Nov 26 2006
a(n) = ceiling(n/2), n >= 1. - Mohammad K. Azarian, May 22 2007
INVERT transformation yields A006054 without leading zeros. INVERTi transformation yields negative of A124745 with the first 5 terms there dropped. - R. J. Mathar, Sep 11 2008
a(n) = A026820(n,2) for n > 1. - Reinhard Zumkeller, Jan 21 2010
a(n) = n - a(n-1) + 1 (with a(0)=1). - Vincenzo Librandi, Nov 19 2010
a(n) = A000217(n) / A110654(n). - Reinhard Zumkeller, Aug 24 2011
a(n+1) = A181971(n,n). - Reinhard Zumkeller, Jul 09 2012
1/(1+2/(2+3/(3+4/(4+5/(5+...(continued fraction))))) = 1/(e-1), see A073333. - Philippe Deléham, Mar 09 2013
a(n) = floor(A000217(n)/n), n > 0. - L. Edson Jeffery, Jul 26 2013
a(n) = n*a(n-1) mod (n+1) = -a(n-1) mod (n+1), the least positive residue modulo n+1 for each expression for n > 0, with a(0) = 1 (basically restatements of Vincenzo Librandi's formula). - Rick L. Shepherd, Apr 02 2014
a(n) = (a(0) + a(1) + ... + a(n-1))/a(n-1), where a(0) = 1. - Melvin Peralta, Jun 16 2015
a(n) = Sum_{k=0..n} (-1)^(n-k) * (k+1). - Rick L. Shepherd, Sep 18 2020
a(n) = a(n-2) + 1 for n >= 2. - Vladimír Modrák, Sep 29 2020
a(n) = A004526(n)+1. - Chai Wah Wu, Jul 07 2022

Extensions

Additional remarks from Daniele Parisse
Edited by N. J. A. Sloane, Sep 06 2009
Partially edited by Joerg Arndt, Mar 11 2010

A211795 Number of (w,x,y,z) with all terms in {1,...,n} and w*x < 2*y*z.

Original entry on oeis.org

0, 1, 11, 58, 177, 437, 894, 1659, 2813, 4502, 6836, 10008, 14121, 19449, 26117, 34372, 44422, 56597, 71044, 88160, 108115, 131328, 158074, 188773, 223604, 263172, 307719, 357715, 413493, 475690, 544480, 620632, 704381, 796413
Offset: 0

Views

Author

Clark Kimberling, Apr 27 2012

Keywords

Comments

Each sequence in the following guide counts 4-tuples
(w,x,y,z) such that the indicated relation holds and the four numbers w,x,y,z are in {1,...,n}. The notation "m div" means that m divides every term of the sequence.
A211058 ... wx <= yz
A211787 ... wx <= 2yz
A211795 ... wx < 2yz
A211797 ... wx > 2yz
A211809 ... wx >= 2yz
A211812 ... wx <= 3yz
A211917 ... wx < 3yz
A211918 ... wx > 3yz
A211919 ... wx >= 3yz
A211920 ... 2wx < 3yz
A211921 ... 2wx <= 3yz
A211922 ... 2wx > 3yz
A211923 ... 2wx >= 3yz
A212019 ... wx = 2yz ..... 2 div
A212020 ... wx = 3yz ..... 2 div
A212021 ... 2wx = 3yz .... 2 div
A212047 ... wx = 4yz
A212048 ... 3wx = 4yz .... 2 div
A212049 ... wx = 5yz ..... 2 div
A212050 ... 2wx = 5yz .... 2 div
A212051 ... 3wx = 5yz .... 2 div
A212052 ... 4wx = 5yz .... 2 div
A209978 ... wx = yz + 1 .. 2 div
A212053 ... wx <= yz + 1
A212054 ... wx > yz + 1
A212055 ... wx <= yz + 2
A212056 ... wx > yz + 2
A197168 ... wx = yz + 2 .. 2 div
A061201 ... w = xyz
A212057 ... w < xyz
A212058 ... w >= xyz
A212059 ... w = xyz - 1
A212060 ... w = xyz - 2
A212061 ... wx = (yz)^2
A212062 ... w^2 = xyz
A212063 ... w^2 < xyz
A212064 ... w^2 >= xyz
A212065 ... w^2 <= xyz
A212066 ... w^2 > xyz
A212067 ... w^3 = xyz
A002623 ... w = 2x + y + z
A006918 ... w = 2x + 2y + z
A000601 ... w = x + 2y + 3z (except for initial 0's)
A212068 ... 2w = x + y + z
A212069 ... 3w = x + y + z (w = average{x,y,z})
A212088 ... 3w < x + y + z
A212089 ... 3w >= x + y + z
A212090 ... w < x + y + z
A000332 ... w >= x + y + z
A212145 ... w < 2x + y + z
A001752 ... w >= 2x + y + z
A001400 ... w = 2x +3y + 4z
A005900 ... w = -x + y + z
A192023 ... w = -x + y + z + 2
A212091 ... w^2 = x^2 + y^2 + z^2 ... 3 div
A212087 ... w^2 + x^2 = y^2 + z^2
A212092 ... w^2 < x^2 + y^2 + z^2
A212093 ... w^2 <= x^2 + y^2 + z^2
A212094 ... w^2 > x^2 + y^2 + z^2
A212095 ... w^2 >= x^2 + y^2 + z^2
A212096 ... w^3 = x^3 + y^3 + z^3 ... 6 div
A212097 ... w^3 < x^3 + y^3 + z^3
A212098 ... w^3 <= x^3 + y^3 + z^3
A212099 ... w^3 > x^3 + y^3 + z^3
A212100 ... w^3 >= x^3 + y^3 + z^3
A212101 ... wx^2 = yz^2
A212102 ... 1/w = 1/x + 1/y + 1/z
A212103 ... 3/w = 1/x + 1/y + 1/z; w = h.m. of {x,y,z}
A212104 ... 3/w >= 1/x + 1/y + 1/z; w >= h.m.
A212105 ... 3/w < 1/x + 1/y + 1/z; w < h.m.
A212106 ... 3/w > 1/x + 1/y + 1/z; w > h.m.
A212107 ... 3/w <= 1/x + 1/y + 1/z; w <= h.m.
A212133 ... median(w,x,y,z) = mean(w,x,y,z)
A212134 ... median(w,x,y,z) <= mean(w,x,y,z)
A212135 ... median(w,x,y,z) > mean(w,x,y,z)
A212241 ... wx + yz > n
A212243 ... 2wx + yz = n
A212244 ... w = xyz - n
A212245 ... w = xyz - 2n
A212246 ... 2w = x + y + z - n
A212247 ... 3w = x + y + z + n
A212249 ... 3w < x + y + z + n
A212250 ... 3w >= x + y + z + n
A212251 ... 3w = x + y + z + n + 1
A212252 ... 3w = x + y + z + n + 2
A212254 ... w = x + 2y + 3z - n
A212255 ... w^2 = mean(x^2, y^2, z^2)
A212256 ... 4/w = 1/x + 1/y +1/z + 1/n
In the list above, if the relation in the second column is of the form "w rel ax + by + cz" then the sequence is linearly recurrent. In the list below, the same is true for expressions involving more than one relation.
A000332 ... w < x <= y < z .... C(n,4)
A000914 ... w < x <= y < z .... Stirling 1st kind
A000914 ... w < x <= y >= z ... Stirling 1st kind
A050534 ... w < x < y >= z .... tritriangular
A001296 ... w <= x <= y >= z .. 4-dim pyramidal
A006322 ... x < x > y >= z
A002418 ... w < x >= y < z
A050534 ... w < x >=y >= z
A212415 ... w < x >= y <= z
A001296 ... w < x >= y <= z
A212246 ... w <= x > y <= z
A006322 ... w <= x >= y <= z
A212501 ... w > x < y >= z
A212503 ... w < 2x and y < 2z ..... A (note below)
A212504 ... w < 2x and y > 2z ..... A
A212505 ... w < 2x and y >= 2z .... A
A212506 ... w <= 2x and y <= 2z ... A
A212507 ... w < 2x and y <= 2z .... B
A212508 ... w < 2x and y < 3z ..... C
A212509 ... w < 2x and y <= 3z .... C
A212510 ... w < 2x and y > 3z ..... C
A212511 ... w < 2x and y >= 3z .... C
A212512 ... w <= 2x and y < 3z .... C
A212513 ... w <= 2x and y <= 3z ... C
A212514 ... w <= 2x and y > 3z .... C
A212515 ... w <= 2x and y >= 3z ... C
A212516 ... w > 2x and y < 3z ..... C
A212517 ... w > 2x and y <= 3z .... C
A212518 ... w > 2x and y > 3z ..... C
A212519 ... w > 2x and y >= 3z .... C
A212520 ... w >= 2x and y < 3z .... C
A212521 ... w >= 2x and y <= 3z ... C
A212522 ... w >= 2x and y > 3z .... C
A212523 ... w + x < y + z
A212560 ... w + x <= y + z
A212561 ... w + x = 2y + 2z
A212562 ... w + x < 2y + 2z ....... B
A212563 ... w + x <= 2y + 2z ...... B
A212564 ... w + x > 2y + 2z ....... B
A212565 ... w + x >= 2y + 2z ...... B
A212566 ... w + x = 3y + 3z
A212567 ... 2w + 2x = 3y + 3z
A212570 ... |w - x| = |x - y| + |y - z|
A212571 ... |w - x| < |x - y| + |y - z| ... B ... 4 div
A212572 ... |w - x| <= |x - y| + |y - z| .. B
A212573 ... |w - x| > |x - y| + |y - z| ... B ... 2 div
A212574 ... |w - x| >= |x - y| + |y - z| .. B
A212575 ... 2|w - x| = |x - y| + |y - z|
A212576 ... |w - x| = 2|x - y| + 2|y - z|
A212577 ... |w - x| = 2|x - y| - |y - z|
A212578 ... 2|w - x| = |x - y| - |y - z|
A212579 ... min{|w-x|,|w-y|} = min{|x-y|,|x-z|}
A212692 ... w = |x - y| + |y - z| ............... 2 div
A212568 ... w < |x - y| + |y - z| ............... 2 div
A212573 ... w <= |x - y| + |y - z| .............. 2 div
A212574 ... w > |x - y| + |y - z|
A212575 ... w >= |x - y| + |y - z|
A212676 ... w + x = |x - y| + |y - z| ......... H
A212677 ... w + y = |x - y| + |y - z|
A212678 ... w + x + y = |x - y| + |y - z|
A006918 ... w + x + y + z = |x - y| + |y - z| . H
A212679 ... |x - y| = |y - z| ................. H
A212680 ... |x - y| = |y - z| + 1 ..............H 2 div
A212681 ... |x - y| < |y - z| ................... 2 div
A212682 ... |x - y| >= |y - z|
A212683 ... |x - y| = w + |y - z| ............... 2 div
A212684 ... |x - y| = n - w + |y - z|
A212685 ... |w - x| = w + |y - z|
A186707 ... |w - x| < w + |y - z| ... (Note D)
A212714 ... |w - x| >= w + |y - z| .......... H . 2 div
A212686 ... 2*|w - x| = n + |y - z| ............. 4 div
A212687 ... 2*|w - x| < n + |y - z| ......... B
A212688 ... 2*|w - x| < n + |y - z| ......... B . 2 div
A212689 ... 2*|w - x| > n + |y - z| ......... B . 2 div
A212690 ... 2*|w - x| <= n + |y - z| ........ B
A212691 ... w + |x - y| = |x - z| + |y - z| . E . 2 div
...
In the above lists, all the terms of (w,x,y,z) are in {1,...,n}, but in the next lists they are all in {0,...,n}, and sequences are all linearly recurrent.
R=range{w,x,y,z}=max{w,x,y,z}-min{w,x,y,z}.
A212740 ... max{w,x,y,z} < 2*min{w,x,y,z} .... A
A212741 ... max{w,x,y,z} >= 2*min{w,x,y,z} ... A
A212742 ... max{w,x,y,z} <= 2*min{w,x,y,z} ... A
A212743 ... max{w,x,y,z} > 2*min{w,x,y,z} .... A . 2 div
A212744 ... w=range (=max-min) ............... E
A212745 ... w=max{w,x,y,z} - 2*min{w,x,y,z}
A212746 ... R is in {w,x,y,z} ................ E
A212569 ... R is not in {w,x,y,z} ............ E
A212749 ... w=R or x
A212750 ... w=R or x=R or y
A212751 ... w=R or x=R or y
A212752 ... wR ......... A
A212753 ... wR or z>R ......... D
A212754 ... wR or y>R or z>R ......... D
A002415 ... w = x + R ........................ D
A212755 ... |w - x| = R ...................... D
A212756 ... 2w = x + R
A212757 ... 2w = R
A212758 ... w = floor(R/2)
A002413 ... w = floor((x+y+z/2))
A212759 ... w, x, y are even
A212760 ... w is even and x = y + z .......... E
A212761 ... w is odd and x and y are even .... F . 2 div
A212762 ... w and x are odd y is even ........ F . 2 div
A212763 ... w, x, y are odd .................. F
A212764 ... w, x, y are even and z is odd .... F
A030179 ... w and x are even and y and z odd
A212765 ... w is even and x,y,z are odd ...... F
A212766 ... w is even and x is odd ........... A . 2 div
A212767 ... w and x are even and w+x=y+z ..... E
A212889 ... R is even ........................ A
A212890 ... R is odd ......................... A . 2 div
A212742 ... w-x, x-y, y-z are all even ....... A
A212892 ... w-x, x-y, y-z are all odd ........ A
A212893 ... w-x, x-y, y-z have same parity ... A
A005915 ... min{|w-x|, |x-y|, |y-z|} = 0
A212894 ... min{|w-x|, |x-y|, |y-z|} = 1
A212895 ... min{|w-x|, |x-y|, |y-z|} = 2
A179824 ... min{|w-x|, |x-y|, |y-z|} > 0
A212896 ... min{|w-x|, |x-y|, |y-z|} <= 1
A212897 ... min{|w-x|, |x-y|, |y-z|} > 1
A212898 ... min{|w-x|, |x-y|, |y-z|} <= 2
A212899 ... min{|w-x|, |x-y|, |y-z|} > 2
A212901 ... |w-x| = |x-y| = |y-z|
A212900 ... |w-x|, |x-y|, |y-z| are distinct . G
A212902 ... |w-x| < |x-y| < |y-z| ............ G
A212903 ... |w-x| <= |x-y| <= |y-z| .......... G
A212904 ... |w-x| + |x-y| + |y-z| = n ........ H
A212905 ... |w-x| + |x-y| + |y-z| = 2n ....... H
...
Note A: A212503-A212506 (and others) have these recurrence coefficients: 2,2,-6,0,6,-2,-2,1.
B: 3,-1,-5,5,1,-3,1
C: 0,2,2,-1,-4,0,2,0,-2,0,4,1,-2,-2,0,1
D: 4,-5,0,5,-4,1
E: 1,3,-3,-3,3,1,-1
F: 1,4,-4,-6,6,4,-4,-1,1
G: 2,1,-3,-1,1,3,-1,-2,1
H: 2,1,-4,1,2,-1

Examples

			a(2)=11 counts these (w,x,y,z): (1,1,1,1), (1,1,1,2), (1,1,2,1), (2,1,2,1), (2,1,1,2), (1,2,2,1), (1,2,1,2), (1,1,2,2), (1,2,2,2), (2,1,2,2), (2,2,2,2).
		

References

  • A. Barvinok, Lattice Points and Lattice Polytopes, Chapter 7 in Handbook of Discrete and Computational Geometry, CRC Press, 1997, 133-152.
  • P. Gritzmann and J. M. Wills, Lattice Points, Chapter 3.2 in Handbook of Convex Geometry, vol. B, North-Holland, 1993, 765-797.

Crossrefs

Programs

  • Mathematica
    t = Compile[{{n, _Integer}}, Module[{s = 0},
        (Do[If[w*x < 2 y*z, s = s + 1], {w, 1, #},
          {x, 1, #}, {y, 1, #}, {z, 1, #}] &[n]; s)]];
    Map[t[#] &, Range[0, 40]] (* A211795 *)
    (* Peter J. C. Moses, Apr 13 2012 *)

Formula

a(n) = n^4 - A211809(n).

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

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

A026820 Euler's table: triangular array T read by rows, where T(n,k) = number of partitions in which every part is <= k for 1 <= k <= n. Also number of partitions of n into at most k parts.

Original entry on oeis.org

1, 1, 2, 1, 2, 3, 1, 3, 4, 5, 1, 3, 5, 6, 7, 1, 4, 7, 9, 10, 11, 1, 4, 8, 11, 13, 14, 15, 1, 5, 10, 15, 18, 20, 21, 22, 1, 5, 12, 18, 23, 26, 28, 29, 30, 1, 6, 14, 23, 30, 35, 38, 40, 41, 42, 1, 6, 16, 27, 37, 44, 49, 52, 54, 55, 56, 1, 7, 19, 34, 47, 58, 65, 70, 73, 75, 76, 77
Offset: 1

Keywords

Examples

			Triangle starts:
  1;
  1, 2;
  1, 2,  3;
  1, 3,  4,  5;
  1, 3,  5,  6,  7;
  1, 4,  7,  9, 10, 11;
  1, 4,  8, 11, 13, 14, 15;
  1, 5, 10, 15, 18, 20, 21, 22;
  1, 5, 12, 18, 23, 26, 28, 29, 30;
  1, 6, 14, 23, 30, 35, 38, 40, 41, 42;
  1, 6, 16, 27, 37, 44, 49, 52, 54, 55, 56;
  ...
		

References

  • G. Chrystal, Algebra, Vol. II, p. 558.
  • D. S. Mitrinovic et al., Handbook of Number Theory, Kluwer, Section XIV.2, p. 493.

Crossrefs

Partial sums of rows of A008284, row sums give A058397, central terms give A171985, mirror is A058400.
T(n,n) = A000041(n), T(n,1) = A000012(n), T(n,2) = A008619(n) for n>1, T(n,3) = A001399(n) for n>2, T(n,4) = A001400(n) for n>3, T(n,5) = A001401(n) for n>4, T(n,6) = A001402(n) for n>5, T(n,7) = A008636(n) for n>6, T(n,8) = A008637(n) for n>7, T(n,9) = A008638(n) for n>8, T(n,10) = A008639(n) for n>9, T(n,11) = A008640(n) for n>10, T(n,12) = A008641(n) for n>11, T(n,n-2) = A007042(n-1) for n>2, T(n,n-1) = A000065(n) for n>1.

Programs

  • Haskell
    import Data.List (inits)
    a026820 n k = a026820_tabl !! (n-1) !! (k-1)
    a026820_row n = a026820_tabl !! (n-1)
    a026820_tabl = zipWith
       (\x -> map (p x) . tail . inits) [1..] $ tail $ inits [1..] where
       p 0 _ = 1
       p _ [] = 0
       p m ks'@(k:ks) = if m < k then 0 else p (m - k) ks' + p m ks
    -- Reinhard Zumkeller, Dec 18 2013
    
  • Maple
    T:= proc(n, k) option remember;
          `if`(n=0 or k=1, 1, T(n, k-1) + `if`(k>n, 0, T(n-k, k)))
        end:
    seq(seq(T(n, k), k=1..n), n=1..12); # Alois P. Heinz, Apr 21 2012
  • Mathematica
    t[n_, k_] := Length@ IntegerPartitions[n, k]; Table[ t[n, k], {n, 12}, {k, n}] // Flatten
    (* Second program: *)
    T[n_, k_] := T[n, k] = If[n==0 || k==1, 1, T[n, k-1] + If[k>n, 0, T[n-k, k]]]; Table[T[n, k], {n, 1, 12}, {k, 1, n}] // Flatten (* Jean-François Alcover, Sep 22 2015, after Alois P. Heinz *)
  • PARI
    T(n,k)=my(s); forpart(v=n,s++,,k); s \\ Charles R Greathouse IV, Feb 27 2018
    
  • SageMath
    from sage.combinat.partition import number_of_partitions_length
    from itertools import accumulate
    for n in (1..11):
        print(list(accumulate([number_of_partitions_length(n, k) for k in (1..n)])))
    # Peter Luschny, Jul 28 2022

Formula

T(T(n,n),n) = A134737(n). - Reinhard Zumkeller, Nov 07 2007
T(A000217(n),n) = A173519(n). - Reinhard Zumkeller, Feb 20 2010
T(n,k) = T(n,k-1) + T(n-k,k). - Thomas Dybdahl Ahle, Jun 13 2011
T(n,k) = Sum_{i=1..min(k,floor(n/2))} T(n-i,i) + Sum_{j=1+floor(n/2)..k} A000041(n-j). - Bob Selcoe, Aug 22 2014 [corrected by Álvar Ibeas, Mar 15 2018]
O.g.f.: Product_{i>=0} 1/(1-y*x^i). - Geoffrey Critzer, Mar 11 2012
T(n,k) = A008284(n+k,k). - Álvar Ibeas, Jan 06 2015

A026810 Number of partitions of n in which the greatest part is 4.

Original entry on oeis.org

0, 0, 0, 0, 1, 1, 2, 3, 5, 6, 9, 11, 15, 18, 23, 27, 34, 39, 47, 54, 64, 72, 84, 94, 108, 120, 136, 150, 169, 185, 206, 225, 249, 270, 297, 321, 351, 378, 411, 441, 478, 511, 551, 588, 632, 672, 720, 764, 816, 864, 920, 972, 1033, 1089, 1154, 1215, 1285, 1350
Offset: 0

Keywords

Comments

Also number of partitions of n into exactly 4 parts.
Also the number of weighted cubic graphs on 4 nodes (=the tetrahedron) with weight n. - R. J. Mathar, Nov 03 2018
From Gus Wiseman, Jun 27 2021: (Start)
Also the number of strict integer partitions of 2n with alternating sum 4, or (by conjugation) partitions of 2n covering an initial interval of positive integers with exactly 4 odd parts. The strict partitions with alternating sum 4 are:
(4) (5,1) (6,2) (7,3) (8,4) (9,5) (10,6)
(5,2,1) (5,3,2) (5,4,3) (6,5,3) (7,6,3)
(6,3,1) (6,4,2) (7,5,2) (8,6,2)
(7,4,1) (8,5,1) (9,6,1)
(6,3,2,1) (6,4,3,1) (6,5,4,1)
(7,4,2,1) (7,4,3,2)
(7,5,3,1)
(8,5,2,1)
(6,4,3,2,1)
(End)

Examples

			From _Gus Wiseman_, Jun 27 2021: (Start)
The a(4) = 1 through a(10) = 9 partitions of length 4:
  (1111)  (2111)  (2211)  (2221)  (2222)  (3222)  (3322)
                  (3111)  (3211)  (3221)  (3321)  (3331)
                          (4111)  (3311)  (4221)  (4222)
                                  (4211)  (4311)  (4321)
                                  (5111)  (5211)  (4411)
                                          (6111)  (5221)
                                                  (5311)
                                                  (6211)
                                                  (7111)
(End)
		

References

  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 275.
  • D. E. Knuth, The Art of Computer Programming, vol. 4,fascicle 3, Generating All Combinations and Partitions, Section 7.2.1.4., p. 56, exercise 31.

Crossrefs

Cf. A001400, A026811, A026812, A026813, A026814, A026815, A026816, A069905 (3 positive parts), A002621 (partial sums), A005044 (first differences).
A non-strict version is A000710 or A088218.
This is column k = 2 of A152146.
A reverse version is A343941.

Programs

  • Magma
    [Round((n^3+3*n^2-9*n*(n mod 2))/144): n in [0..60]]; // Vincenzo Librandi, Oct 14 2015
  • Maple
    A049347 := proc(n)
            op(1+(n mod 3),[1,-1,0]) ;
    end proc:
    A056594 := proc(n)
            op(1+(n mod 4),[1,0,-1,0]) ;
    end proc:
    A026810 := proc(n)
            1/288*(n+1)*(2*n^2+4*n-13+9*(-1)^n) ;
            %-A049347(n)/9 ;
            %+A056594(n)/8 ;
    end proc: # R. J. Mathar, Jul 03 2012
  • Mathematica
    Table[Count[IntegerPartitions[n], {4, _}], {n, 0, 60}]
    LinearRecurrence[{1, 1, 0, 0, -2, 0, 0, 1, 1, -1}, {0, 0, 0, 0, 1, 1, 2, 3, 5, 6}, 60] (* Vincenzo Librandi, Oct 14 2015 *)
    Table[Length[IntegerPartitions[n, {4}]], {n, 0, 60}] (* Eric Rowland, Mar 02 2017 *)
    CoefficientList[Series[x^4/Product[1 - x^k, {k, 1, 4}], {x, 0, 60}], x] (* Robert A. Russell, May 13 2018 *)
  • PARI
    for(n=0, 60, print(n, " ", round((n^3 + 3*n^2 -9*n*(n % 2))/144))); \\ Washington Bomfim, Jul 03 2012
    
  • PARI
    x='x+O('x^60); concat([0, 0, 0, 0], Vec(x^4/((1-x)*(1-x^2)*(1-x^3)*(1-x^4)))) \\ Altug Alkan, Oct 14 2015
    
  • PARI
    vector(60, n, n--; (n+1)*(2*n^2+4*n-13+9*(-1)^n)/288 + real(I^n)/8 - ((n+2)%3-1)/9) \\ Altug Alkan, Oct 26 2015
    
  • PARI
    print1(0,", "); for(n=1,60,j=0;forpart(v=n,j++,,[4,4]); print1(j,", ")) \\ Hugo Pfoertner, Oct 01 2018
    

Formula

G.f.: x^4/((1-x)*(1-x^2)*(1-x^3)*(1-x^4)) = x^4/((1-x)^4*(1+x)^2*(1+x+x^2)*(1+x^2)).
a(n+4) = A001400(n). - Michael Somos, Apr 07 2012
a(n) = round( (n^3 + 3*n^2 -9*n*(n mod 2))/144 ). - Washington Bomfim, Jan 06 2021 and Jul 03 2012
a(n) = (n+1)*(2*n^2+4*n-13+9*(-1)^n)/288 -A049347(n)/9 +A056594(n)/8. - R. J. Mathar, Jul 03 2012
From Gregory L. Simay, Oct 13 2015: (Start)
a(n) = (n^3 + 3*n^2 - 9*n)/144 + a(m) - (m^3 + 3*m^2 - 9*m)/144 if n = 12k + m and m is odd. For example, a(23) = a(12*1 + 11) = (23^3 + 3*23^2 - 9*23)/144 + a(11) - (11^3 + 3*11^2 - 9*11)/144 = 94.
a(n) = (n^3 + 3*n^2)/144 + a(m) - (m^3 + 3*m^2)/144 if n = 12k + m and m is even. For example, a(22) = a(12*1 + 10) = (22^3 + 3*22^2)/144 + a(10) - (10^3 + 3*10^2)/144 = 84. (End)
a(n) = A008284(n,4). - Robert A. Russell, May 13 2018
From Gregory L. Simay, Jul 28 2019: (Start)
a(2n+1) = a(2n) + a(n+1) - a(n-3) and
a(2n) = a(2n-1) + a(n+2) - a(n-2). (End)

A001401 Number of partitions of n into at most 5 parts.

Original entry on oeis.org

1, 1, 2, 3, 5, 7, 10, 13, 18, 23, 30, 37, 47, 57, 70, 84, 101, 119, 141, 164, 192, 221, 255, 291, 333, 377, 427, 480, 540, 603, 674, 748, 831, 918, 1014, 1115, 1226, 1342, 1469, 1602, 1747, 1898, 2062, 2233, 2418, 2611, 2818, 3034, 3266, 3507, 3765, 4033, 4319
Offset: 0

Keywords

Comments

a(n) = T_{r}(n) for r large, where T_{r}(n) = number of outcomes in which r indistinguishable dice yield a sum r+n-1.
a(n) = coefficient of q^n in the expansion of (m choose 5)_q as m goes to infinity. - Y. Kelly Itakura (yitkr(AT)mta.ca), Aug 21 2002
For n > 4: also number of partitions of n into parts <= 5: a(n) = A026820(n,5). - Reinhard Zumkeller, Jan 21 2010
Number of different distributions of n+15 identical balls in 5 boxes as x,y,z,p,q where 0 < x < y < z < p < q. - Ece Uslu and Esin Becenen, Jan 11 2016 [i.e., a(n) is the number of partitions of n+15 into 5 distinct parts. - R. J. Mathar, Feb 28 2021]
Tengely and Ulas prove that a(n) is a square only for n=1 and 2027. - Michel Marcus, Feb 11 2021

Examples

			(5 choose 5)_q = 1;
(6 choose 5)_q = q^5 + q^4 + q^3 + q^2 + q + 1;
(7 choose 5)_q = q^10 + q^9 + 2*q^8 + 2*q^7 + 3*q^6 + 3*q^5 + 3*q^4 + 2*q^3 + 2*q^2 + q + 1;
(8 choose 5)_q = q^15 + q^14 + 2*q^13 + 3*q^12 + 4*q^11 + 5*q^10 + 6*q^9 + 6*q^8 + 6*q^7 + 6*q^6 + 5*q^5 + 4*q^4 + 3*q^3 + 2*q^2 + q + 1;
so the coefficient of q^0 converges to 1, q^1 to 1, q^2 to 2 and so on.
a(3) = 3, i.e., {1,2,3,4,8}, {1,2,3,5,7}, {1,2,4,5,6}. Number of different distributions of 18 identical balls in 5 boxes as x,y,z,p,q where 0 < x < y < z < p < q. - _Ece Uslu_, Esin Becenen, Jan 11 2016
		

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 115, row m=5 of Q(m,n) table.
  • H. Gupta et al., Tables of Partitions. Royal Society Mathematical Tables, Vol. 4, Cambridge Univ. Press, 1958, p. 2.
  • D. E. Knuth, The Art of Computer Programming, vol. 4, fascicle 3, Generating All Combinations and Partitions, Section 7.2.1.4., p. 56, exercise 31.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

a(n) = A008284(n+5, 5), n >= 0.
Cf. A008619, A001400, A001399, A008667 (first differences), A008804.
First differences of A002622.

Programs

  • Maple
    with(combstruct):ZL6:=[S,{S=Set(Cycle(Z,card<6))}, unlabeled]:seq(count(ZL6,size=n),n=0..52); # Zerinvary Lajos, Sep 24 2007
    a:= n-> (Matrix(15, (i,j)-> if (i=j-1) then 1 elif j=1 then [1, 1, 0, 0, -1, -1, -1, 1, 1, 1, 0, 0, -1, -1, 1][i] else 0 fi)^n)[1,1]: seq(a(n), n=0..60); # Alois P. Heinz, Jul 31 2008
    B:=[S,{S = Set(Sequence(Z,1 <= card),card <=5)},unlabelled]: seq(combstruct[count](B, size=n), n=0..52); # Zerinvary Lajos, Mar 21 2009
  • Mathematica
    CoefficientList[ Series[ 1/((1 - x)*(1 - x^2)*(1 - x^3)*(1 - x^4)*(1 - x^5)), {x, 0, 60} ], x ]
    a[n_] := IntegerPartitions[n, 5] // Length; Table[a[n], {n, 0, 52}] (* Jean-François Alcover, Jul 13 2012 *)
    LinearRecurrence[{1,1,0,0,-1,-1,-1,1,1,1,0,0,-1,-1,1},{1,1,2,3,5,7,10,13,18,23,30,37,47,57,70},60] (* Harvey P. Dale, Jan 05 2019 *)
  • PARI
    a(n)=#partitions(n,,5) \\ Charles R Greathouse IV, Sep 15 2014
    
  • PARI
    a(n) = (n^4 + 30*n^3 + 310*n^2 + 1320*n - 90*n*(n%2) + 2880)\2880 \\ Hoang Xuan Thanh, Aug 12 2025

Formula

G.f.: 1/((1-x)*(1-x^2)*(1-x^3)*(1-x^4)*(1-x^5)).
a(n) = 1 + (a(n-2) + a(n-3) + a(n-4)) - (a(n-6) + (2*a(n-7)) + a(n-8)) + (a(n-10) + a(n-11) + a(n-12)) - a(n-14). - Norman J. Meluch (norm(AT)iss.gm.com), Mar 09 2000
Let a1(n) = Sum_{i=0..floor(n/3)} (1 + ceiling((n-3*i-1)/2)), a2(n) = Sum_{i=0..floor(n/4)} (1 + ceiling((n-4*i-1)/2) + a1(n-4*i-3)), then a(n) = Sum_{i=0..floor(n/5)} (1 + ceiling((n-5*i-1)/2) + a1(n-5*i-3) + a2(n-5*i-4)). - Jon Perry, Jun 27 2003
(n choose 5)_q=(q^n-1)*(q^(n-1)-1)*(q^(n-2)-1)*(q^(n-3)-1)*(q^(n-4)-1)/((q^5-1)*(q^4-1)*(q^3-1)*(q^2-1)*(q-1)).
a(n) = round(((n+5)^4 + 10*((n+5)^3 + (n+5)^2) - 75*(n+5) - 45*(n+5)*(-1)^(n+5))/2880). - Washington Bomfim, Jul 03 2012
a(n) = a(n-1) + a(n-2) - a(n-5) - a(n-6) - a(n-7) + a(n-8) + a(n-9) + a(n-10) - a(n-13) - a(n-14) + a(n+15). - David Neil McGrath, Sep 13 2014
a(n+5) = a(n) + A001400(n) = A001400(n)+A026811(n). - Ece Uslu, Esin Becenen, Jan 11 2016
From Vladimír Modrák, Jul 13 2022: (Start)
a(n) = Sum_{k=0..floor(n/5)} Sum_{j=0..floor(n/4)} Sum_{i=0..floor(n/3)} ceiling((max(0, n + 1 - 3*i - 4*j - 5*k))/2).
a(n) = Sum_{j=0..floor(n/5)} Sum_{i=0..floor(n/4)} floor(((max(0, n + 3 - 4*i - 5*j))^2+4)/12). (End)
a(2n) = a(2n-1) + a(n) - a(n-8) = a(n) + Sum_{k=0..n-1} A008804(k). - David García Herrero, Aug 26 2024
a(n) = floor((n^4 + 30*n^3 + 310*n^2 + 1275*n + 45*n*(-1)^n+2880)/2880). - Hoang Xuan Thanh, Aug 12 2025

Extensions

Additional comments from Michael Somos and Branislav Kisacanin (branislav.kisacanin(AT)delphiauto.com)

A266755 Expansion of 1/((1-x^2)*(1-x^3)*(1-x^4)).

Original entry on oeis.org

1, 0, 1, 1, 2, 1, 3, 2, 4, 3, 5, 4, 7, 5, 8, 7, 10, 8, 12, 10, 14, 12, 16, 14, 19, 16, 21, 19, 24, 21, 27, 24, 30, 27, 33, 30, 37, 33, 40, 37, 44, 40, 48, 44, 52, 48, 56, 52, 61, 56, 65, 61, 70, 65, 75, 70, 80, 75, 85, 80, 91, 85, 96, 91, 102, 96, 108, 102, 114, 108, 120, 114, 127, 120, 133, 127, 140, 133, 147, 140, 154, 147, 161, 154, 169
Offset: 0

Author

N. J. A. Sloane, Jan 10 2016

Keywords

Comments

This is the same as A005044 but without the three leading zeros. There are so many situations where one wants this sequence rather than A005044 that it seems appropriate for it to have its own entry.
But see A005044 (still the main entry) for numerous applications and references.
Also, Molien series for invariants of finite Coxeter group D_3.
The Molien series for the finite Coxeter group of type D_k (k >= 3) has g.f. = 1/Product_i (1-x^(1+m_i)) where the m_i are [1,3,5,...,2k-3,k-1]. If k is even only even powers of x appear, and we bisect the sequence.
Also, Molien series for invariants of finite Coxeter group A_3. The Molien series for the finite Coxeter group of type A_k (k >= 1) has g.f. = 1/Product_{i=2..k+1} (1-x^i). Note that this is the root system A_k not the alternating group Alt_k.
a(n) is the number of partitions of n into parts 2, 3, and 4. - Joerg Arndt, Apr 16 2017
From Gus Wiseman, May 23 2021: (Start)
Also the number of integer partitions of n into at most n/2 parts, none greater than 3. The case of any maximum is A110618. The case of any length is A001399. The Heinz numbers of these partitions are given by A344293.
For example, the a(2) = 1 through a(13) = 5 partitions are:
2 3 22 32 33 322 332 333 3322 3332 3333 33322
31 222 331 2222 3222 3331 32222 33222 33331
321 3221 3321 22222 33221 33321 322222
3311 32221 33311 222222 332221
33211 322221 333211
332211
333111
(End)

Examples

			G.f. = 1 + x^2 + x^3 + 2*x^4 + x^5 + 3*x^6 + 2*x^7 + 4*x^8 + ... - _Michael Somos_, Jan 29 2022
		

References

  • J. E. Humphreys, Reflection Groups and Coxeter Groups, Cambridge, 1990. See Table 3.1, page 59.

Crossrefs

Molien series for finite Coxeter groups A_1 through A_12 are A059841, A103221, A266755, A008667, A037145, A001996, and A266776-A266781.
Molien series for finite Coxeter groups D_3 through D_12 are A266755, A266769, A266768, A003402, and A266770-A266775.
A variant of A005044.
Cf. A001400 (partial sums).
Cf. A308065.
Number of partitions of n whose Heinz number is in A344293.
A001399 counts partitions with all parts <= 3, ranked by A051037.
A025065 counts partitions of n with >= n/2 parts, ranked by A344296.
A035363 counts partitions of n with n/2 parts, ranked by A340387.
A110618 counts partitions of n into at most n/2 parts, ranked by A344291.

Programs

  • Magma
    I:=[1,0,1,1,2,1,3,2,4]; [n le 9 select I[n] else Self(n-2)+ Self(n-3)+Self(n-4)-Self(n-5)-Self(n-6)-Self(n-7)+Self(n-9): n in [1..100]]; // Vincenzo Librandi, Jan 11 2016
    
  • Mathematica
    CoefficientList[Series[1/((1-x^2)(1-x^3)(1-x^4)), {x, 0, 100}], x] (* JungHwan Min, Jan 10 2016 *)
    LinearRecurrence[{0,1,1,1,-1,-1,-1,0,1}, {1,0,1,1,2,1,3,2,4}, 100] (* Vincenzo Librandi, Jan 11 2016 *)
    Table[Length[Select[IntegerPartitions[n],Length[#]<=n/2&&Max@@#<=3&]],{n,0,30}] (* Gus Wiseman, May 23 2021 *)
    a[ n_] := Round[(n + 3*(2 - Mod[n,2]))^2/48]; (* Michael Somos, Jan 29 2022 *)
  • PARI
    Vec(1/((1-x^2)*(1-x^3)*(1-x^4)) + O(x^100)) \\ Michel Marcus, Jan 11 2016
    
  • PARI
    {a(n) = round((n + 3*(2-n%2))^2/48)}; /* Michael Somos, Jan 29 2022 */
    
  • Sage
    (1/((1-x^2)*(1-x^3)*(1-x^4))).series(x, 100).coefficients(x, sparse=False) # G. C. Greubel, Jun 13 2019

Formula

a(n) = a(n-2) + a(n-3) + a(n-4) - a(n-5) - a(n-6) - a(n-7) + a(n-9) for n>8. - Vincenzo Librandi, Jan 11 2016
a(n) = a(-9-n) for all n in Z. a(n) = a(n+3) for all n in 2Z. - Michael Somos, Jan 29 2022
E.g.f.: exp(-x)*(81 - 18*x + exp(2*x)*(107 + 60*x + 6*x^2) + 64*exp(x/2)*cos(sqrt(3)*x/2) + 36*exp(x)*(cos(x) - sin(x)))/288. - Stefano Spezia, Mar 05 2023
For n >= 3, if n is even, a(n) = a(n-3) + floor(n/4) + 1, otherwise a(n) = a(n-3). - Robert FERREOL, Feb 05 2024
a(n) = floor((n^2+9*n+(3*n+9)*(-1)^n+39)/48). - Hoang Xuan Thanh, Jun 03 2025

A026811 Number of partitions of n in which the greatest part is 5.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 1, 2, 3, 5, 7, 10, 13, 18, 23, 30, 37, 47, 57, 70, 84, 101, 119, 141, 164, 192, 221, 255, 291, 333, 377, 427, 480, 540, 603, 674, 748, 831, 918, 1014, 1115, 1226, 1342, 1469, 1602, 1747, 1898, 2062, 2233, 2418, 2611, 2818, 3034, 3266, 3507, 3765
Offset: 0

Keywords

Comments

Essentially same as A001401: five zeros followed by A001401.
Also number of partitions of n into exactly 5 parts.

References

  • D. E. Knuth, The Art of Computer Programming, vol. 4, fascicle 3, Generating All Combinations and Partitions, Section 7.2.1.4., p. 56, exercise 31.

Crossrefs

Cf. A026810, A026812, A026813, A026814, A026815, A026816, A002622 (partial sums), A008667 (first differences).

Programs

  • GAP
    List([0..70],n->NrPartitions(n,5)); # Muniru A Asiru, May 17 2018
  • Mathematica
    Table[Count[IntegerPartitions[n], {5, _}], {n, 0, 55}] (* corrected by Harvey P. Dale, Oct 24 2011 *)
    Table[Length[IntegerPartitions[n, {5}]], {n, 0, 55}] (* Eric Rowland, Mar 02 2017 *)
    CoefficientList[Series[x^5/Product[1 - x^k, {k, 1, 5}], {x, 0, 65}], x] (* Robert A. Russell, May 13 2018 *)
    Drop[LinearRecurrence[{1,1,0,0,-1,-1,-1,1,1,1,0,0,-1,-1,1}, Append[Table[0,{14}],1],110],9] (* Robert A. Russell, May 17 2018 *)
  • PARI
    a(n)=round((n^4+10*(n^3+n^2)-75*n-45*(-1)^n*n)/2880);
    for(n=0,10000,print(n," ",a(n))); /* b-file format */
    /* Washington Bomfim, Jul 03 2012 */
    
  • PARI
    x='x+O('x^99); concat(vector(5), Vec(x^5/prod(k=1, 5, 1-x^k))) \\ Altug Alkan, May 17 2018
    

Formula

a(n) = round( ((n^4+10*(n^3+n^2)-75*n -45*n*(-1)^n)) / 2880 ). - Washington Bomfim, Jul 03 2012
G.f.: x^5/((1-x)*(1-x^2)*(1-x^3)*(1-x^4)*(1-x^5)). - Joerg Arndt, Jul 04 2012
a(n) = A008284(n,5). - Robert A. Russell, May 13 2018
From Gregory L. Simay, Jul 28 2019: (Start)
a(2n) = a(2n-1) + a(n+1) + a(n) - a(n-3) - a(n-4);
a(2n+1) = a(2n) + a(n+3) - a(n-5). (End)
From R. J. Mathar, Jun 23 2021: (Start)
a(n) - a(n-5) = A001400(n-5).
a(n) - a(n-4) = A008669(n-5).
a(n) - a(n-3) = A029007(n-5).
a(n) - a(n-2) = A029032(n-5).
a(n) = +a(n-1) +a(n-2) -a(n-5) -a(n-6) -a(n-7) +a(n-8) +a(n-9) +a(n-10) -a(n-13) -a(n-14) +a(n-15). (End)

Extensions

More terms from Robert G. Wilson v, Jan 11 2002
a(0)=0 inserted by Joerg Arndt, Jul 04 2012
Showing 1-10 of 50 results. Next