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

A000009 Expansion of Product_{m >= 1} (1 + x^m); number of partitions of n into distinct parts; number of partitions of n into odd parts.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 4, 5, 6, 8, 10, 12, 15, 18, 22, 27, 32, 38, 46, 54, 64, 76, 89, 104, 122, 142, 165, 192, 222, 256, 296, 340, 390, 448, 512, 585, 668, 760, 864, 982, 1113, 1260, 1426, 1610, 1816, 2048, 2304, 2590, 2910, 3264, 3658, 4097, 4582, 5120, 5718, 6378
Offset: 0

Views

Author

Keywords

Comments

Partitions into distinct parts are sometimes called "strict partitions".
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
The result that number of partitions of n into distinct parts = number of partitions of n into odd parts is due to Euler.
Bijection: given n = L1* 1 + L2*3 + L3*5 + L7*7 + ..., a partition into odd parts, write each Li in binary, Li = 2^a1 + 2^a2 + 2^a3 + ... where the aj's are all different, then expand n = (2^a1 * 1 + ...)*1 + ... by removing the brackets and we get a partition into distinct parts. For the reverse operation, just keep splitting any even number into halves until no evens remain.
Euler transform of period 2 sequence [1,0,1,0,...]. - Michael Somos, Dec 16 2002
Number of different partial sums 1+[1,2]+[1,3]+[1,4]+..., where [1,x] indicates a choice. E.g., a(6)=4, as we can write 1+1+1+1+1+1, 1+2+3, 1+2+1+1+1, 1+1+3+1. - Jon Perry, Dec 31 2003
a(n) is the sum of the number of partitions of x_j into at most j parts, where j is the index for the j-th triangular number and n-T(j)=x_j. For example; a(12)=partitions into <= 4 parts of 12-T(4)=2 + partitions into <= 3 parts of 12-T(3)=6 + partitions into <= 2 parts of 12-T(2)=9 + partitions into 1 part of 12-T(1)=11 = (2)(11) + (6)(51)(42)(411)(33)(321)(222) + (9)(81)(72)(63)(54)+(11) = 2+7+5+1 = 15. - Jon Perry, Jan 13 2004
Number of partitions of n where if k is the largest part, all parts 1..k are present. - Jon Perry, Sep 21 2005
Jack Grahl and Franklin T. Adams-Watters prove this claim of Jon Perry's by observing that the Ferrers dual of a "gapless" partition is guaranteed to have distinct parts; since the Ferrers dual is an involution, this establishes a bijection between the two sets of partitions. - Allan C. Wechsler, Sep 28 2021
The number of connected threshold graphs having n edges. - Michael D. Barrus (mbarrus2(AT)uiuc.edu), Jul 12 2007
Starting with offset 1 = row sums of triangle A146061 and the INVERT transform of A000700 starting: (1, 0, 1, -1, 1, -1, 1, -2, 2, -2, 2, -3, 3, -3, 4, -5, ...). - Gary W. Adamson, Oct 26 2008
Number of partitions of n in which the largest part occurs an odd number of times and all other parts occur an even number of times. (Such partitions are the duals of the partitions with odd parts.) - David Wasserman, Mar 04 2009
Equals A035363 convolved with A010054. The convolution square of A000009 = A022567 = A000041 convolved with A010054. A000041 = A000009 convolved with A035363. - Gary W. Adamson, Jun 11 2009
Considering all partitions of n into distinct parts: there are A140207(n) partitions of maximal size which is A003056(n), and A051162(n) is the greatest number occurring in these partitions. - Reinhard Zumkeller, Jun 13 2009
Equals left border of triangle A091602 starting with offset 1. - Gary W. Adamson, Mar 13 2010
Number of symmetric unimodal compositions of n+1 where the maximal part appears once. Also number of symmetric unimodal compositions of n where the maximal part appears an odd number of times. - Joerg Arndt, Jun 11 2013
Because for these partitions the exponents of the parts 1, 2, ... are either 0 or 1 (j^0 meaning that part j is absent) one could call these partitions also 'fermionic partitions'. The parts are the levels, that is the positive integers, and the occupation number is either 0 or 1 (like Pauli's exclusion principle). The 'fermionic states' are denoted by these partitions of n. - Wolfdieter Lang, May 14 2014
The set of partitions containing only odd parts forms a monoid under the product described in comments to A047993. - Richard Locke Peterson, Aug 16 2018
Ewell (1973) gives a number of recurrences. - N. J. A. Sloane, Jan 14 2020
a(n) equals the number of permutations p of the set {1,2,...,n+1}, written in one line notation as p = p_1p_2...p_(n+1), satisfying p_(i+1) - p_i <= 1 for 1 <= i <= n, (i.e., those permutations that, when read from left to right, never increase by more than 1) whose major index maj(p) := Sum_{p_i > p_(i+1)} i equals n. For example, of the 16 permutations on 5 letters satisfying p_(i+1) - p_i <= 1, 1 <= i <= 4, there are exactly two permutations whose major index is 4, namely, 5 3 4 1 2 and 2 3 4 5 1. Hence a(4) = 2. See the Bala link in A007318 for a proof. - Peter Bala, Mar 30 2022
Conjecture: Each positive integer n can be written as a_1 + ... + a_k, where a_1,...,a_k are strict partition numbers (i.e., terms of the current sequence) with no one dividing another. This has been verified for n = 1..1350. - Zhi-Wei Sun, Apr 14 2023
Conjecture: For each integer n > 7, a(n) divides none of p(n), p(n) - 1 and p(n) + 1, where p(n) is the number of partitions of n given by A000041. This has been verified for n up to 10^5. - Zhi-Wei Sun, May 20 2023 [Verified for n <= 2*10^6. - Vaclav Kotesovec, May 23 2023]
The g.f. Product_{k >= 0} 1 + x^k = Product_{k >= 0} 1 - x^k + 2*x^k == Product_{k >= 0} 1 - x^k == Sum_{k in Z} (-1)^k*x^(k*(3*k-1)/2) (mod 2) by Euler's pentagonal number theorem. It follows that a(n) is odd iff n = k*(3*k - 1)/2 for some integer k, i.e., iff n is a generalized pentagonal number A001318. - Peter Bala, Jan 07 2025

Examples

			G.f. = 1 + x + x^2 + 2*x^3 + 2*x^4 + 3*x^5 + 4*x^6 + 5*x^7 + 6*x^8 + 8*x^9 + ...
G.f. = q + q^25 + q^49 + 2*q^73 + 2*q^97 + 3*q^121 + 4*q^145 + 5*q^169 + ...
The partitions of n into distinct parts (see A118457) for small n are:
  1: 1
  2: 2
  3: 3, 21
  4: 4, 31
  5: 5, 41, 32
  6: 6, 51, 42, 321
  7: 7, 61, 52, 43, 421
  8: 8, 71, 62, 53, 521, 431
  ...
From _Reinhard Zumkeller_, Jun 13 2009: (Start)
a(8)=6, A140207(8)=#{5+2+1,4+3+1}=2, A003056(8)=3, A051162(8)=5;
a(9)=8, A140207(9)=#{6+2+1,5+3+1,4+3+2}=3, A003056(9)=3, A051162(9)=6;
a(10)=10, A140207(10)=#{4+3+2+1}=1, A003056(10)=4, A051162(10)=4. (End)
		

References

  • Mohammad K. Azarian, A Generalization of the Climbing Stairs Problem, Mathematics and Computer Education, Vol. 31, No. 1, pp. 24-28, Winter 1997. MathEduc Database (Zentralblatt MATH, 1997c.01891).
  • Mohammad K. Azarian, A Generalization of the Climbing Stairs Problem II, Missouri Journal of Mathematical Sciences, Vol. 16, No. 1, Winter 2004, pp. 12-17. Zentralblatt MATH, Zbl 1071.05501.
  • George E. Andrews, The Theory of Partitions, Cambridge University Press, 1998, p. 19.
  • George E. Andrews, Number Theory, Dover Publications, 1994, Theorem 12-3, pp. 154-5, and (13-1-1) p. 163.
  • Raymond Ayoub, An Introduction to the Analytic Theory of Numbers, Amer. Math. Soc., 1963; see p. 196.
  • T. J. I'a. Bromwich, Introduction to the Theory of Infinite Series, Macmillan, 2nd. ed. 1949, p. 116, Problem 18.
  • Louis Comtet, Advanced Combinatorics, Reidel, 1974, p. 99.
  • William Dunham, The Mathematical Universe, pp. 57-62, J. Wiley, 1994.
  • Leonhard Euler, De partitione numerorum, Novi commentarii academiae scientiarum Petropolitanae 3 (1750/1), 1753, reprinted in: Commentationes Arithmeticae. (Opera Omnia. Series Prima: Opera Mathematica, Volumen Secundum), 1915, Lipsiae et Berolini, 254-294.
  • Ian P. Goulden and David M. Jackson, Combinatorial Enumeration, Wiley, N.Y., 1983, (2.5.1).
  • G. H. Hardy, Ramanujan: twelve lectures on subjects suggested by his life and work, Cambridge, University Press, 1940, p. 86.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 277, Theorems 344, 346.
  • Carlos J. Moreno and Samuel S. Wagstaff, Jr., Sums of Squares of Integers, Chapman and Hall, 2006, p. 253.
  • Srinivasa Ramanujan, Collected Papers, Ed. G. H. Hardy et al., Cambridge 1927; Chelsea, NY, 1962. See Table V on page 309.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 288-290.

Crossrefs

Apart from the first term, equals A052839-1. The rows of A053632 converge to this sequence. When reduced modulo 2 equals the absolute values of A010815. The positions of odd terms given by A001318.
a(n) = Sum_{n=1..m} A097306(n, m), row sums of triangle of number of partitions of n into m odd parts.
Cf. A001318, A000041, A000700, A003724, A004111, A007837, A010815, A035294, A068049, A078408, A081360, A088670, A109950, A109968, A132312, A146061, A035363, A010054, A057077, A089806, A091602, A237515, A118457 (the partitions), A118459 (partition lengths), A015723 (total number of parts), A230957 (boustrophedon transform).
Cf. A167377 (complement).
Cf. A067659 (odd number of parts), A067661 (even number of parts).
Number of r-regular partitions for r = 2 through 12: A000009, A000726, A001935, A035959, A219601, A035985, A261775, A104502, A261776, A328545, A328546.

Programs

  • Haskell
    import Data.MemoCombinators (memo2, integral)
    a000009 n = a000009_list !! n
    a000009_list = map (pM 1) [0..] where
       pM = memo2 integral integral p
       p _ 0 = 1
       p k m | m < k     = 0
             | otherwise = pM (k + 1) (m - k) + pM (k + 1) m
    -- Reinhard Zumkeller, Sep 09 2015, Nov 05 2013
    
  • Julia
    # uses A010815
    using Memoize
    @memoize function A000009(n)
        n == 0 && return 1
        s = sum((-1)^k*A000009(n - k^2) for k in 1:isqrt(n))
        A010815(n) - 2*s
    end # Peter Luschny, Sep 09 2021
  • Magma
    Coefficients(&*[1+x^m:m in [1..100]])[1..100] where x is PolynomialRing(Integers()).1; // Sergei Haller (sergei(AT)sergei-haller.de), Dec 21 2006
    
  • Maple
    N := 100; t1 := series(mul(1+x^k,k=1..N),x,N); A000009 := proc(n) coeff(t1,x,n); end;
    spec := [ P, {P=PowerSet(N), N=Sequence(Z,card>=1)} ]: [ seq(combstruct[count](spec, size=n), n=0..58) ];
    spec := [ P, {P=PowerSet(N), N=Sequence(Z,card>=1)} ]: combstruct[allstructs](spec, size=10); # to get the actual partitions for n=10
    A000009 := proc(n)
        local x,m;
        product(1+x^m,m=1..n+1) ;
        expand(%) ;
        coeff(%,x,n) ;
    end proc: # R. J. Mathar, Jun 18 2016
    lim := 99; # Enlarge if more terms are needed.
    simplify(expand(QDifferenceEquations:-QPochhammer(-1, x, lim)/2, x)):
    seq(coeff(%, x, n), n=0..55); # Peter Luschny, Nov 17 2016
    # Alternative:
    a:= proc(n) option remember; `if`(n=0, 1, add(a(n-j)*add(
         `if`(d::odd, d, 0), d=numtheory[divisors](j)), j=1..n)/n)
        end:
    seq(a(n), n=0..55);  # Alois P. Heinz, Jun 24 2025
  • Mathematica
    PartitionsQ[Range[0, 60]] (* Harvey Dale, Jul 27 2009 *)
    a[ n_] := SeriesCoefficient[ Product[ 1 + x^k, {k, n}], {x, 0, n}]; (* Michael Somos, Jul 06 2011 *)
    a[ n_] := SeriesCoefficient[ 1 / Product[ 1 - x^k, {k, 1, n, 2}], {x, 0, n}]; (* Michael Somos, Jul 06 2011 *)
    a[ n_] := With[ {t = Log[q] / (2 Pi I)}, SeriesCoefficient[ q^(-1/24) DedekindEta[2 t] / DedekindEta[ t], {q, 0, n}]]; (* Michael Somos, Jul 06 2011 *)
    a[ n_] := SeriesCoefficient[ 1 / QPochhammer[ x, x^2], {x, 0, n}]; (* Michael Somos, May 24 2013 *)
    a[ n_] := SeriesCoefficient[ Series[ QHypergeometricPFQ[ {q}, {q x}, q, - q x], {q, 0, n}] /. x -> 1, {q, 0, n}]; (* Michael Somos, Mar 04 2014 *)
    a[ n_] := SeriesCoefficient[ QHypergeometricPFQ[{}, {}, q, -1] / 2, {q, 0, n}]; (* Michael Somos, Mar 04 2014 *)
    nmax = 60; CoefficientList[Series[Exp[Sum[(-1)^(k+1)/k*x^k/(1-x^k), {k, 1, nmax}]], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 25 2015 *)
    nmax = 100; poly = ConstantArray[0, nmax + 1]; poly[[1]] = 1; poly[[2]] = 1; Do[Do[poly[[j + 1]] += poly[[j - k + 1]], {j, nmax, k, -1}];, {k, 2, nmax}]; poly (* Vaclav Kotesovec, Jan 14 2017 *)
  • Maxima
    num_distinct_partitions(60,list); /* Emanuele Munarini, Feb 24 2014 */
    
  • Maxima
    h(n):=if oddp(n)=true then 1 else 0;
    S(n,m):=if n=0 then 1 else if nVladimir Kruchinin, Sep 07 2014 */
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff( prod( k=1, n, 1 + x^k, 1 + x * O(x^n)), n))}; /* Michael Somos, Nov 17 1999 */
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A) / eta(x + A), n))};
    
  • PARI
    {a(n) = my(c); forpart(p=n, if( n<1 || p[1]<2, c++; for(i=1, #p-1, if( p[i+1] > p[i]+1, c--; break)))); c}; /* Michael Somos, Aug 13 2017 */
    
  • PARI
    lista(nn) = {q='q+O('q^nn); Vec(eta(q^2)/eta(q))} \\ Altug Alkan, Mar 20 2018
    
  • Python
    # uses A010815
    from functools import lru_cache
    from math import isqrt
    @lru_cache(maxsize=None)
    def A000009(n): return 1 if n == 0 else A010815(n)+2*sum((-1)**(k+1)*A000009(n-k**2) for k in range(1,isqrt(n)+1)) # Chai Wah Wu, Sep 08 2021
    
  • Python
    import numpy as np
    n = 1000
    arr = np.zeros(n,dtype=object)
    arr[0] = 1
    for i in range(1,n):
        arr[i:] += arr[:n-i]
    print(arr) # Yigit Oktar, Jul 12 2025
    
  • SageMath
    # uses[EulerTransform from A166861]
    a = BinaryRecurrenceSequence(0, 1)
    b = EulerTransform(a)
    print([b(n) for n in range(56)]) # Peter Luschny, Nov 11 2020
    

Formula

G.f.: Product_{m>=1} (1 + x^m) = 1/Product_{m>=0} (1-x^(2m+1)) = Sum_{k>=0} Product_{i=1..k} x^i/(1-x^i) = Sum_{n>=0} x^(n*(n+1)/2) / Product_{k=1..n} (1-x^k).
G.f.: Sum_{n>=0} x^n*Product_{k=1..n-1} (1+x^k) = 1 + Sum_{n>=1} x^n*Product_{k>=n+1} (1+x^k). - Joerg Arndt, Jan 29 2011
Product_{k>=1} (1+x^(2k)) = Sum_{k>=0} x^(k*(k+1))/Product_{i=1..k} (1-x^(2i)) - Euler (Hardy and Wright, Theorem 346).
Asymptotics: a(n) ~ exp(Pi l_n / sqrt(3)) / ( 4 3^(1/4) l_n^(3/2) ) where l_n = (n-1/24)^(1/2) (Ayoub).
For n > 1, a(n) = (1/n)*Sum_{k=1..n} b(k)*a(n-k), with a(0)=1, b(n) = A000593(n) = sum of odd divisors of n; cf. A000700. - Vladeta Jovovic, Jan 21 2002
a(n) = t(n, 0), t as defined in A079211.
a(n) = Sum_{k=0..n-1} A117195(n,k) = A117192(n) + A117193(n) for n>0. - Reinhard Zumkeller, Mar 03 2006
a(n) = A026837(n) + A026838(n) = A118301(n) + A118302(n); a(A001318(n)) = A051044(n); a(A090864(n)) = A118303(n). - Reinhard Zumkeller, Apr 22 2006
Expansion of 1 / chi(-x) = chi(x) / chi(-x^2) = f(-x) / phi(x) = f(x) / phi(-x^2) = psi(x) / f(-x^2) = f(-x^2) / f(-x) = f(-x^4) / psi(-x) in powers of x where phi(), psi(), chi(), f() are Ramanujan theta functions. - Michael Somos, Mar 12 2011
G.f. is a period 1 Fourier series which satisfies f(-1 / (1152 t)) = 2^(-1/2) / f(t) where q = exp(2 Pi i t). - Michael Somos, Aug 16 2007
Expansion of q^(-1/24) * eta(q^2) / eta(q) in powers of q.
Expansion of q^(-1/24) 2^(-1/2) f2(t) in powers of q = exp(2 Pi i t) where f2() is a Weber function. - Michael Somos, Oct 18 2007
Given g.f. A(x), then B(x) = x * A(x^3)^8 satisfies 0 = f(B(x), B(x^2)) where f(u, v) = v - u^2 + 16*u*v^2 . - Michael Somos, May 31 2005
Given g.f. A(x), then B(x) = x * A(x^8)^3 satisfies 0 = f(B(x), B(x^3)) where f(u, v) = (u^3 - v) * (u + v^3) - 9 * u^3 * v^3. - Michael Somos, Mar 25 2008
From Evangelos Georgiadis, Andrew V. Sutherland, Kiran S. Kedlaya (egeorg(AT)mit.edu), Mar 03 2009: (Start)
a(0)=1; a(n) = 2*(Sum_{k=1..floor(sqrt(n))} (-1)^(k+1) a(n-k^2)) + sigma(n) where sigma(n) = (-1)^j if (n=(j*(3*j+1))/2 OR n=(j*(3*j-1))/2) otherwise sigma(n)=0 (simpler: sigma = A010815). (End)
From Gary W. Adamson, Jun 13 2009: (Start)
The product g.f. = (1/(1-x))*(1/(1-x^3))*(1/(1-x^5))*...; = (1,1,1,...)*
(1,0,0,1,0,0,1,0,0,1,...)*(1,0,0,0,0,1,0,0,0,0,1,0,0,0,0,1,...) * ...; =
a*b*c*... where a, a*b, a*b*c, ... converge to A000009:
1, 1, 1, 2, 2, 2, 3, 3, 3, 4, ... = a*b
1, 1, 1, 2, 2, 3, 4, 4, 5, 6, ... = a*b*c
1, 1, 1, 2, 2, 3, 4, 5, 6, 7, ... = a*b*c*d
1, 1, 1, 2, 2, 3, 4, 5, 6, 8, ... = a*b*c*d*e
1, 1, 1, 2, 2, 3, 4, 5, 6, 8, ... = a*b*c*d*e*f
... (cf. analogous example in A000041). (End)
a(A004526(n)) = A172033(n). - Reinhard Zumkeller, Jan 23 2010
a(n) = P(n) - P(n-2) - P(n-4) + P(n-10) + P(n-14) + ... + (-1)^m P(n-2p_m) + ..., where P(n) is the partition function (A000041) and p_m = m(3m-1)/2 is the m-th generalized pentagonal number (A001318). - Jerome Malenfant, Feb 16 2011
a(n) = A054242(n,0) = A201377(n,0). - Reinhard Zumkeller, Dec 02 2011
More precise asymptotics: a(n) ~ exp(Pi*sqrt((n-1/24)/3)) / (4*3^(1/4)*(n-1/24)^(3/4)) * (1 + (Pi^2-27)/(24*Pi*sqrt(3*(n-1/24))) + (Pi^4-270*Pi^2-1215)/(3456*Pi^2*(n-1/24))). - Vaclav Kotesovec, Nov 30 2015
a(n) = A067661(n) + A067659(n). Wolfdieter Lang, Jan 18 2016
From Vaclav Kotesovec, May 29 2016: (Start)
a(n) ~ exp(Pi*sqrt(n/3))/(4*3^(1/4)*n^(3/4)) * (1 + (Pi/(48*sqrt(3)) - (3*sqrt(3))/(8*Pi))/sqrt(n) + (Pi^2/13824 - 5/128 - 45/(128*Pi^2))/n).
a(n) ~ exp(Pi*sqrt(n/3) + (Pi/(48*sqrt(3)) - 3*sqrt(3)/(8*Pi))/sqrt(n) - (1/32 + 9/(16*Pi^2))/n) / (4*3^(1/4)*n^(3/4)).
(End)
a(n) = A089806(n)*A010815(floor(n/2)) + a(n-1) + a(n-2) - a(n-5) - a(n-7) + a(n-12) + ... + A057077(m-1)*a(n-A001318(m)) + ..., where n > A001318(m). - Gevorg Hmayakyan, Jul 07 2016
a(n) ~ Pi*BesselI(1, Pi*sqrt((n+1/24)/3)) / sqrt(24*n+1). - Vaclav Kotesovec, Nov 08 2016
a(n) = A000041(n) - A047967(n). - R. J. Mathar, Nov 20 2017
Sum_{n>=1} 1/a(n) = A237515. - Amiram Eldar, Nov 15 2020
From Peter Bala, Jan 15 2021: (Start)
G.f.: (1 + x)*Sum_{n >= 0} x^(n*(n+3)/2)/Product_{k = 1..n} (1 - x^k) =
(1 + x)*(1 + x^2)*Sum_{n >= 0} x^(n*(n+5)/2)/Product_{k = 1..n} (1 - x^k) = (1 + x)*(1 + x^2)*(1 + x^3)*Sum_{n >= 0} x^(n*(n+7)/2)/Product_{k = 1..n} (1 - x^k) = ....
G.f.: (1/2)*Sum_{n >= 0} x^(n*(n-1)/2)/Product_{k = 1..n} (1 - x^k) =
(1/2)*(1/(1 + x))*Sum_{n >= 0} x^((n-1)*(n-2)/2)/Product_{k = 1..n} (1 - x^k) = (1/2)*(1/((1 + x)*(1 + x^2)))*Sum_{n >= 0} x^((n-2)*(n-3)/2)/Product_{k = 1..n} (1 - x^k) = ....
G.f.: Sum_{n >= 0} x^n/Product_{k = 1..n} (1 - x^(2*k)) = (1/(1 - x)) * Sum_{n >= 0} x^(3*n)/Product_{k = 1..n} (1 - x^(2*k)) = (1/((1 - x)*(1 - x^3))) * Sum_{n >= 0} x^(5*n)/Product_{k = 1..n} (1 - x^(2*k)) = (1/((1 - x)*(1 - x^3)*(1 - x^5))) * Sum_{n >= 0} x^(7*n)/Product_{k = 1..n} (1 - x^(2*k)) = .... (End)
From Peter Bala, Feb 02 2021: (Start)
G.f.: A(x) = Sum_{n >= 0} x^(n*(2*n-1))/Product_{k = 1..2*n} (1 - x^k). (Set z = x and q = x^2 in Mc Laughlin et al. (2019 ArXiv version), Section 1.3, Identity 7.)
Similarly, A(x) = Sum_{n >= 0} x^(n*(2*n+1))/Product_{k = 1..2*n+1} (1 - x^k). (End)
a(n) = A001227(n) + A238005(n) + A238006(n). - R. J. Mathar, Sep 08 2021
G.f.: A(x) = exp ( Sum_{n >= 1} x^n/(n*(1 - x^(2*n))) ) = exp ( Sum_{n >= 1} (-1)^(n+1)*x^n/(n*(1 - x^n)) ). - Peter Bala, Dec 23 2021
Sum_{n>=0} a(n)/exp(Pi*n) = exp(Pi/24)/2^(1/8) = A292820. - Simon Plouffe, May 12 2023 [Proof: Sum_{n>=0} a(n)/exp(Pi*n) = phi(exp(-2*Pi)) / phi(exp(-Pi)), where phi(q) is the Euler modular function. We have phi(exp(-2*Pi)) = exp(Pi/12) * Gamma(1/4) / (2 * Pi^(3/4)) and phi(exp(-Pi)) = exp(Pi/24) * Gamma(1/4) / (2^(7/8) * Pi^(3/4)), see formulas (14) and (13) in I. Mező, 2013. - Vaclav Kotesovec, May 12 2023]
a(2*n) = Sum_{j=1..n} p(n+j, 2*j) and a(2*n+1) = Sum_{j=1..n+1} p(n+j,2*j-1), where p(n, s) is the number of partitions of n having exactly s parts. - Gregory L. Simay, Aug 30 2023

A027193 Number of partitions of n into an odd number of parts.

Original entry on oeis.org

0, 1, 1, 2, 2, 4, 5, 8, 10, 16, 20, 29, 37, 52, 66, 90, 113, 151, 190, 248, 310, 400, 497, 632, 782, 985, 1212, 1512, 1851, 2291, 2793, 3431, 4163, 5084, 6142, 7456, 8972, 10836, 12989, 15613, 18646, 22316, 26561, 31659, 37556, 44601, 52743, 62416, 73593, 86809, 102064, 120025, 140736
Offset: 0

Views

Author

Keywords

Comments

Number of partitions of n in which greatest part is odd.
Number of partitions of n+1 into an even number of parts, the least being 1. Example: a(5)=4 because we have [5,1], [3,1,1,1], [2,1,1] and [1,1,1,1,1,1].
Also number of partitions of n+1 such that the largest part is even and occurs only once. Example: a(5)=4 because we have [6], [4,2], [4,1,1] and [2,1,1,1,1]. - Emeric Deutsch, Apr 05 2006
Also the number of partitions of n such that the number of odd parts and the number of even parts have opposite parities. Example: a(8)=10 is a count of these partitions: 8, 611, 521, 431, 422, 41111, 332, 32111, 22211, 2111111. - Clark Kimberling, Feb 01 2014, corrected Jan 06 2021
In Chaves 2011 see page 38 equation (3.20). - Michael Somos, Dec 28 2014
Suppose that c(0) = 1, that c(1), c(2), ... are indeterminates, that d(0) = 1, and that d(n) = -c(n) - c(n-1)*d(1) - ... - c(0)*d(n-1). When d(n) is expanded as a polynomial in c(1), c(2),..,c(n), the terms are of the form H*c(i_1)*c(i_2)*...*c(i_k). Let P(n) = [c(i_1), c(i_2), ..., c(i_k)], a partition of n. Then H is negative if P has an odd number of parts, and H is positive if P has an even number of parts. That is, d(n) has A027193(n) negative coefficients, A027187(n) positive coefficients, and A000041 terms. The maximal coefficient in d(n), in absolute value, is A102462(n). - Clark Kimberling, Dec 15 2016

Examples

			G.f. = x + x^2 + 2*x^3 + 2*x^4 + 4*x^5 + 5*x^6 + 8*x^7 + 10*x^8 + 16*x^9 + 20*x^10 + ...
From _Gus Wiseman_, Feb 11 2021: (Start)
The a(1) = 1 through a(8) = 10 partitions into an odd number of parts are the following. The Heinz numbers of these partitions are given by A026424.
  (1)  (2)  (3)    (4)    (5)      (6)      (7)        (8)
            (111)  (211)  (221)    (222)    (322)      (332)
                          (311)    (321)    (331)      (422)
                          (11111)  (411)    (421)      (431)
                                   (21111)  (511)      (521)
                                            (22111)    (611)
                                            (31111)    (22211)
                                            (1111111)  (32111)
                                                       (41111)
                                                       (2111111)
The a(1) = 1 through a(8) = 10 partitions whose greatest part is odd are the following. The Heinz numbers of these partitions are given by A244991.
  (1)  (11)  (3)    (31)    (5)      (33)      (7)        (53)
             (111)  (1111)  (32)     (51)      (52)       (71)
                            (311)    (321)     (322)      (332)
                            (11111)  (3111)    (331)      (521)
                                     (111111)  (511)      (3221)
                                               (3211)     (3311)
                                               (31111)    (5111)
                                               (1111111)  (32111)
                                                          (311111)
                                                          (11111111)
(End)
		

References

  • N. J. Fine, Basic Hypergeometric Series and Applications, Amer. Math. Soc., 1988; p. 39, Example 7.

Crossrefs

The Heinz numbers of these partitions are A026424 or A244991.
The even-length version is A027187.
The case of odd sum as well as length is A160786, ranked by A340931.
The case of odd maximum as well as length is A340385.
Other cases of odd length:
- A024429 counts set partitions of odd length.
- A067659 counts strict partitions of odd length.
- A089677 counts ordered set partitions of odd length.
- A166444 counts compositions of odd length.
- A174726 counts ordered factorizations of odd length.
- A332304 counts strict compositions of odd length.
- A339890 counts factorizations of odd length.
A000009 counts partitions into odd parts, ranked by A066208.
A026804 counts partitions whose least part is odd.
A058695 counts partitions of odd numbers, ranked by A300063.
A072233 counts partitions by sum and length.
A101707 counts partitions of odd positive rank.

Programs

  • Maple
    g:=sum(x^(2*k)/product(1-x^j,j=1..2*k-1),k=1..40): gser:=series(g,x=0,50): seq(coeff(gser,x,n),n=1..45); # Emeric Deutsch, Apr 05 2006
  • Mathematica
    nn=40;CoefficientList[Series[ Sum[x^(2j+1)Product[1/(1- x^i),{i,1,2j+1}],{j,0,nn}],{x,0,nn}],x]  (* Geoffrey Critzer, Dec 01 2012 *)
    a[ n_] := If[ n < 0, 0, Length@Select[ IntegerPartitions[ n], OddQ[ Length@#] &]]; (* Michael Somos, Dec 28 2014 *)
    a[ n_] := If[ n < 1, 0, Length@Select[ IntegerPartitions[ n], OddQ[ First@#] &]]; (* Michael Somos, Dec 28 2014 *)
    a[ n_] := If[ n < 0, 0, Length@Select[ IntegerPartitions[ n + 1], #[[-1]] == 1 && EvenQ[ Length@#] &]]; (* Michael Somos, Dec 28 2014 *)
    a[ n_] := If[ n < 1, 0, Length@Select[ IntegerPartitions[ n + 1], EvenQ[ First@#] && (Length[#] < 2 || #[[1]] != #[[2]]) &]]; (* Michael Somos, Dec 28 2014 *)
  • PARI
    {a(n) = if( n<1, 0, polcoeff( sum( k=1, n, if( k%2, x^k / prod( j=1, k, 1 - x^j, 1 + x * O(x^(n-k)) ))), n))}; /* Michael Somos, Jul 24 2012 */
    
  • PARI
    q='q+O('q^66); concat([0], Vec( (1/eta(q)-eta(q)/eta(q^2))/2 ) ) \\ Joerg Arndt, Mar 23 2014

Formula

a(n) = (A000041(n) - (-1)^n*A000700(n)) / 2.
For g.f. see under A027187.
G.f.: Sum(k>=1, x^(2*k-1)/Product(j=1..2*k-1, 1-x^j ) ). - Emeric Deutsch, Apr 05 2006
G.f.: - Sum(k>=1, (-x)^(k^2)) / Product(k>=1, 1-x^k ). - Joerg Arndt, Feb 02 2014
G.f.: Sum(k>=1, x^(k*(2*k-1)) / Product(j=1..2*k, 1-x^j)). - Michael Somos, Dec 28 2014
a(2*n) = A000701(2*n), a(2*n-1) = A046682(2*n-1); a(n) = A000041(n)-A027187(n). - Reinhard Zumkeller, Apr 22 2006

A339890 Number of odd-length factorizations of n into factors > 1.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 3, 1, 1, 2, 2, 1, 2, 1, 4, 1, 1, 1, 4, 1, 1, 1, 3, 1, 2, 1, 2, 2, 1, 1, 6, 1, 2, 1, 2, 1, 3, 1, 3, 1, 1, 1, 5, 1, 1, 2, 5, 1, 2, 1, 2, 1, 2, 1, 8, 1, 1, 2, 2, 1, 2, 1, 6, 2, 1, 1, 5, 1, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Dec 28 2020

Keywords

Examples

			The a(n) factorizations for n = 24, 48, 60, 72, 96, 120:
  24      48          60       72          96          120
  2*2*6   2*3*8       2*5*6    2*4*9       2*6*8       3*5*8
  2*3*4   2*4*6       3*4*5    2*6*6       3*4*8       4*5*6
          3*4*4       2*2*15   3*3*8       4*4*6       2*2*30
          2*2*12      2*3*10   3*4*6       2*2*24      2*3*20
          2*2*2*2*3            2*2*18      2*3*16      2*4*15
                               2*3*12      2*4*12      2*5*12
                               2*2*2*3*3   2*2*2*2*6   2*6*10
                                           2*2*2*3*4   3*4*10
                                                       2*2*2*3*5
		

Crossrefs

The case of set partitions (or n squarefree) is A024429.
The case of partitions (or prime powers) is A027193.
The ordered version is A174726 (even: A174725).
The remaining (even-length) factorizations are counted by A339846.
A000009 counts partitions into odd parts, ranked by A066208.
A001055 counts factorizations, with strict case A045778.
A027193 counts partitions of odd length, ranked by A026424.
A058695 counts partitions of odd numbers, ranked by A300063.
A160786 counts odd-length partitions of odd numbers, ranked by A300272.
A316439 counts factorizations by product and length.
A340101 counts factorizations into odd factors.
A340102 counts odd-length factorizations into odd factors.

Programs

  • Maple
    g:= proc(n, k, t) option remember; `if`(n>k, 0, t)+
          `if`(isprime(n), 0, add(`if`(d>k, 0, g(n/d, d, 1-t)),
              d=numtheory[divisors](n) minus {1, n}))
        end:
    a:= n-> `if`(n<2, 0, g(n$2, 1)):
    seq(a(n), n=1..100);  # Alois P. Heinz, Dec 30 2020
  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n],OddQ@Length[#]&]],{n,100}]

Formula

a(n) + A339846(n) = A001055(n).

A058695 Number of ways to partition 2n+1 into positive integers.

Original entry on oeis.org

1, 3, 7, 15, 30, 56, 101, 176, 297, 490, 792, 1255, 1958, 3010, 4565, 6842, 10143, 14883, 21637, 31185, 44583, 63261, 89134, 124754, 173525, 239943, 329931, 451276, 614154, 831820, 1121505, 1505499, 2012558, 2679689, 3554345, 4697205, 6185689, 8118264, 10619863
Offset: 0

Views

Author

N. J. A. Sloane, Dec 31 2000

Keywords

Comments

A bisection of A000041, the other one is A058696.
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700). - Michael Somos, Feb 16 2014
a(n) is the number of partitions of 3n-1 having n as a part, for n >=1. Also, a(n+1) is the number of partitions of 3n having n as a part, for n >= 1. - Clark Kimberling, Mar 02 2014

Examples

			G.f. = 1 + 3*x + 7*x^2 + 15*x^3 + 30*x^4 + 56*x^5 + 101*x^6 + 176*x^7 + 297*x^8 + ...
G.f. = q^23 + 3*q^71 + 7*q^119 + 15*q^167 + 30*q^215 + 56*q^263 + 101*q^311 + ...
		

Crossrefs

Programs

  • Maple
    a:= n-> combinat[numbpart](2*n+1):
    seq(a(n), n=0..42);  # Alois P. Heinz, Jan 29 2020
  • Mathematica
    nn=100;Table[CoefficientList[Series[Product[1/(1-x^i),{i,1,nn}],{x,0,nn}],x][[2i]],{i,1,nn/2}] (* Geoffrey Critzer, Sep 28 2013 *)
    (* also *)
    Table[PartitionsP[2 n + 1], {n, 0, 40}] (* Clark Kimberling, Mar 02 2014 *)
    (* also *)
    Table[Count[IntegerPartitions[3 n - 1], p_ /; MemberQ[p, n]], {n, 20}]   (* Clark Kimberling, Mar 02 2014 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( 1 / eta(x + O(x^(2*n + 2))), 2*n + 1))}; /* Michael Somos, Apr 25 2003 */
    
  • PARI
    a(n) = numbpart(2*n+1); \\ Michel Marcus, Sep 28 2013

Formula

a(n) = A000041(2*n + 1).
Euler transform of period 16 sequence [ 3, 1, 2, 2, 2, 2, 3, 1, 3, 2, 2, 2, 2, 1, 3, 1, ...]. - Michael Somos, Apr 25 2003
G.f.: (Sum_{k>=0} x^A074377(k)) / (Product_{k>0} (1 - x^k))^2. - Michael Somos, Apr 25 2003
Expansion of f(x^1, x^7) / f(-x)^2 in powers of x where f() is a Ramanujan theta function. - Michael Somos, Feb 16 2014
Convolution of A000041 and A078408. - Michael Somos, Feb 16 2014

A006827 Number of partitions of 2n with all subsums different from n.

Original entry on oeis.org

1, 2, 5, 8, 17, 24, 46, 64, 107, 147, 242, 302, 488, 629, 922, 1172, 1745, 2108, 3104, 3737, 5232, 6419, 8988, 10390, 14552, 17292, 23160, 27206, 36975, 41945, 57058, 65291, 85895, 99384, 130443, 145283, 193554, 218947, 281860, 316326, 413322, 454229, 594048
Offset: 1

Views

Author

Keywords

Comments

Partitions of this type are also called non-biquanimous partitions. - Gus Wiseman, Apr 19 2024

Examples

			From _Gus Wiseman_, Apr 19 2024: (Start)
The a(1) = 1 through a(5) = 17 partitions (A = 10):
  (2)  (4)   (6)    (8)     (A)
       (31)  (42)   (53)    (64)
             (51)   (62)    (73)
             (222)  (71)    (82)
             (411)  (332)   (91)
                    (521)   (433)
                    (611)   (442)
                    (5111)  (622)
                            (631)
                            (721)
                            (811)
                            (3331)
                            (4222)
                            (6211)
                            (7111)
                            (22222)
                            (61111)
(End)
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

The complement is counted by A002219, ranks A357976.
Central diagonal of A046663.
The strict case is A321142, even bisection of A371794 (odd A078408).
This is the "bi-" version of A321451, ranks A321453.
Column k = 0 of A367094.
These partitions have Heinz numbers A371731.
Even bisection of A371795 (odd A058695).
A371783 counts k-quanimous partitions.

Programs

  • Maple
    b:= proc(n, i, s) option remember;
          `if`(0 in s or n in s, 0, `if`(n=0, 1, `if`(i<1, 0, b(n, i-1, s)+
          `if`(i<=n, b(n-i, i, select(y-> 0<=y and y<=n-i,
                     map(x-> [x, x-i][], s))), 0))))
        end:
    a:= n-> b(2*n, 2*n, {n}):
    seq(a(n), n=1..25);  # Alois P. Heinz, Jul 10 2012
  • Mathematica
    b[n_, i_, s_] := b[n, i, s] = If[MemberQ[s, 0 | n], 0, If[n == 0, 1, If[i < 1, 0, b[n, i-1, s] + If[i <= n, b[n-i, i, Select[Flatten[Transpose[{s, s-i}]], 0 <= # <= n-i &]], 0]]]]; a[n_] := b[2*n, 2*n, {n}]; Table[Print[an = a[n]]; an, {n, 1, 25}] (* Jean-François Alcover, Nov 12 2013, after Alois P. Heinz *)
  • Python
    from itertools import combinations_with_replacement
    from collections import Counter
    from sympy import npartitions
    from sympy.utilities.iterables import partitions
    def A006827(n): return npartitions(n<<1)-len({tuple(sorted((p+q).items())) for p, q in combinations_with_replacement(tuple(Counter(p) for p in partitions(n)),2)}) # Chai Wah Wu, Sep 20 2023

Formula

a(n) = A000041(2*n) - A002219(n).
a(n) = A046663(2*n,n).

Extensions

More terms from Don Reble, Nov 03 2001
More terms from Alois P. Heinz, Jul 10 2012

A035457 Number of partitions of n into parts of the form 4*k + 2.

Original entry on oeis.org

1, 0, 1, 0, 1, 0, 2, 0, 2, 0, 3, 0, 4, 0, 5, 0, 6, 0, 8, 0, 10, 0, 12, 0, 15, 0, 18, 0, 22, 0, 27, 0, 32, 0, 38, 0, 46, 0, 54, 0, 64, 0, 76, 0, 89, 0, 104, 0, 122, 0, 142, 0, 165, 0, 192, 0, 222, 0, 256, 0, 296, 0, 340, 0, 390, 0, 448, 0, 512, 0, 585, 0, 668, 0, 760, 0, 864, 0, 982, 0
Offset: 0

Views

Author

Keywords

Comments

Also number of partitions of n into distinct even parts. Example: a(10)=3 because we have [10],[8,2] and [6,4]. - Emeric Deutsch, Feb 22 2006
Also number of partitions of n into odd parts, each part occurring an even number of times. Example: a(10)=3 because we have [5,5], [3,3,1,1,1,1] and [1,1,1,1,1,1,1,1,1,1]. - Emeric Deutsch, Apr 08 2006

Examples

			a(10)=3 because we have [10], [6,2,2] and [2,2,2,2,2].
		

References

  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 288-289.

Crossrefs

Programs

  • Maple
    g:=product(1+x^(2*j),j=1..45): gser:=series(g,x=0,85): seq(coeff(gser,x,n),n=0..79); # Emeric Deutsch, Feb 22 2006; a(0) added by Georg Fischer, Dec 10 2020
  • Mathematica
    nn=80;CoefficientList[Series[Product[1+ x^(2i),{i,1,nn}],{x,0,nn}],x] (* Geoffrey Critzer, Jun 20 2014 *)
    nmax = 50; kmax = nmax/4; s = Range[0, kmax]*4 + 2;
    Table[Count[IntegerPartitions@n, x_ /; SubsetQ[s, x]], {n, 0, nmax}] (* Robert Price, Aug 03 2020 *)
  • PARI
    N=166;  S=2+sqrtint(N);  x='x+O('x^N);
    gf=sum(n=0, S, x^(n^2+n)/prod(k=1,n, 1-x^(2*k)) );
    Vec(gf)
    \\ Joerg Arndt, Feb 18 2014

Formula

G.f.: 1/Product_{n>=0} (1 - x^(4*n+2)).
G.f.: 1/Product_{j>=0} (1 - x^(8*j+2))*(1 - x^(8*j+6)).
G.f.: Product_{j>=1} (1 + x^(2*j)). - Emeric Deutsch, Feb 22 2006
a(2*n-1) = 0, a(2*n) = A000009(n). a(n) = A116675(n,0). - Emeric Deutsch, Feb 22 2006
G.f.: Sum_{n>=1} (x^(n*(n+1)) / Product_{k=1..n} (1 - x^(2*k))). - Joerg Arndt, Mar 10 2011
If n is even, a(n) ~ exp(Pi*sqrt(n/6)) / (2^(5/4) * 3^(1/4) * n^(3/4)). - Vaclav Kotesovec, Feb 26 2015
a(4*n) = A035294(n) and a(4*n+2) = A078408(n). - George Beck, Aug 19 2017

Extensions

More terms from Emeric Deutsch, Feb 22 2006
Description simplified by Joerg Arndt, Jun 24 2009
a(0)=1 prepended by Joerg Arndt, Mar 11 2011

A160786 The number of odd partitions of consecutive odd integers.

Original entry on oeis.org

1, 2, 4, 8, 16, 29, 52, 90, 151, 248, 400, 632, 985, 1512, 2291, 3431, 5084, 7456, 10836, 15613, 22316, 31659, 44601, 62416, 86809, 120025, 165028, 225710, 307161, 416006, 560864, 752877, 1006426, 1340012, 1777365, 2348821, 3093095, 4059416, 5310255, 6924691
Offset: 0

Views

Author

Utpal Sarkar (doetoe(AT)gmail.com), May 26 2009

Keywords

Comments

It seems that these are partitions of odd length and sum, ranked by A340931. The parts do not have to be odd. - Gus Wiseman, Apr 06 2021

Examples

			From _Gus Wiseman_, Apr 06 2021: (Start)
The a(0) = 1 through a(4) = 16 partitions:
  (1)  (3)    (5)      (7)        (9)
       (111)  (221)    (322)      (333)
              (311)    (331)      (432)
              (11111)  (421)      (441)
                       (511)      (522)
                       (22111)    (531)
                       (31111)    (621)
                       (1111111)  (711)
                                  (22221)
                                  (32211)
                                  (33111)
                                  (42111)
                                  (51111)
                                  (2211111)
                                  (3111111)
                                  (111111111)
(End)
		

Crossrefs

Partitions with all odd parts are counted by A000009 and ranked by A066208.
This is a bisection of A027193 (odd-length partitions), which is ranked by A026424.
The case of all odd parts is counted by A078408 and ranked by A300272.
The even version is A236913, ranked by A340784.
A multiplicative version is A340102.
These partitions are ranked by A340931.
A047993 counts balanced partitions, ranked by A106529.
A058695 counts partitions of odd numbers, ranked by A300063.
A072233 counts partitions by sum and length.
A236914 counts partition of type OO, ranked by A341448.
A340385 counts partitions with odd length and maximum, ranked by A340386.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, [1, 0$3],
          `if`(i<1, [0$4], b(n, i-1)+`if`(i>n, [0$4], (p->
          `if`(irem(i, 2)=0, [p[3], p[4], p[1], p[2]],
              [p[2], p[1], p[4], p[3]]))(b(n-i, i)))))
        end:
    a:= n-> b(2*n+1$2)[2]:
    seq(a(n), n=0..40);  # Alois P. Heinz, Feb 16 2014
  • Mathematica
    b[n_, i_] := b[n, i] = If[n==0, {1, 0, 0, 0}, If[i<1, {0, 0, 0, 0}, b[n, i-1] + If[i>n, {0, 0, 0, 0}, Function[{p}, If[Mod[i, 2]==0, p[[{3, 4, 1, 2}]], p[[{2, 1, 4, 3}]]]][b[n-i, i]]]]]; a[n_] := b[2*n+1, 2*n+1][[2]]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Jul 01 2015, after Alois P. Heinz *)
    (* Slow but easy to read *)
    a[n_] := Length@IntegerPartitions[2 n + 1, {1, 2 n + 1, 2}]
    a /@ Range[0, 25]
    (* Leo C. Stein, Nov 11 2020 *)
    (* Faster, don't build the partitions themselves *)
    (* Number of partitions of n into exactly k parts *)
    P[0, 0] = 1;
    P[n_, k_] := 0 /; ((k <= 0) || (n <= 0))
    P[n_, k_] := P[n, k] = P[n - k, k] + P[n - 1, k - 1]
    a[n_] := Sum[P[2 n + 1, k], {k, 1, 2 n + 1, 2}]
    a /@ Range[0, 40]
    (* Leo C. Stein, Nov 11 2020 *)
  • Python
    # Could be memoized for speedup
    def numoddpart(n, m=1):
        """The number of partitions of n into an odd number of parts of size at least m"""
        if n < m:
            return 0
        elif n == m:
            return 1
        else:
            # 1 (namely n = n) and all partitions of the form
            # k + even partitions that start with >= k
            return 1 + sum([numevenpart(n - k,  k) for k in range(m, n//3 + 1)])
    def numevenpart(n, m=1):
        """The number of partitions of n into an even number of parts of size at least m"""
        if n < 2*m:
            return 0
        elif n == 2*m:
            return 1
        else:
            return sum([numoddpart(n - k,  k) for k in range(m,  n//2 + 1)])
    [numoddpart(n) for n in range(1, 70, 2)]
    
  • Python
    # dict to memoize
    ps = {(0,0): 1}
    def p(n, k):
        """Number of partitions of n into exactly k parts"""
        if (n,k) in ps: return ps[(n,k)]
        if (n<=0) or (k<=0): return 0
        ps[(n,k)] = p(n-k,k) + p(n-1,k-1)
        return ps[(n,k)]
    def a(n): return sum([p(2*n+1, k) for k in range(1,2*n+3,2)])
    [a(n) for n in range(0,41)]
    # Leo C. Stein, Nov 11 2020

Formula

a(n) = A027193(2n+1).

A300272 Sorted list of Heinz numbers of odd partitions.

Original entry on oeis.org

2, 5, 8, 11, 17, 20, 23, 31, 32, 41, 44, 47, 50, 59, 67, 68, 73, 80, 83, 92, 97, 103, 109, 110, 124, 125, 127, 128, 137, 149, 157, 164, 167, 170, 176, 179, 188, 191, 197, 200, 211, 227, 230, 233, 236, 241, 242, 257, 268, 269, 272, 275, 277, 283, 292, 307, 310
Offset: 1

Views

Author

Gus Wiseman, Mar 01 2018

Keywords

Comments

An odd partition is an integer partition of an odd number into an odd number of parts, all of which are odd.
Any product of three members of this sequence is also in the sequence.

Examples

			Sequence of odd partitions begins: (1), (3), (111), (5), (7), (311), (9), (11), (11111), (13), (511), (15), (331), (17), (19), (711), (21), (31111), (23), (911), (25), (27), (29), (531), (1111), (333), (31), (1111111).
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],OddQ[Total[primeMS[#]]]&&And@@OddQ/@primeMS[#]&]

A096441 Number of palindromic and unimodal compositions of n. Equivalently, the number of orbits under conjugation of even nilpotent n X n matrices.

Original entry on oeis.org

1, 2, 2, 4, 3, 7, 5, 11, 8, 17, 12, 26, 18, 37, 27, 54, 38, 76, 54, 106, 76, 145, 104, 199, 142, 266, 192, 357, 256, 472, 340, 621, 448, 809, 585, 1053, 760, 1354, 982, 1740, 1260, 2218, 1610, 2818, 2048, 3559, 2590, 4485, 3264, 5616, 4097, 7018, 5120, 8728, 6378
Offset: 1

Views

Author

Nolan R. Wallach (nwallach(AT)ucsd.edu), Aug 10 2004

Keywords

Comments

Number of partitions of n such that all differences between successive parts are even, see example. [Joerg Arndt, Dec 27 2012]
Number of partitions of n where either all parts are odd or all parts are even. - Omar E. Pol, Aug 16 2013
From Gus Wiseman, Jan 13 2022: (Start)
Also the number of integer partitions of n with all even multiplicities (or run-lengths) except possibly the first. These are the conjugates of the partitions described by Joerg Arndt above. For example, the a(1) = 1 through a(8) = 11 partitions are:
(1) (2) (3) (4) (5) (6) (7) (8)
(11) (111) (22) (311) (33) (322) (44)
(211) (11111) (222) (511) (422)
(1111) (411) (31111) (611)
(2211) (1111111) (2222)
(21111) (3311)
(111111) (22211)
(41111)
(221111)
(2111111)
(11111111)
(End)

Examples

			From _Joerg Arndt_, Dec 27 2012: (Start)
There are a(10)=17 partitions of 10 where all differences between successive parts are even:
[ 1]  [ 1 1 1 1 1 1 1 1 1 1 ]
[ 2]  [ 2 2 2 2 2 ]
[ 3]  [ 3 1 1 1 1 1 1 1 ]
[ 4]  [ 3 3 1 1 1 1 ]
[ 5]  [ 3 3 3 1 ]
[ 6]  [ 4 2 2 2 ]
[ 7]  [ 4 4 2 ]
[ 8]  [ 5 1 1 1 1 1 ]
[ 9]  [ 5 3 1 1 ]
[10]  [ 5 5 ]
[11]  [ 6 2 2 ]
[12]  [ 6 4 ]
[13]  [ 7 1 1 1 ]
[14]  [ 7 3 ]
[15]  [ 8 2 ]
[16]  [ 9 1 ]
[17]  [ 10 ]
(End)
		

References

  • A. G. Elashvili and V. G. Kac, Classification of good gradings of simple Lie algebras. Lie groups and invariant theory, 85-104, Amer. Math. Soc. Transl. Ser. 2, 213, Amer. Math. Soc., Providence, RI, 2005.

Crossrefs

Bisections are A078408 and A096967.
The complement in partitions is counted by A006477
A version for compositions is A016116.
A pointed version is A035363, ranked by A066207.
A000041 counts integer partitions.
A025065 counts palindromic partitions.
A027187 counts partitions with even length/maximum.
A035377 counts partitions using multiples of 3.
A058696 counts partitions of even numbers, ranked by A300061.
A340785 counts factorizations into even factors.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(i>n, 0,
          `if`(irem(n, i)=0, 1, 0) +add(`if`(irem(j, 2)=0,
           b(n-i*j, i+1), 0), j=0..n/i))
        end:
    a:= n-> b(n, 1):
    seq(a(n), n=1..60);  # Alois P. Heinz, Mar 26 2014
  • Mathematica
    (* The following Mathematica program first generates all of the palindromic, unimodal compositions of n and then counts them. *)
    Pal[n_] := Block[{i, j, k, m, Q, L}, If[n == 1, Return[{{1}}]]; If[n == 2, Return[{{1, 1}, {2}}]]; L = {{n}}; If[Mod[n, 2] == 0, L = Append[L, {n/2, n/2}]]; For[i = 1, i < n, i++, Q = Pal[n - 2i]; m = Length[Q]; For[j = 1, j <= m, j++, If[i <= Q[[j, 1]], L = Append[L, Append[Prepend[Q[[j]], i], i]]]]]; L] NoPal[n_] := Length[Pal[n]]
    a[n_] := PartitionsQ[n] + If[EvenQ[n], PartitionsP[n/2], 0]; Table[a[n], {n, 1, 55}] (* Jean-François Alcover, Mar 17 2014, after Vladeta Jovovic *)
    Table[Length[Select[IntegerPartitions[n],And@@EvenQ/@Rest[Length/@Split[#]]&]],{n,1,30}] (* Gus Wiseman, Jan 13 2022 *)
  • PARI
    my(x='x+O('x^66)); Vec(eta(x^2)/eta(x)+1/eta(x^2)-2) \\ Joerg Arndt, Jan 17 2016

Formula

G.f.: sum(j>=1, q^j * (1-q^j)/prod(i=1..j, 1-q^(2*i) ) ).
G.f.: F + G - 2, where F = Product_{j>=1} 1/(1-q^(2*j)), G = Product_{j>=0} 1/(1-q^(2*j+1)).
a(2*n) = A000041(n) + A000009(2*n); a(2*n-1) = A000009(2*n-1). - Vladeta Jovovic, Aug 11 2004
a(n) = A000009(n) + A035363(n) = A000041(n) - A006477(n). - Omar E. Pol, Aug 16 2013

A340101 Number of factorizations of 2n + 1 into odd factors > 1.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Dec 28 2020

Keywords

Examples

			The factorizations for 2n + 1 = 27, 45, 135, 225, 315, 405, 1155:
  27      45      135       225       315       405         1155
  3*9     5*9     3*45      3*75      5*63      5*81        15*77
  3*3*3   3*15    5*27      5*45      7*45      9*45        21*55
          3*3*5   9*15      9*25      9*35      15*27       33*35
                  3*5*9     15*15     15*21     3*135       3*385
                  3*3*15    5*5*9     3*105     5*9*9       5*231
                  3*3*3*5   3*3*25    5*7*9     3*3*45      7*165
                            3*5*15    3*3*35    3*5*27      11*105
                            3*3*5*5   3*5*21    3*9*15      3*5*77
                                      3*7*15    3*3*5*9     3*7*55
                                      3*3*5*7   3*3*3*15    5*7*33
                                                3*3*3*3*5   3*11*35
                                                            5*11*21
                                                            7*11*15
                                                            3*5*7*11
		

Crossrefs

The version for partitions is A160786, ranked by A300272.
The even version is A340785.
The odd-length case is A340102.
A000009 counts partitions into odd parts, ranked by A066208.
A001055 counts factorizations, with strict case A045778.
A027193 counts partitions of odd length, ranked by A026424.
A058695 counts partitions of odd numbers, ranked by A300063.
A316439 counts factorizations by product and length.
Odd bisection of A001055, and also of A349907.

Programs

  • Maple
    g:= proc(n, k) option remember; `if`(n>k, 0, 1)+
          `if`(isprime(n), 0, add(`if`(d>k, 0, g(n/d, d)),
              d=numtheory[divisors](n) minus {1, n}))
        end:
    a:= n-> g(2*n+1$2):
    seq(a(n), n=0..100);  # Alois P. Heinz, Dec 30 2020
  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n],OddQ[Times@@#]&]],{n,1,100,2}]
  • PARI
    A001055(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A001055(n/d, d))); (s)); \\ After code in A001055
    A340101(n) = A001055(n+n+1); \\ Antti Karttunen, Dec 13 2021

Formula

a(n) = A001055(2n+1).
a(n) = A349907(2n+1). - Antti Karttunen, Dec 13 2021

Extensions

Data section extended up to 105 terms by Antti Karttunen, Dec 13 2021
Showing 1-10 of 59 results. Next