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

A000070 a(n) = Sum_{k=0..n} p(k) where p(k) = number of partitions of k (A000041).

Original entry on oeis.org

1, 2, 4, 7, 12, 19, 30, 45, 67, 97, 139, 195, 272, 373, 508, 684, 915, 1212, 1597, 2087, 2714, 3506, 4508, 5763, 7338, 9296, 11732, 14742, 18460, 23025, 28629, 35471, 43820, 53963, 66273, 81156, 99133, 120770, 146785, 177970, 215308, 259891, 313065, 376326, 451501
Offset: 0

Views

Author

Keywords

Comments

Also the total number of all different integers in all partitions of n + 1. E.g., a(3) = 7 because the partitions of 4 comprise the sets {1},{1, 2},{2},{1, 3},{4} of different integers and their total number is 7. - Thomas Wieder, Apr 10 2004
With offset 1, also the number of 1's in all partitions of n. For example, 3 = 2+1 = 1+1+1, a(3) = (zero 1's) + (one 1's) + (three 1's), so a(3) = 4. - Naohiro Nomoto, Jan 09 2002. See the Riordan reference p. 184, last formula, first term, for a proof based on Fine's identity given in Riordan, p. 182 (20).
Also, number of partitions of n into parts when there are two kinds of parts of size one.
Also number of graphical forest partitions of 2n+2.
a(n) = count 2 for each partition of n and 1 for each decrement. E.g., the partitions of 4 are 4 (2), 31 (3), 22 (2), 211 (3) and 1111 (2). 2 + 3 + 2 + 3 + 2 = 12. This is related to the Ferrers representation. We can see that taking the Ferrers diagram for each partition of n and adding a new * to all available columns, we generate each partition of n+1, but with repeats (A058884). - Jon Perry, Feb 06 2004
Also the number of 1-transitions among all integer partitions of n. A 1-transition is the removal of a digit "1" from a partition containing at least one "1" and subsequent addition of that "1" to another digit in that partition. This other digit may be a "1" also, but all digits of equal amount are considered as undistinquishable (unlabeled). E.g., for n=6 one has the partition [1113] for which the following two 1-transitions are possible: [1113] --> [123] and [1113] --> [114]. The 1-transitions of n form a partial order (poset). For n=6 one has 12 1-transitions: [111111] --> [11112], [11112] --> [1113], [11112] --> [1122], [1113] --> [114], [1113] --> [123], [1122] --> [123], [1122] --> [222], [123] --> [33], [123] --> [24], [114] --> [15], [114] --> [24], [15] --> [6]. - Thomas Wieder, Mar 08 2005
Also number of partitions of 2n+1 where one of the parts is greater than n (also where there are more than n parts) and of 2n+2 where one of the parts is greater than n+1 (or with more than n+1 parts). - Henry Bottomley, Aug 01 2005
Equals left border of triangle A137633 - Gary W. Adamson, Jan 31 2008
Equals row sums of triangle A027293. - Gary W. Adamson, Oct 26 2008
Convolved with A010815 = [1,1,1,...]. n-th partial sum of A000041 convolved with A010815 = the binomial sequence starting (1, n, ...). - Gary W. Adamson, Nov 09 2008
Equals A036469 convolved with A035363. - Gary W. Adamson, Jun 09 2009
a(A004526(n)) = A025065(n). - Reinhard Zumkeller, Jan 23 2010
a(n) = if n <= 1 then A054225(1,n) else A054225(n,1). - Reinhard Zumkeller, Nov 30 2011
Also the total number of 1's among all hook-lengths in all partitions of n. E.g., a(4)=7 because hooks of the partitions of n = 4 comprise the multisets {4,3,2,1}, {4,2,1,1}, {3,2,2,1}, {4,1,2,1}, {4,3,2,1} and their total number of 1's is 7. - T. Amdeberhan, Jun 03 2012
With offset 1, a(n) is also the difference between the sum of largest and the sum of second largest elements in all partitions of n. More generally, the number of occurrences of k in all partitions of n equals the difference between the sum of k-th largest and the sum of (k+1)st largest elements in all partitions of n. And more generally, the sum of the number of occurrences of k, k+1, k+2..k+m in all partitions of n equals the difference between the sum of k-th largest and the sum of (k+m+1)st largest elements in all partitions of n. - Omar E. Pol, Oct 25 2012
a(0) = 1 and 2*a(n-1) >= a(n) for all n > 0. Hence a(n) is a complete sequence. - Frank M Jackson, Apr 08 2013
a(n) is the number of conjugacy classes in the order-preserving, order-decreasing and (order-preserving and order-decreasing) injective transformation semigroups. - Ugbene Ifeanyichukwu, Jun 03 2015
a(n) is also the number of unlabeled subgraphs of the n-cycle C_n. For example, for n = 3, there are 3 unlabeled subgraphs of the triangle C_3 with 0 edges, 2 with 1 edge, 1 with 2 edges, and 1 with 3 edges (C_3 itself), so a(3) = 3 + 2 + 1 + 1 = 7. - John P. McSorley, Nov 21 2016
a(n) is also the number of partitions of 2n with all parts either even or equal to 1. Proof: the number of such partitions of 2n with exactly 2k 1's is p(n-k), for k = 0,..,n. Summing over k gives the formula. - Leonard Chastkofsky, Jul 24 2018
a(n) is the total number of polygamma functions that appear in the expansion of the (n+1)st derivative of x! with respect to x. More specifically, a(n) is the number of times the string "PolyGamma" appears in the expansion of D[x!, {x, n + 1}] in Mathematica. For example, D[x!, {x, 3 + 1}] = Gamma[1 + x] PolyGamma[0, 1 + x]^4 + 6 Gamma[1 + x] PolyGamma[0, 1 + x]^2 PolyGamma[1, 1 + x] + 3 Gamma[1 + x] PolyGamma[1, 1 + x]^2 + 4 Gamma[1 + x] PolyGamma[0, 1 + x] PolyGamma[2, 1 + x] + Gamma[1 + x] PolyGamma[3, 1 + x], and we see that the string "PolyGamma" appears a total of a(3) = 7 times in this expansion. - John M. Campbell, Aug 11 2018
With offset 1, also the number of integer partitions of 2n that do not comprise the multiset of vertex-degrees of any multigraph (i.e., non-multigraphical partitions); see A209816 for multigraphical partitions. - Gus Wiseman, Oct 26 2018
Also a(n) is the number of partitions of 2n+1 with exactly one odd part.
Delete the odd part 2k+1, k=0, ..., n, to get a partition of 2n-2k into even parts. There are as many unrestricted partitions of n-k; now sum those numbers from 0 to n to get a(n). - George Beck, Jul 22 2019
In the Young's lattice, a(n) is the number of branches that connect the (n-1)-th layer to the n-th layer. - Shouvik Datta, Sep 19 2021
a(n) is the number of multiset partitions of the multiset {r^n, s^1}, equivalently, factorization patterns of any number m=p^n*q^1 where p and q are primes. - Joerg Arndt, Jan 01 2024
a(n) is the number of positive integers whose divisors are the parts of the partitions of n + 1. - Omar E. Pol, Nov 07 2024

Examples

			G.f. = 1 + 2*x + 4*x^2 + 7*x^3 + 12*x^4 + 19*x^5 + 30*x^6 + 45*x^7 + 67*x^8 + ...
From _Omar E. Pol_, Oct 25 2012: (Start)
For n = 5 consider the partitions of n+1:
--------------------------------------
.                         Number
Partitions of 6           of 1's
--------------------------------------
6 .......................... 0
3 + 3 ...................... 0
4 + 2 ...................... 0
2 + 2 + 2 .................. 0
5 + 1 ...................... 1
3 + 2 + 1 .................. 1
4 + 1 + 1 .................. 2
2 + 2 + 1 + 1 .............. 2
3 + 1 + 1 + 1 .............. 3
2 + 1 + 1 + 1 + 1 .......... 4
1 + 1 + 1 + 1 + 1 + 1 ...... 6
------------------------------------
35-16 =                     19
.
The difference between the sum of the first column and the sum of the second column of the set of partitions of 6 is 35 - 16 = 19 and equals the number of 1's in all partitions of 6, so the 6th term of this sequence is a(5) = 19.
(End)
From _Gus Wiseman_, Oct 26 2018: (Start)
With offset 1, the a(1) = 1 through a(6) = 19 partitions of 2*n whose greatest part is > n:
  (2)  (4)   (6)    (8)     (A)      (C)
       (31)  (42)   (53)    (64)     (75)
             (51)   (62)    (73)     (84)
             (411)  (71)    (82)     (93)
                    (521)   (91)     (A2)
                    (611)   (622)    (B1)
                    (5111)  (631)    (732)
                            (721)    (741)
                            (811)    (822)
                            (6211)   (831)
                            (7111)   (921)
                            (61111)  (A11)
                                     (7221)
                                     (7311)
                                     (8211)
                                     (9111)
                                     (72111)
                                     (81111)
                                     (711111)
With offset 1, the a(1) = 1 through a(6) = 19 partitions of 2*n whose number of parts is > n:
  (11)  (211)   (2211)    (22211)     (222211)      (2222211)
        (1111)  (3111)    (32111)     (322111)      (3222111)
                (21111)   (41111)     (331111)      (3321111)
                (111111)  (221111)    (421111)      (4221111)
                          (311111)    (511111)      (4311111)
                          (2111111)   (2221111)     (5211111)
                          (11111111)  (3211111)     (6111111)
                                      (4111111)     (22221111)
                                      (22111111)    (32211111)
                                      (31111111)    (33111111)
                                      (211111111)   (42111111)
                                      (1111111111)  (51111111)
                                                    (222111111)
                                                    (321111111)
                                                    (411111111)
                                                    (2211111111)
                                                    (3111111111)
                                                    (21111111111)
                                                    (111111111111)
(End)
From _Joerg Arndt_, Jan 01 2024: (Start)
The a(5) = 19 multiset partitions of the multiset {1^5, 2^1} are:
   1:  {{1, 1, 1, 1, 1, 2}}
   2:  {{1, 1, 1, 1, 1}, {2}}
   3:  {{1, 1, 1, 1, 2}, {1}}
   4:  {{1, 1, 1, 1}, {1, 2}}
   5:  {{1, 1, 1, 1}, {1}, {2}}
   6:  {{1, 1, 1, 2}, {1, 1}}
   7:  {{1, 1, 1, 2}, {1}, {1}}
   8:  {{1, 1, 1}, {1, 1, 2}}
   9:  {{1, 1, 1}, {1, 1}, {2}}
  10:  {{1, 1, 1}, {1, 2}, {1}}
  11:  {{1, 1, 1}, {1}, {1}, {2}}
  12:  {{1, 1, 2}, {1, 1}, {1}}
  13:  {{1, 1, 2}, {1}, {1}, {1}}
  14:  {{1, 1}, {1, 1}, {1, 2}}
  15:  {{1, 1}, {1, 1}, {1}, {2}}
  16:  {{1, 1}, {1, 2}, {1}, {1}}
  17:  {{1, 1}, {1}, {1}, {1}, {2}}
  18:  {{1, 2}, {1}, {1}, {1}, {1}}
  19:  {{1}, {1}, {1}, {1}, {1}, {2}}
(End)
		

References

  • H. Gupta, An asymptotic formula in partitions. J. Indian Math. Soc., (N. S.) 10 (1946), 73-76.
  • H. Gupta et al., Tables of Partitions. Royal Society Mathematical Tables, Vol. 4, Cambridge Univ. Press, 1958, p. 90.
  • R. Honsberger, Mathematical Gems III, M.A.A., 1985, p. 6.
  • D. E. Knuth, The Art of Computer Programming, Vol. 4A, Table A-1, page 778. - N. J. A. Sloane, Dec 30 2018
  • A. M. Odlyzko, Asymptotic Enumeration Methods, p. 19
  • J. Riordan, Combinatorial Identities, Wiley, 1968, p. 199.
  • 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).
  • Stanley, R. P., Exercise 1.26 in Enumerative Combinatorics, Vol. 1. Cambridge, England: Cambridge University Press, p. 59, 1999.

Crossrefs

A diagonal of A066633.
Also second column of A126442. - George Beck, May 07 2011
Row sums of triangle A092905.
Also row sums of triangle A261555. - Omar E. Pol, Sep 14 2016
Also row sums of triangle A278427. - John P. McSorley, Nov 25 2016
Column k=2 of A292508.

Programs

  • GAP
    List([0..45],n->Sum([0..n],k->NrPartitions(k))); # Muniru A Asiru, Jul 25 2018
    
  • Haskell
    a000070 = p a028310_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, Nov 06 2012
    
  • Maple
    with(combinat): a:=n->add(numbpart(j),j=0..n): seq(a(n), n=0..44); # Zerinvary Lajos, Aug 26 2008
  • Mathematica
    CoefficientList[ Series[1/(1 - x)*Product[1/(1 - x^k), {k, 75}], {x, 0, 45}], x] (* Robert G. Wilson v, Jul 13 2004 *)
    Table[ Count[ Flatten@ IntegerPartitions@ n, 1], {n, 45}] (* Robert G. Wilson v, Aug 06 2008 *)
    Join[{1}, Accumulate[PartitionsP[Range[50]]]+1] (* _Harvey P. Dale, Mar 12 2013 *)
    a[ n_] := SeriesCoefficient[ 1 / (1 - x) / QPochhammer[ x], {x, 0, n}]; (* Michael Somos, Nov 09 2013 *)
    Accumulate[PartitionsP[Range[0,49]]] (* George Beck, Oct 23 2014; typo fixed by Virgile Andreani, Jul 10 2016 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( 1 / prod(m=1, n, 1 - x^m, 1 + x * O(x^n)) / (1 - x), n))}; /* Michael Somos, Nov 08 2002 */
    
  • PARI
    x='x+O('x^66); Vec(1/((1-x)*eta(x))) /* Joerg Arndt, May 15 2011 */
    
  • PARI
    a(n) = sum(k=0, n, numbpart(k)); \\ Michel Marcus, Sep 16 2016
    
  • Python
    from itertools import accumulate
    def A000070iter(n):
        L = [0]*n; L[0] = 1
        def numpart(n):
            S = 0; J = n-1; k = 2
            while 0 <= J:
                T = L[J]
                S = S+T if (k//2)%2 else S-T
                J -= k  if (k)%2 else k//2
                k += 1
            return S
        for j in range(1, n): L[j] = numpart(j)
        return accumulate(L)
    print(list(A000070iter(100))) # Peter Luschny, Aug 30 2019
    
  • Python
    # Using function A365676Row. Compare also A365675.
    from itertools import accumulate
    def A000070List(size: int) -> list[int]:
        return [sum(accumulate(reversed(A365676Row(n)))) for n in range(size)]
    print(A000070List(45))  # Peter Luschny, Sep 16 2023
  • Sage
    def A000070_list(leng):
        p = [number_of_partitions(n) for n in range(leng)]
        return [add(p[:k+1]) for k in range(leng)]
    A000070_list(45) # Peter Luschny, Sep 15 2014
    

Formula

Euler transform of [ 2, 1, 1, 1, 1, 1, 1, ...].
log(a(n)) ~ -3.3959 + 2.44613*sqrt(n). - Robert G. Wilson v, Jan 11 2002
a(n) = (1/n)*Sum_{k=1..n} (sigma(k)+1)*a(n-k), n > 1, a(0) = 1. - Vladeta Jovovic, Aug 22 2002
G.f.: (1/(1 - x))*Product_{m >= 1} 1/(1 - x^m).
a(n) seems to have the same parity as A027349(n+1). Comment from James Sellers, Mar 08 2006: that is true.
a(n) = A000041(2n+1) - A110618(2n+1) = A000041(2n+2) - A110618(2n+2). - Henry Bottomley, Aug 01 2005
Row sums of triangle A133735. - Gary W. Adamson, Sep 22 2007
a(n) = A092269(n+1) - A195820(n+1). - Omar E. Pol, Oct 20 2011
a(n) = A181187(n+1,1) - A181187(n+1,2). - Omar E. Pol, Oct 25 2012
From Peter Bala, Dec 23 2013: (Start)
Gupta gives the asymptotic result a(n-1) ~ sqrt(6/Pi^2)* sqrt(n)*p(n), where p(n) is the partition function A000041(n).
Let P(2,n) denote the set of partitions of n into parts k >= 2.
a(n-2) = Sum_{parts k in all partitions in P(2,n)} phi(k), where phi(k) is the Euler totient function (see A000010). Using this result and Mertens's theorem on the average order of the phi function, leads to the asymptotic result
a(n-2) ~ (6/Pi^2)*n*(p(n) - p(n-1)) = (6/Pi^2)*A138880(n) as n -> infinity. (End)
a(n) ~ exp(Pi*sqrt(2*n/3)) / (2^(3/2)*Pi*sqrt(n)) * (1 + 11*Pi/(24*sqrt(6*n)) + (73*Pi^2 - 1584)/(6912*n)). - Vaclav Kotesovec, Oct 26 2016
a(n) = A024786(n+2) + A024786(n+1). - Vaclav Kotesovec, Nov 05 2016
G.f.: exp(Sum_{k>=1} (sigma_1(k) + 1)*x^k/k). - Ilya Gutkovskiy, Aug 21 2018
a(n) = A025065(2n). - Gus Wiseman, Oct 26 2018
a(n - 1) = A000041(2n) - A209816(n). - Gus Wiseman, Oct 26 2018

A045778 Number of factorizations of n into distinct factors greater than 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 3, 1, 2, 2, 2, 1, 3, 1, 3, 2, 2, 1, 5, 1, 2, 2, 3, 1, 5, 1, 3, 2, 2, 2, 5, 1, 2, 2, 5, 1, 5, 1, 3, 3, 2, 1, 7, 1, 3, 2, 3, 1, 5, 2, 5, 2, 2, 1, 9, 1, 2, 3, 4, 2, 5, 1, 3, 2, 5, 1, 9, 1, 2, 3, 3, 2, 5, 1, 7, 2, 2, 1, 9, 2, 2, 2, 5, 1, 9, 2, 3, 2, 2, 2, 10, 1, 3, 3, 5, 1, 5, 1, 5
Offset: 1

Views

Author

Keywords

Comments

This sequence depends only on the prime signature of n and not on the actual value of n.
Also the number of strict multiset partitions (sets of multisets) of the prime factors of n. - Gus Wiseman, Dec 03 2016
Number of sets of integers greater than 1 whose product is n. - Antti Karttunen, Feb 20 2024

Examples

			24 can be factored as 24, 2*12, 3*8, 4*6, or 2*3*4, so a(24) = 5. The factorization 2*2*6 is not permitted because the factor 2 is present twice. a(1) = 1 represents the empty factorization.
		

Crossrefs

Cf. A036469, A114591, A114592, A316441 (Dirichlet inverse).
Cf. A156648 (2*Dgf at s=2), A073017 (2*Dgf at s=3), A258870 (2*Dgf at s=4).
Cf. also A069626 (Number of sets of integers > 1 whose least common multiple is n).
Cf. A287549 (partial sums).

Programs

  • APL
    ⍝ Dyalog dialect
    divisors ← {ð←⍵{(0=⍵|⍺)/⍵}⍳⌊⍵*÷2 ⋄ 1=⍵:ð ⋄ ð, (⍵∘÷)¨(⍵=(⌊⍵*÷2)*2)↓⌽ð}
    A045778 ← { D←1↓divisors(⍵) ⋄ T←(⍴D)⍴2 ⋄ +/⍵⍷{×/D/⍨T⊤⍵}¨(-∘1)⍳2*⍴D } ⍝ (simple, but a memory hog)
    A045778 ← { ⍺←⌽divisors(⍵) ⋄ 1=⍵:1 ⋄ 0=≢⍺:0 ⋄ R←⍺↓⍨⍺⍳⍵∘÷ ⋄ Ð←{⍺/⍨0=⍺|⍵} ⋄ +/(((R)Ð⊢)∇⊢)¨(⍵∘÷)¨⍺ } ⍝ (more efficient) - Antti Karttunen, Feb 20 2024
  • Maple
    with(numtheory):
    b:= proc(n, k) option remember;
          `if`(n>k, 0, 1) +`if`(isprime(n), 0,
          add(`if`(d>k, 0, b(n/d, d-1)), d=divisors(n) minus {1, n}))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=1..120);  # Alois P. Heinz, May 26 2013
  • Mathematica
    gd[m_, 1] := 1; gd[1, n_] := 0; gd[1, 1] := 1; gd[0, n_] := 0; gd[m_, n_] := gd[m, n] = Total[gd[# - 1, n/#] & /@ Select[Divisors[n], # <= m &]]; Array[ gd[#, #] &, 100]  (* Alexander Adam, Dec 28 2012 *)
  • PARI
    v=vector(100,k,k==1); for(n=2,#v, v+=dirmul(v,vector(#v,k,k==n)) ); v /* Max Alekseyev, Jul 16 2014 */
    
  • PARI
    A045778(n, k=n) = ((n<=k) + sumdiv(n, d, if(d > 1 && d <= k && d < n, A045778(n/d, d-1)))); \\ After Alois P. Heinz's Maple-code by Antti Karttunen, Jul 23 2017, edited Feb 20 2024
    
  • PARI
    A045778(n, m=n) = if(1==n, 1, sumdiv(n,d,if((d>1)&&(d<=m),A045778(n/d,d-1)))); \\ Antti Karttunen, Feb 20 2024
    
  • Python
    from sympy.core.cache import cacheit
    from sympy import divisors, isprime
    @cacheit
    def b(n, k): return (0 if n>k else 1) + (0 if isprime(n) else sum(0 if d>k else b(n//d, d - 1) for d in divisors(n)[1:-1]))
    def a(n): return b(n, n)
    print([a(n) for n in range(1, 121)]) # Indranil Ghosh, Aug 19 2017, after Maple code
    

Formula

Dirichlet g.f.: Product_{n>=2} (1 + 1/n^s).
Let p and q be two distinct prime numbers and k a natural number. Then a(p^k) = A000009(k) and a(p^k*q) = A036469(k). - Alexander Adam, Dec 28 2012
Let p_i with 1<=i<=k k distinct prime numbers. Then a(Product_{i=1..k} p_i) = A000110(k). - Alexander Adam, Dec 28 2012

Extensions

Edited by Franklin T. Adams-Watters, Jun 04 2009

A035363 Number of partitions of n into even parts.

Original entry on oeis.org

1, 0, 1, 0, 2, 0, 3, 0, 5, 0, 7, 0, 11, 0, 15, 0, 22, 0, 30, 0, 42, 0, 56, 0, 77, 0, 101, 0, 135, 0, 176, 0, 231, 0, 297, 0, 385, 0, 490, 0, 627, 0, 792, 0, 1002, 0, 1255, 0, 1575, 0, 1958, 0, 2436, 0, 3010, 0, 3718, 0, 4565, 0, 5604, 0, 6842, 0, 8349, 0, 10143, 0, 12310, 0
Offset: 0

Views

Author

Keywords

Comments

Convolved with A036469 = A000070. - Gary W. Adamson, Jun 09 2009
Note that these partitions are located in the head of the last section of the set of partitions of n (see A135010). - Omar E. Pol, Nov 20 2009
Number of symmetric unimodal compositions of n+2 where the maximal part appears twice, see example. Also number of symmetric unimodal compositions of n where the maximal part appears an even number of times. - Joerg Arndt, Jun 11 2013
Number of partitions of n having parts of even multiplicity. These are the conjugates of the partitions from the definition. Example: a(8)=5 because we have [4,4],[3,3,1,1],[2,2,2,2],[2,2,1,1,1,1], and [1,1,1,1,1,1,1,1]. - Emeric Deutsch, Jan 27 2016
From Gus Wiseman, May 22 2021: (Start)
The Heinz numbers of the conjugate partitions described in Emeric Deutsch's comment above are given by A000290.
For n > 1, also the number of integer partitions of n-1 whose only odd part is the smallest. The Heinz numbers of these partitions are given by A341446. For example, the a(2) = 1 through a(14) = 15 partitions (empty columns shown as dots, A..D = 10..13) are:
1 . 3 . 5 . 7 . 9 . B . D
21 41 43 63 65 85
221 61 81 83 A3
421 441 A1 C1
2221 621 443 643
4221 641 661
22221 821 841
4421 A21
6221 4441
42221 6421
222221 8221
44221
62221
422221
2222221
Also the number of integer partitions of n whose greatest part is the sum of all the other parts. The Heinz numbers of these partitions are given by A344415. For example, the a(2) = 1 through a(12) = 11 partitions (empty columns not shown) are:
(11) (22) (33) (44) (55) (66)
(211) (321) (422) (532) (633)
(3111) (431) (541) (642)
(4211) (5221) (651)
(41111) (5311) (6222)
(52111) (6321)
(511111) (6411)
(62211)
(63111)
(621111)
(6111111)
Also the number of integer partitions of n of length n/2. The Heinz numbers of these partitions are given by A340387. For example, the a(2) = 1 through a(14) = 15 partitions (empty columns not shown) are:
(2) (22) (222) (2222) (22222) (222222) (2222222)
(31) (321) (3221) (32221) (322221) (3222221)
(411) (3311) (33211) (332211) (3322211)
(4211) (42211) (333111) (3332111)
(5111) (43111) (422211) (4222211)
(52111) (432111) (4322111)
(61111) (441111) (4331111)
(522111) (4421111)
(531111) (5222111)
(621111) (5321111)
(711111) (5411111)
(6221111)
(6311111)
(7211111)
(8111111)
(End)

Examples

			From _Joerg Arndt_, Jun 11 2013: (Start)
There are a(12)=11 symmetric unimodal compositions of 12+2=14 where the maximal part appears twice:
01:  [ 1 1 1 1 1 2 2 1 1 1 1 1 ]
02:  [ 1 1 1 1 3 3 1 1 1 1 ]
03:  [ 1 1 1 4 4 1 1 1 ]
04:  [ 1 1 2 3 3 2 1 1 ]
05:  [ 1 1 5 5 1 1 ]
06:  [ 1 2 4 4 2 1 ]
07:  [ 1 6 6 1 ]
08:  [ 2 2 3 3 2 2 ]
09:  [ 2 5 5 2 ]
10:  [ 3 4 4 3 ]
11:  [ 7 7 ]
There are a(14)=15 symmetric unimodal compositions of 14 where the maximal part appears an even number of times:
01:  [ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ]
02:  [ 1 1 1 1 1 2 2 1 1 1 1 1 ]
03:  [ 1 1 1 1 3 3 1 1 1 1 ]
04:  [ 1 1 1 2 2 2 2 1 1 1 ]
05:  [ 1 1 1 4 4 1 1 1 ]
06:  [ 1 1 2 3 3 2 1 1 ]
07:  [ 1 1 5 5 1 1 ]
08:  [ 1 2 2 2 2 2 2 1 ]
09:  [ 1 2 4 4 2 1 ]
10:  [ 1 3 3 3 3 1 ]
11:  [ 1 6 6 1 ]
12:  [ 2 2 3 3 2 2 ]
13:  [ 2 5 5 2 ]
14:  [ 3 4 4 3 ]
15:  [ 7 7 ]
(End)
a(8)=5 because we  have [8], [6,2], [4,4], [4,2,2], and [2,2,2,2]. - _Emeric Deutsch_, Jan 27 2016
From _Gus Wiseman_, May 22 2021: (Start)
The a(0) = 1 through a(12) = 11 partitions into even parts are the following (empty columns shown as dots, A = 10, C = 12). The Heinz numbers of these partitions are given by A066207.
  ()  .  (2)  .  (4)   .  (6)    .  (8)     .  (A)      .  (C)
                 (22)     (42)      (44)       (64)        (66)
                          (222)     (62)       (82)        (84)
                                    (422)      (442)       (A2)
                                    (2222)     (622)       (444)
                                               (4222)      (642)
                                               (22222)     (822)
                                                           (4422)
                                                           (6222)
                                                           (42222)
                                                           (222222)
(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.

Crossrefs

Bisection (even part) gives the partition numbers A000041.
Column k=0 of A103919, A264398.
Note: A-numbers of ranking sequences are in parentheses below.
The version for odd instead of even parts is A000009 (A066208).
The version for parts divisible by 3 instead of 2 is A035377.
The strict case is A035457.
The Heinz numbers of these partitions are given by A066207.
The ordered version (compositions) is A077957 prepended by (1,0).
This is column k = 2 of A168021.
The multiplicative version (factorizations) is A340785.
A000569 counts graphical partitions (A320922).
A004526 counts partitions of length 2 (A001358).
A025065 counts palindromic partitions (A265640).
A027187 counts partitions with even length/maximum (A028260/A244990).
A058696 counts partitions of even numbers (A300061).
A067661 counts strict partitions of even length (A030229).
A236913 counts partitions of even length and sum (A340784).
A340601 counts partitions of even rank (A340602).
The following count partitions of even length:
- A096373 cannot be partitioned into strict pairs (A320891).
- A338914 can be partitioned into strict pairs (A320911).
- A338915 cannot be partitioned into distinct pairs (A320892).
- A338916 can be partitioned into distinct pairs (A320912).
- A339559 cannot be partitioned into distinct strict pairs (A320894).
- A339560 can be partitioned into distinct strict pairs (A339561).

Programs

  • Maple
    ZL:= [S, {C = Cycle(B), S = Set(C), E = Set(B), B = Prod(Z,Z)}, unlabelled]: seq(combstruct[count](ZL, size=n), n=0..69); # Zerinvary Lajos, Mar 26 2008
    g := 1/mul(1-x^(2*k), k = 1 .. 100): gser := series(g, x = 0, 80): seq(coeff(gser, x, n), n = 0 .. 78); # Emeric Deutsch, Jan 27 2016
    # Using the function EULER from Transforms (see link at the bottom of the page).
    [1,op(EULER([0,1,seq(irem(n,2),n=0..66)]))]; # Peter Luschny, Aug 19 2020
    # next Maple program:
    a:= n-> `if`(n::odd, 0, combinat[numbpart](n/2)):
    seq(a(n), n=0..84);  # Alois P. Heinz, Jun 22 2021
  • Mathematica
    nmax = 50; s = Range[2, nmax, 2];
    Table[Count[IntegerPartitions@n, x_ /; SubsetQ[s, x]], {n, 0, nmax}] (* Robert Price, Aug 05 2020 *)
  • Python
    from sympy import npartitions
    def A035363(n): return 0 if n&1 else npartitions(n>>1) # Chai Wah Wu, Sep 23 2023

Formula

G.f.: Product_{k even} 1/(1 - x^k).
Convolution with the number of partitions into distinct parts (A000009, which is also number of partitions into odd parts) gives the number of partitions (A000041). - Franklin T. Adams-Watters, Jan 06 2006
If n is even then a(n)=A000041(n/2) otherwise a(n)=0. - Omar E. Pol, Nov 20 2009
G.f.: 1 + x^2*(1 - G(0))/(1-x^2) where G(k) = 1 - 1/(1-x^(2*k+2))/(1-x^2/(x^2-1/G(k+1) )); (recursively defined continued fraction). - Sergei N. Gladkovskii, Jan 23 2013
a(n) = A096441(n) - A000009(n), n >= 1. - Omar E. Pol, Aug 16 2013
G.f.: exp(Sum_{k>=1} x^(2*k)/(k*(1 - x^(2*k)))). - Ilya Gutkovskiy, Aug 13 2018

A038348 Expansion of (1/(1-x^2))*Product_{m>=0} 1/(1-x^(2m+1)).

Original entry on oeis.org

1, 1, 2, 3, 4, 6, 8, 11, 14, 19, 24, 31, 39, 49, 61, 76, 93, 114, 139, 168, 203, 244, 292, 348, 414, 490, 579, 682, 801, 938, 1097, 1278, 1487, 1726, 1999, 2311, 2667, 3071, 3531, 4053, 4644, 5313, 6070, 6923, 7886, 8971, 10190, 11561
Offset: 0

Views

Author

Keywords

Comments

Number of partitions of n+2 with exactly one even part. - Vladeta Jovovic, Sep 10 2003
Also, number of partitions of n with at most one even part. - Vladeta Jovovic, Sep 10 2003
Also total number of parts, counted without multiplicity, in all partitions of n into odd parts, offset 1. - Vladeta Jovovic, Mar 27 2005
a(n) = Sum_{k>=1} k*A116674(n+1,k). - Emeric Deutsch, Feb 22 2006
Equals row sums of triangle A173305. - Gary W. Adamson, Feb 15 2010
Equals partial sums of A025147 (observed by Jonathan Vos Post, proved by several correspondents).
Conjecture: The n-th derivative of Gamma(x+1) at x = 0 has a(n+1) terms. For example, d^4/dx^4_(x = 0) Gamma(x+1) = 8*eulergamma*zeta(3) + eulergamma^4 + eulergamma^2*Pi^2 + 3*Pi^4/20 which has a(5) = 4 terms. - David Ulgenes, Dec 05 2023

Examples

			From _Gus Wiseman_, Sep 23 2019: (Start)
Also the number of integer partitions of n that are strict except possibly for any number of 1's. For example, the a(1) = 1 through a(7) = 11 partitions are:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)
       (11)  (21)   (31)    (32)     (42)      (43)
             (111)  (211)   (41)     (51)      (52)
                    (1111)  (311)    (321)     (61)
                            (2111)   (411)     (421)
                            (11111)  (3111)    (511)
                                     (21111)   (3211)
                                     (111111)  (4111)
                                               (31111)
                                               (211111)
                                               (1111111)
(End)
		

Crossrefs

Programs

  • Maple
    f:=1/(1-x^2)/product(1-x^(2*j-1),j=1..32): fser:=series(f,x=0,62): seq(coeff(fser,x,n),n=0..58); # Emeric Deutsch, Feb 22 2006
  • Mathematica
    mmax = 47; CoefficientList[ Series[ (1/(1-x^2))*Product[1/(1-x^(2m+1)), {m, 0, mmax}], {x, 0, mmax}], x] (* Jean-François Alcover, Jun 21 2011 *)
  • SageMath
    # uses[EulerTransform from A166861]
    def g(n): return n % 2 if n > 2 else 1
    a = EulerTransform(g)
    print([a(n) for n in range(48)]) # Peter Luschny, Dec 04 2020

Formula

a(n) = A036469(n) - a(n-1) = Sum_{k=0..n} (-1)^k*A036469(n-k). - Vladeta Jovovic, Sep 10 2003
a(n) = A000009(n) + a(n-2). - Vladeta Jovovic, Feb 10 2004
G.f.: 1/((1-x^2)*Product_{j>=1} (1 - x^(2*j-1))). - Emeric Deutsch, Feb 22 2006
From Vaclav Kotesovec, Aug 16 2015: (Start)
a(n) ~ (1/2) * A036469(n).
a(n) ~ 3^(1/4) * exp(Pi*sqrt(n/3)) / (4*Pi*n^(1/4)). (End)
Euler transform of the sequence [1, 1, period(1, 0)] (A266591). - Georg Fischer, Dec 04 2020

A090858 Number of partitions of n such that there is exactly one part which occurs twice, while all other parts occur only once.

Original entry on oeis.org

0, 0, 1, 0, 2, 2, 2, 4, 6, 7, 8, 13, 15, 21, 25, 30, 39, 50, 58, 74, 89, 105, 129, 156, 185, 221, 264, 309, 366, 433, 505, 593, 696, 805, 941, 1090, 1258, 1458, 1684, 1933, 2225, 2555, 2922, 3346, 3823, 4349, 4961, 5644, 6402, 7267, 8234, 9309, 10525, 11886, 13393
Offset: 0

Views

Author

Vladeta Jovovic, Feb 12 2004

Keywords

Comments

Number of solutions (p(1),p(2),...,p(n)), p(i)>=0,i=1..n, to p(1)+2*p(2)+...+n*p(n)=n such that |{i: p(i)<>0}| = p(1)+p(2)+...+p(n)-1.
Also number of partitions of n such that if k is the largest part, then, with exactly one exception, all the integers 1,2,...,k occur as parts. Example: a(7)=4 because we have [4,2,1], [3,3,1], [3,2,2] and [3,1,1,1,1]. - Emeric Deutsch, Apr 18 2006

Examples

			a(7) = 4 because we have 4 such partitions of 7: [1,1,2,3], [1,1,5], [2,2,3], [1,3,3].
From _Gus Wiseman_, Apr 19 2019: (Start)
The a(2) = 1 through a(11) = 13 partitions described in the name are the following (empty columns not shown). The Heinz numbers of these partitions are given by A060687.
  (11)  (22)   (221)  (33)   (322)   (44)    (441)   (55)    (443)
        (211)  (311)  (411)  (331)   (332)   (522)   (433)   (533)
                             (511)   (422)   (711)   (442)   (551)
                             (3211)  (611)   (3321)  (622)   (722)
                                     (3221)  (4221)  (811)   (911)
                                     (4211)  (4311)  (5221)  (4322)
                                             (5211)  (5311)  (4331)
                                                     (6211)  (4421)
                                                             (5411)
                                                             (6221)
                                                             (6311)
                                                             (7211)
                                                             (43211)
The a(2) = 1 through a(10) = 8 partitions described in Emeric Deutsch's comment are the following (empty columns not shown). The Heinz numbers of these partitions are given by A325284.
  (2)  (22)  (32)   (222)   (322)    (332)     (432)      (3322)
       (31)  (311)  (3111)  (331)    (431)     (3222)     (3331)
                            (421)    (2222)    (4221)     (22222)
                            (31111)  (3311)    (4311)     (42211)
                                     (4211)    (33111)    (43111)
                                     (311111)  (42111)    (331111)
                                               (3111111)  (421111)
                                                          (31111111)
(End)
		

Crossrefs

Programs

  • Maple
    g:=sum(x^(k*(k+1)/2)*((1-x^k)/x^(k-1)/(1-x)-k)/product(1-x^i,i=1..k),k=1..15): gser:=series(g,x=0,64): seq(coeff(gser,x,n),n=1..54); # Emeric Deutsch, Apr 18 2006
    # second Maple program:
    b:= proc(n, i, t) option remember; `if`(n>i*(i+3-2*t)/2, 0,
         `if`(n=0, t, b(n, i-1, t)+`if`(i>n, 0, b(n-i, i-1, t)+
         `if`(t=1 or 2*i>n, 0, b(n-2*i, i-1, 1)))))
        end:
    a:= n-> b(n$2, 0):
    seq(a(n), n=0..100);  # Alois P. Heinz, Dec 28 2015
  • Mathematica
    b[n_, i_, t_] := b[n, i, t] = If[n > i*(i + 3 - 2*t)/2, 0, If[n == 0, t, b[n, i - 1, t] + If[i > n, 0,  b[n - i, i - 1, t] + If[t == 1 || 2*i > n, 0, b[n - 2*i, i - 1, 1]]]]]; a[n_] := b[n, n, 0]; Table[a[n], {n, 0, 100} ] (* Jean-François Alcover, Jan 20 2016, after Alois P. Heinz *)
    Table[Length[Select[IntegerPartitions[n],Length[#]-Length[Union[#]]==1&]],{n,0,30}] (* Gus Wiseman, Apr 19 2019 *)
  • PARI
    alist(n)=concat([0,0],Vec(sum(k=1,n\2,(x^(2*k)+x*O(x^n))/(1+x^k)*prod(j=1,n-2*k,1+x^j+x*O(x^n))))) \\ Franklin T. Adams-Watters, Nov 02 2015

Formula

G.f.: Sum_{k>0} x^(2*k)/(1+x^k) * Product_{k>0} (1+x^k). Convolution of 1-A048272(n) and A000009(n). a(n) = A036469(n) - A015723(n).
G.f.: sum(x^(k(k+1)/2)[(1-x^k)/x^(k-1)/(1-x)-k]/product(1-x^i,i=1..k), k=1..infinity). - Emeric Deutsch, Apr 18 2006
a(n) ~ c * exp(Pi*sqrt(n/3)) / n^(1/4), where c = 3^(1/4) * (1 - log(2)) / (2*Pi) = 0.064273294789... - Vaclav Kotesovec, May 24 2018

Extensions

More terms from Pab Ter (pabrlos(AT)yahoo.com), May 26 2004
a(0) added by Franklin T. Adams-Watters, Nov 02 2015

A346520 Number A(n,k) of partitions of the (n+k)-multiset {0,...,0,1,2,...,k} with n 0's into distinct multisets; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 2, 2, 1, 5, 5, 3, 2, 15, 15, 9, 5, 2, 52, 52, 31, 16, 7, 3, 203, 203, 120, 59, 25, 10, 4, 877, 877, 514, 244, 100, 38, 14, 5, 4140, 4140, 2407, 1112, 442, 161, 56, 19, 6, 21147, 21147, 12205, 5516, 2134, 750, 249, 80, 25, 8, 115975, 115975, 66491, 29505, 11147, 3799, 1213, 372, 111, 33, 10
Offset: 0

Views

Author

Alois P. Heinz, Jul 21 2021

Keywords

Comments

Also number A(n,k) of factorizations of 2^n * Product_{i=1..k} prime(i+1) into distinct factors; A(3,1) = 5: 2*3*4, 4*6, 3*8, 2*12, 24; A(1,2) = 5: 2*3*5, 5*6, 3*10, 2*15, 30.

Examples

			A(2,2) = 9: 00|1|2, 001|2, 1|002, 0|01|2, 0|1|02, 01|02, 00|12, 0|012, 0012.
Square array A(n,k) begins:
  1,  1,   2,   5,   15,    52,   203,    877,    4140, ...
  1,  2,   5,  15,   52,   203,   877,   4140,   21147, ...
  1,  3,   9,  31,  120,   514,  2407,  12205,   66491, ...
  2,  5,  16,  59,  244,  1112,  5516,  29505,  168938, ...
  2,  7,  25, 100,  442,  2134, 11147,  62505,  373832, ...
  3, 10,  38, 161,  750,  3799, 20739, 121141,  752681, ...
  4, 14,  56, 249, 1213,  6404, 36332, 220000, 1413937, ...
  5, 19,  80, 372, 1887, 10340, 60727, 379831, 2516880, ...
  6, 25, 111, 539, 2840, 16108, 97666, 629346, 4288933, ...
  ...
		

Crossrefs

Main diagonal gives A346519.
Antidiagonal sums give A346521.

Programs

  • Maple
    g:= proc(n) option remember; `if`(n=0, 1, add(g(n-j)*add(
         `if`(d::odd, d, 0), d=numtheory[divisors](j)), j=1..n)/n)
        end:
    s:= proc(n) option remember; expand(`if`(n=0, 1,
          x*add(s(n-j)*binomial(n-1, j-1), j=1..n)))
        end:
    S:= proc(n, k) option remember; coeff(s(n), x, k) end:
    b:= proc(n, i) option remember; `if`(n=0, 1,
         `if`(i=0, g(n), add(b(n-j, i-1), j=0..n)))
        end:
    A:= (n, k)-> add(S(k, j)*b(n, j), j=0..k):
    seq(seq(A(n, d-n), n=0..d), d=0..12);
  • Mathematica
    g[n_] := g[n] = If[n == 0, 1, Sum[g[n - j]*Sum[If[OddQ[d], d, 0], {d, Divisors[j]}], {j, 1, n}]/n];
    s[n_] := s[n] = Expand[If[n == 0, 1, x*Sum[s[n - j]*Binomial[n - 1, j - 1], {j, 1, n}]]];
    S[n_, k_] := S[n, k] = Coefficient[s[n], x, k];
    b[n_, i_] := b[n, i] = If[n == 0, 1, If[i == 0, g[n], Sum[b[n - j, i - 1], {j, 0, n}]]];
    A[n_, k_] := Sum[S[k, j]*b[n, j], {j, 0, k}];
    Table[Table[A[n, d - n], {n, 0, d}], {d, 0, 12}] // Flatten (* Jean-François Alcover, Jul 31 2021, after Alois P. Heinz *)

Formula

A(n,k) = A045778(A000079(n)*A070826(k+1)).
A(n,k) = Sum_{j=0..k} Stirling2(k,j)*Sum_{i=0..n} binomial(j+i-1,i)*A000009(n-i).

A087787 a(n) = Sum_{k=0..n} (-1)^(n-k)*A000041(k).

Original entry on oeis.org

1, 0, 2, 1, 4, 3, 8, 7, 15, 15, 27, 29, 48, 53, 82, 94, 137, 160, 225, 265, 362, 430, 572, 683, 892, 1066, 1370, 1640, 2078, 2487, 3117, 3725, 4624, 5519, 6791, 8092, 9885, 11752, 14263, 16922, 20416, 24167, 29007, 34254, 40921, 48213, 57345, 67409
Offset: 0

Views

Author

Vladeta Jovovic, Oct 07 2003

Keywords

Comments

Essentially first differences of A024786 (see the formula). Also, a(n) is the number of 2's in the last section of the set of partitions of n+2 (see A135010). - Omar E. Pol, Sep 10 2008
From Gus Wiseman, May 20 2024: (Start)
Also the number of integer partitions of n containing an even number of ones, ranked by A003159. The a(0) = 1 through a(8) = 15 partitions are:
() . (2) (3) (4) (5) (6) (7) (8)
(11) (22) (32) (33) (43) (44)
(211) (311) (42) (52) (53)
(1111) (222) (322) (62)
(411) (511) (332)
(2211) (3211) (422)
(21111) (31111) (611)
(111111) (2222)
(3311)
(4211)
(22211)
(41111)
(221111)
(2111111)
(11111111)
Also the number of integer partitions of n + 1 containing an odd number of ones, ranked by A036554.
(End)

Crossrefs

The unsigned version is A000070, strict A036469.
For powers of 2 instead number of partitions we have A001045.
The strict or odd version is A025147 or A096765.
The ordered version (compositions instead of partitions) is A078008.
For powers of 2 instead of -1 we have A259401, cf. A259400.
A002865 counts partitions with no ones, column k=0 of A116598.
A072233 counts partitions by sum and length.

Programs

  • Mathematica
    Table[Sum[(-1)^(n-k)*PartitionsP[k], {k,0,n}], {n,0,50}] (* Vaclav Kotesovec, Aug 16 2015 *)
    (* more efficient program *) sig = 1; su = 1; Flatten[{1, Table[sig = -sig; su = su + sig*PartitionsP[n]; Abs[su], {n, 1, 50}]}] (* Vaclav Kotesovec, Nov 06 2016 *)
    Table[Length[Select[IntegerPartitions[n], EvenQ[Count[#,1]]&]],{n,0,30}] (* Gus Wiseman, May 20 2024 *)
  • Python
    from sympy import npartitions
    def A087787(n): return sum(-npartitions(k) if n-k&1 else npartitions(k) for k in range(n+1)) # Chai Wah Wu, Oct 25 2023

Formula

G.f.: 1/(1+x)*1/Product_{k>0} (1-x^k).
a(n) = 1/n*Sum_{k=1..n} (sigma(k)+(-1)^k)*a(n-k).
a(n) = A024786(n+2)-A024786(n+1). - Omar E. Pol, Sep 10 2008
a(n) ~ exp(Pi*sqrt(2*n/3)) / (8*sqrt(3)*n) * (1 + (11*Pi/(24*sqrt(6)) - sqrt(3/2)/Pi)/sqrt(n) - (11/16 + (23*Pi^2)/6912)/n). - Vaclav Kotesovec, Nov 05 2016
a(n) = A000041(n) - a(n-1). - Jon Maiga, Aug 29 2019
Alternating partial sums of A000041. - Gus Wiseman, May 20 2024

A054242 Triangle read by rows: row n (n>=0) gives the number of partitions of (n,0), (n-1,1), (n-2,2), ..., (0,n) respectively into sums of distinct pairs.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 2, 3, 3, 2, 2, 5, 5, 5, 2, 3, 7, 9, 9, 7, 3, 4, 10, 14, 17, 14, 10, 4, 5, 14, 21, 27, 27, 21, 14, 5, 6, 19, 31, 42, 46, 42, 31, 19, 6, 8, 25, 44, 64, 74, 74, 64, 44, 25, 8, 10, 33, 61, 93, 116, 123, 116, 93, 61, 33, 10
Offset: 0

Views

Author

Marc LeBrun, Feb 08 2000 and Jul 01 2003

Keywords

Comments

By analogy with ordinary partitions into distinct parts (A000009). The empty partition gives T(0,0)=1 by definition. A054225 and A201376 give pair partitions with repeats allowed.
Also number of partitions into pairs which are not both even.
In the paper by S. M. Luthra: "Partitions of bipartite numbers when the summands are unequal", the square table on page 370 contains an errors. In the formula (6, p. 372) for fixed m there should be factor 1/m!. The correct asymptotic formula is q(m, n) ~ (sqrt(12*n)/Pi)^m * exp(Pi*sqrt(n/3)) / (4*3^(1/4)*m!*n^(3/4)). The same error is also in article by F. C. Auluck (see A054225). - Vaclav Kotesovec, Feb 02 2016

Examples

			The second row (n=1) is 1,1 since (1,0) and (0,1) each have a single partition.
The third row (n=2) is 1, 2, 1 from (2,0), (1,1) or (1,0)+(0,1), (0,2).
In the fourth row, T(1,3)=5 from (1,3), (0,3)+(1,0), (0,2)+(1,1), (0,2)+(0,1)+(1,0), (0,1)+(1,2).
The triangle begins:
  1;
  1,  1;
  1,  2,  1;
  2,  3,  3,  2;
  2,  5,  5,  5,  2;
  3,  7,  9,  9,  7,  3;
  4, 10, 14, 17, 14, 10,  4;
  5, 14, 21, 27, 27, 21, 14,  5;
  6, 19, 31, 42, 46, 42, 31, 19,  6;
  8, 25, 44, 64, 74, 74, 64, 44, 25, 8;
  ...
		

Crossrefs

See A201377 for the same triangle formatted in a different way.
The outer diagonals are T(n,0) = T(n,n) = A000009(n).
Cf. A054225.
T(2*n,n) = A219554(n). Row sums give A219555. - Alois P. Heinz, Nov 22 2012

Programs

  • Haskell
    see Zumkeller link.
  • Mathematica
    max = 10; f[x_, y_] := Product[1 + x^n*y^k, {n, 0, max}, {k, 0, max}]/2; se = Series[f[x, y], {x, 0, max}, {y, 0, max}] ; coes = CoefficientList[ se, {x, y}]; t[n_, k_] := coes[[n-k+1, k+1]]; Flatten[ Table[ t[n, k], {n, 0, max}, {k, 0, n}]] (* Jean-François Alcover, Dec 06 2011 *)

Formula

G.f.: (1/2)*Product(1+x^i*y^j), i, j>=0.

Extensions

Entry revised by N. J. A. Sloane, Nov 30 2011, to incorporate corrections provided by Reinhard Zumkeller, who also contributed the alternative version A201377.

A342086 Number of strict factorizations of divisors of n.

Original entry on oeis.org

1, 2, 2, 3, 2, 5, 2, 5, 3, 5, 2, 9, 2, 5, 5, 7, 2, 9, 2, 9, 5, 5, 2, 16, 3, 5, 5, 9, 2, 15, 2, 10, 5, 5, 5, 18, 2, 5, 5, 16, 2, 15, 2, 9, 9, 5, 2, 25, 3, 9, 5, 9, 2, 16, 5, 16, 5, 5, 2, 31, 2, 5, 9, 14, 5, 15, 2, 9, 5, 15, 2, 34, 2, 5, 9, 9, 5, 15, 2, 25, 7, 5
Offset: 1

Views

Author

Gus Wiseman, Mar 05 2021

Keywords

Comments

A strict factorization of n is a set of distinct positive integers > 1 with product n.

Examples

			The a(1) = 1 through a(12) = 9 factorizations:
  ()  ()   ()   ()   ()   ()     ()   ()     ()   ()     ()    ()
      (2)  (3)  (2)  (5)  (2)    (7)  (2)    (3)  (2)    (11)  (2)
                (4)       (3)         (4)    (9)  (5)          (3)
                          (6)         (8)         (10)         (4)
                          (2*3)       (2*4)       (2*5)        (6)
                                                               (12)
                                                               (2*3)
                                                               (2*6)
                                                               (3*4)
		

Crossrefs

A version for partitions is A026906 (strict partitions of 1..n).
A version for partitions is A036469 (strict partitions of 0..n).
A version for partitions is A047966 (strict partitions of divisors).
The non-strict version is A057567.
A000005 counts divisors, with sum A000203.
A000009 counts strict partitions.
A001055 counts factorizations, with strict case A045778.
A001221 counts prime divisors, with sum A001414.
A001222 counts prime-power divisors.
A005117 lists squarefree numbers.

Programs

  • Maple
    sf1:= proc(n,m)
      local D,d;
      if n = 1 then return 1 fi;
      D:= select(`<`,numtheory:-divisors(n) minus {1},m);
      add( procname(n/d,d), d= D)
    end proc:
    sf:= proc(n) option remember; sf1(n,n+1) end proc:f:= proc(n) local d; add(sf(d),d=numtheory:-divisors(n)) end proc:map(f, [$1..100]); # Robert Israel, Mar 10 2021
  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Sum[Length[Select[facs[k],UnsameQ@@#&]],{k,Divisors[n]}],{n,30}]

A340793 Sequence whose partial sums give A000203.

Original entry on oeis.org

1, 2, 1, 3, -1, 6, -4, 7, -2, 5, -6, 16, -14, 10, 0, 7, -13, 21, -19, 22, -10, 4, -12, 36, -29, 11, -2, 16, -26, 42, -40, 31, -15, 6, -6, 43, -53, 22, -4, 34, -48, 54, -52, 40, -6, -6, -24, 76, -67, 36, -21, 26, -44, 66, -48, 48, -40, 10, -30, 108, -106, 34, 8
Offset: 1

Views

Author

Omar E. Pol, Jan 21 2021

Keywords

Comments

Essentially a duplicate of A053222.
Convolved with the nonzero terms of A000217 gives A175254, the volume of the stepped pyramid described in A245092.
Convolved with the nonzero terms of A046092 gives A244050, the volume of the stepped pyramid described in A244050.
Convolved with A000027 gives A024916.
Convolved with A000041 gives A138879.
Convolved with A000070 gives the nonzero terms of A066186.
Convolved with the nonzero terms of A002088 gives A086733.
Convolved with A014153 gives A182738.
Convolved with A024916 gives A000385.
Convolved with A036469 gives the nonzero terms of A277029.
Convolved with A091360 gives A276432.
Convolved with A143128 gives the nonzero terms of A000441.
For the correspondence between divisors and partitions see A336811.

Crossrefs

Programs

  • Maple
    a:= n-> (s-> s(n)-s(n-1))(numtheory[sigma]):
    seq(a(n), n=1..77);  # Alois P. Heinz, Jan 21 2021
  • Mathematica
    Join[{1}, Differences @ Table[DivisorSigma[1, n], {n, 1, 100}]] (* Amiram Eldar, Jan 21 2021 *)
  • PARI
    a(n) = if (n==1, 1, sigma(n)-sigma(n-1)); \\ Michel Marcus, Jan 22 2021

Formula

a(n) = A053222(n-1) for n>1. - Michel Marcus, Jan 22 2021
Showing 1-10 of 49 results. Next