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

A330469 Number of series-reduced rooted trees whose leaves are multisets with a total of n elements covering an initial interval of positive integers.

Original entry on oeis.org

1, 1, 4, 24, 250, 3744, 73408, 1768088, 50468854, 1664844040, 62304622320, 2607765903568, 120696071556230, 6120415124163512, 337440974546042416, 20096905939846645064, 1285779618228281270718, 87947859243850506008984, 6404472598196204610148232
Offset: 0

Views

Author

Gus Wiseman, Dec 22 2019

Keywords

Comments

Also the number of different colorings of phylogenetic trees with n labels using a multiset of colors covering an initial interval of positive integers. A phylogenetic tree is a series-reduced rooted tree whose leaves are (usually disjoint) sets.

Examples

			The a(3) = 24 trees:
  (123)          (122)          (112)          (111)
  ((1)(23))      ((1)(22))      ((1)(12))      ((1)(11))
  ((2)(13))      ((2)(12))      ((2)(11))      ((1)(1)(1))
  ((3)(12))      ((1)(2)(2))    ((1)(1)(2))    ((1)((1)(1)))
  ((1)(2)(3))    ((1)((2)(2)))  ((1)((1)(2)))
  ((1)((2)(3)))  ((2)((1)(2)))  ((2)((1)(1)))
  ((2)((1)(3)))
  ((3)((1)(2)))
		

Crossrefs

The singleton-reduced version is A316651.
The unlabeled version is A330465.
The strongly normal case is A330467.
The case when leaves are sets is A330764.
Row sums of A330762.

Programs

  • Mathematica
    allnorm[n_]:=If[n<=0,{{}},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[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    multing[t_,n_]:=Array[(t+#-1)/#&,n,1,Times];
    amemo[m_]:=amemo[m]=1+Sum[Product[multing[amemo[s[[1]]],Length[s]],{s,Split[c]}],{c,Select[mps[m],Length[#]>1&]}];
    Table[Sum[amemo[m],{m,allnorm[n]}],{n,0,5}]
  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
    R(n, k)={my(v=[]); for(n=1, n, v=concat(v, EulerT(concat(v, [binomial(n+k-1, k-1)]))[n])); v}
    seq(n)={concat([1], sum(k=1, n, R(n,k)*sum(r=k, n, binomial(r,k)*(-1)^(r-k))))} \\ Andrew Howroyd, Dec 29 2019

Extensions

Terms a(9) and beyond from Andrew Howroyd, Dec 29 2019

A330675 Number of balanced reduced multisystems of maximum depth whose atoms constitute a strongly normal multiset of size n.

Original entry on oeis.org

1, 1, 2, 6, 43, 440, 7158, 151414
Offset: 0

Views

Author

Gus Wiseman, Dec 30 2019

Keywords

Comments

A balanced reduced multisystem is either a finite multiset, or a multiset partition with at least two parts, not all of which are singletons, of a balanced reduced multisystem.
A finite multiset is strongly normal if it covers an initial interval of positive integers with weakly decreasing multiplicities.

Examples

			The a(2) = 2 and a(3) = 6 multisystems:
  {1,1}  {{1},{1,1}}
  {1,2}  {{1},{1,2}}
         {{1},{2,3}}
         {{2},{1,1}}
         {{2},{1,3}}
         {{3},{1,2}}
The a(4) = 43 multisystems (commas and outer brackets elided):
  {{1}}{{1}{11}} {{1}}{{1}{12}} {{1}}{{1}{22}} {{1}}{{1}{23}} {{1}}{{2}{34}}
  {{11}}{{1}{1}} {{11}}{{1}{2}} {{11}}{{2}{2}} {{11}}{{2}{3}} {{12}}{{3}{4}}
                 {{1}}{{2}{11}} {{1}}{{2}{12}} {{1}}{{2}{13}} {{1}}{{3}{24}}
                 {{12}}{{1}{1}} {{12}}{{1}{2}} {{12}}{{1}{3}} {{13}}{{2}{4}}
                 {{2}}{{1}{11}} {{2}}{{1}{12}} {{1}}{{3}{12}} {{1}}{{4}{23}}
                                {{2}}{{2}{11}} {{13}}{{1}{2}} {{14}}{{2}{3}}
                                {{22}}{{1}{1}} {{2}}{{1}{13}} {{2}}{{1}{34}}
                                               {{2}}{{3}{11}} {{2}}{{3}{14}}
                                               {{23}}{{1}{1}} {{23}}{{1}{4}}
                                               {{3}}{{1}{12}} {{2}}{{4}{13}}
                                               {{3}}{{2}{11}} {{24}}{{1}{3}}
                                                              {{3}}{{1}{24}}
                                                              {{3}}{{2}{14}}
                                                              {{3}}{{4}{12}}
                                                              {{34}}{{1}{2}}
                                                              {{4}}{{1}{23}}
                                                              {{4}}{{2}{13}}
                                                              {{4}}{{3}{12}}
		

Crossrefs

The case with all atoms equal is A000111.
The case with all atoms different is A006472.
The version allowing all depths is A330475.
The unlabeled version is A330663.
The version where the atoms are the prime indices of n is A330665.
The (weakly) normal version is A330676.
The version where the degrees are the prime indices of n is A330728.
Multiset partitions of strongly normal multisets are A035310.
Series-reduced rooted trees with strongly normal leaves are A316652.

Programs

  • Mathematica
    strnorm[n_]:=Flatten[MapIndexed[Table[#2,{#1}]&,#]]&/@IntegerPartitions[n];
    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]]]];
    totm[m_]:=Prepend[Join@@Table[totm[p],{p,Select[mps[m],1
    				

A330935 Irregular triangle read by rows where T(n,k) is the number of length-k chains from minimum to maximum in the poset of factorizations of n into factors > 1, ordered by refinement.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jan 04 2020

Keywords

Comments

This poset is equivalent to the poset of multiset partitions of the prime indices of n, ordered by refinement.

Examples

			Triangle begins:
   1:          16: 0 1 3 2    31: 1            46: 0 1
   2: 1        17: 1          32: 0 1 5 8 4    47: 1
   3: 1        18: 0 1 2      33: 0 1          48: 0 1 10 23 15
   4: 0 1      19: 1          34: 0 1          49: 0 1
   5: 1        20: 0 1 2      35: 0 1          50: 0 1 2
   6: 0 1      21: 0 1        36: 0 1 7 7      51: 0 1
   7: 1        22: 0 1        37: 1            52: 0 1 2
   8: 0 1 1    23: 1          38: 0 1          53: 1
   9: 0 1      24: 0 1 5 5    39: 0 1          54: 0 1 5 5
  10: 0 1      25: 0 1        40: 0 1 5 5      55: 0 1
  11: 1        26: 0 1        41: 1            56: 0 1 5 5
  12: 0 1 2    27: 0 1 1      42: 0 1 3        57: 0 1
  13: 1        28: 0 1 2      43: 1            58: 0 1
  14: 0 1      29: 1          44: 0 1 2        59: 1
  15: 0 1      30: 0 1 3      45: 0 1 2        60: 0 1 9 11
Row n = 48 counts the following chains (minimum and maximum not shown):
  ()  (6*8)      (2*3*8)->(6*8)       (2*2*2*6)->(2*4*6)->(6*8)
      (2*24)     (2*4*6)->(6*8)       (2*2*3*4)->(2*3*8)->(6*8)
      (3*16)     (2*3*8)->(2*24)      (2*2*3*4)->(2*4*6)->(6*8)
      (4*12)     (2*3*8)->(3*16)      (2*2*2*6)->(2*4*6)->(2*24)
      (2*3*8)    (2*4*6)->(2*24)      (2*2*2*6)->(2*4*6)->(4*12)
      (2*4*6)    (2*4*6)->(4*12)      (2*2*3*4)->(2*3*8)->(2*24)
      (3*4*4)    (3*4*4)->(3*16)      (2*2*3*4)->(2*3*8)->(3*16)
      (2*2*12)   (3*4*4)->(4*12)      (2*2*3*4)->(2*4*6)->(2*24)
      (2*2*2*6)  (2*2*12)->(2*24)     (2*2*3*4)->(2*4*6)->(4*12)
      (2*2*3*4)  (2*2*12)->(4*12)     (2*2*3*4)->(3*4*4)->(3*16)
                 (2*2*2*6)->(6*8)     (2*2*3*4)->(3*4*4)->(4*12)
                 (2*2*3*4)->(6*8)     (2*2*2*6)->(2*2*12)->(2*24)
                 (2*2*2*6)->(2*24)    (2*2*2*6)->(2*2*12)->(4*12)
                 (2*2*2*6)->(4*12)    (2*2*3*4)->(2*2*12)->(2*24)
                 (2*2*3*4)->(2*24)    (2*2*3*4)->(2*2*12)->(4*12)
                 (2*2*3*4)->(3*16)
                 (2*2*3*4)->(4*12)
                 (2*2*2*6)->(2*4*6)
                 (2*2*3*4)->(2*3*8)
                 (2*2*3*4)->(2*4*6)
                 (2*2*3*4)->(3*4*4)
                 (2*2*2*6)->(2*2*12)
                 (2*2*3*4)->(2*2*12)
		

Crossrefs

Row lengths are A001222.
Row sums are A317176.
Column k = 1 is A010051.
Column k = 2 is A066247.
Column k = 3 is A330936.
Final terms of each row are A317145.
The version for set partitions is A008826, with row sums A005121.
The version for integer partitions is A330785, with row sums A213427.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    upfacs[q_]:=Union[Sort/@Join@@@Tuples[facs/@q]];
    paths[eds_,start_,end_]:=If[start==end,Prepend[#,{}],#]&[Join@@Table[Prepend[#,e]&/@paths[eds,Last[e],end],{e,Select[eds,First[#]==start&]}]];
    Table[Length[Select[paths[Join@@Table[{y,#}&/@DeleteCases[upfacs[y],y],{y,facs[n]}],{n},First[facs[n]]],Length[#]==k-1&]],{n,100},{k,PrimeOmega[n]}]

Formula

T(2^n,k) = A330785(n,k).
T(n,1) + T(n,2) = 1.

A318849 Number of orderless tree-partitions of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 2, 2, 4, 6, 11, 8, 27, 20, 30, 38, 96, 74, 114, 58, 308, 234, 1052, 176, 509, 278, 3648, 374, 600, 1076, 1760, 814, 13003, 1306, 47006, 612, 2226, 4200, 3094, 2914, 172605, 16588, 9814, 2168, 640662, 6998, 2402388, 3698, 11496, 65936, 9082538, 4914, 17996
Offset: 1

Views

Author

Gus Wiseman, Sep 04 2018

Keywords

Comments

This multiset is generally not the same as the multiset of prime indices of n. For example, the prime indices of 12 are {1,1,2}, while a multiset whose multiplicities are {1,1,2} is {1,1,2,3}.
A tree-partition of m is either m itself or a multiset of tree-partitions, one of each part of a multiset partition of m with at least two parts.

Examples

			The a(7) = 11 orderless tree-partitions of {1,1,1,1}:
  (1111)
  ((1)(111))
  ((11)(11))
  ((1)(1)(11))
  ((1)((1)(11)))
  ((11)((1)(1)))
  ((1)(1)(1)(1))
  ((1)((1)(1)(1)))
  ((1)(1)((1)(1)))
  ((1)((1)((1)(1))))
  (((1)(1))((1)(1)))
		

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]]]];
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    olmsptrees[m_]:=Prepend[Union@@Table[Sort/@Tuples[olmsptrees/@p],{p,Select[mps[m],Length[#]>1&]}],m];
    Table[Length[olmsptrees[nrmptn[n]]],{n,15}]

Formula

a(n) = A292504(A181821(n)).
a(prime(n)) = A141268(n).
a(2^n) = A005804(n).

Extensions

More terms from Jinyuan Wang, Jun 26 2020

A323788 Number of non-isomorphic weight-n sets of multisets of multisets.

Original entry on oeis.org

1, 1, 5, 19, 88, 391, 1995, 10281
Offset: 0

Views

Author

Gus Wiseman, Jan 27 2019

Keywords

Comments

Also the number of non-isomorphic strict multiset partitions of multiset partitions of weight n.
All sets and multisets must be finite, and only the outermost may be empty.
The weight of an atom is 1, and the weight of a multiset is the sum of weights of its elements, counting multiplicity.

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(3) = 19 multiset partitions:
  {{1}}  {{11}}      {{111}}
         {{12}}      {{112}}
         {{1}{1}}    {{123}}
         {{1}{2}}    {{1}{11}}
         {{1}}{{2}}  {{1}{12}}
                     {{1}{23}}
                     {{2}{11}}
                     {{1}}{{11}}
                     {{1}{1}{1}}
                     {{1}}{{12}}
                     {{1}{1}{2}}
                     {{1}}{{23}}
                     {{1}{2}{3}}
                     {{2}}{{11}}
                     {{1}}{{1}{1}}
                     {{1}}{{1}{2}}
                     {{1}}{{2}{3}}
                     {{2}}{{1}{1}}
                     {{1}}{{2}}{{3}}
		

Crossrefs

A323789 Number of non-isomorphic weight-n sets of sets of multisets.

Original entry on oeis.org

1, 1, 4, 15, 64, 269, 1310, 6460
Offset: 0

Views

Author

Gus Wiseman, Jan 27 2019

Keywords

Comments

Also the number of non-isomorphic strict multiset partitions, with strict parts, of multiset partitions of weight n.
All sets and multisets must be finite, and only the outermost may be empty.
The weight of an atom is 1, and the weight of a multiset is the sum of weights of its elements, counting multiplicity.

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(3) = 15 multiset partition partitions:
  {{1}}  {{11}}      {{111}}
         {{12}}      {{112}}
         {{1}{2}}    {{123}}
         {{1}}{{2}}  {{1}{11}}
                     {{1}{12}}
                     {{1}{23}}
                     {{2}{11}}
                     {{1}}{{11}}
                     {{1}}{{12}}
                     {{1}}{{23}}
                     {{1}{2}{3}}
                     {{2}}{{11}}
                     {{1}}{{1}{2}}
                     {{1}}{{2}{3}}
                     {{1}}{{2}}{{3}}
		

Crossrefs

A323791 Number of non-isomorphic weight-n sets of multisets of sets.

Original entry on oeis.org

1, 1, 4, 13, 52, 196, 877, 3917
Offset: 0

Views

Author

Gus Wiseman, Jan 27 2019

Keywords

Comments

All sets and multisets must be finite, and only the outermost may be empty.
The weight of an atom is 1, and the weight of a multiset is the sum of weights of its elements, counting multiplicity.

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(3) = 13 sets of multisets of sets:
  {{1}}  {{12}}      {{123}}
         {{1}{1}}    {{1}{12}}
         {{1}{2}}    {{1}{23}}
         {{1}}{{2}}  {{1}{1}{1}}
                     {{1}}{{12}}
                     {{1}{1}{2}}
                     {{1}}{{23}}
                     {{1}{2}{3}}
                     {{1}}{{1}{1}}
                     {{1}}{{1}{2}}
                     {{1}}{{2}{3}}
                     {{2}}{{1}{1}}
                     {{1}}{{2}}{{3}}
		

Crossrefs

A330665 Number of balanced reduced multisystems of maximal depth whose atoms are the prime indices of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 5, 1, 1, 1, 2, 1, 3, 1, 5, 1, 1, 1, 7, 1, 1, 1, 5, 1, 3, 1, 2, 2, 1, 1, 16, 1, 2, 1, 2, 1, 5, 1, 5, 1, 1, 1, 11, 1, 1, 2, 16, 1, 3, 1, 2, 1, 3, 1, 27, 1, 1, 2, 2, 1, 3, 1, 16, 2, 1, 1, 11, 1
Offset: 1

Views

Author

Gus Wiseman, Dec 27 2019

Keywords

Comments

First differs from A317145 at a(32) = 5, A317145(32) = 4.
A balanced reduced multisystem is either a finite multiset, or a multiset partition with at least two parts, not all of which are singletons, of a balanced reduced multisystem.
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 series/singleton-reduced factorizations of n with Omega(n) levels of parentheses. See A001055, A050336, A050338, A050340, etc.

Examples

			The a(n) multisystems for n = 2, 6, 12, 24, 48:
  {1}  {1,2}  {{1},{1,2}}  {{{1}},{{1},{1,2}}}  {{{{1}}},{{{1}},{{1},{1,2}}}}
              {{2},{1,1}}  {{{1,1}},{{1},{2}}}  {{{{1}}},{{{1,1}},{{1},{2}}}}
                           {{{1}},{{2},{1,1}}}  {{{{1},{1}}},{{{1}},{{1,2}}}}
                           {{{1,2}},{{1},{1}}}  {{{{1},{1,1}}},{{{1}},{{2}}}}
                           {{{2}},{{1},{1,1}}}  {{{{1,1}}},{{{1}},{{1},{2}}}}
                                                {{{{1}}},{{{1}},{{2},{1,1}}}}
                                                {{{{1}}},{{{1,2}},{{1},{1}}}}
                                                {{{{1},{1}}},{{{2}},{{1,1}}}}
                                                {{{{1},{1,2}}},{{{1}},{{1}}}}
                                                {{{{1,1}}},{{{2}},{{1},{1}}}}
                                                {{{{1}}},{{{2}},{{1},{1,1}}}}
                                                {{{{1},{2}}},{{{1}},{{1,1}}}}
                                                {{{{1,2}}},{{{1}},{{1},{1}}}}
                                                {{{{2}}},{{{1}},{{1},{1,1}}}}
                                                {{{{2}}},{{{1,1}},{{1},{1}}}}
                                                {{{{2},{1,1}}},{{{1}},{{1}}}}
		

Crossrefs

The last nonzero term in row n of A330667 is a(n).
The chain version is A317145.
The non-maximal version is A318812.
Unlabeled versions are A330664 and A330663.
Other labeled versions are A330675 (strongly normal) and A330676 (normal).

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    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]]]];
    totm[m_]:=Prepend[Join@@Table[totm[p],{p,Select[mps[m],1
    				

Formula

a(2^n) = A000111(n - 1).
a(product of n distinct primes) = A006472(n).

A330676 Number of balanced reduced multisystems of weight n and maximum depth whose atoms cover an initial interval of positive integers.

Original entry on oeis.org

1, 1, 2, 8, 70, 1012, 21944, 665708, 26917492, 1399033348, 90878863352, 7214384973908, 687197223963640, 77354805301801012, 10158257981179981304, 1539156284259756811748, 266517060496258245459352, 52301515332984084095078308, 11546416513975694879642736152
Offset: 0

Views

Author

Gus Wiseman, Dec 30 2019

Keywords

Comments

A balanced reduced multisystem is either a finite multiset, or a multiset partition with at least two parts, not all of which are singletons, of a balanced reduced multisystem. The weight of an atom is 1, while the weight of a multiset is the sum of weights of its elements.
A finite multiset is normal if it covers an initial interval of positive integers.

Examples

			The a(0) = 1 through a(3) = 8 multisystems:
  {}  {1}  {1,1}  {{1},{1,1}}
           {1,2}  {{1},{1,2}}
                  {{1},{2,2}}
                  {{1},{2,3}}
                  {{2},{1,1}}
                  {{2},{1,2}}
                  {{2},{1,3}}
                  {{3},{1,2}}
		

Crossrefs

Row sums of A330778.
The case with all atoms equal is A000111.
The case with all atoms different is A006472.
The version allowing all depths is A330655.
The unlabeled version is A330663.
The version where the atoms are the prime indices of n is A330665.
The strongly normal version is A330675.
The version where the degrees are the prime indices of n is A330728.
Multiset partitions of normal multisets are A255906.
Series-reduced rooted trees with normal leaves are A316651.

Programs

  • Mathematica
    allnorm[n_]:=If[n<=0,{{}},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[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    totm[m_]:=Prepend[Join@@Table[totm[p],{p,Select[mps[m],1
    				
  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
    R(n, k)={my(v=vector(n), u=vector(n)); v[1]=k; for(n=1, #v, for(i=n, #v, u[i] += v[i]*(-1)^(i-n)*binomial(i-1, n-1)); v=EulerT(v)); u}
    seq(n)={concat([1], sum(k=1, n, R(n, k)*sum(r=k, n, binomial(r, k)*(-1)^(r-k))))} \\ Andrew Howroyd, Dec 30 2020

Extensions

Terms a(8) and beyond from Andrew Howroyd, Dec 30 2019

A318846 Number of balanced reduced multisystems whose atoms cover an initial interval of positive integers with multiplicities equal to the prime indices of n.

Original entry on oeis.org

1, 1, 1, 1, 2, 3, 6, 4, 15, 11, 20, 21, 90, 51, 80, 32, 468, 166, 2910, 124, 521, 277, 20644, 266, 621, 1761, 1866, 841, 165874, 1374, 1484344, 436, 3797, 12741, 5383, 3108, 14653890, 103783, 31323, 2294, 158136988, 12419, 1852077284, 6382, 20786, 939131, 23394406084
Offset: 1

Views

Author

Gus Wiseman, Sep 04 2018

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. A multiset whose multiplicities are the prime indices of n (such as row n of A305936) is generally not the same as the multiset of prime indices of n. For example, the prime indices of 12 are {1,1,2}, while a multiset whose multiplicities are {1,1,2} is {1,1,2,3}.
A balanced reduced multisystem is either a finite multiset, or a multiset partition with at least two parts, not all of which are singletons, of a balanced reduced multisystem.

Examples

			The a(12) = 21 multisystems on {1,1,2,3} (commas elided):
  {1123}  {{1}{123}}  {{1}{1}{23}}  {{{1}}{{1}{23}}}
          {{2}{113}}  {{1}{2}{13}}  {{{23}}{{1}{1}}}
          {{3}{112}}  {{1}{3}{12}}  {{{1}}{{2}{13}}}
          {{11}{23}}  {{2}{3}{11}}  {{{2}}{{1}{13}}}
          {{12}{13}}                {{{13}}{{1}{2}}}
                                    {{{1}}{{3}{12}}}
                                    {{{3}}{{1}{12}}}
                                    {{{12}}{{1}{3}}}
                                    {{{2}}{{3}{11}}}
                                    {{{3}}{{2}{11}}}
                                    {{{11}}{{2}{3}}}
		

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]]]];
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    tmsp[m_]:=Prepend[Join@@Table[tmsp[c],{c,Select[mps[m],1
    				

Formula

a(n) = A318812(A181821(n)).
a(prime(n)) = A318813(n).
a(2^n) = A005121(n).

Extensions

Terminology corrected by Gus Wiseman, Jan 04 2020
More terms from Jinyuan Wang, Jun 26 2020
Previous Showing 11-20 of 30 results. Next