cp's OEIS Frontend

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

Previous Showing 21-30 of 198 results. Next

A003983 Array read by antidiagonals with T(n,k) = min(n,k).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 2, 3, 2, 1, 1, 2, 3, 3, 2, 1, 1, 2, 3, 4, 3, 2, 1, 1, 2, 3, 4, 4, 3, 2, 1, 1, 2, 3, 4, 5, 4, 3, 2, 1, 1, 2, 3, 4, 5, 5, 4, 3, 2, 1, 1, 2, 3, 4, 5, 6, 5, 4, 3, 2, 1, 1, 2, 3, 4, 5, 6, 6, 5, 4, 3, 2, 1, 1, 2, 3, 4, 5, 6, 7, 6, 5, 4, 3, 2, 1, 1, 2, 3, 4, 5, 6, 7, 7, 6, 5, 4, 3, 2, 1
Offset: 1

Views

Author

Keywords

Comments

Also, "correlation triangle" for the constant sequence 1. - Paul Barry, Jan 16 2006
Antidiagonal sums are in A002620.
As a triangle, row sums are A002620. T(2n,n)=n+1. Diagonal sums are A001399. Construction: Take antidiagonal triangle of MM^T where M is the sequence array for the constant sequence 1 (lower triangular matrix with all 1's). - Paul Barry, Jan 16 2006
From Franklin T. Adams-Watters, Sep 25 2011: (Start)
As a triangle, count up to ceiling(n/2) and back down again (repeating the central term when n is even).
When the first two instances of each number are removed from the sequence, the original sequence is recovered.
(End)

Examples

			Triangle version begins
  1;
  1, 1;
  1, 2, 1;
  1, 2, 2, 1;
  1, 2, 3, 2, 1;
  1, 2, 3, 3, 2, 1;
  1, 2, 3, 4, 3, 2, 1;
  1, 2, 3, 4, 4, 3, 2, 1;
  1, 2, 3, 4, 5, 4, 3, 2, 1;
  ...
		

Crossrefs

Programs

  • Haskell
    a003983 n k = a003983_tabl !! (n-1) !! (k-1)
    a003983_tabl = map a003983_row [1..]
    a003983_row n = hs ++ drop m (reverse hs)
       where hs = [1..n' + m]
             (n',m) = divMod n 2
    -- Reinhard Zumkeller, Aug 14 2011
    
  • Maple
    a(n) = min(floor(1/2 + sqrt(2*n)) - (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2+1, (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2) # Leonid Bedratyuk, Dec 13 2009
  • Mathematica
    Flatten[Table[Min[n-k+1, k], {n, 1, 14}, {k, 1, n}]] (* Jean-François Alcover, Feb 23 2012 *)
  • PARI
    T(n,k) = min(n,k) \\ Charles R Greathouse IV, Feb 06 2017
    
  • Python
    from math import isqrt
    def A003983(n):
        a = (m:=isqrt(k:=n<<1))+(k>m*(m+1))
        x = n-(a*(a-1)>>1)
        return min(x,a-x+1) # Chai Wah Wu, Jun 14 2025

Formula

Number triangle T(n, k) = Sum_{j=0..n} [j<=k][j<=n-k]. - Paul Barry, Jan 16 2006
G.f.: 1/((1-x)*(1-x*y)*(1-x^2*y)). - Christian G. Bower, Jan 17 2006
a(n) = min(floor( 1/2 + sqrt(2*n)) - (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2+1, (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2). - Leonid Bedratyuk, Dec 13 2009

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Nov 08 2000
Entry revised by N. J. A. Sloane, Dec 05 2006

A220377 Number of partitions of n into three distinct and mutually relatively prime parts.

Original entry on oeis.org

1, 0, 2, 1, 3, 1, 6, 1, 7, 3, 7, 3, 14, 3, 15, 6, 14, 6, 25, 6, 22, 10, 25, 9, 42, 8, 34, 15, 37, 15, 53, 13, 48, 22, 53, 17, 78, 17, 65, 30, 63, 24, 99, 24, 88, 35, 84, 30, 126, 34, 103, 45, 103, 38, 166, 35, 124, 57, 128, 51, 184, 44, 150, 67, 172, 52, 218
Offset: 6

Views

Author

Carl Najafi, Dec 13 2012

Keywords

Comments

The Heinz numbers of these partitions are the intersection of A005117 (strict), A014612 (triples), and A302696 (coprime). - Gus Wiseman, Oct 14 2020

Examples

			For n=10 we have three such partitions: 1+2+7, 1+4+5 and 2+3+5.
From _Gus Wiseman_, Oct 14 2020: (Start)
The a(6) = 1 through a(20) = 15 triples (empty column indicated by dot, A..H = 10..17):
321  .  431  531  532  731  543  751  743  753  754  971  765  B53  875
        521       541       651       752  951  853  B51  873  B71  974
                  721       732       761  B31  871  D31  954  D51  A73
                            741       851       952       972       A91
                            831       941       B32       981       B54
                            921       A31       B41       A71       B72
                                      B21       D21       B43       B81
                                                          B52       C71
                                                          B61       D43
                                                          C51       D52
                                                          D32       D61
                                                          D41       E51
                                                          E31       F41
                                                          F21       G31
                                                                    H21
(End)
		

Crossrefs

A023022 is the 2-part version.
A101271 is the relative prime instead of pairwise coprime version.
A220377*6 is the ordered version.
A305713 counts these partitions of any length, with Heinz numbers A302797.
A307719 is the non-strict version.
A337461 is the non-strict ordered version.
A337563 is the case with no 1's.
A337605 is the pairwise non-coprime instead of pairwise coprime version.
A001399(n-6) counts strict 3-part partitions, with Heinz numbers A007304.
A008284 counts partitions by sum and length, with strict case A008289.
A318717 counts pairwise non-coprime strict partitions.
A326675 ranks pairwise coprime sets.
A327516 counts pairwise coprime partitions.
A337601 counts 3-part partitions whose distinct parts are pairwise coprime.

Programs

  • Mathematica
    Table[Length@Select[ IntegerPartitions[ n, {3}], #[[1]] != #[[2]] != #[[3]] && GCD[#[[1]], #[[2]]] == 1 && GCD[#[[1]], #[[3]]] == 1 && GCD[#[[2]], #[[3]]] == 1 &], {n, 6, 100}]
    Table[Count[IntegerPartitions[n,{3}],?(CoprimeQ@@#&&Length[ Union[#]] == 3&)],{n,6,100}] (* _Harvey P. Dale, May 22 2020 *)
  • PARI
    a(n)=my(P=partitions(n));sum(i=1,#P,#P[i]==3&&P[i][1]Charles R Greathouse IV, Dec 14 2012

Formula

a(n > 2) = A307719(n) - 1. - Gus Wiseman, Oct 15 2020

A103221 Number of partitions of n into parts 2 and 3.

Original entry on oeis.org

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

Views

Author

Michael Somos, Jan 25 2005

Keywords

Comments

Essentially the same as A008615.
Poincaré series [or Poincare series] for modular forms of weight w for the full modular group. As generators one may take the Eisenstein series E_4 (A004009) and E_6 (A013973).
Dimension of the space of weight 2n+12 cusp forms for Gamma_0( 1 ).
Dimension of the space of weight 2n cuspidal newforms for Gamma_0( 5 ).
a(n) is the number of partitions of n into two nonnegative parts congruent modulo 3. - Andrew Baxter, Jun 28 2006
Also number of equivalence classes of period 2n billiards on an equilateral triangle. - Andrew Baxter, Jun 06 2008
a(n) is also the number of 2-regular multigraphs on n vertices, where each component is either a pair of parallel edges, or a triangle. - Jason Kimberley, Oct 14 2011
For n>1, a(n) is the number of partitions of 2n into positive parts x,y, and z such that x>=y and y=z. This sequence is used in calculating the probability of the need for a run-off election when n voters randomly cast ballots for two of three candidates running for two empty slots on a county commission. - Dennis P. Walsh, Apr 25 2013
Also, Molien series for invariants of finite Coxeter group A_2. 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. - N. J. A. Sloane, Jan 11 2016
The coefficient of x^(2*n+1) in the power series expansion of the Weierstrass sigma function is a polynomial in the invariants g2 and g3 with a(n) terms. - Michael Somos, Jun 14 2016
a(n) is also the dimension of the complex vector space of modular forms M_{2*n} of weight 2*n and level 1 (full modular group). See Apostol p. 119, eq. (9) for k=2*n, and Ash and Gross, p. 178, Table 13.1. For a(6*k+1) = a(6*k+j)-1 for j = 0,2,3,4,5 and k >= 0 see A016921 (so-called dips, cf. Ash and Gross, p. 178.). - Wolfdieter Lang, Sep 16 2016
In an hexagonal tiling of the plane where the base tile is (0,0)--(2,1)--(3,3)--(1,4)--(-1,3)--(-2,1)--(0,0), a(n) is the number of vertices on the (n,0)--(n,n) closed line segment. - Luc Rousseau, Mar 22 2018

Examples

			For n=8, a(n)=2 since there are two partitions of 16 into 3 positive parts x, y, and z such that x >= y and y=z, namely, 16 = 8+4+4 and 16 = 6+5+5. - _Dennis P. Walsh_, Apr 25 2013
G.f. = 1 + x^2 + x^3 + x^4 + x^5 + 2*x^6 + x^7 + 2*x^8 + 2*x^9 + 2*x^10 + 2*x^11 + ...
		

References

  • T. M. Apostol, Modular Functions and Dirichlet Series in Number Theory, Springer-Verlag, 1990, page 119.
  • Avner Ash and Robert Gross, Summing it up, Princeton University Press, 2016, p. 178.
  • D. J. Benson, Polynomial Invariants of Finite Groups, Cambridge, 1993, p. 100.
  • E. Freitag, Siegelsche Modulfunktionen, Springer-Verlag, Berlin, 1983; p. 141, Th. 1.1.
  • R. C. Gunning, Lectures on Modular Forms. Princeton Univ. Press, Princeton, NJ, 1962.
  • J. E. Humphreys, Reflection Groups and Coxeter Groups, Cambridge, 1990. See Table 3.1, page 59.
  • J.-M. Kantor, Ou en sont les mathématiques, La formule de Molien-Weyl, SMF, Vuibert, p. 79
  • S. Mukai, An Introduction to Invariants and Moduli, Cambridge, 2003; see p. 26. - N. J. A. Sloane, Aug 28 2010.

Crossrefs

Cf. A008615, A001399 (partial sums), A128115, A171386, A081753.
Molien series for finite Coxeter groups A_1 through A_12 are A059841, A103221, A266755, A008667, A037145, A001996, and A266776, A266777, A266778, A266779, A266780, A266781.

Programs

  • Magma
    [Floor((n+2)/2)-Floor((n+2)/3): n in [0..100]]; // Vincenzo Librandi, Sep 18 2016
    
  • Maple
    A103221:=n->floor((n+2)/2)-floor((n+2)/3): # Andrew Baxter, Jun 06 2008
  • Mathematica
    a=b=c=d=0;Table[e=a+b-d+1;a=b;b=c;c=d;d=e,{n,100}] (* Vladimir Joseph Stephan Orlovsky, Feb 26 2011 *)
    LinearRecurrence[{0, 1, 1, 0, -1},{1, 0, 1, 1, 1},88] (* Ray Chandler, Sep 23 2015 *)
    a[ n_] := With[{m = Max[-5 - n, n]}, (-1)^Boole[n < 0] SeriesCoefficient[ 1 / ((1 - x^2) (1 - x^3)), {x, 0, m}]]; (* Michael Somos, Jun 02 2019 *)
  • PARI
    {a(n) = if( n<-4, -a(-5-n), polcoeff( 1 / ((1 - x^2) * (1 - x^3)) + x * O(x^n), n))};
    
  • PARI
    a(n)=n+=2; n\2 - n\3 \\ Charles R Greathouse IV, Jul 31 2017
    
  • Python
    def A103221(n): return (n>>1)+1-(n+2)//3 # Chai Wah Wu, Apr 15 2025
  • Sage
    def a(n) : return( len( CuspForms( Gamma0( 1), 2*n + 12, prec=1). basis())); # Michael Somos, May 29 2013
    

Formula

Euler transform of finite sequence [0, 1, 1] with offset 1, which is A171386.
a(n) = A008615(n+2). First differences of A001399.
a(n) = a(n-6) + 1 = a(n-2) + a(n-3) - a(n-5). - Henry Bottomley, Sep 02 2000
G.f.: 1/((1-x^2)*(1-x^3)).
a(n) = floor((n+2)/2) - floor((n+2)/3). - Andrew Baxter, Jun 06 2008
For odd n, a(n)=floor((n+3)/6). For even n, a(n)=floor((n+6)/6). - Dennis P. Walsh, Apr 25 2013
a(n) = floor(n/6)+1 unless n == 1 (mod 6); if n == 1 (mod 6), a(n) = floor(n/6). - Bob Selcoe, Sep 27 2014
a(n) = A081753(2*n); see the Dennis P. Walsh formula. - Wolfdieter Lang, Sep 16 2016
a(n)-a(n-2) = A079978(n). - R. J. Mathar, Jun 23 2021
E.g.f.: (3*(4 + x)*cosh(x) + exp(-x/2)*(6*cos(sqrt(3)*x/2) - 2*sqrt(3)*sin(sqrt(3)*x/2)) + 3*(1 + x)*sinh(x))/18. - Stefano Spezia, Mar 05 2023
a(n) = A008615(n-1)+A059841(n). - R. J. Mathar, May 03 2023

Extensions

Name changed by Wolfdieter Lang, Sep 16 2016

A337461 Number of pairwise coprime ordered triples of positive integers summing to n.

Original entry on oeis.org

0, 0, 0, 1, 3, 3, 9, 3, 15, 9, 21, 9, 39, 9, 45, 21, 45, 21, 87, 21, 93, 39, 87, 39, 153, 39, 135, 63, 153, 57, 255, 51, 207, 93, 225, 93, 321, 81, 291, 135, 321, 105, 471, 105, 393, 183, 381, 147, 597, 147, 531, 213, 507, 183, 759, 207, 621, 273, 621, 231
Offset: 0

Views

Author

Gus Wiseman, Sep 02 2020

Keywords

Examples

			The a(3) = 1 through a(9) = 9 triples:
  (1,1,1)  (1,1,2)  (1,1,3)  (1,1,4)  (1,1,5)  (1,1,6)  (1,1,7)
           (1,2,1)  (1,3,1)  (1,2,3)  (1,5,1)  (1,2,5)  (1,3,5)
           (2,1,1)  (3,1,1)  (1,3,2)  (5,1,1)  (1,3,4)  (1,5,3)
                             (1,4,1)           (1,4,3)  (1,7,1)
                             (2,1,3)           (1,5,2)  (3,1,5)
                             (2,3,1)           (1,6,1)  (3,5,1)
                             (3,1,2)           (2,1,5)  (5,1,3)
                             (3,2,1)           (2,5,1)  (5,3,1)
                             (4,1,1)           (3,1,4)  (7,1,1)
                                               (3,4,1)
                                               (4,1,3)
                                               (4,3,1)
                                               (5,1,2)
                                               (5,2,1)
                                               (6,1,1)
		

Crossrefs

A000212 counts the unimodal instead of coprime version.
A220377*6 is the strict case.
A307719 is the unordered version.
A337462 counts these compositions of any length.
A337563 counts the case of partitions with no 1's.
A337603 only requires the *distinct* parts to be pairwise coprime.
A337604 is the intersecting instead of coprime version.
A014612 ranks 3-part partitions.
A302696 ranks pairwise coprime partitions.
A327516 counts pairwise coprime partitions.
A333228 ranks compositions whose distinct parts are pairwise coprime.

Programs

  • Mathematica
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n,{3}],CoprimeQ@@#&]],{n,0,30}]

A307719 Number of partitions of n into 3 mutually coprime parts.

Original entry on oeis.org

0, 0, 0, 1, 1, 1, 2, 1, 3, 2, 4, 2, 7, 2, 8, 4, 8, 4, 15, 4, 16, 7, 15, 7, 26, 7, 23, 11, 26, 10, 43, 9, 35, 16, 38, 16, 54, 14, 49, 23, 54, 18, 79, 18, 66, 31, 64, 25, 100, 25, 89, 36, 85, 31, 127, 35, 104, 46, 104, 39, 167, 36, 125, 58, 129, 52, 185, 45
Offset: 0

Views

Author

Wesley Ivan Hurt, Apr 24 2019

Keywords

Comments

The Heinz numbers of these partitions are the intersection of A014612 (triples) and A302696 (pairwise coprime). - Gus Wiseman, Oct 16 2020

Examples

			There are 2 partitions of 9 into 3 mutually coprime parts: 7+1+1 = 5+3+1, so a(9) = 2.
There are 4 partitions of 10 into 3 mutually coprime parts: 8+1+1 = 7+2+1 = 5+4+1 = 5+3+2, so a(10) = 4.
There are 2 partitions of 11 into 3 mutually coprime parts: 9+1+1 = 7+3+1, so a(11) = 2.
There are 7 partitions of 12 into 3 mutually coprime parts: 10+1+1 = 9+2+1 = 8+3+1 = 7+4+1 = 6+5+1 = 7+3+2 = 5+4+3, so a(12) = 7.
		

Crossrefs

A023022 is the version for pairs.
A220377 is the strict case, with ordered version A220377*6.
A327516 counts these partitions of any length, with strict version A305713 and Heinz numbers A302696.
A337461 is the ordered version.
A337563 is the case with no 1's.
A337599 is the pairwise non-coprime instead of pairwise coprime version.
A337601 only requires the distinct parts to be pairwise coprime.
A001399(n-3) = A069905(n) = A211540(n+2) counts 3-part partitions.
A002865 counts partitions with no 1's, with strict case A025147.
A007359 and A337485 count pairwise coprime partitions with no 1's.
A200976 and A328673 count pairwise non-coprime partitions.

Programs

  • Maple
    N:= 200: # to get a(0)..a(N)
    A:= Array(0..N):
    for a from 1 to N/3 do
      for b from a to (N-a)/2 do
        if igcd(a,b) > 1 then next fi;
        ab:= a*b;
        for c from b to N-a-b do
           if igcd(ab,c)=1 then A[a+b+c]:= A[a+b+c]+1 fi
    od od od:
    convert(A,list); # Robert Israel, May 09 2019
  • Mathematica
    Table[Sum[Sum[Floor[1/(GCD[i, j] GCD[j, n - i - j] GCD[i, n - i - j])], {i, j, Floor[(n - j)/2]}], {j, Floor[n/3]}], {n, 0, 100}]
    Table[Length[Select[IntegerPartitions[n,{3}],CoprimeQ@@#&]],{n,0,100}] (* Gus Wiseman, Oct 15 2020 *)

Formula

a(n) = Sum_{j=1..floor(n/3)} Sum_{i=j..floor((n-j)/2)} [gcd(i,j) * gcd(j,n-i-j) * gcd(i,n-i-j) = 1], where [] is the Iverson bracket.
a(n > 2) = A220377(n) + 1. - Gus Wiseman, Oct 15 2020

A049451 Twice second pentagonal numbers.

Original entry on oeis.org

0, 4, 14, 30, 52, 80, 114, 154, 200, 252, 310, 374, 444, 520, 602, 690, 784, 884, 990, 1102, 1220, 1344, 1474, 1610, 1752, 1900, 2054, 2214, 2380, 2552, 2730, 2914, 3104, 3300, 3502, 3710, 3924, 4144, 4370, 4602, 4840, 5084, 5334, 5590, 5852, 6120, 6394, 6674, 6960, 7252, 7550, 7854
Offset: 0

Views

Author

Joe Keane (jgk(AT)jgk.org)

Keywords

Comments

From Floor van Lamoen, Jul 21 2001: (Start)
Write 1,2,3,4,... in a hexagonal spiral around 0, then a(n) is the sequence found by reading the line from 0 in the direction 0,4,... . The spiral begins:
.
52
. \
33--32--31--30 51
/ . \ \
34 16--15--14 29 50
/ / . \ \ \
35 17 5---4 13 28 49
/ / / . \ \ \ \
36 18 6 0 3 12 27 48
/ / / / / / / /
37 19 7 1---2 11 26 47
\ \ \ / / /
38 20 8---9--10 25 46
\ \ / /
39 21--22--23--24 45
\ /
40--41--42--43--44
(End)
Number of edges in the join of the complete bipartite graph of order 2n and the cycle graph of order n, K_n,n * C_n. - Roberto E. Martinez II, Jan 07 2002
The average of the first n elements starting from a(1) is equal to (n+1)^2. - Mario Catalani (mario.catalani(AT)unito.it), Apr 10 2003
If Y is a 4-subset of an n-set X then, for n >= 4, a(n-4) is the number of (n-4)-subsets of X having either one element or two elements in common with Y. - Milan Janjic, Dec 28 2007
With offset 1: the maximum possible sum of numbers in an N x N standard Minesweeper grid. - Dmitry Kamenetsky, Dec 14 2008
a(n) = A001399(6*n-2), number of partitions of 6*n-2 into parts < 4. For example a(2)=14 where the partitions of 6*2-2=10 into parts < 4 are [1,1,1,1,1,1,1,1,1,1], [1,1,1,1,1,1,1,1,2], [1,1,1,1,1,1,1,3], [1,1,1,1,1,1,2,2], [1,1,1,1,1,2,3], [1,1,1,1,2,2,2], [1,1,1,1,3,3], [1,1,1,2,2,3], [1,1,2,2,2,2], [1,1,2,3,3], [1,2,2,2,3], [2,2,2,2,2], [1,3,3,3], [2,2,3,3]. - Adi Dani, Jun 07 2011
A003056 is the following array A read by antidiagonals:
0, 1, 2, 3, 4, 5, ...
1, 2, 3, 4, 5, 6, ...
2, 3, 4, 5, 6, 7, ...
3, 4, 5, 6, 7, 8, ...
4, 5, 6, 7, 8, 9, ...
5, 6, 7, 8, 9, 10, ...
and a(n) is the hook sum Sum_{k=0..n} A(n,k) + Sum_{r=0..n-1} A(r,n). - R. J. Mathar, Jun 30 2013
a(n)*Pi is the total length of 3 points circle center spiral after n rotations. The spiral length at each rotation (L(n)) is A016957. The spiral length ratio rounded down [floor(L(n)/L(1))] is A001651. See illustration in links. - Kival Ngaokrajang, Dec 27 2013
Partial sums give A114364. - Leo Tavares, Feb 25 2022
For n >= 1, the continued fraction expansion of sqrt(27*a(n)) is [9n+1; {2, 2n-1, 1, 4, 1, 2n-1, 2, 18n+2}]. - Magus K. Chu, Oct 13 2022

Examples

			From _Dmitry Kamenetsky_, Dec 14 2008, with slight rewording by Raymond Martineau (mart0258(AT)yahoo.com), Dec 16 2008: (Start)
For an N x N Minesweeper grid the highest sum of numbers is (N-1)(3*N-2). This is achieved by filling every second row with mines (shown as 'X'). For example, when N=5 the best grids are:
.
  X X X X X
  4 6 6 6 4
  X X X X X
  4 6 6 6 4
  X X X X X
.
  and
.
  2 3 3 3 2
  X X X X X
  4 6 6 6 4
  X X X X X
  2 3 3 3 2
.
each giving a total of 52. (End)
		

References

  • L. B. W. Jolley, Summation of Series, Dover Publications, 1961, p. 12.

Crossrefs

Similar sequences are listed in A316466.

Programs

Formula

a(n) = n*(3*n+1).
G.f.: 2*x*(2+x)/(1-x)^3.
Sum_{i=1..n} a(i) = A045991(n+1). - Gary W. Adamson, Dec 20 2006
a(n) = 2*A005449(n). - Omar E. Pol, Dec 18 2008
a(n) = a(n-1) + 6*n -2, n > 0. - Vincenzo Librandi, Aug 06 2010
a(n) = A100104(n+1) - A100104(n). - Reinhard Zumkeller, Jul 07 2012
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) with a(0) = 0, a(1) = 4, a(2) = 14. - Philippe Deléham, Mar 26 2013
a(n) = A174709(6*n+3). - Philippe Deléham, Mar 26 2013
a(n) = (24/(n+2)!)*Sum_{j=0..n} (-1)^(n-j)*binomial(n,j)*j^(n+2). - Bruno Berselli, Jun 04 2013 - after the similar formula of Vladimir Kruchinin in A002411
a(n) = A002061(n+1) + A056220(n). - Bruce J. Nicholson, Sep 21 2017
a(n) = Sum_{i = 2..5} P(i,n), where P(i,m) = m*((i-2)*m-(i-4))/2. - Bruno Berselli, Jul 04 2018
E.g.f.: x*(4 + 3*x)*exp(x). - G. C. Greubel, Sep 01 2019
a(n) = A003215(n) - A005408(n). - Leo Tavares, Feb 25 2022
From Amiram Eldar, Feb 27 2022: (Start)
Sum_{n>=1} 1/a(n) = 3 - Pi/(2*sqrt(3)) - 3*log(3)/2.
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/sqrt(3) + 2*log(2) - 3. (End)
a(n) = A001105(n) + A002378(n). - Torlach Rush, Jul 11 2022

A007997 a(n) = ceiling((n-3)(n-4)/6).

Original entry on oeis.org

0, 0, 1, 1, 2, 4, 5, 7, 10, 12, 15, 19, 22, 26, 31, 35, 40, 46, 51, 57, 64, 70, 77, 85, 92, 100, 109, 117, 126, 136, 145, 155, 166, 176, 187, 199, 210, 222, 235, 247, 260, 274, 287, 301, 316, 330, 345, 361, 376, 392, 409, 425, 442, 460, 477, 495, 514, 532, 551, 571, 590, 610
Offset: 3

Views

Author

Keywords

Comments

Number of solutions to x+y+z=0 (mod m) with 0<=x<=y<=z
Nonorientable genus of complete graph on n nodes.
Also (with different offset) Molien series for alternating group A_3.
(1+x^3 ) / ((1-x)*(1-x^2)*(1-x^3)) is the Poincaré series [or Poincare series] (or Molien series) for H^*(S_6, F_2).
a(n+5) is the number of necklaces with 3 black beads and n white beads.
The g.f./x^5 is Z(C_3,x), the 3-variate cycle index polynomial for the cyclic group C_3, with substitution x[i]->1/(1-x^i), i=1,2,3. Therefore by Polya enumeration a(n+5) is the number of cyclically inequivalent 3-necklaces whose 3 beads are labeled with nonnegative integers such that the sum of labels is n, for n=0,1,2,... . See A102190 for Z(C_3,x). - Wolfdieter Lang, Feb 15 2005
a(n+1) is the number of pairs (x,y) with x and y in {0,...,n}, x = (y mod 3), and x+y < n. - Clark Kimberling, Jul 02 2012
From Gus Wiseman, Oct 17 2020: (Start)
Also the number of 3-part integer compositions of n - 2 that are either weakly increasing or strictly decreasing. For example, the a(5) = 1 through a(13) = 15 compositions are:
(111) (112) (113) (114) (115) (116) (117) (118) (119)
(122) (123) (124) (125) (126) (127) (128)
(222) (133) (134) (135) (136) (137)
(321) (223) (224) (144) (145) (146)
(421) (233) (225) (226) (155)
(431) (234) (235) (227)
(521) (333) (244) (236)
(432) (334) (245)
(531) (532) (335)
(621) (541) (344)
(631) (542)
(721) (632)
(641)
(731)
(821)
(End)

Examples

			For m=7 (n=12), the 12 solutions are xyz = 000 610 520 511 430 421 331 322 662 653 644 554.
		

References

  • A. Adem and R. J. Milgram, Cohomology of Finite Groups, Springer-Verlag, 2nd. ed., 2004, p. 204.
  • D. J. Benson, Polynomial Invariants of Finite Groups, Cambridge, 1993, p. 105.
  • J. L. Gross and T. W. Tucker, Topological Graph Theory, Wiley, 1987; see \bar{I}(n) p. 221.
  • J. L. Gross and J. Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 740.
  • E. V. McLaughlin, Numbers of factorizations in non-unique factorial domains, Senior Thesis, Allegeny College, Meadville, PA, 2004.

Crossrefs

Apart from initial term, same as A058212.
A001399(n-6)*2 = A069905(n-3)*2 = A211540(n-1)*2 counts the strict case.
A014311 intersected with A225620 U A333256 ranks these compositions.
A218004 counts these compositions of any length.
A000009 counts strictly decreasing compositions.
A000041 counts weakly increasing compositions.
A001523 counts unimodal compositions, with complement counted by A115981.
A007318 and A097805 count compositions by length.
A032020 counts strict compositions, ranked by A233564.
A333149 counts neither increasing nor decreasing strict compositions.

Programs

  • Haskell
    a007997 n = ceiling $ (fromIntegral $ (n - 3) * (n - 4)) / 6
    a007997_list = 0 : 0 : 1 : zipWith (+) a007997_list [1..]
    -- Reinhard Zumkeller, Dec 18 2013
    
  • Maple
    x^5*(1+x^3)/((1-x)*(1-x^2)*(1-x^3));
    seq(ceil(binomial(n,2)/3), n=0..63); # Zerinvary Lajos, Jan 12 2009
    a := n -> (n*(n-7)-2*([1,1,-1][n mod 3 +1]-7))/6;
    seq(a(n), n=3..64); # Peter Luschny, Jan 13 2015
  • Mathematica
    k = 3; Table[Apply[Plus, Map[EulerPhi[ # ]Binomial[n/#, k/# ] &, Divisors[GCD[n, k]]]]/n, {n, k, 30}] (* Robert A. Russell, Sep 27 2004 *)
    Table[Ceiling[((n-3)(n-4))/6],{n,3,100}] (* or *) LinearRecurrence[ {2,-1,1,-2,1},{0,0,1,1,2},100] (* Harvey P. Dale, Jan 21 2014 *)
  • PARI
    a(n)=(n^2-7*n+16)\6 \\ Charles R Greathouse IV, Sep 24 2015

Formula

a(n) = a(n-3) + n - 2, a(0)=0, a(1)=0, a(2)=1 [Offset 0]. - Paul Barry, Jul 14 2004
G.f.: x^5*(1+x^3)/((1-x)*(1-x^2)*(1-x^3)) = x^5*(1-x+x^2)/((1-x)^2*(1-x^3)).
a(n+5) = Sum_{k=0..floor(n/2)} C(n-k,L(k/3)), where L(j/p) is the Legendre symbol of j and p. - Paul Barry, Mar 16 2006
a(3)=0, a(4)=0, a(5)=1, a(6)=1, a(7)=2, a(n) = 2*a(n-1) - a(n-2) + a(n-3) - 2*a(n-4) + a(n-5). - Harvey P. Dale, Jan 21 2014
a(n) = (n^2 - 7*n + 14 - 2*(-1)^(2^(n + 1 - 3*floor((n+1)/3))))/6. - Luce ETIENNE, Dec 27 2014
a(n) = A001399(n-3) + A001399(n-6). Compare to A140106(n) = A001399(n-3) - A001399(n-6). - Gus Wiseman, Oct 17 2020
a(n) = (40 + 3*(n - 7)*n - 4*cos(2*n*Pi/3) - 4*sqrt(3)*sin(2*n*Pi/3))/18. - Stefano Spezia, Dec 14 2021
Sum_{n>=5} 1/a(n) = 6 - 2*Pi/sqrt(3) + 2*Pi*tanh(sqrt(5/3)*Pi/2)/sqrt(15). - Amiram Eldar, Oct 01 2022

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)

A052307 Triangle read by rows: T(n,k) = number of bracelets (reversible necklaces) with n beads, k of which are black and n - k are white.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 3, 3, 3, 1, 1, 1, 1, 3, 4, 4, 3, 1, 1, 1, 1, 4, 5, 8, 5, 4, 1, 1, 1, 1, 4, 7, 10, 10, 7, 4, 1, 1, 1, 1, 5, 8, 16, 16, 16, 8, 5, 1, 1, 1, 1, 5, 10, 20, 26, 26, 20, 10, 5, 1, 1, 1, 1, 6, 12, 29, 38, 50, 38, 29, 12, 6, 1, 1, 1, 1, 6, 14, 35, 57, 76, 76, 57, 35, 14, 6, 1, 1
Offset: 0

Author

Christian G. Bower, Nov 15 1999

Keywords

Comments

Equivalently, T(n,k) is the number of orbits of k-element subsets of the vertices of a regular n-gon under the usual action of the dihedral group D_n, or under the action of Euclidean plane isometries. Note that each row of the table is symmetric and unimodal. - Austin Shapiro, Apr 20 2009
Also, the number of k-chords in n-tone equal temperament, up to (musical) transposition and inversion. Example: there are 29 tetrachords, 38 pentachords, 50 hexachords in the familiar 12-tone equal temperament. Called "Forte set-classes," after Allen Forte who first cataloged them. - Jon Wild, May 21 2004

Examples

			Triangle T(n,k) (with rows n >= 0 and columns k = 0..n) begins:
   1;
   1,  1;
   1,  1,  1;
   1,  1,  1,  1;
   1,  1,  2,  1,  1;
   1,  1,  2,  2,  1,  1;
   1,  1,  3,  3,  3,  1,  1;
   1,  1,  3,  4,  4,  3,  1,  1;
   1,  1,  4,  5,  8,  5,  4,  1,  1;
   1,  1,  4,  7, 10, 10,  7,  4,  1,  1;
   1,  1,  5,  8, 16, 16, 16,  8,  5,  1,  1;
   1,  1,  5, 10, 20, 26, 26, 20, 10,  5,  1,  1;
   1,  1,  6, 12, 29, 38, 50, 38, 29, 12,  6,  1,  1;
   ...
		

References

  • Martin Gardner, "New Mathematical Diversions from Scientific American" (Simon and Schuster, New York, 1966), pages 245-246.
  • N. Zagaglia Salvi, Ordered partitions and colourings of cycles and necklaces, Bull. Inst. Combin. Appl., 27 (1999), 37-40.

Programs

  • Maple
    A052307 := proc(n,k)
            local hk,a,d;
            if k = 0 then
                    return 1 ;
            end if;
            hk := k mod 2 ;
            a := 0 ;
            for d in numtheory[divisors](igcd(k,n)) do
                    a := a+ numtheory[phi](d)*binomial(n/d-1,k/d-1) ;
            end do:
            %/k + binomial(floor((n-hk)/2),floor(k/2)) ;
            %/2 ;
    end proc: # R. J. Mathar, Sep 04 2011
  • Mathematica
    Table[If[m*n===0,1,1/2*If[EvenQ[n], If[EvenQ[m], Binomial[n/2, m/2], Binomial[(n-2)/2, (m-1)/2 ]], If[EvenQ[m], Binomial[(n-1)/2, m/2], Binomial[(n-1)/2, (m-1)/2]]] + 1/2*Fold[ #1 +(EulerPhi[ #2]*Binomial[n/#2, m/#2])/n &, 0, Intersection[Divisors[n], Divisors[m]]]], {n,0,12}, {m,0,n}] (* Wouter Meeussen, Aug 05 2002, Jan 19 2009 *)
  • PARI
    B(n,k)={ if(n==0, return(1)); GCD = gcd(n, k); S = 0;
    for(d = 1, GCD, if((k%d==0)&&(n%d==0), S+=eulerphi(d)*binomial(n/d,k/d)));
    return (binomial(floor(n/2)- k%2*(1-n%2), floor(k/2))/2 + S/(2*n)); }
    n=0;k=0; for(L=0, 8645, print(L, " ", B(n,k)); k++; if(k>n, k=0; n++))
    /* Washington Bomfim, Jun 30 2012 */
    
  • Python
    from sympy import binomial as C, totient, divisors, gcd
    def T(n, k): return 1 if n==0 else C((n//2) - k%2 * (1 - n%2), (k//2))/2 + sum(totient(d)*C(n//d, k//d) for d in divisors(gcd(n, k)))/(2*n)
    for n in range(11): print([T(n, k) for k in range(n + 1)]) # Indranil Ghosh, Apr 23 2017

Formula

T(0,0) = 1. If n > 0, T(n,k) = binomial(floor(n/2) - (k mod 2) * (1 - (n mod 2)), floor(k/2)) / 2 + Sum_{d|n, d|k} (phi(d)*binomial(n/d, k/d)) / (2*n). - Washington Bomfim, Jun 30 2012 [edited by Petros Hadjicostas, May 29 2019]
From Freddy Barrera, Apr 21 2019: (Start)
T(n,k) = (1/2) * (A119963(n,k) + A047996(n,k)).
T(n,k) = T(n, n-k) for each k < n (Theorem 2 of H. Gupta). (End)
G.f. for column k >= 1: (x^k/2) * ((1/k) * Sum_{m|k} phi(m)/(1 - x^m)^(k/m) + (1 + x)/(1 - x^2)^floor((k/2) + 1)). (This formula is due to Herbert Kociemba.) - Petros Hadjicostas, May 25 2019
Bivariate o.g.f.: Sum_{n,k >= 0} T(n, k)*x^n*y^k = (1/2) * ((x + 1) * (x*y + 1) / (1 - x^2 * (y^2 + 1)) + 1 - Sum_{d >= 1} (phi(d)/d) * log(1 - x^d * (1 + y^d))). - Petros Hadjicostas, Jun 13 2019

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
Previous Showing 21-30 of 198 results. Next