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-4 of 4 results.

A001935 Number of partitions with no even part repeated; partitions of n in which no parts are multiples of 4.

Original entry on oeis.org

1, 1, 2, 3, 4, 6, 9, 12, 16, 22, 29, 38, 50, 64, 82, 105, 132, 166, 208, 258, 320, 395, 484, 592, 722, 876, 1060, 1280, 1539, 1846, 2210, 2636, 3138, 3728, 4416, 5222, 6163, 7256, 8528, 10006, 11716, 13696, 15986, 18624, 21666, 25169, 29190, 33808, 39104, 45164
Offset: 0

Views

Author

Keywords

Comments

Also number of partitions of n where no part appears more than three times.
a(n) satisfies Euler's pentagonal number (A001318) theorem, unless n is in A062717 (see Fink et al.).
Also number of partitions of n in which the least part and the differences between consecutive parts is at most 3. Example: a(5)=6 because we have [4,1], [3,2], [3,1,1], [2,2,1], [2,1,1,1] and [1,1,1,1,1]. - Emeric Deutsch, Apr 19 2006
Equals A000009 convolved with its aerated variant, = polcoeff A000009 * A000041 * A010054 (with alternate signs). - Gary W. Adamson, Mar 16 2010
Equals left border of triangle A174715. - Gary W. Adamson, Mar 27 2010
The Cayley reference is actually to A083365. - Michael Somos, Feb 24 2011
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Convolution of A000009 and A035457. - Vaclav Kotesovec, Aug 23 2015
Convolution inverse is A082303. - Michael Somos, Sep 30 2017
The g.f. in the form Sum_{n >= 0} x^(n*(n+1)/2) * Product_{k = 1..n} (1+x^k)/(1-x^k) = Sum_{n >= 0} x^(n*(n+1)/2) * Product_{k = 1..n} (1+x^k)/(1+x^k-2*x^k) == Sum_{n >= 0} x^(n*(n+1)/2) (mod 2). It follows that a(n) is odd iff n = k*(k + 1)/2 for some nonnegative integer k. Cf. A333374. - Peter Bala, Jan 08 2025

Examples

			G.f. = 1 + x + 2*x^2 + 3*x^3 + 4*x^4 + 6*x^5 + 9*x^6 + 12*x^7 + 16*x^8 + 22*x^9 + ...
G.f. = q + q^9 + 2*q^17 + 3*q^25 + 4*q^33 + 6*q^41 + 9*q^49 + 12*q^57 + 16*q^65 + 22*q^73 + ...
a(5)=6 because we have [5], [4,1], [3,2], [3,1,1], [2,1,1,1] and [1,1,1,1,1].
		

References

  • A. Cayley, A memoir on the transformation of elliptic functions, Collected Mathematical Papers. Vols. 1-13, Cambridge Univ. Press, London, 1889-1897, Vol. 9, p. 128.
  • I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, Wiley, N.Y., 1983, (2.5.2).
  • M. D. Hirschhorn, The Power of q, Springer, 2017. See ped page 303ff.
  • R. Honsberger, Mathematical Gems III, M.A.A., 1985, p. 241.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A000041, A010054. - Gary W. Adamson, Mar 16 2010
Cf. A174715. - Gary W. Adamson, Mar 27 2010
Cf. A082303.
Number of r-regular partitions for r = 2 through 12: A000009, A000726, A001935, A035959, A219601, A035985, A261775, A104502, A261776, A328545, A328546.

Programs

  • Haskell
    a001935 = p a042968_list where
       p _          0 = 1
       p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Sep 02 2012
  • Maple
    g:=product((1+x^j)*(1+x^(2*j)),j=1..50): gser:=series(g,x=0,55): seq(coeff(gser,x,n),n=0..48); # Emeric Deutsch, Apr 19 2006
    # second Maple program:
    with(numtheory):
    a:= proc(n) option remember; `if`(n=0, 1, add(a(n-j)*add(
         `if`(irem(d, 4)=0, 0, d), d=divisors(j)), j=1..n)/n)
        end:
    seq(a(n), n=0..50);  # Alois P. Heinz, Nov 24 2015
  • Mathematica
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 2, 0, q] / EllipticTheta[ 2, Pi/4, q^(1/2)] / (16 q)^(1/8), {q, 0, n}]; (* Michael Somos, Jul 11 2011 *)
    a[ n_] := SeriesCoefficient[ Product[ 1 - x^k, {k, 4, n, 4}] / Product[ 1 - x^k, {k, n}], {x, 0, n}]; (* Michael Somos, Jul 08 2011 *)
    CoefficientList[Series[Product[1+x^j+x^(2j)+x^(3j), {j,1,48}], {x,0,48}],x] (* Jean-François Alcover, May 26 2011, after Jon Perry *)
    QP = QPochhammer; CoefficientList[QP[q^4]/QP[q] + O[q]^50, q] (* Jean-François Alcover, Nov 24 2015 *)
    a[0] = 1; a[n_] := a[n] = Sum[a[n-j] DivisorSum[j, If[Divisible[#, 4], 0, #]&], {j, 1, n}]/n; Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Feb 19 2016, after Alois P. Heinz *)
    Table[Count[IntegerPartitions@n, x_ /; ! MemberQ [Mod[x, 4], 0, 2] ], {n, 0, 49}] (* Robert Price, Jul 28 2020 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( eta(x^4 + x * O(x^n)) / eta(x + x * O(x^n)), n))};
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff( sum(k=0, (sqrtint( 8*n + 1) - 1)\2, prod(i=1, k, (1 + x^i) / (x^-i - 1), 1 + x * O(x^n))), n))}; /* Michael Somos, Jun 01 2004 */
    
  • PARI
    {a(n)=polcoeff(exp(sum(m=1, n+1, x^m/(1+(-x)^m+x*O(x^n))/m)),n)} \\ Paul D. Hanna, Jul 24 2013
    

Formula

Euler transform of period 4 sequence [ 1, 1, 1, 0, ...].
Expansion of q^(-1/8) * eta(q^4) / eta(q) in powers of q. - Michael Somos, Mar 19 2004
Expansion of psi(-x) / phi(-x) = psi(x) / phi(-x^2) = psi(x^2) / psi(-x) = chi(x) / chi(-x^2)^2 = 1 / (chi(x) * chi(-x)^2) = 1 / (chi(-x) * chi(-x^2)) = f(-x^4) / f(-x) in powers of x where phi(), psi(), chi(), f() are Ramanujan theta functions. - Michael Somos, Jul 08 2011
G.f.: Product(j>=1, 1 + x^j + x^(2*j) + x^(3*j)). - Jon Perry, Mar 30 2004
G.f.: Product_{k>=1} (1+x^k)^(2-k%2). - Jon Perry, May 05 2005
G.f.: Product_{k>0} (1 + x^(2*k)) / (1 - x^(2*k-1)) = 1 + Sum_{k>0}(Product_{i=1..k} (x^i + 1) / (x^-i - 1)).
G.f.: Sum_{n>=0} ( x^(n*(n+1)/2) * Product_{k=1..n} (1+x^k)/(1-x^k) ). - Joerg Arndt, Apr 07 2011
G.f.: P(x^4)/P(x) where P(x) = Product_{k>=1} 1-x^k. - Joerg Arndt, Jun 21 2011
A083365(n) = (-1)^n a(n). Convolution square is A001936. a(n) = A098491(n) + A098492(n). a(2*n) = A081055(n). a(2*n + 1) = A081056(n).
G.f.: (1+ 1/G(0))/2, where G(k) = 1 - x^(2*k+1) - x^(2*k+1)/(1 + x^(2*k+2) + x^(2*k+2)/G(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Jul 03 2013
G.f.: exp( Sum_{n>=1} (x^n/n) / (1 + (-x)^n) ). - Paul D. Hanna, Jul 24 2013
a(n) ~ Pi * BesselI(1, sqrt(8*n + 1)*Pi/4) / (2*sqrt(8*n + 1)) ~ exp(Pi*sqrt(n/2)) / (4 * (2*n)^(3/4)) * (1 + (Pi/(16*sqrt(2)) - 3/(4*Pi*sqrt(2))) / sqrt(n) + (Pi^2/1024 - 15/(64*Pi^2) - 15/128) / n). - Vaclav Kotesovec, Aug 23 2015, extended Jan 14 2017
a(n) = (1/n)*Sum_{k=1..n} A046897(k)*a(n-k), a(0) = 1. - Seiichi Manyama, Mar 25 2017
G.f. is a period 1 Fourier series which satisfies f(-1 / (256 t)) = 1/2 g(t) where q = exp(2 Pi i t) and g() is the g.f. for A082303. - Michael Somos, Sep 30 2017

Extensions

More terms from James Sellers

A035959 Number of partitions of n in which no parts are multiples of 5.

Original entry on oeis.org

1, 1, 2, 3, 5, 6, 10, 13, 19, 25, 34, 44, 60, 76, 100, 127, 164, 205, 262, 325, 409, 505, 628, 769, 950, 1156, 1414, 1713, 2081, 2505, 3026, 3625, 4352, 5192, 6200, 7364, 8756, 10357, 12258, 14450, 17034, 20006, 23500, 27510, 32200, 37582, 43846, 51022
Offset: 0

Views

Author

Keywords

Comments

Also number of partitions with at most 4 parts of size 1 and differences between parts at distance 6 are greater than 1.
Also number of partitions of n where no part appears more than four times.
Case k=7, i=5 of Gordon Theorem.

Examples

			G.f. = 1 + x + 2*x^2 + 3*x^3 + 5*x^4 + 6*x^5 + 10*x^6 + 13*x^7 + 19*x^8 + ...
G.f. = q + q^7 + 2*q^13 + 3*q^19 + 5*q^25 + 6*q^31 + 10*q^37 + 13*q^43 + ...
a(6) counts these partitions: 6, 42, 411, 33, 321, 3111, 2211, 21111, 111111. - _Clark Kimberling_, Mar 09 2014
		

References

  • G. E. Andrews, The Theory of Partitions, Addison-Wesley, 1976, p. 109.

Crossrefs

Cf. A000009 (m=2), A000726 (m=3), A001935 (m=4), A219601 (m=6), A035985 (m=7), A261775 (m=8), A104502 (m=9), A261776 (m=10).
Number of r-regular partitions for r = 2 through 12: A000009, A000726, A001935, A035959, A219601, A035985, A261775, A104502, A261776, A328545, A328546.

Programs

  • Haskell
    a035959 = p a047201_list where
       p _      0 = 1
       p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Dec 17 2011
  • Mathematica
    max = 47; f[x_] := (x^5-1)/(x-1); g[x_] := Product[f[x^k], {k, 1, max}]; CoefficientList[ Series[g[x], {x, 0, max}], x] (* Jean-François Alcover, Nov 29 2011, after Michael Somos *)
    t = Flatten[Table[5 n + r, {n, 0, 60}, {r, 1, 4}]]; p[n_] := IntegerPartitions[n, All, t]; Table[p[n], {n, 0, 8}] (* shows partitions *)
    a[n_] := Length@p@n; a /@ Range[0, 50] (* Clark Kimberling, Mar 09 2014 *)
    nmax = 50; CoefficientList[Series[Product[(1 - x^(5*k))/(1 - x^k), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 31 2015 *)
    QP = QPochhammer; s = QP[q^5]/QP[q] + O[q]^50; CoefficientList[s, q] (* Jean-François Alcover, Nov 25 2015, after Michael Somos *)
    Table[Count[IntegerPartitions@n, x_ /; ! MemberQ [Mod[x, 5], 0, 2] ], {n, 0, 47}] (* Robert Price, Jul 28 2020 *)
    Table[Count[IntegerPartitions[n],?(NoneTrue[Mod[#,5]==0&])],{n,0,50}] (* _Harvey P. Dale, Dec 25 2021 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( eta(x^5 + x * O(x^n)) / eta(x + x * O(x^n)), n))}; /* Michael Somos, May 28 2006 */
    

Formula

G.f.: Product_{j>=1} (1 + x^j + x^2j + x^3j + x^4j). - Jon Perry, Mar 30 2004
G.f.: Product_{n>0, n==1, 2, 3, 4 mod 5} 1/(1-q^n).
Given g.f. A(x) then B(x) = x * A(x^3)^2 satisfies 0 = f(B(x), B(x^2)) where f(u, v) = u^3 + v^3 - u*v - 5*u^2*v^2. - Michael Somos, May 28 2006
Given g.f. A(x) then B(x) = x * A(x^3)^2 satisfies 0 = f(B(x), B(x^2), B(x^4)) where f(u, v, w) = v + 5*v^2*(u + w) - (u^2 + u*w + w^2). - Michael Somos, May 28 2006
Euler transform of period 5 sequence [ 1, 1, 1, 1, 0, ...]. - Michael Somos, May 28 2006
G.f.: Product_{k > 0} P5(x^k) where P5 is 5th cyclotomic polynomial.
Convolution inverse is A145466. - Michael Somos, Jun 26 2014
a(n) ~ 2*Pi * BesselI(1, 2*sqrt((6*n + 1)/5) * Pi/3) / (5*sqrt(6*n + 1)) ~ exp(2*Pi*sqrt(2*n/15)) / (3^(1/4) * 10^(3/4) * n^(3/4)) * (1 + (Pi/(3*sqrt(15)) - 3*sqrt(15)/(16*Pi)) / sqrt(2*n) + (Pi^2/540 - 225/(1024*Pi^2) - 5/32) / n). - Vaclav Kotesovec, Aug 31 2015, extended Jan 14 2017
a(n) = (1/n)*Sum_{k=1..n} A116073(k)*a(n-k), a(0) = 1. - Seiichi Manyama, Mar 25 2017
G.f.: exp(Sum_{k>=1} x^k*(1 + x^k + x^(2*k) + x^(3*k))/(k*(1 - x^(5*k)))). - Ilya Gutkovskiy, Aug 15 2018

A286653 Square array A(n,k), n>=0, k>=1, read by antidiagonals, where column k is the expansion of Product_{j>=1} (1 - x^(k*j))/(1 - x^j).

Original entry on oeis.org

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

Views

Author

Ilya Gutkovskiy, May 11 2017

Keywords

Comments

A(n,k) is the number of partitions of n in which no parts are multiples of k.
A(n,k) is also the number of partitions of n into at most k-1 copies of each part.

Examples

			Square array begins:
  1,  1,  1,  1,  1,  1,  ...
  0,  1,  1,  1,  1,  1,  ...
  0,  1,  2,  2,  2,  2,  ...
  0,  2,  2,  3,  3,  3,  ...
  0,  2,  4,  4,  5,  5,  ...
  0,  3,  5,  6,  6,  7,  ...
		

Crossrefs

Main diagonal gives A000041.
Mirror of A061198.

Programs

  • Maple
    b:= proc(n, i, k) option remember; `if`(n=0, [1, 0], `if`(k*i*(i+1)/2[0, l[1]*j]+l)(b(n-i*j, i-1, k)), j=0..min(n/i, k))))
        end:
    A:= (n, k)-> b(n$2, k-1)[1]:
    seq(seq(A(n, 1+d-n), n=0..d), d=0..16);  # Alois P. Heinz, Oct 17 2018
  • Mathematica
    Table[Function[k, SeriesCoefficient[Product[(1 - x^(i k))/(1 - x^i), {i, Infinity}], {x, 0, n}]][j - n + 1], {j, 0, 12}, {n, 0, j}] // Flatten
    Table[Function[k, SeriesCoefficient[QPochhammer[x^k, x^k]/QPochhammer[x, x], {x, 0, n}]][j - n + 1], {j, 0, 12}, {n, 0, j}] // Flatten

Formula

G.f. of column k: Product_{j>=1} (1 - x^(k*j))/(1 - x^j).

A061199 Upper right triangle read by columns where T(n,k), with k >= n, is the number of partitions of k where no part appears more than n times; also partitions of k where no parts are multiples of (n+1).

Original entry on oeis.org

1, 0, 1, 0, 1, 2, 0, 2, 2, 3, 0, 2, 4, 4, 5, 0, 3, 5, 6, 6, 7, 0, 4, 7, 9, 10, 10, 11, 0, 5, 9, 12, 13, 14, 14, 15, 0, 6, 13, 16, 19, 20, 21, 21, 22, 0, 8, 16, 22, 25, 27, 28, 29, 29, 30, 0, 10, 22, 29, 34, 37, 39, 40, 41, 41, 42, 0, 12, 27, 38, 44, 49, 51, 53, 54, 55, 55, 56, 0, 15, 36
Offset: 0

Views

Author

Henry Bottomley, Apr 20 2001

Keywords

Examples

			T(2,4) = 4 since the possible partitions of 4 are on the first definition (no term more than twice) 1+1+2, 2+2, 1+3, or 4 and on the second definition (no term a multiple of 3) 1+1+1+1, 1+1+2, 2+2, or 4.
Triangle T(n,k) begins:
1, 0, 0, 0, 0, 0,  0,  0,  0,  0, ...
   1, 1, 2, 2, 3,  4,  5,  6,  8, ...
      2, 2, 4, 5,  7,  9, 13, 16, ...
         3, 4, 6,  9, 12, 16, 22, ...
            5, 6, 10, 13, 19, 25, ...
               7, 10, 14, 20, 27, ...
                  11, 14, 21, 28, ...
                      15, 21, 29, ...
                          22, 29, ...
                              30, ...
		

Crossrefs

Rows effectively include A000007, A000009, A000726, A001935, A035959.
Main diagonal is A000041.
A061198 is the same table but includes cases where n>k.
T(n,2*n) gives: A232623.

Programs

  • Maple
    b:= proc(n, i, k) option remember;
          `if`(n=0, 1, `if`(i<1, 0,
          add(b(n-i*j, i-1, k), j=0..min(n/i, k))))
        end:
    T:= (n, k)-> b(k$2, n):
    seq(seq(T(n, k), n=0..k), k=0..12);  # Alois P. Heinz, Nov 27 2013
  • Mathematica
    b[n_, i_, k_] := b[n, i, k] = If[n == 0, 1, If[i<1, 0, Sum[b[n - i*j, i-1, k], {j, 0, Min[n/i, k]}]]]; T[n_, k_] := b[k, k, n]; Table[Table[T[n, k], {n, 0, k}], {k, 0, 12}] // Flatten (* Jean-François Alcover, Jan 28 2015, after Alois P. Heinz *)
Showing 1-4 of 4 results.