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.

Previous Showing 31-40 of 332 results. Next

A045782 Number of factorizations of n for some n (image of A001055).

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 9, 11, 12, 15, 16, 19, 21, 22, 26, 29, 30, 31, 36, 38, 42, 45, 47, 52, 56, 57, 64, 66, 67, 74, 77, 92, 97, 98, 101, 105, 109, 118, 135, 137, 139, 141, 162, 165, 171, 176, 181, 189, 195, 198, 203, 212, 231, 249, 250, 254, 257, 267, 269, 272, 289
Offset: 1

Views

Author

Keywords

Comments

Also the image of A318284. - Gus Wiseman, Jan 11 2020

Crossrefs

Factorizations are A001055 with image this sequence and complement A330976.
Strict factorizations are A045778 with image A045779 and complement A330975.
The least number with exactly a(n) factorizations is A045783(n).
The least number with exactly n factorizations is A330973(n).

Programs

  • Mathematica
    terms = 61; m0 = 10^5; dm = 10^4;
    f[1, ] = 1; f[n, k_] := f[n, k] = Sum[f[n/d, d], {d, Select[Divisors[n], 1 < # <= k &]}];
    Clear[seq]; seq[m_] := seq[m] = Sort[Tally[Table[f[n, n], {n, 1, m}]][[All, 1]]][[1 ;; terms]]; seq[m = m0]; seq[m += dm]; While[Print[m]; seq[m] != seq[m - dm], m += dm];
    seq[m] (* Jean-François Alcover, Oct 04 2018 *)

Formula

The Luca et al. paper shows that the number of terms with a(n) <= x is x^{ O( log log log x / log log x )}. - N. J. A. Sloane, Jun 12 2009

Extensions

Name edited by Gus Wiseman, Jan 11 2020

A379722 Numbers whose prime indices do not have the same sum as product.

Original entry on oeis.org

1, 4, 6, 8, 10, 12, 14, 15, 16, 18, 20, 21, 22, 24, 25, 26, 27, 28, 32, 33, 34, 35, 36, 38, 39, 40, 42, 44, 45, 46, 48, 49, 50, 51, 52, 54, 55, 56, 57, 58, 60, 62, 63, 64, 65, 66, 68, 69, 70, 72, 74, 75, 76, 77, 78, 80, 81, 82, 85, 86, 87, 88, 90, 91, 92, 93
Offset: 1

Views

Author

Gus Wiseman, Jan 08 2025

Keywords

Comments

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.
Partitions of this type are counted by A379736.
The complement is A301987, counted by A001055.

Examples

			The terms together with their prime indices begin:
    1: {}
    4: {1,1}
    6: {1,2}
    8: {1,1,1}
   10: {1,3}
   12: {1,1,2}
   14: {1,4}
   15: {2,3}
   16: {1,1,1,1}
   18: {1,2,2}
   20: {1,1,3}
   21: {2,4}
   22: {1,5}
   24: {1,1,1,2}
   25: {3,3}
   26: {1,6}
   27: {2,2,2}
   28: {1,1,4}
   32: {1,1,1,1,1}
		

Crossrefs

Nonzeros of A325036.
A000040 lists the primes, differences A001223.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A324851 finds numbers > 1 divisible by the sum of their prime indices.
A379666 counts partitions by sum and product, without 1's A379668.
A379681 gives sum plus product of prime indices, firsts A379682.
Counting and ranking multisets by comparing sum and product:
- same: A001055 (strict A045778), ranks A301987
- divisible: A057567, ranks A326155
- divisor: A057568, ranks A326149, see A379733
- greater: A096276 shifted right, ranks A325038
- greater or equal: A096276, ranks A325044
- less: A114324, ranks A325037, see A318029
- less or equal: A319005, ranks A379721
- different: A379736, ranks A379722 (this)

Programs

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

A292886 Number of knapsack factorizations of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 2, 2, 1, 4, 1, 2, 2, 4, 1, 4, 1, 4, 2, 2, 1, 7, 2, 2, 3, 4, 1, 5, 1, 6, 2, 2, 2, 8, 1, 2, 2, 7, 1, 5, 1, 4, 4, 2, 1, 11, 2, 4, 2, 4, 1, 7, 2, 7, 2, 2, 1, 11, 1, 2, 4, 7, 2, 5, 1, 4, 2, 5, 1, 15, 1, 2, 4, 4, 2, 5, 1, 11, 4, 2, 1, 11, 2
Offset: 1

Views

Author

Gus Wiseman, Sep 26 2017

Keywords

Comments

A knapsack factorization is a finite multiset of positive integers greater than one such that every distinct submultiset has a different product.
The sequence giving the number of factorizations of n is described as "the multiplicative partition function" (see A001055), so knapsack factorizations are a multiplicative generalization of knapsack partitions. - Gus Wiseman, Oct 24 2017

Examples

			The a(36) = 8 factorizations are 2*2*3*3, 2*2*9, 2*18, 3*3*4, 3*12, 4*9, 6*6, 36. The factorization 2*3*6 is not knapsack.
		

Crossrefs

Programs

  • Mathematica
    postfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[postfacs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[postfacs[n],UnsameQ@@Times@@@Union[Subsets[#]]&]],{n,100}]

A303707 Number of factorizations of n using elements of A007916 (numbers that are not perfect powers).

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 3, 1, 2, 2, 1, 1, 3, 1, 3, 2, 2, 1, 4, 1, 2, 1, 3, 1, 5, 1, 1, 2, 2, 2, 5, 1, 2, 2, 4, 1, 5, 1, 3, 3, 2, 1, 5, 1, 3, 2, 3, 1, 4, 2, 4, 2, 2, 1, 9, 1, 2, 3, 1, 2, 5, 1, 3, 2, 5, 1, 8, 1, 2, 3, 3, 2, 5, 1, 5, 1, 2, 1, 9, 2, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Apr 29 2018

Keywords

Comments

First differs from A081707 at a(60) = 9, A081707(60) = 8.

Examples

			The a(60) = 9 factorizations are (2*2*3*5), (2*2*15), (2*3*10), (2*5*6), (2*30), (3*20), (5*12), (6*10), (60).
		

Crossrefs

Programs

  • Mathematica
    radQ[n_]:=Or[n===1,GCD@@FactorInteger[n][[All,2]]===1];
    facsr[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facsr[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],radQ]}]];
    Table[Length[facsr[n]],{n,100}]

Formula

Dirichlet g.f.: Product_{n in A007916} 1/(1 - n^s).

A370813 Number of non-condensed integer factorizations of n into unordered factors > 1.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 2, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 2, 1, 0, 0, 0, 0, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Mar 04 2024

Keywords

Comments

A multiset is condensed iff it is possible to choose a different divisor of each element.

Examples

			The a(96) = 4 factorizations: (2*2*2*2*2*3), (2*2*2*2*6), (2*2*2*3*4), (2*2*2*12).
		

Crossrefs

Partitions not of this type are counted by A239312, ranks A368110.
Factors instead of divisors: A368413, complement A368414, unique A370645.
Partitions of this type are counted by A370320, ranks A355740.
Subsets of this type: A370583 and A370637, complement A370582 and A370636.
The complement is counted by A370814, partitions A370592, ranks A368100.
For a unique choice we have A370815, partitions A370595, ranks A370810.
A000005 counts divisors.
A001055 counts factorizations, strict A045778.
A355731 counts choices of a divisor of each prime index, firsts A355732.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join @@ Table[Map[Prepend[#,d]&,Select[facs[n/d],Min @@ #>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n],Length[Select[Tuples[Divisors /@ #],UnsameQ@@#&]]==0&]],{n,100}]

A317791 Number of non-isomorphic multiset partitions of the multiset of prime indices of n (row n of A112798).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 2, 2, 1, 4, 1, 2, 2, 5, 1, 4, 1, 4, 2, 2, 1, 7, 2, 2, 3, 4, 1, 3, 1, 7, 2, 2, 2, 7, 1, 2, 2, 7, 1, 3, 1, 4, 4, 2, 1, 12, 2, 4, 2, 4, 1, 7, 2, 7, 2, 2, 1, 9, 1, 2, 4, 11, 2, 3, 1, 4, 2, 3, 1, 16, 1, 2, 4, 4, 2, 3, 1, 12, 5, 2, 1, 9, 2, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Aug 07 2018

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n.
a(n) depends only on prime signature of n (cf. A025487). - Antti Karttunen, Dec 03 2018
Are any terms of the complement known? In particular, does this sequence contain 6? - Gus Wiseman, Oct 21 2022

Examples

			Non-isomorphic representatives of the a(42) = 3 multiset partitions are {{1,2,4}}, {{1},{2,4}}, {{1},{2},{4}}.
Non-isomorphic representatives of the a(60) = 9 multiset partitions:
  {1123},
  {1}{123}, {2}{113}, {11}{23}, {12}{13},
  {1}{1}{23}, {1}{2}{13}, {2}{3}{11},
  {1}{1}{2}{3}.
Missing from this list are {3}{112} and {1}{3}{12}, which are isomorphic to {2}{113} and {1}{2}{13} respectively.
For n = 180 = 2^2 * 3^2 * 5, there are A001055(180) = 26 different factorizations to one or more factors larger than 1. Of these 18 are such that by swapping 2 and 3 in each factor of that factorization the result is another, different factorization of 180, while the other 8 cases are such that 2 <-> 3 swap doesn't change the factorization. Thus a(180) = 18/2 + 8 = 17. - _Antti Karttunen_, Dec 03 2018
		

Crossrefs

Programs

  • Mathematica
    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]]]];
    sysnorm[{}] := {};sysnorm[m_]:=If[Union@@m!=Range[Max@@Flatten[m]],sysnorm[m/.Rule@@@Table[{(Union@@m)[[i]],i},{i,Length[Union@@m]}]],First[Sort[sysnorm[m,1]]]];sysnorm[m_,aft_]:=If[Length[Union@@m]<=aft,{m},With[{mx=Table[Count[m,i,{2}],{i,Select[Union@@m,#>=aft&]}]},Union@@(sysnorm[#,aft+1]&/@Union[Table[Map[Sort,m/.{par+aft-1->aft,aft->par+aft-1},{0,1}],{par,First/@Position[mx,Max[mx]]}]])]];
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[Union[sysnorm/@mps[primeMS[n]]]],{n,100}]

Formula

For all n, a(n) <= A001055(n). - Antti Karttunen, Dec 01 2018
If n is squarefree with k prime factors, or if n = p^k for p prime, we have a(n) = A000041(k).
a(n) = A318285(A181819(n)). - Andrew Howroyd, Jan 17 2023

Extensions

Terms corrected by Gus Wiseman, Dec 04 2018

A340596 Number of co-balanced factorizations of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jan 15 2021

Keywords

Comments

We define a factorization of n into factors > 1 to be co-balanced if it has exactly A001221(n) factors.

Examples

			The a(n) co-balanced factorizations for n = 12, 24, 36, 72, 120, 144, 180:
  2*6    3*8     4*9     8*9     3*5*8     2*72     4*5*9
  3*4    4*6     6*6     2*36    4*5*6     3*48     5*6*6
         2*12    2*18    3*24    2*2*30    4*36     2*2*45
                 3*12    4*18    2*3*20    6*24     2*3*30
                         6*12    2*4*15    8*18     2*5*18
                                 2*5*12    9*16     2*6*15
                                 2*6*10    12*12    2*9*10
                                 3*4*10             3*3*20
                                                    3*4*15
                                                    3*5*12
                                                    3*6*10
		

Crossrefs

Positions of terms > 1 are A126706.
Positions of 1's are A303554.
The version for unlabeled multiset partitions is A319616.
The alt-balanced version is A340599.
The balanced version is A340653.
The cross-balanced version is A340654.
The twice-balanced version is A340655.
A001055 counts factorizations.
A045778 counts strict factorizations.
A316439 counts factorizations by product and length.
Other balance-related sequences:
- A010054 counts balanced strict partitions.
- A047993 counts balanced partitions.
- A098124 counts balanced compositions.
- A106529 lists Heinz numbers of balanced partitions.
- A340597 lists numbers with an alt-balanced factorization.
- A340598 counts balanced set partitions.
- A340600 counts unlabeled balanced multiset partitions.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n],Length[#]==PrimeNu[n]&]],{n,100}]
  • PARI
    A340596(n, m=n, om=omega(n)) = if(1==n,(0==om), sumdiv(n, d, if((d>1)&&(d<=m), A340596(n/d, d, om-1)))); \\ Antti Karttunen, Jun 10 2024

Extensions

Data section extended up to a(120) by Antti Karttunen, Jun 10 2024

A379666 Array read by antidiagonals downward where A(n,k) is the number of integer partitions of n with product k.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jan 01 2025

Keywords

Comments

Counts finite multisets of positive integers by sum and product.

Examples

			Array begins:
        k=1 k=2 k=3 k=4 k=5 k=6 k=7 k=8 k=9 k10 k11 k12
        -----------------------------------------------
   n=0:  1   0   0   0   0   0   0   0   0   0   0   0
   n=1:  1   0   0   0   0   0   0   0   0   0   0   0
   n=2:  1   1   0   0   0   0   0   0   0   0   0   0
   n=3:  1   1   1   0   0   0   0   0   0   0   0   0
   n=4:  1   1   1   2   0   0   0   0   0   0   0   0
   n=5:  1   1   1   2   1   1   0   0   0   0   0   0
   n=6:  1   1   1   2   1   2   0   2   1   0   0   0
   n=7:  1   1   1   2   1   2   1   2   1   1   0   2
   n=8:  1   1   1   2   1   2   1   3   1   1   0   3
   n=9:  1   1   1   2   1   2   1   3   2   1   0   3
  n=10:  1   1   1   2   1   2   1   3   2   2   0   3
  n=11:  1   1   1   2   1   2   1   3   2   2   1   3
  n=12:  1   1   1   2   1   2   1   3   2   2   1   4
For example, the A(9,12) = 3 partitions are: (6,2,1), (4,3,1,1), (3,2,2,1,1).
Antidiagonals begin:
   n+k=1: 1
   n+k=2: 0 1
   n+k=3: 0 0 1
   n+k=4: 0 0 1 1
   n+k=5: 0 0 0 1 1
   n+k=6: 0 0 0 1 1 1
   n+k=7: 0 0 0 0 1 1 1
   n+k=8: 0 0 0 0 2 1 1 1
   n+k=9: 0 0 0 0 0 2 1 1 1
  n+k=10: 0 0 0 0 0 1 2 1 1 1
  n+k=11: 0 0 0 0 0 1 1 2 1 1 1
  n+k=12: 0 0 0 0 0 0 2 1 2 1 1 1
  n+k=13: 0 0 0 0 0 0 0 2 1 2 1 1 1
  n+k=14: 0 0 0 0 0 0 2 1 2 1 2 1 1 1
  n+k=15: 0 0 0 0 0 0 1 2 1 2 1 2 1 1 1
  n+k=16: 0 0 0 0 0 0 0 1 3 1 2 1 2 1 1 1
For example, antidiagonal n+k=10 counts the following partitions:
  n=5: (5)
  n=6: (411), (2211)
  n=7: (31111)
  n=8: (2111111)
  n=9: (111111111)
so the 10th antidiagonal is: (0,0,0,0,0,1,2,1,1,1).
		

Crossrefs

Row sums are A000041 = partitions of n, strict A000009, no ones A002865.
Diagonal A(n,n) is A001055(n) = factorizations of n, strict A045778.
Antidiagonal sums are A379667.
The case without ones is A379668, antidiagonal sums A379669 (zeros A379670).
The strict case is A379671, antidiagonal sums A379672.
The strict case without ones is A379678, antidiagonal sums A379679 (zeros A379680).
A316439 counts factorizations by length, partitions A008284.
A326622 counts factorizations with integer mean, strict A328966.
Counting and ranking multisets by comparing sum and product:
- same: A001055, ranks A301987
- divisible: A057567, ranks A326155
- divisor: A057568, ranks A326149, see A379733
- greater than: A096276 shifted right, ranks A325038
- greater or equal: A096276, ranks A325044
- less than: A114324, ranks A325037, see A318029
- less or equal: A319005, ranks A379721, see A025147
- different: A379736, ranks A379722, see A111133

Programs

  • Mathematica
    nn=12;
    tt=Table[Length[Select[IntegerPartitions[n],Times@@#==k&]],{n,0,nn},{k,1,nn}] (* array *)
    tr=Table[tt[[j,i-j]],{i,2,nn},{j,i-1}] (* antidiagonals *)
    Join@@tr (* sequence *)

A381441 Number of multisets that can be obtained by partitioning the prime indices of n into a set of sets (set system) and taking their sums.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Mar 06 2025

Keywords

Comments

First differs from A050326 at a(210) = 13, A050326(210) = 15. This comes from the set systems {{3},{1,2,4}} and {{1,2},{3,4}}, and from {{4},{1,2,3}} and {{1,3},{2,4}}.
Also the number of multisets that can be obtained by taking the sums of prime indices of each factor in a strict factorization of n into squarefree numbers > 1.
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.
A multiset partition can be regarded as an arrow in the ranked poset of integer partitions. For example, we have {{1},{1,2},{1,3},{1,2,3}}: {1,1,1,1,2,2,3,3} -> {1,3,4,6}, or (33221111) -> (6431) (depending on notation).
Sets of sets are generally not transitive. For example, we have arrows: {{1},{1,2}}: {1,1,2} -> {1,3} and {{1,3}}: {1,3} -> {4}, but there is no set of sets {1,1,2} -> {4}.

Examples

			The prime indices of 60 are {1,1,2,3}, with partitions into sets of sets:
  {{1},{1,2,3}}
  {{1,2},{1,3}}
  {{1},{2},{1,3}}
  {{1},{3},{1,2}}
with block-sums: {1,6}, {3,4}, {1,2,4}, {1,3,3}, which are all different, so a(60) = 4.
		

Crossrefs

Before taking sums we had A050326, non-strict A050320.
Positions of 0 are A293243.
Positions of 1 are A293511.
This is the strict version of A381078 (lower A381454).
For distinct block-sums (instead of blocks) we have A381634, before sums A381633.
Other multiset partitions of prime indices:
- For multisets of constant multisets (A000688) see A381455 (upper), A381453 (lower).
- For multiset partitions (A001055) see A317141 (upper), A300383 (lower).
- For strict multiset partitions (A045778) see A381452.
- For sets of constant multisets (A050361) see A381715.
- For strict multiset partitions with distinct sums (A321469) see A381637.
- For sets of constant multisets with distinct sums (A381635) see A381716, A381636.
More on set systems: A050342, A116539, A279785, A296120, A318361.
A000041 counts integer partitions, strict A000009.
A000040 lists the primes, differences A001223.
A003963 gives product of prime indices.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A122111 represents conjugation in terms of Heinz numbers.
A265947 counts refinement-ordered pairs of integer partitions.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Union[Sort[Total/@prix/@#]&/@Select[facs[n],UnsameQ@@#&&And@@SquareFreeQ/@#&]]],{n,100}]

Formula

a(A002110(n)) = A066723(n).

A381718 Number of normal multiset partitions of weight n into sets with distinct sums.

Original entry on oeis.org

1, 1, 2, 6, 23, 106, 549, 3184, 20353, 141615, 1063399, 8554800, 73281988, 665141182, 6369920854, 64133095134, 676690490875, 7462023572238, 85786458777923, 1025956348473929, 12739037494941490
Offset: 0

Views

Author

Gus Wiseman, Mar 26 2025

Keywords

Comments

We call a multiset or multiset partition normal iff it covers an initial interval of positive integers. The weight of a multiset partition is the sum of sizes of its blocks.

Examples

			The a(1) = 1 through a(3) = 6 multiset partitions:
  {{1}}  {{1,2}}    {{1,2,3}}
         {{1},{2}}  {{1},{1,2}}
                    {{1},{2,3}}
                    {{2},{1,2}}
                    {{2},{1,3}}
                    {{1},{2},{3}}
The a(4) = 23 factorizations:
  2*3*6  5*30    3*30    2*30    210
         10*15   6*15    6*10    2*105
         2*5*15  2*3*15  2*3*10  3*70
         3*5*10                  5*42
                                 7*30
                                 6*35
                                 10*21
                                 2*3*35
                                 2*5*21
                                 2*7*15
                                 3*5*14
                                 2*3*5*7
		

Crossrefs

For distinct blocks instead of sums we have A116539, see A050326.
Without distinct sums we have A116540 (normal set multipartitions).
Twice-partitions of this type are counted by A279785.
Without strict blocks we have A326519.
Factorizations of this type are counted by A381633.
For constant instead of strict blocks we have A382203.
For distinct sizes instead of sums we have A382428, non-strict blocks A326517.
For equal instead of distinct block-sums we have A382429, non-strict blocks A326518.
A000670 counts patterns, ranked by A055932 and A333217, necklace A019536.
A001055 count factorizations, strict A045778.
Normal multiset partitions: A034691, A035310, A255906.
Set multipartitions: A089259, A270995, A296119, A318360.

Programs

  • Mathematica
    allnorm[n_Integer]:=Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[mset_]:=Union[Sort[Sort/@(#/.x_Integer:>mset[[x]])]&/@sps[Range[Length[mset]]]];
    Table[Length[Join@@(Select[mps[#],UnsameQ@@Total/@#&&And@@UnsameQ@@@#&]&/@allnorm[n])],{n,0,5}]

Extensions

a(10)-a(11) from Robert Price, Mar 31 2025
a(12)-a(20) from Christian Sievers, Apr 05 2025
Previous Showing 31-40 of 332 results. Next