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

A073470 Trisection of A007294.

Original entry on oeis.org

1, 2, 4, 6, 10, 15, 22, 32, 44, 60, 81, 107, 141, 181, 233, 295, 372, 465, 575, 710, 869, 1059, 1286, 1549, 1863, 2225, 2654, 3149, 3724, 4391, 5157, 6043, 7059, 8220, 9553, 11066, 12798, 14760, 16988, 19509, 22357, 25576, 29200, 33274, 37856, 42985
Offset: 0

Views

Author

Randall L Rathbun, Aug 26 2002

Keywords

Crossrefs

Programs

  • Magma
    function a007294(n, k) return n gt k*(k+1)/2 select a007294(n-k*(k+1)/2, k) + a007294(n, k+1) else n eq k*(k+1)/2 select 1 else 0; end function; [ a007294(n, 1): n in [0..135 by 3] ]; // Klaus Brockhaus, Jan 23 2009

Extensions

More terms from Klaus Brockhaus, Jan 23 2009

A073471 Trisection of A007294.

Original entry on oeis.org

1, 2, 4, 7, 11, 17, 24, 35, 48, 66, 89, 117, 153, 197, 252, 320, 401, 501, 619, 763, 933, 1134, 1374, 1654, 1985, 2372, 2821, 3345, 3949, 4652, 5460, 6390, 7457, 8676, 10071, 11662, 13468, 15526, 17850, 20487, 23463, 26818, 30598, 34842, 39614, 44962
Offset: 0

Views

Author

Randall L Rathbun, Aug 26 2002

Keywords

Crossrefs

Programs

  • Magma
    function a007294(n, k) return n gt k*(k+1)/2 select a007294(n-k*(k+1)/2, k) + a007294(n, k+1) else n eq k*(k+1)/2 select 1 else 0; end function; [ a007294(n, 1): n in [1..136 by 3] ]; // Klaus Brockhaus, Jan 23 2009

Extensions

More terms from Klaus Brockhaus, Jan 23 2009

A073472 Trisection of A007294.

Original entry on oeis.org

1, 2, 4, 7, 11, 17, 25, 36, 50, 68, 92, 121, 159, 205, 262, 332, 417, 520, 645, 793, 971, 1180, 1430, 1722, 2065, 2467, 2934, 3478, 4109, 4835, 5676, 6638, 7747, 9013, 10459, 12108, 13980, 16110, 18522, 21247, 24333, 27800, 31715, 36106, 41039, 46567
Offset: 0

Views

Author

Randall L Rathbun, Aug 26 2002

Keywords

Crossrefs

Programs

  • Magma
    function a007294(n, k) return n gt k*(k+1)/2 select a007294(n-k*(k+1)/2, k) + a007294(n, k+1) else n eq k*(k+1)/2 select 1 else 0; end function; [ a007294(n, 2): n in [2..137 by 3] ]; // Klaus Brockhaus, Jan 23 2009

Extensions

More terms from Klaus Brockhaus, Jan 23 2009

A325325 Number of integer partitions of n with distinct differences between successive parts.

Original entry on oeis.org

1, 1, 2, 2, 4, 5, 5, 8, 11, 12, 16, 22, 21, 30, 34, 42, 49, 64, 67, 87, 95, 117, 132, 160, 169, 207, 230, 274, 301, 360, 395, 463, 506, 602, 656, 762, 834, 960, 1042, 1220, 1311, 1505, 1643, 1859, 2000, 2341, 2491, 2827, 3083, 3464, 3747, 4302, 4561, 5154
Offset: 0

Views

Author

Gus Wiseman, Apr 23 2019

Keywords

Comments

The Heinz numbers of these partitions are given by A325368.

Examples

			The a(0) = 1 through a(9) = 12 partitions:
  ()  (1)  (2)   (3)   (4)    (5)    (6)    (7)    (8)     (9)
           (11)  (21)  (22)   (32)   (33)   (43)   (44)    (54)
                       (31)   (41)   (42)   (52)   (53)    (63)
                       (211)  (221)  (51)   (61)   (62)    (72)
                              (311)  (411)  (322)  (71)    (81)
                                            (331)  (332)   (441)
                                            (421)  (422)   (522)
                                            (511)  (431)   (621)
                                                   (521)   (711)
                                                   (611)   (4221)
                                                   (4211)  (4311)
                                                           (5211)
For example, (5,2,1,1) has differences (-3,-1,0), which are distinct, so (5,2,1,1) is counted under a(9).
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@Differences[#]&]],{n,0,30}]

A007862 Number of triangular numbers that divide n.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Also a(n) is the total number of ways to represent n+1 as a centered polygonal number of the form km(m+1)/2+1 for k>1. - Alexander Adamchuk, Apr 26 2007
Number of oblong numbers that divide 2n. - Ray Chandler, Jun 24 2008
The number of divisors d of 2n such that d+1 is also a divisor of 2n, see first formula. - Michel Marcus, Jun 18 2015
From Gus Wiseman, May 03 2019: (Start)
Also the number of integer partitions of n forming a finite arithmetic progression with offset 0, i.e. the differences are all equal (with the last part taken to be 0). The Heinz numbers of these partitions are given by A325327. For example, the a(1) = 1 through a(12) = 3 partitions are (A = 10, B = 11, C = 12):
1 2 3 4 5 6 7 8 9 A B C
21 42 63 4321 84
321 642
(End)

Crossrefs

Programs

  • Haskell
    a007862 = sum . map a010054 . a027750_row
    -- Reinhard Zumkeller, Jul 05 2014
    
  • Mathematica
    sup=90; TriN=Array[ (#+1)(#+2)/2&, Floor[ N[ Sqrt[ sup*2 ] ] ]-1 ]; Array[ Function[n, 1+Count[ Map[ Mod[ n, # ]&, TriN ], 0 ] ], sup ]
    Table[Count[Divisors[k], ?(IntegerQ[Sqrt[8 # + 1]] &)], {k, 105}] (* _Jayanta Basu, Aug 12 2013 *)
    Table[Length[Select[IntegerPartitions[n],SameQ@@Differences[Append[#,0]]&]],{n,0,30}] (* Gus Wiseman, May 03 2019 *)
  • PARI
    a(n) = sumdiv(n, d, ispolygonal(d, 3)); \\ Michel Marcus, Jun 18 2015
    
  • Python
    from itertools import pairwise
    from sympy import divisors
    def A007862(n): return sum(1 for a, b in pairwise(divisors(n<<1)) if a+1==b)  # Chai Wah Wu, Jun 09 2025

Formula

a(n) = Sum_{d|2*n,d+1|2*n} 1.
G.f.: Sum_{k>=1} x^A000217(k)/(1-x^A000217(k)). - Jon Perry, Jul 03 2004
a(A130317(n)) = n and a(m) <> n for m < A130317(n). - Reinhard Zumkeller, May 23 2007
a(n) = A129308(2n). - Ray Chandler, Jun 24 2008
a(n) = Sum_{k=1..A000005(n)} A010054(A027750(n,k)). - Reinhard Zumkeller, Jul 05 2014
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 2. - Amiram Eldar, Dec 31 2023

Extensions

Extended by Ray Chandler, Jun 24 2008

A107742 G.f.: Product_{j>=1} Product_{i>=1} (1 + x^(i*j)).

Original entry on oeis.org

1, 1, 2, 4, 6, 10, 17, 25, 38, 59, 86, 125, 184, 260, 369, 524, 726, 1005, 1391, 1894, 2576, 3493, 4687, 6272, 8373, 11090, 14647, 19294, 25265, 32991, 42974, 55705, 72025, 92895, 119349, 152965, 195592, 249280, 316991, 402215, 508932, 642598, 809739, 1017850, 1276959, 1599015, 1997943, 2491874, 3102477, 3855165, 4782408, 5922954
Offset: 0

Views

Author

Vladeta Jovovic, Jun 11 2005

Keywords

Comments

From Gus Wiseman, Sep 13 2022: (Start)
Also the number of multiset partitions of integer partitions of n into intervals, where an interval is a set of positive integers with all differences of adjacent elements equal to 1. For example, the a(1) = 1 through a(4) = 6 multiset partitions are:
{{1}} {{2}} {{3}} {{4}}
{{1},{1}} {{1,2}} {{1},{3}}
{{1},{2}} {{2},{2}}
{{1},{1},{1}} {{1},{1,2}}
{{1},{1},{2}}
{{1},{1},{1},{1}}
Intervals are counted by A001227, ranked by A073485.
The initial version is A007294.
The strict version is A327731.
The version for gapless multisets instead of intervals is A356941.
The case of strict partitions is A356957.
Also the number of multiset partitions of integer partitions of n into distinct constant blocks. For example, the a(1) = 1 through a(4) = 6 multiset partitions are:
{{1}} {{2}} {{3}} {{4}}
{{1,1}} {{1,1,1}} {{2,2}}
{{1},{2}} {{1},{3}}
{{1},{1,1}} {{1,1,1,1}}
{{2},{1,1}}
{{1},{1,1,1}}
Constant multisets are counted by A000005, ranked by A000961.
The non-strict version is A006171.
The unlabeled version is A089259.
The non-constant block version is A261049.
The version for twice-partitions is A279786, factorizations A296131.
Also the number of multiset partitions of integer partitions of n into constant blocks of odd length. For example, a(1) = 1 through a(4) = 6 multiset partitions are:
{{1}} {{2}} {{3}} {{4}}
{{1},{1}} {{1,1,1}} {{1},{3}}
{{1},{2}} {{2},{2}}
{{1},{1},{1}} {{1},{1,1,1}}
{{1},{1},{2}}
{{1},{1},{1},{1}}
The strict version is A327731 (also).
(End)

Crossrefs

Product_{k>=1} (1 + x^k)^sigma_m(k): this sequence (m=0), A192065 (m=1), A288414 (m=2), A288415 (m=3), A301548 (m=4), A301549 (m=5), A301550 (m=6), A301551 (m=7), A301552 (m=8).
A000041 counts integer partitions, strict A000009.
A000110 counts set partitions.
A072233 counts partitions by sum and length.

Programs

  • Mathematica
    nmax = 50; CoefficientList[Series[Product[(1+x^(i*j)), {i, 1, nmax}, {j, 1, nmax/i}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Jan 04 2017 *)
    nmax = 50; CoefficientList[Series[Product[(1+x^k)^DivisorSigma[0, k], {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Mar 23 2018 *)
    nmax = 50; s = 1 + x; Do[s *= Sum[Binomial[DivisorSigma[0, k], j]*x^(j*k), {j, 0, nmax/k}]; s = Expand[s]; s = Take[s, Min[nmax + 1, Exponent[s, x] + 1, Length[s]]];, {k, 2, nmax}]; Take[CoefficientList[s, x], nmax + 1] (* Vaclav Kotesovec, Aug 28 2018 *)
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    chQ[y_]:=Length[y]<=1||Union[Differences[y]]=={1};
    Table[Length[Select[Join@@mps/@IntegerPartitions[n],And@@chQ/@#&]],{n,0,5}] (* Gus Wiseman, Sep 13 2022 *)
  • PARI
    a(n)=polcoeff(prod(k=1,n,prod(j=1,n\k,1+x^(j*k)+x*O(x^n))),n) /* Paul D. Hanna */
    
  • PARI
    N=66;  x='x+O('x^N); gf=1/prod(j=0,N, eta(x^(2*j+1))); gf=prod(j=1,N,(1+x^j)^numdiv(j)); Vec(gf) /* Joerg Arndt, May 03 2008 */
    
  • PARI
    {a(n)=if(n==0,1,polcoeff(exp(sum(m=1,n,sigma(m)*x^m/(1-x^(2*m)+x*O(x^n))/m)),n))} /* Paul D. Hanna, Mar 28 2009 */

Formula

Euler transform of A001227.
Weigh transform of A000005.
G.f. satisfies: log(A(x)) = Sum_{n>=1} A109386(n)/n*x^n, where A109386(n) = Sum_{d|n} d*Sum_{m|d} (m mod 2). - Paul D. Hanna, Jun 26 2005
G.f.: A(x) = exp( Sum_{n>=1} sigma(n)*x^n/(1-x^(2n)) /n ). - Paul D. Hanna, Mar 28 2009
G.f.: Product_{n>=1} Q(x^n) where Q(x) is the g.f. of A000009. - Joerg Arndt, Feb 27 2014
a(0) = 1, a(n) = (1/n)*Sum_{k=1..n} A109386(k)*a(n-k) for n > 0. - Seiichi Manyama, Jun 04 2017
Conjecture: log(a(n)) ~ Pi*sqrt(n*log(n)/6). - Vaclav Kotesovec, Aug 29 2018

Extensions

More terms from Paul D. Hanna, Jun 26 2005

A000294 Expansion of g.f. Product_{k >= 1} (1 - x^k)^(-k*(k+1)/2).

Original entry on oeis.org

1, 1, 4, 10, 26, 59, 141, 310, 692, 1483, 3162, 6583, 13602, 27613, 55579, 110445, 217554, 424148, 820294, 1572647, 2992892, 5652954, 10605608, 19765082, 36609945, 67405569, 123412204, 224728451, 407119735, 733878402, 1316631730, 2351322765, 4180714647, 7401898452, 13051476707, 22922301583, 40105025130, 69909106888, 121427077241, 210179991927, 362583131144
Offset: 0

Views

Author

Keywords

Comments

Number of partitions of n if there are k(k+1)/2 kinds of k (k=1,2,...). E.g., a(3)=10 because we have six kinds of 3, three kinds of 2+1 because there are three kinds of 2 and 1+1+1+1. - Emeric Deutsch, Mar 23 2005
Euler transform of the triangular numbers 1,3,6,10,...
Equals A028377: [1, 1, 3, 9, 19, 46, 100, ...] convolved with the aerated version of A000294: [1, 0, 1, 0, 4, 0, 10, 0, 26, 0, 59, ...]. - Gary W. Adamson, Jun 13 2009
The formula for p3(n) in the article by S. Finch (page 2) is incomplete, terms with n^(1/2) and n^(1/4) are also needed. These terms are in the article by Mustonen and Rajesh (page 2) and agree with my results, but in both articles the multiplicative constant is marked only as C, resp. c3(m). The following is a closed form of this constant: Pi^(1/24) * exp(1/24 - Zeta(3) / (8*Pi^2) + 75*Zeta(3)^3 / (2*Pi^8)) / (A^(1/2) * 2^(157/96) * 15^(13/96)) = A255939 = 0.213595160470..., where A = A074962 is the Glaisher-Kinkelin constant and Zeta(3) = A002117. - Vaclav Kotesovec, Mar 11 2015 [The new version of "Integer Partitions" by S. Finch contains the missing terms, see pages 2 and 5. - Vaclav Kotesovec, May 12 2015]
Number of solid partitions of corner-hook volume n (see arXiv:2009.00592 among links for definition). E.g., a(2) = 1 because there is only one solid partition [[[2]]] with cohook volume 2; a(3) = 4 because we have three solid partitions with two 1's (entry at (1,1,1) contributes 1, another entry at (2,1,1) or (1,2,1) or (1,1,2) contributes 2 to corner-hook volume) and one solid partition with single entry 3 (which contributes 3 to the corner-hook volume). Namely as 3D arrays [[[1],[1]]],[[[1]],[[1]]],[[[1]],[[1]]], [[[3]]]. - Alimzhan Amanov, Jul 13 2021

References

  • R. Chandra, Tables of solid partitions, Proceedings of the Indian National Science Academy, 26 (1960), 134-139.
  • V. S. Nanda, Tables of solid partitions, Proceedings of the Indian National Science Academy, 19 (1953), 313-314.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. also A278403 (log of o.g.f.).

Programs

  • Maple
    with(numtheory): etr:= proc(p) local b; b:=proc(n) option remember; local d,j; if n=0 then 1 else add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n fi end end: a:= etr(n-> n*(n+1)/2): seq(a(n), n=0..30);  # Alois P. Heinz, Sep 08 2008
  • Mathematica
    a[0] = 1; a[n_] := a[n] = 1/(2*n)*Sum[(DivisorSigma[2, k]+DivisorSigma[3, k])*a[n-k], {k, 1, n}]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Mar 05 2014, after Vladeta Jovovic *)
    nmax=50; CoefficientList[Series[Product[1/(1-x^k)^(k*(k+1)/2),{k,1,nmax}],{x,0,nmax}],x] (* Vaclav Kotesovec, Mar 11 2015 *)
  • PARI
    a(n)=if(n<0, 0, polcoeff(exp(sum(k=1, n, x^k/(1-x^k)^3/k, x*O(x^n))), n)) \\ Joerg Arndt, Apr 16 2010
    
  • SageMath
    # uses[EulerTransform from A166861]
    b = EulerTransform(lambda n: binomial(n+1, 2))
    print([b(n) for n in range(37)]) # Peter Luschny, Nov 11 2020

Formula

a(n) = (1/(2*n))*Sum_{k=1..n} (sigma[2](k)+sigma[3](k))*a(n-k). - Vladeta Jovovic, Sep 17 2002
a(n) ~ Pi^(1/24) * exp(1/24 - Zeta(3) / (8*Pi^2) + 75*Zeta(3)^3 / (2*Pi^8) - 15^(5/4) * Zeta(3)^2 * n^(1/4) / (2^(7/4)*Pi^5) + 15^(1/2) * Zeta(3) * n^(1/2) / (2^(1/2)*Pi^2) + 2^(7/4) * Pi * n^(3/4) / (3*15^(1/4))) / (A^(1/2) * 2^(157/96) * 15^(13/96) * n^(61/96)), where A = A074962 = 1.2824271291... is the Glaisher-Kinkelin constant and Zeta(3) = A002117 = 1.202056903... . - Vaclav Kotesovec, Mar 11 2015
G.f.: exp(Sum_{k>=1} (sigma_2(k) + sigma_3(k))*x^k/(2*k)). - Ilya Gutkovskiy, Aug 21 2018

Extensions

More terms from Sascha Kurz, Aug 15 2002

A024940 Number of partitions of n into distinct triangular numbers.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, 1, 1, 2, 1, 1, 2, 1, 2, 2, 0, 2, 3, 1, 1, 3, 2, 1, 4, 3, 0, 3, 3, 2, 4, 3, 3, 3, 2, 3, 3, 2, 4, 6, 4, 2, 5, 4, 2, 6, 5, 3, 7, 6, 3, 5, 5, 5, 6, 5, 4, 7, 7, 6, 8, 6, 5, 9, 7, 4, 9, 9, 6, 10, 9, 4, 9, 10, 8, 11, 11, 9, 10, 10, 9, 10, 10, 9, 14, 14, 7, 14, 14, 7, 15, 15, 8, 15, 17, 13
Offset: 0

Views

Author

Keywords

Examples

			a(31) counts these partitions:  [28,3], [21,10], [21,6,3,1], [15,10,6] _Clark Kimberling_, Mar 09 2014
		

Crossrefs

Programs

  • Haskell
    a024940 = p $ tail a000217_list where
       p _  0 = 1
       p (k:ks) m = if m < k then 0 else p ks (m - k) + p ks m
    -- Reinhard Zumkeller, Jun 28 2013
  • Mathematica
    Drop[ CoefficientList[ Series[ Product[(1 + x^(k*(k + 1)/2)), {k, 1, 15}], {x, 0, 102}], x], 1]
    (* also *)
    t = Table[n (n + 1)/2, {n, 1, 200}] ; p[n_] := IntegerPartitions[n, All, t]; Table[p[n], {n, 0, 12}] (*shows unrestricted partitions*)
    d[n_] := Select[p[n], Max[Length /@ Split@#] == 1 &]; Table[d[n], {n, 1, 31}] (*shows strict partitions*)
    Table[Length[d[n]], {n, 1, 70}] (* Clark Kimberling, Mar 09 2014 *)
    nmax = 100; nn = Floor[Sqrt[8*nmax + 1]/2] + 1; poly = ConstantArray[0, nn*(nn+1)/2 + 1]; poly[[1]] = 1; poly[[2]] = 1; Do[Do[poly[[j + 1]] += poly[[j - k*(k+1)/2 + 1]], {j, nn*(nn+1)/2, k*(k+1)/2, -1}];, {k, 2, nn}]; Take[poly, nmax + 1] (* Vaclav Kotesovec, Dec 10 2016 *)

Formula

For n>0: a(n) = b(n, 1) where b(n, k) = if n>k*(k+1)/2 then b(n-k*(k+1)/2, k+1) + b(n, k+1) else (if n=k*(k+1)/2 then 1 else 0). - Reinhard Zumkeller, Aug 26 2003
a(n) ~ exp(3*Pi^(1/3) * ((sqrt(2)-1)*Zeta(3/2))^(2/3) * n^(1/3) / 2^(4/3)) * ((sqrt(2)-1)*Zeta(3/2))^(1/3) / (2^(5/3) * sqrt(3) * Pi^(1/3) * n^(5/6)). - Vaclav Kotesovec, Jan 02 2017
G.f.: prod_{i>=1} (1+x^A000217(i)). - R. J. Mathar, Sep 20 2020

A320348 Number of partition into distinct parts (a_1, a_2, ... , a_m) (a_1 > a_2 > ... > a_m and Sum_{k=1..m} a_k = n) such that a1 - a2, a2 - a_3, ... , a_{m-1} - a_m, a_m are different.

Original entry on oeis.org

1, 1, 1, 2, 3, 2, 4, 4, 4, 6, 9, 7, 13, 12, 13, 16, 22, 17, 28, 28, 31, 36, 50, 45, 63, 62, 74, 78, 102, 92, 123, 123, 146, 148, 191, 181, 228, 233, 280, 283, 348, 350, 420, 437, 518, 523, 616, 641, 727, 774, 884, 911, 1038, 1102, 1240, 1292, 1463, 1530, 1715, 1861, 2002
Offset: 1

Views

Author

Seiichi Manyama, Oct 11 2018

Keywords

Comments

Also the number of integer partitions of n whose parts cover an initial interval of positive integers with distinct multiplicities. Also the number of integer partitions of n whose multiplicities cover an initial interval of positive integers and are distinct (see A048767 for a bijection). - Gus Wiseman, May 04 2019

Examples

			n = 9
[9]        *********  a_1 = 9.
           ooooooooo
------------------------------------
[8, 1]             *        a_2 = 1.
            *******o  a_1 - a_2 = 7.
            oooooooo
------------------------------------
[7, 2]            **        a_2 = 2.
             *****oo  a_1 - a_2 = 5.
             ooooooo
------------------------------------
[5, 4]          ****        a_2 = 4.
               *oooo  a_1 - a_2 = 1.
               ooooo
------------------------------------
a(9) = 4.
From _Gus Wiseman_, May 04 2019: (Start)
The a(1) = 1 through a(11) = 9 strict partitions with distinct differences (where the last part is taken to be 0) are the following (A = 10, B = 11). The Heinz numbers of these partitions are given by A325388.
  (1)  (2)  (3)  (4)   (5)   (6)   (7)   (8)   (9)   (A)    (B)
                 (31)  (32)  (51)  (43)  (53)  (54)  (64)   (65)
                       (41)        (52)  (62)  (72)  (73)   (74)
                                   (61)  (71)  (81)  (82)   (83)
                                                     (91)   (92)
                                                     (631)  (A1)
                                                            (632)
                                                            (641)
                                                            (731)
The a(1) = 1 through a(10) = 6 partitions covering an initial interval of positive integers with distinct multiplicities are the following. The Heinz numbers of these partitions are given by A325326.
  1  11  111  211   221    21111   2221     22211     22221      222211
              1111  2111   111111  22111    221111    2211111    322111
                    11111          211111   2111111   21111111   2221111
                                   1111111  11111111  111111111  22111111
                                                                 211111111
                                                                 1111111111
The a(1) = 1 through a(10) = 6 partitions whose multiplicities cover an initial interval of positive integers and are distinct are the following (A = 10). The Heinz numbers of these partitions are given by A325337.
  (1)  (2)  (3)  (4)    (5)    (6)    (7)    (8)    (9)    (A)
                 (211)  (221)  (411)  (322)  (332)  (441)  (433)
                        (311)         (331)  (422)  (522)  (442)
                                      (511)  (611)  (711)  (622)
                                                           (811)
                                                           (322111)
(End)
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&UnsameQ@@Differences[Append[#,0]]&]],{n,30}] (* Gus Wiseman, May 04 2019 *)

A320387 Number of partitions of n into distinct parts such that the successive differences of consecutive parts are nonincreasing, and first difference <= first part.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 3, 2, 2, 4, 3, 4, 5, 3, 5, 7, 4, 7, 8, 6, 8, 11, 7, 9, 13, 9, 11, 16, 12, 15, 18, 13, 17, 20, 17, 21, 24, 19, 24, 30, 22, 28, 34, 26, 34, 38, 30, 37, 43, 37, 42, 48, 41, 50, 58, 48, 55, 64, 53, 64, 71, 59, 73, 81, 69, 79, 89, 79, 90, 101, 87, 100, 111
Offset: 0

Views

Author

Seiichi Manyama, Oct 12 2018

Keywords

Comments

Partitions are usually written with parts in descending order, but the conditions are easier to check "visually" if written in ascending order.
Generating function of the "second integrals" of partitions: given a partition (p_1, ..., p_s) written in weakly decreasing order, write the sequence B = (b_1, b_2, ..., b_s) = (p_1, p_1 + p_2, ..., p_1 + ... + p_s). The sequence gives the coefficients of the generating function summing q^(b_1 + ... + b_s) over all partitions of all nonnegative integers. - William J. Keith, Apr 23 2022
From Gus Wiseman, Jan 17 2023: (Start)
Equivalently, a(n) is the number of multisets (weakly increasing sequences of positive integers) with weighted sum n. For example, the Heinz numbers of the a(0) = 1 through a(15) = 7 multisets are:
1 2 3 4 7 6 8 10 15 12 16 18 20 26 24 28
5 11 9 17 19 14 21 22 27 41 30 32
13 23 29 31 33 55 39 34
25 35 37 43 45
49 77 47
65
121
These multisets are counted by A264034. The reverse version is A007294. The zero-based version is A359678.
(End)

Examples

			There are a(29) = 15 such partitions of 29:
  01: [29]
  02: [10, 19]
  03: [11, 18]
  04: [12, 17]
  05: [13, 16]
  06: [14, 15]
  07: [5, 10, 14]
  08: [6, 10, 13]
  09: [6, 11, 12]
  10: [7, 10, 12]
  11: [8, 10, 11]
  12: [3, 6, 9, 11]
  13: [5, 7, 8, 9]
  14: [2, 4, 6, 8, 9]
  15: [3, 5, 6, 7, 8]
There are a(30) = 18 such partitions of 30:
  01: [30]
  02: [10, 20]
  03: [11, 19]
  04: [12, 18]
  05: [13, 17]
  06: [14, 16]
  07: [5, 10, 15]
  08: [6, 10, 14]
  09: [6, 11, 13]
  10: [7, 10, 13]
  11: [7, 11, 12]
  12: [8, 10, 12]
  13: [3, 6, 9, 12]
  14: [9, 10, 11]
  15: [4, 7, 9, 10]
  16: [2, 4, 6, 8, 10]
  17: [6, 7, 8, 9]
  18: [4, 5, 6, 7, 8]
		

Crossrefs

Number of appearances of n > 0 in A304818, reverse A318283.
A053632 counts compositions by weighted sum.
A358194 counts partitions by weighted sum, reverse A264034.
Weighted sum of prime indices: A359497, A359676, A359682, A359754, A359755.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    ots[y_]:=Sum[i*y[[i]],{i,Length[y]}];
    Table[Length[Select[Range[2^n],ots[prix[#]]==n&]],{n,10}] (* Gus Wiseman, Jan 17 2023 *)
  • PARI
    seq(n)={Vec(sum(k=1, (sqrtint(8*n+1)+1)\2, my(t=binomial(k,2)); x^t/prod(j=1, k-1, 1 - x^(t-binomial(j,2)) + O(x^(n-t+1)))))} \\ Andrew Howroyd, Jan 22 2023
  • Ruby
    def partition(n, min, max)
      return [[]] if n == 0
      [max, n].min.downto(min).flat_map{|i| partition(n - i, min, i - 1).map{|rest| [i, *rest]}}
    end
    def f(n)
      return 1 if n == 0
      cnt = 0
      partition(n, 1, n).each{|ary|
        ary << 0
        ary0 = (1..ary.size - 1).map{|i| ary[i - 1] - ary[i]}
        cnt += 1 if ary0.sort == ary0
      }
      cnt
    end
    def A320387(n)
      (0..n).map{|i| f(i)}
    end
    p A320387(50)
    

Formula

G.f.: Sum_{k>=1} x^binomial(k,2)/Product_{j=1..k-1} (1 - x^(binomial(k,2)-binomial(j,2))). - Andrew Howroyd, Jan 22 2023
Showing 1-10 of 115 results. Next