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

A002865 Number of partitions of n that do not contain 1 as a part.

Original entry on oeis.org

1, 0, 1, 1, 2, 2, 4, 4, 7, 8, 12, 14, 21, 24, 34, 41, 55, 66, 88, 105, 137, 165, 210, 253, 320, 383, 478, 574, 708, 847, 1039, 1238, 1507, 1794, 2167, 2573, 3094, 3660, 4378, 5170, 6153, 7245, 8591, 10087, 11914, 13959, 16424, 19196, 22519, 26252, 30701
Offset: 0

Views

Author

Keywords

Comments

Also the number of partitions of n-1, n >= 2, such that the least part occurs exactly once. See A096373, A097091, A097092, A097093. - Robert G. Wilson v, Jul 24 2004 [Corrected by Wolfdieter Lang, Feb 18 2009]
Number of partitions of n+1 where the number of parts is itself a part. Take a partition of n (with k parts) which does not contain 1, remove 1 from each part and add a new part of size k+1. - Franklin T. Adams-Watters, May 01 2006
Number of partitions where the largest part occurs at least twice. - Joerg Arndt, Apr 17 2011
Row sums of triangle A147768. - Gary W. Adamson, Nov 11 2008
From Lewis Mammel (l_mammel(AT)att.net), Oct 06 2009: (Start)
a(n) is the number of sets of n disjoint pairs of 2n things, called a pairing, disjoint with a given pairing (A053871), that are unique under permutations preserving the given pairing.
Can be seen immediately from a graphical representation which must decompose into even numbered cycles of 4 or more things, as connected by pairs alternating between the pairings. Each thing is in a single cycle, so this is a partition of 2n into even parts greater than 2, equivalent to a partition of n into parts greater than 1. (End)
Convolution product (1, 1, 2, 2, 4, 4, ...) * (1, 2, 3, ...) = A058682 starting (1, 3, 7, 13, 23, 37, ...); with row sums of triangle A171239 = A058682. - Gary W. Adamson, Dec 05 2009
Also the number of 2-regular multigraphs with loops forbidden. - Jason Kimberley, Jan 05 2011
Number of appearances of the multiplicity n, n-1, ..., n-k in all partitions of n, for k < n/2. (Only populated by multiplicities of large numbers of 1's.) - William Keith, Nov 20 2011
Also the number of equivalence classes of n X n binary matrices with exactly 2 1's in each row and column, up to permutations of rows and columns (cf. A133687). - N. J. A. Sloane, Sep 16 2013
Starting at a(2) this sequence gives the number of vertices on a nim tree created in the game of edge removal for a path P_{n} where n is the number of vertices on the path. This is the number of nonisomorphic graphs that can result from the path when the game of edge removal is played. - Lyndsey Wong, Jul 09 2016
The number of different ways to climb a staircase taking at least two stairs at a time. - Mohammad K. Azarian, Nov 20 2016
Let 1,0,1,1,1,... (offset 0) count unlabeled, connected, loopless 1-regular digraphs. This here is the Euler transform of that sequence, counting unlabeled loopless 1-regular digraphs. A145574 is the associated multiset transformation. A000166 are the labeled loopless 1-regular digraphs. - R. J. Mathar, Mar 25 2019
For n > 1, also the number of partitions with no part greater than the number of ones. - George Beck, May 09 2019 [See A187219 which is the correct sequence for this interpretation for n >= 1. - Spencer Miller, Jan 30 2023]
From Gus Wiseman, May 19 2019: (Start)
Conjecture: Also the number of integer partitions of n - 1 that have a consecutive subsequence summing to each positive integer from 1 to n - 1. For example, (32211) is such a partition because we have consecutive subsequences:
1: (1)
2: (2)
3: (3) or (21)
4: (22) or (211)
5: (32) or (221)
6: (2211)
7: (322)
8: (3221)
9: (32211)
(End)
There is a sufficient and necessary condition to characterize the partitions defined by Gus Wiseman. It is that the largest part must be less than or equal to the number of ones plus one. Hence, the number of partitions of n with no part greater than the number of ones is the same as the number of partitions of n-1 that have a consecutive subsequence summing to each integer from 1 to n-1. Gus Wiseman's conjecture can be proved bijectively. - Andrew Yezhou Wang, Dec 14 2019
From Peter Bala, Dec 01 2024: (Start)
Let P(2, n) denote the set of partitions of n into parts k > 1. Then A000041(n) = - Sum_{parts k in all partitions in P(2, n+2)} mu(k). For example, with n = 5, there are 4 partitions of n + 2 = 7 into parts greater than 1, namely, 7, 5 + 2, 4 + 3, 3 + 2 + 2, and mu(7) + (mu(5) + mu(2)) + (mu(4 ) + mu(3)) + (mu(3) + mu(2) + mu(2)) = -7 = - A000041(5). (End)

Examples

			a(6) = 4 from 6 = 4+2 = 3+3 = 2+2+2.
G.f. = 1 + x^2 + x^3 + 2*x^4 + 2*x^5 + 4*x^6 + 4*x^7 + 7*x^8 + 8*x^9 + ...
From _Gus Wiseman_, May 19 2019: (Start)
The a(2) = 1 through a(9) = 8 partitions not containing 1 are the following. The Heinz numbers of these partitions are given by A005408.
  (2)  (3)  (4)   (5)   (6)    (7)    (8)     (9)
            (22)  (32)  (33)   (43)   (44)    (54)
                        (42)   (52)   (53)    (63)
                        (222)  (322)  (62)    (72)
                                      (332)   (333)
                                      (422)   (432)
                                      (2222)  (522)
                                              (3222)
The a(2) = 1 through a(9) = 8 partitions of n - 1 whose least part appears exactly once are the following. The Heinz numbers of these partitions are given by A247180.
  (1)  (2)  (3)   (4)   (5)    (6)    (7)     (8)
            (21)  (31)  (32)   (42)   (43)    (53)
                        (41)   (51)   (52)    (62)
                        (221)  (321)  (61)    (71)
                                      (331)   (332)
                                      (421)   (431)
                                      (2221)  (521)
                                              (3221)
The a(2) = 1 through a(9) = 8 partitions of n + 1 where the number of parts is itself a part are the following. The Heinz numbers of these partitions are given by A325761.
  (21)  (22)  (32)   (42)   (52)    (62)    (72)     (82)
              (311)  (321)  (322)   (332)   (333)    (433)
                            (331)   (431)   (432)    (532)
                            (4111)  (4211)  (531)    (631)
                                            (4221)   (4222)
                                            (4311)   (4321)
                                            (51111)  (4411)
                                                     (52111)
The a(2) = 1 through a(8) = 7 partitions of n whose greatest part appears at least twice are the following. The Heinz numbers of these partitions are given by A070003.
  (11)  (111)  (22)    (221)    (33)      (331)      (44)
               (1111)  (11111)  (222)     (2221)     (332)
                                (2211)    (22111)    (2222)
                                (111111)  (1111111)  (3311)
                                                     (22211)
                                                     (221111)
                                                     (11111111)
Nonisomorphic representatives of the a(2) = 1 through a(6) = 4 2-regular multigraphs with n edges and n vertices are the following.
  {12,12}  {12,13,23}  {12,12,34,34}  {12,12,34,35,45}  {12,12,34,34,56,56}
                       {12,13,24,34}  {12,13,24,35,45}  {12,12,34,35,46,56}
                                                        {12,13,23,45,46,56}
                                                        {12,13,24,35,46,56}
The a(2) = 1 through a(9) = 8 partitions of n with no part greater than the number of ones are the following. The Heinz numbers of these partitions are given by A325762.
  (11)  (111)  (211)   (2111)   (2211)    (22111)    (22211)     (33111)
               (1111)  (11111)  (3111)    (31111)    (32111)     (222111)
                                (21111)   (211111)   (41111)     (321111)
                                (111111)  (1111111)  (221111)    (411111)
                                                     (311111)    (2211111)
                                                     (2111111)   (3111111)
                                                     (11111111)  (21111111)
                                                                 (111111111)
(End)
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 836.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 115, p*(n).
  • H. P. Robinson, Letter to N. J. A. Sloane, Jan 04 1974.
  • 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).
  • P. G. Tait, Scientific Papers, Cambridge Univ. Press, Vol. 1, 1898, Vol. 2, 1900, see Vol. 1, p. 334.

Crossrefs

First differences of partition numbers A000041. Cf. A053445, A072380, A081094, A081095, A232697.
Pairwise sums seem to be in A027336.
Essentially the same as A085811.
A column of A090824 and of A133687 and of A292508 and of A292622. Cf. A229161.
2-regular not necessarily connected graphs: A008483 (simple graphs), A000041 (multigraphs with loops allowed), this sequence (multigraphs with loops forbidden), A027336 (graphs with loops allowed but no multiple edges). - Jason Kimberley, Jan 05 2011
See also A098743 (parts that do not divide n).
Numbers n such that in the edge-delete game on the path P_{n} the first player does not have a winning strategy: A274161. - Lyndsey Wong, Jul 09 2016
Row sums of characteristic array A145573.
Number of partitions of n into parts >= m: A008483 (m = 3), A008484 (m = 4), A185325 - A185329 (m = 5 through 9).

Programs

  • GAP
    Concatenation([1],List([1..41],n->NrPartitions(n)-NrPartitions(n-1))); # Muniru A Asiru, Aug 20 2018
    
  • Magma
    A41 := func; [A41(n)-A41(n-1):n in [0..50]]; // Jason Kimberley, Jan 05 2011
    
  • Maple
    with(combstruct): ZL1:=[S, {S=Set(Cycle(Z,card>1))}, unlabeled]: seq(count(ZL1,size=n), n=0..50);  # Zerinvary Lajos, Sep 24 2007
    G:= {P=Set (Set (Atom, card>1))}: combstruct[gfsolve](G, unlabeled, x): seq  (combstruct[count] ([P, G, unlabeled], size=i), i=0..50);  # Zerinvary Lajos, Dec 16 2007
    with(combstruct):a:=proc(m) [ZL, {ZL=Set(Cycle(Z, card>=m))}, unlabeled]; end: A:=a(2):seq(count(A, size=n), n=0..50);  # Zerinvary Lajos, Jun 11 2008
    # alternative Maple program:
    A002865:= proc(n) option remember; `if`(n=0, 1, add(
          (numtheory[sigma](j)-1)*A002865(n-j), j=1..n)/n)
        end:
    seq(A002865(n), n=0..60);  # Alois P. Heinz, Sep 17 2017
  • Mathematica
    Table[ PartitionsP[n + 1] - PartitionsP[n], {n, -1, 50}] (* Robert G. Wilson v, Jul 24 2004 *)
    f[1, 1] = 1; f[n_, k_] := f[n, k] = If[n < 0, 0, If[k > n, 0, If[k == n, 1, f[n, k + 1] + f[n - k, k]]]]; Table[ f[n, 2], {n, 50}] (* Robert G. Wilson v *)
    Table[SeriesCoefficient[Exp[Sum[x^(2*k)/(k*(1 - x^k)), {k, 1, n}]], {x, 0, n}], {n, 0, 50}] (* Vaclav Kotesovec, Aug 18 2018 *)
    CoefficientList[Series[1/QPochhammer[x^2, x], {x,0,50}], x] (* G. C. Greubel, Nov 03 2019 *)
    Table[Count[IntegerPartitions[n],?(FreeQ[#,1]&)],{n,0,50}] (* _Harvey P. Dale, Feb 12 2023 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( (1 - x) / eta(x + x * O(x^n)), n))};
    
  • PARI
    a(n)=if(n,numbpart(n)-numbpart(n-1),1) \\ Charles R Greathouse IV, Nov 26 2012
    
  • Python
    from sympy import npartitions
    def A002865(n): return npartitions(n)-npartitions(n-1) if n else 1 # Chai Wah Wu, Mar 30 2023
  • SageMath
    def A002865_list(prec):
        P. = PowerSeriesRing(ZZ, prec)
        return P( 1/product((1-x^(m+2)) for m in (0..60)) ).list()
    A002865_list(50) # G. C. Greubel, Nov 03 2019
    

Formula

G.f.: Product_{m>1} 1/(1-x^m).
a(0)=1, a(n) = p(n) - p(n-1), n >= 1, with the partition numbers p(n) := A000041(n).
a(n) = A085811(n+3). - James Sellers, Dec 06 2005 [Corrected by Gionata Neri, Jun 14 2015]
a(n) = A116449(n) + A116450(n). - Reinhard Zumkeller, Feb 16 2006
a(n) = Sum_{k=2..floor((n+2)/2)} A008284(n-k+1,k-1) for n > 0. - Reinhard Zumkeller, Nov 04 2007
G.f.: 1 + Sum_{n>=2} x^n / Product_{k>=n} (1 - x^k). - Joerg Arndt, Apr 13 2011
G.f.: Sum_{n>=0} x^(2*n) / Product_{k=1..n} (1 - x^k). - Joerg Arndt, Apr 17 2011
a(n) = A090824(n,1) for n > 0. - Reinhard Zumkeller, Oct 10 2012
a(n) ~ Pi * exp(sqrt(2*n/3)*Pi) / (12*sqrt(2)*n^(3/2)) * (1 - (3*sqrt(3/2)/Pi + 13*Pi/(24*sqrt(6)))/sqrt(n) + (217*Pi^2/6912 + 9/(2*Pi^2) + 13/8)/n). - Vaclav Kotesovec, Feb 26 2015, extended Nov 04 2016
G.f.: exp(Sum_{k>=1} (sigma_1(k) - 1)*x^k/k). - Ilya Gutkovskiy, Aug 21 2018
a(0) = 1, a(n) = A232697(n) - 1. - George Beck, May 09 2019
From Peter Bala, Feb 19 2021: (Start)
G.f.: A(q) = Sum_{n >= 0} q^(n^2)/( (1 - q)*Product_{k = 2..n} (1 - q^k)^2 ).
More generally, A(q) = Sum_{n >= 0} q^(n*(n+r))/( (1 - q) * Product_{k = 2..n} (1 - q^k)^2 * Product_{i = 1..r} (1 - q^(n+i)) ) for r = 0,1,2,.... (End)
G.f.: 1 + Sum_{n >= 1} x^(n+1)/Product_{k = 1..n-1} 1 - x^(k+2). - Peter Bala, Dec 01 2024

A257993 Least gap in the partition having Heinz number n; index of the least prime not dividing n.

Original entry on oeis.org

1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3
Offset: 1

Views

Author

Emeric Deutsch, May 18 2015

Keywords

Comments

The "least gap" of a partition is the least positive integer that is not a part of the partition. For example, the least gap of the partition [7,4,2,2,1] is 3.
We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product(p_j-th prime, j=1...r) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] we get 2*2*3*7*29 = 2436.
In the Maple program the subprogram B yields the partition with Heinz number n.
Sum of least gaps of all partitions of m = A022567(m).
From Antti Karttunen, Aug 22 2016: (Start)
Index of the least prime not dividing n. (After a formula given by Heinz.)
Least k such that A002110(k) does not divide n.
One more than the number of trailing zeros in primorial base representation of n, A049345.
(End)
The least gap is also called the mex (minimal excludant) of the partition. - Gus Wiseman, Apr 20 2021

Examples

			a(18) = 3 because the partition having Heinz number 18 = 2*3*3 is [1,2,2], having least gap equal to 3.
		

References

  • G. E. Andrews and K. Eriksson, Integer Partitions, Cambridge Univ. Press, 2004, Cambridge.
  • Miklós Bóna, A Walk Through Combinatorics, World Scientific Publishing Co., 2002.

Crossrefs

Positions of 1's are A005408.
Positions of 2's are A047235.
The number of gaps is A079067.
The version for crank is A257989.
The triangle counting partitions by this statistic is A264401.
One more than A276084.
The version for greatest difference is A286469 or A286470.
A maximal instead of minimal version is A339662.
Positions of even terms are A342050.
Positions of odd terms are A342051.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A056239 adds up prime indices, row sums of A112798.
A073491 lists numbers with gap-free prime indices.
A238709 counts partitions by sum and least difference.
A333214 lists positions of adjacent unequal prime gaps.
A339737 counts partitions by sum and greatest gap.

Programs

  • Maple
    with(numtheory): a := proc (n) local B, q: B := proc (n) local nn, j, m: nn := op(2, ifactors(n)): for j to nops(nn) do m[j] := op(j, nn) end do: [seq(seq(pi(op(1, m[i])), q = 1 .. op(2, m[i])), i = 1 .. nops(nn))] end proc: for q while member(q, B(n)) = true do  end do: q end proc: seq(a(n), n = 1 .. 150);
    # second Maple program:
    a:= n-> `if`(n=1, 1, (s-> min({$1..(max(s)+1)} minus s))(
            {map(x-> numtheory[pi](x[1]), ifactors(n)[2])[]})):
    seq(a(n), n=1..100);  # Alois P. Heinz, May 09 2016
    # faster:
    A257993 := proc(n) local p, c; c := 1; p := 2;
    while n mod p = 0 do p := nextprime(p); c := c + 1 od: c end:
    seq(A257993(n), n=1..100); # Peter Luschny, Jun 04 2017
  • Mathematica
    A053669[n_] := For[p = 2, True, p = NextPrime[p], If[CoprimeQ[p, n], Return[p]]]; a[n_] := PrimePi[A053669[n]]; Array[a, 100] (* Jean-François Alcover, Nov 28 2016 *)
    Table[k = 1; While[! CoprimeQ[Prime@ k, n], k++]; k, {n, 100}] (* Michael De Vlieger, Jun 22 2017 *)
  • PARI
    a(n) = forprime(p=2,, if (n % p, return(primepi(p)))); \\ Michel Marcus, Jun 22 2017
  • Python
    from sympy import nextprime, primepi
    def a053669(n):
        p = 2
        while True:
            if n%p!=0: return p
            else: p=nextprime(p)
    def a(n): return primepi(a053669(n)) # Indranil Ghosh, May 12 2017
    
  • Scheme
    (define (A257993 n) (let loop ((n n) (i 1)) (let* ((p (A000040 i)) (d (modulo n p))) (if (not (zero? d)) i (loop (/ (- n d) p) (+ 1 i))))))
    ;; Antti Karttunen, Aug 22 2016
    

Formula

a(n) = A000720(A053669(n)). - Alois P. Heinz, May 18 2015
From Antti Karttunen, Aug 22-30 2016: (Start)
a(n) = 1 + A276084(n).
a(n) = A055396(A276086(n)).
A276152(n) = A002110(a(n)).
(End)
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1 + Sum_{k>=1} 1/A002110(k) = 1.705230... (1 + A064648). - Amiram Eldar, Jul 23 2022
a(n) << log n/log log n. - Charles R Greathouse IV, Dec 03 2022

Extensions

A simpler description added to the name by Antti Karttunen, Aug 22 2016

A349156 Number of integer partitions of n whose mean is not an integer.

Original entry on oeis.org

1, 0, 0, 1, 1, 5, 3, 13, 11, 21, 28, 54, 31, 99, 111, 125, 165, 295, 259, 488, 425, 648, 933, 1253, 943, 1764, 2320, 2629, 2962, 4563, 3897, 6840, 6932, 9187, 11994, 12840, 12682, 21635, 25504, 28892, 28187, 44581, 42896, 63259, 66766, 74463, 104278, 124752
Offset: 0

Views

Author

Gus Wiseman, Nov 14 2021

Keywords

Comments

Equivalently, partitions whose length does not divide their sum.
By conjugation, also the number of integer partitions of n with greatest part not dividing n.

Examples

			The a(3) = 1 through a(8) = 11 partitions:
  (21)  (211)  (32)    (2211)   (43)      (332)
               (41)    (3111)   (52)      (422)
               (221)   (21111)  (61)      (431)
               (311)            (322)     (521)
               (2111)           (331)     (611)
                                (421)     (22211)
                                (511)     (32111)
                                (2221)    (41111)
                                (3211)    (221111)
                                (4111)    (311111)
                                (22111)   (2111111)
                                (31111)
                                (211111)
		

Crossrefs

Below, "!" means either enumerative or set theoretical complement.
The version for nonempty subsets is !A051293.
The complement is counted by A067538, ranked by A316413.
The geometric version is !A067539, strict !A326625, ranked by !A326623.
The strict case is !A102627.
The version for prime factors is A175352, complement A078175.
The version for distinct prime factors is A176587, complement A078174.
The ordered version (compositions) is !A271654, ranked by !A096199.
The multiplicative version (factorizations) is !A326622, geometric !A326028.
The conjugate is ranked by !A326836.
The conjugate strict version is !A326850.
These partitions are ranked by A348551.
A000041 counts integer partitions.
A326567/A326568 give the mean of prime indices, conjugate A326839/A326840.
A236634 counts unbalanced partitions, complement of A047993.
A327472 counts partitions not containing their mean, complement of A237984.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],!IntegerQ[Mean[#]]&]],{n,0,30}]

Formula

a(n > 0) = A000041(n) - A067538(n).

A264401 Triangle read by rows: T(n,k) is the number of partitions of n having least gap k.

Original entry on oeis.org

1, 0, 1, 1, 1, 1, 1, 1, 2, 2, 1, 2, 3, 2, 4, 4, 2, 1, 4, 6, 4, 1, 7, 8, 5, 2, 8, 11, 8, 3, 12, 15, 10, 4, 1, 14, 20, 15, 6, 1, 21, 26, 19, 9, 2, 24, 35, 27, 12, 3, 34, 45, 34, 17, 5, 41, 58, 47, 23, 6, 1, 55, 75, 59, 31, 10, 1, 66, 96, 79, 41, 13, 2
Offset: 0

Views

Author

Emeric Deutsch, Nov 21 2015

Keywords

Comments

The "least gap" or "mex" of a partition is the least positive integer that is not a part of the partition. For example, the least gap of the partition [7,4,2,2,1] is 3.
Sum of entries in row n is A000041(n).
T(n,1) = A002865(n).
Sum_{k>=1} k*T(n,k) = A022567(n).

Examples

			Row n=5 is 2,3,2; indeed, the least gaps of [5], [4,1], [3,2], [3,1,1], [2,2,1], [2,1,1,1], and [1,1,1,1,1] are 1, 2, 1, 2, 3, 3, and 2, respectively (i.e., two 1s, three 2s, and two 3s).
Triangle begins:
   1
   0   1
   1   1
   1   1   1
   2   2   1
   2   3   2
   4   4   2   1
   4   6   4   1
   7   8   5   2
   8  11   8   3
  12  15  10   4   1
  14  20  15   6   1
  21  26  19   9   2
		

Crossrefs

Row sums are A000041.
Row lengths are A002024.
Column k = 1 is A002865.
Column k = 2 is A027336.
The strict case is A343348.
A000009 counts strict partitions.
A000041 counts partitions.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A015723 counts strict partitions with a selected part.
A257993 gives the least gap of the partition with Heinz number n.
A339564 counts factorizations with a selected factor.
A342050 ranks partitions with even least gap.
A342051 ranks partitions with odd least gap.

Programs

  • Maple
    g := (sum(t^j*x^((1/2)*j*(j-1))*(1-x^j), j = 1 .. 80))/(product(1-x^i, i = 1 .. 80)): gser := simplify(series(g, x = 0, 23)): for n from 0 to 30 do P[n] := sort(coeff(gser, x, n)) end do: for n from 0 to 25 do seq(coeff(P[n], t, j), j = 1 .. degree(P[n])) end do; # yields sequence in triangular form
    # second Maple program:
    b:= proc(n, i) option remember; `if`(n=0, `if`(i=0, [1, 0],
          [0, x]), `if`(i<1, 0, (p-> [0, p[2] +p[1]*x^i])(
          b(n, i-1)) +add(b(n-i*j, i-1), j=1..n/i)))
        end:
    T:= n->(p->seq(coeff(p, x, i), i=1..degree(p)))(b(n, n+1)[2]):
    seq(T(n), n=0..20);  # Alois P. Heinz, Nov 29 2015
  • Mathematica
    Needs["Combinatorica`"]; {1, 0}~Join~Flatten[Table[Count[Map[If[# == {}, 0, First@ #] &@ Complement[Range@ n, #] &, Combinatorica`Partitions@ n], n_ /; n == k], {n, 17}, {k, n}] /. 0 -> Nothing] (* Michael De Vlieger, Nov 21 2015 *)
    mingap[q_]:=Min@@Complement[Range[If[q=={},0,Max[q]]+1],q];Table[Length[Select[IntegerPartitions[n],mingap[#]==k&]],{n,0,15},{k,Round[Sqrt[2*(n+1)]]}] (* Gus Wiseman, Apr 19 2021 *)
    b[n_, i_] := b[n, i] = If[n == 0, If[i == 0, {1, 0}, {0, x}], If[i<1, {0, 0}, {0, #[[2]] + #[[1]]*x^i}&[b[n, i-1]] + Sum[b[n-i*j, i - 1], {j, 1, n/i}]]];
    T[n_] := CoefficientList[b[n, n + 1], x][[2]] // Rest;
    T /@ Range[0, 20] // Flatten (* Jean-François Alcover, May 21 2021, after Alois P. Heinz *)

Formula

G.f.: G(t,x) = Sum_{j>=1} (t^j*x^{j(j-1)/2}*(1-x^j))/Product_{i>=1}(1-x^i).

A348551 Heinz numbers of integer partitions whose mean is not an integer.

Original entry on oeis.org

1, 6, 12, 14, 15, 18, 20, 24, 26, 33, 35, 36, 38, 40, 42, 44, 45, 48, 50, 51, 52, 54, 56, 58, 60, 63, 65, 66, 69, 70, 72, 74, 75, 76, 77, 80, 86, 92, 93, 95, 96, 102, 104, 106, 108, 112, 114, 117, 119, 120, 122, 123, 124, 126, 130, 132, 135, 136, 140, 141, 142
Offset: 1

Views

Author

Gus Wiseman, Nov 14 2021

Keywords

Comments

Equivalently, partitions whose length does not divide their sum.
The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.

Examples

			The terms and their prime indices begin:
   1: {}
   6: {1,2}
  12: {1,1,2}
  14: {1,4}
  15: {2,3}
  18: {1,2,2}
  20: {1,1,3}
  24: {1,1,1,2}
  26: {1,6}
  33: {2,5}
  35: {3,4}
  36: {1,1,2,2}
  38: {1,8}
  40: {1,1,1,3}
  42: {1,2,4}
  44: {1,1,5}
  45: {2,2,3}
  48: {1,1,1,1,2}
		

Crossrefs

A version counting nonempty subsets is A000079 - A051293.
A version counting factorizations is A001055 - A326622.
A version counting compositions is A011782 - A271654.
A version for prime factors is A175352, complement A078175.
A version for distinct prime factors A176587, complement A078174.
The complement is A316413, counted by A067538, strict A102627.
The geometric version is the complement of A326623.
The conjugate version is the complement of A326836.
These partitions are counted by A349156.
A000041 counts partitions.
A001222 counts prime factors with multiplicity.
A018818 counts partitions into divisors, ranked by A326841.
A143773 counts partitions into multiples of the length, ranked by A316428.
A236634 counts unbalanced partitions.
A047993 counts balanced partitions, ranked by A106529.
A056239 adds up prime indices, row sums of A112798.
A326567/A326568 gives the mean of prime indices, conjugate A326839/A326840.
A327472 counts partitions not containing their mean, complement A237984.

Programs

  • Maple
    q:= n-> (l-> nops(l)=0 or irem(add(i, i=l), nops(l))>0)(map
            (i-> numtheory[pi](i[1])$i[2], ifactors(n)[2])):
    select(q, [$1..142])[];  # Alois P. Heinz, Nov 19 2021
  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],!IntegerQ[Mean[primeMS[#]]]&]

A097986 Number of strict integer partitions of n with a part dividing all the other parts.

Original entry on oeis.org

1, 1, 2, 2, 2, 4, 3, 5, 5, 7, 6, 12, 9, 13, 15, 20, 18, 28, 26, 37, 39, 47, 49, 71, 68, 85, 94, 117, 120, 159, 160, 201, 216, 257, 277, 348, 357, 430, 470, 562, 592, 720, 758, 901, 981, 1134, 1220, 1457, 1542, 1798, 1952, 2250, 2419, 2819, 3023, 3482, 3773, 4291
Offset: 1

Views

Author

Vladeta Jovovic, Oct 23 2004

Keywords

Comments

If n > 0, we can assume such a part is the smallest. - Gus Wiseman, Apr 23 2021
Also the number of uniform (constant multiplicity) partitions of n containing 1, ranked by A367586. The strict case is A096765. The version without 1 is A329436. - Gus Wiseman, Dec 01 2023

Examples

			From _Gus Wiseman_, Dec 01 2023: (Start)
The a(1) = 1 through a(8) = 5 strict partitions with a part dividing all the other parts:
  (1)  (2)  (3)    (4)    (5)    (6)      (7)      (8)
            (2,1)  (3,1)  (4,1)  (4,2)    (6,1)    (6,2)
                                 (5,1)    (4,2,1)  (7,1)
                                 (3,2,1)           (4,3,1)
                                                   (5,2,1)
The a(1) = 1 through a(8) = 5 uniform partitions containing 1:
  (1)  (11)  (21)   (31)    (41)     (51)      (61)       (71)
             (111)  (1111)  (11111)  (321)     (421)      (431)
                                     (2211)    (1111111)  (521)
                                     (111111)             (3311)
                                                          (11111111)
(End)
		

Crossrefs

The non-strict version is A083710.
The case with no 1's is A098965.
The Heinz numbers of these partitions are A339563.
The strict complement is counted by A341450.
The version for "divisible by" instead of "dividing" is A343347.
The case where there is also a part divisible by all the others is A343378.
The case where there is no part divisible by all the others is A343381.
A000005 counts divisors.
A000009 counts strict partitions.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A015723 counts strict partitions with a selected part.
A018818 counts partitions into divisors (strict: A033630).
A167865 counts strict chains of divisors > 1 summing to n.

Programs

  • Mathematica
    Take[ CoefficientList[ Expand[ Sum[x^k*Product[1 + x^(k*i), {i, 2, 62}], {k, 62}]], x], {2, 60}] (* Robert G. Wilson v, Nov 01 2004 *)
    Table[Length[Select[IntegerPartitions[n], UnsameQ@@#&&Or@@Table[And@@IntegerQ/@(#/x), {x,#}]&]], {n,0,30}] (* Gus Wiseman, Apr 23 2021 *)
  • PARI
    A_x(N) = {my(x='x+O('x^N)); Vec(sum(k=1,N,x^k*prod(i=2,N-k, (1+x^(k*i)))))}
    A_x(50) \\ John Tyler Rascoe, Nov 19 2024

Formula

a(n) = Sum_{d|n} A025147(d-1).
G.f.: Sum_{k>=1} (x^k*Product_{i>=2} (1+x^(k*i))).
a(n) ~ exp(Pi*sqrt(n/3)) / (8*3^(1/4)*n^(3/4)). - Vaclav Kotesovec, Jul 06 2025

Extensions

More terms from Robert G. Wilson v, Nov 01 2004
Name shortened by Gus Wiseman, Apr 23 2021

A365830 Heinz numbers of incomplete integer partitions, meaning not every number from 0 to A056239(n) is the sum of some submultiset.

Original entry on oeis.org

3, 5, 7, 9, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 25, 26, 27, 28, 29, 31, 33, 34, 35, 37, 38, 39, 41, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 55, 57, 58, 59, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79, 81, 82, 83, 85, 86, 87, 88, 89
Offset: 1

Views

Author

Gus Wiseman, Sep 26 2023

Keywords

Comments

First differs from A325798 in lacking 156.
The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.
The complement (complete partitions) is A325781.

Examples

			The terms together with their prime indices begin:
   3: {2}
   5: {3}
   7: {4}
   9: {2,2}
  10: {1,3}
  11: {5}
  13: {6}
  14: {1,4}
  15: {2,3}
  17: {7}
  19: {8}
  21: {2,4}
  22: {1,5}
  23: {9}
  25: {3,3}
  26: {1,6}
  27: {2,2,2}
  28: {1,1,4}
For example, the submultisets of (1,1,2,6) (right column) and their sums (left column) are:
   0: ()
   1: (1)
   2: (2)  or (11)
   3: (12)
   4: (112)
   6: (6)
   7: (16)
   8: (26) or (116)
   9: (126)
  10: (1126)
But 5 is missing, so 156 is in the sequence.
		

Crossrefs

For prime indices instead of sums we have A080259, complement of A055932.
The complement is A325781, counted by A126796, strict A188431.
Positions of nonzero terms in A325799, complement A304793.
These partitions are counted by A365924, strict A365831.
A056239 adds up prime indices, row sums of A112798.
A276024 counts positive subset-sums of partitions, strict A284640
A299701 counts distinct subset-sums of prime indices.
A365918 counts distinct non-subset-sums of partitions, strict A365922.
A365923 counts partitions by distinct non-subset-sums, strict A365545.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    nmz[y_]:=Complement[Range[Total[y]],Total/@Subsets[y]];
    Select[Range[100],Length[nmz[prix[#]]]>0&]

A200745 Number of partitions of n into distinct non-divisors of n.

Original entry on oeis.org

1, 0, 0, 0, 0, 1, 0, 2, 1, 2, 2, 6, 1, 9, 5, 6, 5, 20, 4, 28, 7, 19, 24, 55, 6, 51, 45, 49, 27, 136, 16, 180, 50, 117, 143, 146, 28, 403, 242, 260, 68, 668, 91, 852, 246, 260, 649, 1370, 90, 1191, 493, 1110, 634, 2701, 386, 1635, 462, 2160, 2486, 5154, 167
Offset: 0

Views

Author

Reinhard Zumkeller, Nov 22 2011

Keywords

Examples

			a(10) = #{7+3, 6+4} = 2;
a(11) = #{9+2, 8+3, 7+4, 6+5, 6+3+2, 5+4+2} = 6;
a(12) = #{7+5} = 1;
a(13) = #{11+2, 10+3, 9+4, 8+5, 8+3+2, 7+6, 7+4+2, 6+5+2, 6+4+3} = 9;
a(14) = #{11+3, 10+4, 9+5, 8+6, 6+5+3} = 5;
a(15) = #{13+2, 11+5, 9+6, 9+4+2, 8+7, 8+5+2} = 6.
		

Crossrefs

Programs

  • Haskell
    a200745 n = p [nd | nd <- [1..n], mod n nd /= 0] n where
       p _  0 = 1
       p [] _ = 0
       p (k:ks) m | m < k = 0 | otherwise = p ks (m - k) + p ks m
  • Maple
    a:= proc(n) option remember; local b, l;
          l:= sort([({$1..n} minus numtheory[divisors](n))[]]);
          b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
                 b(n, i-1) +`if`(l[i]>n, 0, b(n-l[i], i-1))))
              end: forget(b):
          b(n, nops(l))
        end:
    seq(a(n), n=0..80);  # Alois P. Heinz, Jan 18 2013
  • Mathematica
    a[0] = 1; a[n_] := a[n] = Module[{b, l}, l = Sort[Range[n] ~Complement~ Divisors[n]]; b[m_, i_] := b[m, i] = If[m == 0, 1, If[i < 1, 0, b[m, i - 1] + If[l[[i]] > m, 0, b[m - l[[i]], i - 1]]]]; b[n, Length[l]]];
    Table[a[n], {n, 0, 80}] (* Jean-François Alcover, Feb 06 2017, after Alois P. Heinz *)

A342193 Numbers with no prime index dividing all the other prime indices.

Original entry on oeis.org

1, 15, 33, 35, 45, 51, 55, 69, 75, 77, 85, 91, 93, 95, 99, 105, 119, 123, 135, 141, 143, 145, 153, 155, 161, 165, 175, 177, 187, 195, 201, 203, 205, 207, 209, 215, 217, 219, 221, 225, 231, 245, 247, 249, 253, 255, 265, 275, 279, 285, 287, 291, 295, 297, 299
Offset: 1

Views

Author

Gus Wiseman, Apr 11 2021

Keywords

Comments

Alternative name: 1 and numbers with smallest prime index not dividing all the other prime indices.
First differs from A339562 in having 45.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
Also 1 and Heinz numbers of integer partitions with smallest part not dividing all the others (counted by A338470). The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.

Examples

			The sequence of terms together with their prime indices begins:
      1: {}         105: {2,3,4}      201: {2,19}
     15: {2,3}      119: {4,7}        203: {4,10}
     33: {2,5}      123: {2,13}       205: {3,13}
     35: {3,4}      135: {2,2,2,3}    207: {2,2,9}
     45: {2,2,3}    141: {2,15}       209: {5,8}
     51: {2,7}      143: {5,6}        215: {3,14}
     55: {3,5}      145: {3,10}       217: {4,11}
     69: {2,9}      153: {2,2,7}      219: {2,21}
     75: {2,3,3}    155: {3,11}       221: {6,7}
     77: {4,5}      161: {4,9}        225: {2,2,3,3}
     85: {3,7}      165: {2,3,5}      231: {2,4,5}
     91: {4,6}      175: {3,3,4}      245: {3,4,4}
     93: {2,11}     177: {2,17}       247: {6,8}
     95: {3,8}      187: {5,7}        249: {2,23}
     99: {2,2,5}    195: {2,3,6}      253: {5,9}
		

Crossrefs

The complement is counted by A083710 (strict: A097986).
The complement with no 1's is A083711 (strict: A098965).
These partitions are counted by A338470 (strict: A341450).
The squarefree case is A339562, with squarefree complement A339563.
The case with maximum prime index not divisible by all others is A343338.
The case with maximum prime index divisible by all others is A343339.
A000005 counts divisors.
A000070 counts partitions with a selected part.
A001221 counts distinct prime factors.
A006128 counts partitions with a selected position (strict: A015723).
A056239 adds up prime indices, row sums of A112798.
A299702 lists Heinz numbers of knapsack partitions.
A339564 counts factorizations with a selected factor.

Programs

  • Mathematica
    Select[Range[100],#==1||With[{p=PrimePi/@First/@FactorInteger[#]},!And@@IntegerQ/@(p/Min@@p)]&]

A083711 a(n) = A083710(n) - A000041(n-1).

Original entry on oeis.org

1, 1, 1, 2, 1, 4, 1, 5, 3, 7, 1, 14, 1, 13, 8, 20, 1, 33, 1, 40, 14, 44, 1, 85, 6, 79, 25, 117, 1, 181, 1, 196, 45, 233, 17, 389, 1, 387, 80, 545, 1, 750, 1, 839, 165, 1004, 1, 1516, 12, 1612, 234, 2040, 1, 2766, 48, 3142, 388, 3720, 1, 5295, 1, 5606, 663, 7038, 83, 9194, 1, 10379, 1005
Offset: 1

Views

Author

N. J. A. Sloane, Jun 16 2003

Keywords

Comments

Number of integer partitions of n with no 1's with a part dividing all the others. If n > 0, we can assume such a part is the smallest. - Gus Wiseman, Apr 18 2021

Examples

			From _Gus Wiseman_, Apr 18 2021: (Start)
The a(6) = 4 through a(12) = 13 partitions:
  (6)      (7)  (8)        (9)      (10)         (11)  (12)
  (3,3)         (4,4)      (6,3)    (5,5)              (6,6)
  (4,2)         (6,2)      (3,3,3)  (8,2)              (8,4)
  (2,2,2)       (4,2,2)             (4,4,2)            (9,3)
                (2,2,2,2)           (6,2,2)            (10,2)
                                    (4,2,2,2)          (4,4,4)
                                    (2,2,2,2,2)        (6,3,3)
                                                       (6,4,2)
                                                       (8,2,2)
                                                       (3,3,3,3)
                                                       (4,4,2,2)
                                                       (6,2,2,2)
                                                       (4,2,2,2,2)
                                                       (2,2,2,2,2,2)
(End)
		

References

  • L. M. Chawla, M. O. Levan and J. E. Maxfield, On a restricted partition function and its tables, J. Natur. Sci. and Math., 12 (1972), 95-101.

Crossrefs

Allowing 1's gives A083710.
The strict case is A098965.
The complement (except also without 1's) is counted by A338470.
The dual version is A339619.
A000005 counts divisors.
A000041 counts partitions.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A018818 counts partitions into divisors (strict: A033630).
A167865 counts strict chains of divisors > 1 summing to n.
A339564 counts factorizations with a selected factor.

Programs

  • Maple
    with(combinat): with(numtheory): a := proc(n) c := 0: l := sort(convert(divisors(n), list)): for i from 1 to nops(l)-1 do c := c+numbpart(l[i]-1) od: RETURN(c): end: for j from 2 to 100 do printf(`%d,`,a(j)) od: # James Sellers, Jun 21 2003
    # second Maple program:
    a:= n-> max(1, add(combinat[numbpart](d-1), d=numtheory[divisors](n) minus {n})):
    seq(a(n), n=1..69);  # Alois P. Heinz, Feb 15 2023
  • Mathematica
    a[n_] := If[n==1, 1, Sum[PartitionsP[d-1], {d, Most@Divisors[n]}]];
    Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Feb 15 2023 *)

Formula

a(n) = Sum_{ d|n, dA000041(d-1).

Extensions

More terms from James Sellers, Jun 21 2003
Showing 1-10 of 30 results. Next