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

A161919 Permutation of natural numbers: concatenation of subsequences A161924(A000070(k-1)..A026905(k)), k >= 1, each sorted into ascending order.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 6, 8, 9, 11, 15, 10, 13, 16, 17, 19, 23, 31, 12, 14, 18, 21, 27, 32, 33, 35, 39, 47, 63, 20, 22, 25, 29, 34, 37, 43, 55, 64, 65, 67, 71, 79, 95, 127, 24, 26, 30, 36, 38, 41, 45, 51, 59, 66, 69, 75, 87, 111, 128, 129, 131, 135, 143, 159, 191, 255, 28, 40
Offset: 1

Views

Author

Alford Arnold, Jun 23 2009

Keywords

Comments

This is the lexicographically earliest sequence a_n for which it holds that A161511(a(n)) = A036042(n) for all n.
Triangle T(n,k) read by rows. Row n lists in increasing order the viabin numbers of the integer partitions of n (n >= 1, k >= 1). The viabin number of an integer partition is defined in the following way. Consider the southeast border of the Ferrers board of the integer partition and consider the binary number obtained by replacing each east step with 1 and each north step, except the last one, with 0. The corresponding decimal form is, by definition, the viabin number of the given integer partition. "Viabin" is coined from "via binary". For example, consider the integer partition [3,1,1] of 5. The southeast border of its Ferrers board yields 10011, leading to the viabin number 19 (an entry in the 5th row). - Emeric Deutsch, Sep 06 2017
After specifying the value of n, the first Maple program yields the entries of row n. - Emeric Deutsch, Feb 26 2016
After specifying the value of m, the third Maple program yields the first m rows; the command partovi(p) yields the viabin number of the partition p = [a,b,c,...]. - Emeric Deutsch, Aug 31 2017

Examples

			This can be viewed as an irregular table, where row r (>= 1) has A000041(r) elements, i.e., as 1; 2,3; 4,5,7; 6,8,9,11,15; 10,13,16,17,19,23,31; etc. A125106 illustrates how each number is mapped to a partition.
		

Crossrefs

Inverse: A166277. Sequence A161924 gives the same rows before sorting.

Programs

  • Maple
    n := 11: s := proc (b) local t, i, j: t := 0: for i to nops(b) do for j from i+1 to nops(b) do if b[j]-b[i] = 1 then t := t+1 else  end if end do end do: t end proc: A[n] := {}: for i to 2^n do a[i] := convert(2*i, base, 2) end do: for k to 2^n do if s(a[k]) = n then A[n] := `union`(A[n], {k}) else  end if end do: A[n]; # Emeric Deutsch, Feb 26 2016
    # second Maple program:
    f:= proc(l) local i, r; r:= 0; for i to nops(l)-1 do
           r:= 2*((x-> 2*x+1)@@(l[i+1]-l[i]))(r) od; r/2
        end:
    b:= proc(n, i) `if`(n=0 or i=1, [[0, 1$n]], [b(n, i-1)[],
          `if`(i>n, [], map(x-> [x[], i], b(n-i, i)))[]])
        end:
    T:= n-> sort(map(f, b(n$2)))[]:
    seq(T(n), n=1..10);  # Alois P. Heinz, Jul 25 2017
    # 3rd Maple program:
    m := 10; with(combinat): ff := proc (X) local s: s := [1, seq(0, j = 1 .. X[2])]: s := map(convert, s, string): return cat(op(s)) end proc: partovi := proc (P) local X, n, Y, i: X := convert(P, multiset): n := X[-1][1]: Y := map(proc (t) options operator, arrow: t[1] end proc, X): for i to n do if member(i, Y) = false then X := [op(X), [i, 0]] end if end do: X := sort(X, proc (s, t) options operator, arrow: evalb(s[1] < t[1]) end proc): X := map(ff, X): X := cat(op(X)): n := parse(X): n := convert(n, decimal, binary): (1/2)*n end proc: for n to m do {seq(partovi(partition(n)[q]), q = 1 .. numbpart(n))} end do; # Emeric Deutsch, Aug 31 2017
  • Mathematica
    columns = 10;
    row[n_] := n - 2^Floor[Log2[n]];
    col[0] = 0; col[n_] := If[EvenQ[n], col[n/2] + DigitCount[n/2, 2, 1], col[(n-1)/2] + 1];
    Clear[T]; T[, ] = 0; Do[T[row[k], col[k]] = k, {k, 1, 2^columns}];
    Table[DeleteCases[Sort @ Table[T[n-1, k], {n, 1, 2^(k-1)}], 0], {k, 1, columns}] // Flatten (* Jean-François Alcover, Feb 16 2021 *)

Extensions

Edited and extended by Antti Karttunen, Oct 12 2009

A134286 Characteristic sequence for sequence A026905.

Original entry on oeis.org

1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 1

Views

Author

Wolfdieter Lang, Nov 13 2007

Keywords

Comments

This partition array is the member k=1 in the family M_0(k), with M_0(2)=M_0= A048996, M_0(3)= A134283, etc.
When read as partition array (tabf with sequence of row lengths given by the partition numbers A000041) in Abramowitz-Stegun order (see A117506 for the reference) a(n,k) is the characteristic partition array for the partition (1^n) of n.

Programs

  • Mathematica
    terms = 105; nmax = 10;
    pp = PartitionsP[Range[nmax]] // Accumulate;
    a[n_] := If[n > pp[[-1]], Print["nmax = ", nmax, " too small"], Boole[ MemberQ[ pp, n]]];
    Array[a, terms] (* Jean-François Alcover, Jun 19 2019 *)

Formula

a(n)=1 if n from A026905, else 0.

A085360 Partial sums of A026905; the convolution of the natural numbers with the partition function.

Original entry on oeis.org

1, 4, 10, 21, 39, 68, 112, 178, 274, 412, 606, 877, 1249, 1756, 2439, 3353, 4564, 6160, 8246, 10959, 14464, 18971, 24733, 32070, 41365, 53096, 67837, 86296, 109320, 137948
Offset: 1

Views

Author

Jon Perry, Jun 25 2003

Keywords

Examples

			a(4) = A026905(1) + A026905(2) + A026905(3) + A026905(4) = 1 + 3 + 6 + 11 = 21.
		

Crossrefs

Cf. A026905.

Programs

Formula

a(n) = A086716(n) - A086716(n-1). - Vaclav Kotesovec, Jun 23 2015
a(n) ~ sqrt(3) * exp(Pi*sqrt(2*n/3)) / (2*Pi^2). - Vaclav Kotesovec, Jun 23 2015

A056870 Difference between partial sums of partition numbers (A026905) and partial sums of numbers of partitions into distinct parts (A026906).

Original entry on oeis.org

0, 1, 2, 5, 9, 16, 26, 42, 64, 96, 140, 202, 285, 398, 547, 746, 1005, 1344, 1780, 2343, 3059, 3972, 5123, 6576, 8392, 10663, 13481, 16977, 21286, 26594, 33096, 41055, 50750, 62548, 76846, 94155, 115032, 140183, 170386, 206611, 249934, 301682
Offset: 1

Views

Author

Robert G. Wilson v, Sep 02 2000

Keywords

Crossrefs

Programs

  • Mathematica
    Table[ Sum[ PartitionsP[k], {k, 1, n}] - Sum[PartitionsQ[k], {k, 1, n}], {n, 1, 45}]
    (* second program *)
    Module[{nn=50},#[[1]]-#[[2]]&/@Thread[{Accumulate[PartitionsP[Range[nn]]],Accumulate[PartitionsQ[Range[nn]]]}]] (* Harvey P. Dale, Jun 19 2025 *)

A056871 Sum of partial sums of partition numbers (A026905) and partial sums of numbers of partitions into distinct parts (A026906).

Original entry on oeis.org

2, 5, 10, 17, 27, 42, 62, 90, 128, 180, 248, 340, 459, 616, 819, 1082, 1417, 1848, 2392, 3083, 3951, 5042, 6401, 8098, 10198, 12799, 16001, 19941, 24762, 30662, 37844, 46583, 57174, 69996, 85464, 104109, 126506, 153385, 185552, 224003, 269846
Offset: 1

Views

Author

Robert G. Wilson v, Sep 02 2000

Keywords

Comments

Table[ Sum[ PartitionsP[k], {k, 1, n}] + Sum[PartitionsQ[k], {k, 1, n}], {n, 1, 45}]

Crossrefs

A086741 a(n) = A026905(n) - A014284(n).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 2, 7, 18, 37, 64, 110, 174, 268, 401, 585, 829, 1155, 1584, 2144, 2865, 3794, 4970, 6462, 8331, 10670, 13579, 17194, 21652, 27147, 33876, 42098, 52110, 64283, 79027, 96855, 118341, 144199, 175221, 212392, 256802, 309797, 372877, 447861, 536802, 642163, 766718, 913780
Offset: 1

Views

Author

Jon Perry, Jul 29 2003

Keywords

Comments

Partial sums of partition numbers - partial sums of primes, if 1 is considered as a prime.

Crossrefs

Programs

  • PARI
    a(n) = sum(k=1, n, numbpart(k)) - (1 + vecsum(primes(n-1))); \\ Michel Marcus, Aug 10 2023

Extensions

More terms from Michel Marcus, Aug 10 2023

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

A006128 Total number of parts in all partitions of n. Also, sum of largest parts of all partitions of n.

Original entry on oeis.org

0, 1, 3, 6, 12, 20, 35, 54, 86, 128, 192, 275, 399, 556, 780, 1068, 1463, 1965, 2644, 3498, 4630, 6052, 7899, 10206, 13174, 16851, 21522, 27294, 34545, 43453, 54563, 68135, 84927, 105366, 130462, 160876, 198014, 242812, 297201, 362587, 441546, 536104, 649791, 785437, 947812, 1140945, 1371173, 1644136, 1968379, 2351597, 2805218, 3339869, 3970648, 4712040, 5584141, 6606438, 7805507, 9207637
Offset: 0

Views

Author

Keywords

Comments

a(n) = degree of Kac determinant at level n as polynomial in the conformal weight (called h). (Cf. C. Itzykson and J.-M. Drouffe, Statistical Field Theory, Vol. 2, p. 533, eq.(98); reference p. 643, Cambridge University Press, (1989).) - Wolfdieter Lang
Also the number of one-element transitions from the integer partitions of n to the partitions of n-1 for labeled parts with the assumption that from any part z > 1 one can take an element of amount 1 in one way only. That means z is composed of z unlabeled parts of amount 1, i.e. z = 1 + 1 + ... + 1. E.g., for n=3 to n=2 we have a(3) = 6 and [111] --> [11], [111] --> [11], [111] --> [11], [12] --> [11], [12] --> [2], [3] --> [2]. For the case of z composed by labeled elements, z = 1_1 + 1_2 + ... + 1_z, see A066186. - Thomas Wieder, May 20 2004
Number of times a derivative of any order (not 0 of course) appears when expanding the n-th derivative of 1/f(x). For instance (1/f(x))'' = (2 f'(x)^2-f(x) f''(x)) / f(x)^3 which makes a(2) = 3 (by counting k times the k-th power of a derivative). - Thomas Baruchel, Nov 07 2005
Starting with offset 1, = the partition triangle A008284 * [1, 2, 3, ...]. - Gary W. Adamson, Feb 13 2008
Starting with offset 1 equals A000041: (1, 1, 2, 3, 5, 7, 11, ...) convolved with A000005: (1, 2, 2, 3, 2, 4, ...). - Gary W. Adamson, Jun 16 2009
Apart from initial 0 row sums of triangle A066633, also the Möbius transform is A085410. - Gary W. Adamson, Mar 21 2011
More generally, the total number of parts >= k in all partitions of n equals the sum of k-th largest parts of all partitions of n. In this case k = 1. Apart from initial 0 the first column of A181187. - Omar E. Pol, Feb 14 2012
Row sums of triangle A221530. - Omar E. Pol, Jan 21 2013
From Omar E. Pol, Feb 04 2021: (Start)
a(n) is also the total number of divisors of all positive integers in a sequence with n blocks where the m-th block consists of A000041(n-m) copies of m, with 1 <= m <= n. The mentioned divisors are also all parts of all partitions of n.
Apart from initial zero this is also as follows:
Convolution of A000005 and A000041.
Convolution of A006218 and A002865.
Convolution of A341062 and A000070.
Row sums of triangles A221531, A245095, A339258, A340525, A340529. (End)
Number of ways to choose a part index of an integer partition of n, i.e., partitions of n with a selected position. Selecting a part value instead of index gives A000070. - Gus Wiseman, Apr 19 2021

Examples

			For n = 4 the partitions of 4 are [4], [2, 2], [3, 1], [2, 1, 1], [1, 1, 1, 1]. The total number of parts is 12. On the other hand, the sum of the largest parts of all partitions is 4 + 2 + 3 + 2 + 1 = 12, equaling the total number of parts, so a(4) = 12. - _Omar E. Pol_, Oct 12 2018
		

References

  • S. M. Luthra, On the average number of summands in partitions of n, Proc. Nat. Inst. Sci. India Part. A, 23 (1957), p. 483-498.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Main diagonal of A210485.
Column k=1 of A256193.
The version for normal multisets is A001787.
The unordered version is A001792.
The strict case is A015723.
The version for factorizations is A066637.
A000041 counts partitions.
A000070 counts partitions with a selected part.
A336875 counts compositions with a selected part.
A339564 counts factorizations with a selected factor.

Programs

  • GAP
    List([0..60],n->Length(Flat(Partitions(n)))); # Muniru A Asiru, Oct 12 2018
  • Haskell
    a006128 = length . concat . ps 1 where
       ps _ 0 = [[]]
       ps i j = [t:ts | t <- [i..j], ts <- ps t (j - t)]
    -- Reinhard Zumkeller, Jul 13 2013
    
  • Maple
    g:= add(n*x^n*mul(1/(1-x^k), k=1..n), n=1..61):
    a:= n-> coeff(series(g,x,62),x,n):
    seq(a(n), n=0..61);
    # second Maple program:
    a:= n-> add(combinat[numbpart](n-j)*numtheory[tau](j), j=1..n):
    seq(a(n), n=0..61);  # Alois P. Heinz, Aug 23 2019
  • Mathematica
    a[n_] := Sum[DivisorSigma[0, m] PartitionsP[n - m], {m, 1, n}]; Table[ a[n], {n, 0, 41}]
    CoefficientList[ Series[ Sum[n*x^n*Product[1/(1 - x^k), {k, n}], {n, 100}], {x, 0, 100}], x]
    a[n_] := Plus @@ Max /@ IntegerPartitions@ n; Array[a, 45] (* Robert G. Wilson v, Apr 12 2011 *)
    Join[{0}, ((Log[1 - x] + QPolyGamma[1, x])/(Log[x] QPochhammer[x]) + O[x]^60)[[3]]] (* Vladimir Reshetnikov, Nov 17 2016 *)
    Length /@ Table[IntegerPartitions[n] // Flatten, {n, 50}] (* Shouvik Datta, Sep 12 2021 *)
  • PARI
    f(n)= {local(v,i,k,s,t);v=vector(n,k,0);v[n]=2;t=0;while(v[1]1,i--;s+=i*(v[i]=(n-s)\i));t+=sum(k=1,n,v[k]));t } /* Thomas Baruchel, Nov 07 2005 */
    
  • PARI
    a(n) = sum(m=1, n, numdiv(m)*numbpart(n-m)) \\ Michel Marcus, Jul 13 2013
    
  • Python
    from sympy import divisor_count, npartitions
    def a(n): return sum([divisor_count(m)*npartitions(n - m) for m in range(1, n + 1)]) # Indranil Ghosh, Apr 25 2017
    

Formula

G.f.: Sum_{n>=1} n*x^n / Product_{k=1..n} (1-x^k).
G.f.: Sum_{k>=1} x^k/(1-x^k) / Product_{m>=1} (1-x^m).
a(n) = Sum_{k=1..n} k*A008284(n, k).
a(n) = Sum_{m=1..n} of the number of divisors of m * number of partitions of n-m.
Note that the formula for the above comment is a(n) = Sum_{m=1..n} d(m)*p(n-m) = Sum_{m=1..n} A000005(m)*A000041(n-m), if n >= 1. - Omar E. Pol, Jan 21 2013
Erdős and Lehner show that if u(n) denotes the average largest part in a partition of n, then u(n) ~ constant*sqrt(n)*log n.
a(n) = A066897(n) + A066898(n), n>0. - Reinhard Zumkeller, Mar 09 2012
a(n) = A066186(n) - A196087(n), n >= 1. - Omar E. Pol, Apr 22 2012
a(n) = A194452(n) + A024786(n+1). - Omar E. Pol, May 19 2012
a(n) = A000203(n) + A220477(n). - Omar E. Pol, Jan 17 2013
a(n) = Sum_{m=1..p(n)} A194446(m) = Sum_{m=1..p(n)} A141285(m), where p(n) = A000041(n), n >= 1. - Omar E. Pol, May 12 2013
a(n) = A198381(n) + A026905(n), n >= 1. - Omar E. Pol, Aug 10 2013
a(n) = O(sqrt(n)*log(n)*p(n)), where p(n) is the partition function A000041(n). - Peter Bala, Dec 23 2013
a(n) = Sum_{m=1..n} A006218(m)*A002865(n-m), n >= 1. - Omar E. Pol, Jul 14 2014
From Vaclav Kotesovec, Jun 23 2015: (Start)
Asymptotics (Luthra, 1957): a(n) = p(n) * (C*N^(1/2) + C^2/2) * (log(C*N^(1/2)) + gamma) + (1+C^2)/4 + O(N^(-1/2)*log(N)), where N = n - 1/24, C = sqrt(6)/Pi, gamma is the Euler-Mascheroni constant A001620 and p(n) is the partition function A000041(n).
The formula a(n) = p(n) * (sqrt(3*n/(2*Pi)) * (log(n) + 2*gamma - log(Pi/6)) + O(log(n)^3)) in the abstract of the article by Kessler and Livingston (cited also in the book by Sandor, p. 495) is incorrect!
Right is: a(n) = p(n) * (sqrt(3*n/2)/Pi * (log(n) + 2*gamma - log(Pi^2/6)) + O(log(n)^3))
or a(n) ~ exp(Pi*sqrt(2*n/3)) * (log(6*n/Pi^2) + 2*gamma) / (4*Pi*sqrt(2*n)).
(End)
a(n) = Sum_{m=1..n} A341062(m)*A000070(n-m), n >= 1. - Omar E. Pol, Feb 05 2021 2014

A363225 Number of integer partitions of n containing three parts (a,b,c) (repeats allowed) such that a + b = c. A variation of sum-full partitions.

Original entry on oeis.org

0, 0, 0, 1, 1, 2, 4, 5, 9, 14, 21, 29, 43, 58, 81, 109, 148, 195, 263, 339, 445, 574, 744, 942, 1209, 1515, 1923, 2399, 3005, 3721, 4629, 5693, 7024, 8589, 10530, 12804, 15596, 18876, 22870, 27538, 33204, 39816, 47766, 57061, 68161, 81099, 96510, 114434, 135634
Offset: 0

Views

Author

Gus Wiseman, Jul 19 2023

Keywords

Comments

Note that, by this definition, the partition (2,1) is sum-full, because (1,1,2) is a triple satisfying a + b = c.

Examples

			The a(3) = 1 through a(9) = 14 partitions:
  (21)  (211)  (221)   (42)     (421)     (422)      (63)
               (2111)  (321)    (2221)    (431)      (432)
                       (2211)   (3211)    (521)      (621)
                       (21111)  (22111)   (3221)     (3321)
                                (211111)  (4211)     (4221)
                                          (22211)    (4311)
                                          (32111)    (5211)
                                          (221111)   (22221)
                                          (2111111)  (32211)
                                                     (42111)
                                                     (222111)
                                                     (321111)
                                                     (2211111)
                                                     (21111111)
		

Crossrefs

For subsets of {1..n} we have A093971, A088809 without re-using parts.
The complement for subsets is A007865, A085489 without re-using parts.
Without re-using parts we have A237113, complement A236912.
For sums of any length > 1 (without re-usable parts) we have A237668, complement A237667.
The strict case is A363226.
The complement is counted by A364345, strict A364346.
These partitions have ranks A364348, complement A364347.
The strict linear combination-free version is A364350.
A000041 counts partitions, strict A000009.
A008284 counts partitions by length, strict A008289.
A323092 counts double-free partitions, ranks A320340.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],Select[Tuples[#,3],#[[1]]+#[[2]]==#[[3]]&]!={}&]],{n,0,15}]
  • Python
    from collections import Counter
    from itertools import combinations_with_replacement
    from sympy.utilities.iterables import partitions
    def A363225(n): return sum(1 for p in partitions(n) if any(q[0]+q[1]==q[2] for q in combinations_with_replacement(sorted(Counter(p).elements()),3))) # Chai Wah Wu, Sep 21 2023

Extensions

a(31)-a(48) from Chai Wah Wu, Sep 21 2023

A185974 Partitions in Abramowitz-Stegun order A036036 mapped one-to-one to positive integers.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 8, 7, 10, 9, 12, 16, 11, 14, 15, 20, 18, 24, 32, 13, 22, 21, 25, 28, 30, 27, 40, 36, 48, 64, 17, 26, 33, 35, 44, 42, 50, 45, 56, 60, 54, 80, 72, 96, 128, 19, 34, 39, 55, 49, 52, 66, 70, 63, 75, 88, 84, 100, 90, 81, 112, 120, 108, 160, 144, 192, 256, 23, 38, 51, 65, 77, 68, 78, 110, 98, 99, 105, 125, 104, 132, 140, 126, 150, 135, 176, 168, 200, 180, 162, 224, 240, 216, 320, 288, 384, 512, 29, 46, 57, 85, 91, 121, 76, 102, 130, 154, 117, 165, 147, 175, 136, 156, 220, 196, 198, 210, 250, 189, 225, 208, 264, 280, 252, 300, 270, 243, 352, 336, 400, 360, 324, 448, 480, 432, 640, 576, 768, 1024
Offset: 0

Views

Author

Wolfdieter Lang, Feb 10 2011

Keywords

Comments

First differs from A334438 (shifted left once) at a(75) = 98, A334438(76) = 99. - Gus Wiseman, May 20 2020
This mapping of the set of all partitions of N >= 0 to {1, 2, 3, ...} (set of natural numbers) is one to one (bijective). The empty partition for N = 0 maps to 1.
A129129 seems to be analogous, except that the partition ordering A080577 is used. This ordering, however, does not care about the number of parts: e.g., 1^2,4 = 4,1^2 comes before 3^2, so a(23)=28 and a(22)=25 are interchanged.
Also Heinz numbers of all reversed integer partitions (finite weakly increasing sequences of positive integers), sorted first by sum, then by length, and finally lexicographically, where the Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). The version for non-reversed partitions is A334433. - Gus Wiseman, May 20 2020

Examples

			a(22) = 25 = prime(3)^2 because the 22nd partition in A-St order is the 2-part partition (3,3) of N = 6, because A026905(5) = 18 < 22 <= A026905(6) = 29.
a(23) = 28 = prime(1)^2*prime(4) corresponds to the partition 1+1+4 = 4+1+1 with three parts, also of N = 6.
From _Gus Wiseman_, May 20 2020: (Start)
Triangle begins:
   1
   2
   3   4
   5   6   8
   7  10   9  12  16
  11  14  15  20  18  24  32
  13  22  21  25  28  30  27  40  36  48  64
  17  26  33  35  44  42  50  45  56  60  54  80  72  96 128
As a triangle of reversed partitions we have:
                             0
                            (1)
                          (2)(11)
                        (3)(12)(111)
                   (4)(13)(22)(112)(1111)
             (5)(14)(23)(113)(122)(1112)(11111)
  (6)(15)(24)(33)(114)(123)(222)(1113)(1122)(11112)(111111)
(End)
		

Crossrefs

Row lengths are A000041.
The constructive version is A036036.
Also Heinz numbers of the partitions in A036037.
The generalization to compositions is A124734.
The version for non-reversed partitions is A334433.
The non-reversed length-insensitive version is A334434.
The opposite version (sum/length/revlex) is A334435.
Ignoring length gives A334437.
Sorting reversed partitions by Heinz number gives A112798.
Partitions in lexicographic order are A193073.
Partitions in colexicographic order are A211992.
Graded Heinz numbers are A215366.

Programs

  • Mathematica
    Join@@Table[Times@@Prime/@#&/@Sort[Reverse/@IntegerPartitions[n]],{n,0,8}] (* Gus Wiseman, May 21 2020 *)
  • PARI
    A185974_row(n)=[vecprod([prime(i)|i<-p])|p<-partitions(n)] \\ below a helper function:
    index_of_partition(n)={for(r=0, oo, my(c = numbpart(r)); n >= c || return([r,n+1]); n -= c)}
    /* A185974(n,k), 1 <= k <= A000041(n), gives the k-th partition of n >= 0; if k is omitted, A185974(n) return the term of index n of the flattened sequence a(n >= 0).
      This function is used in other sequences (such as A122172) which need to access the n-th partition as listed in A-S order. */
    A185974(n, k=index_of_partition(n))=A185974_row(iferr(k[1], E, k=[k,k]; n))[k[2]] \\ (End)

Formula

a(n) = Product_{j=1..N(n)} p(j)^e(j), with p(j):=A000040(j) (j-th prime), and the exponent e(j) >= 0 of the part j in the n-th partition written in Abramowitz-Stegun (A-St) order, indicated in A036036. Note that j^0 is not 1 but has to be omitted in the partition. N(n) is the index (argument) of the smallest A026905-number greater than or equal to n (the index of the A026905-ceiling of n).
From Gus Wiseman, May 21 2020: (Start)
A001221(a(n)) = A103921(n).
A001222(a(n)) = A036043(n).
A056239(a(n)) = A036042(n).
A061395(a(n)) = A049085(n).
(End)

Extensions

Examples edited by M. F. Hasler, Jan 07 2024
Showing 1-10 of 60 results. Next