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

A335127 A multiset whose multiplicities are the prime indices of n is separable.

Original entry on oeis.org

1, 2, 4, 6, 8, 9, 12, 15, 16, 18, 20, 24, 25, 27, 30, 32, 35, 36, 40, 42, 45, 48, 49, 50, 54, 56, 60, 63, 64, 70, 72, 75, 77, 80, 81, 84, 90, 96, 98, 99, 100, 105, 108, 110, 112, 120, 121, 125, 126, 128, 132, 135, 140, 143, 144, 147, 150, 154, 160, 162, 165
Offset: 1

Views

Author

Gus Wiseman, Jul 02 2020

Keywords

Comments

A multiset is separable if it has a permutation that is an anti-run, meaning there are no adjacent equal parts.
A multiset whose multiplicities are the prime indices of n (such as row n of A305936) is 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}.

Examples

			The sequence together with the corresponding multisets begins:
   1: {}
   2: {1}
   4: {1,2}
   6: {1,1,2}
   8: {1,2,3}
   9: {1,1,2,2}
  12: {1,1,2,3}
  15: {1,1,1,2,2}
  16: {1,2,3,4}
  18: {1,1,2,2,3}
  20: {1,1,1,2,3}
  24: {1,1,2,3,4}
  25: {1,1,1,2,2,2}
  27: {1,1,2,2,3,3}
  30: {1,1,1,2,2,3}
		

Crossrefs

The complement is A335126.
Anti-run compositions are A003242.
Anti-runs are ranked by A333489.
Separable partitions are A325534.
Inseparable partitions are A325535.
Separable factorizations are A335434.
Inseparable factorizations are A333487.
Separable partitions are ranked by A335433.
Inseparable partitions are ranked by A335448.
Anti-run permutations of prime indices are A335452.
Patterns contiguously matched by compositions are A335457.

Programs

  • Mathematica
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Select[Range[100],Select[Permutations[nrmptn[#]],!MatchQ[#,{_,x_,x_,_}]&]!={}&]

A318286 Number of strict multiset partitions of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 2, 5, 5, 5, 3, 9, 4, 7, 9, 15, 5, 18, 6, 16, 14, 10, 8, 31, 17, 14, 40, 25, 10, 34, 12, 52, 21, 19, 27, 70, 15, 25, 31, 59, 18, 57, 22, 38, 80, 33, 27, 120, 46, 67, 44, 56, 32, 172, 42, 100, 61, 43, 38, 141, 46, 55, 143, 203, 64, 91, 54, 80
Offset: 1

Views

Author

Gus Wiseman, Aug 23 2018

Keywords

Crossrefs

Programs

  • Mathematica
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    strfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strfacs[n/d],Min@@#>d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[strfacs[Times@@Prime/@nrmptn[n]]],{n,60}]
  • PARI
    permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
    sig(n)={my(f=factor(n)); concat(vector(#f~, i, vector(f[i, 2], j, primepi(f[i, 1]))))}
    count(sig)={my(r=0, A=O(x*x^vecmax(sig))); for(n=1, vecsum(sig)+1, my(s=0); forpart(p=n, my(q=1/prod(i=1, #p, 1 - x^p[i] + A)); s+=prod(i=1, #sig, polcoef(q, sig[i]))*(-1)^#p*permcount(p)); r+=(-1)^n*s/n!); r/2}
    a(n)={if(n==1, 1, count(sig(n)))} \\ Andrew Howroyd, Dec 18 2018

Formula

a(n) = A045778(A181821(n)).
a(prime(n)^k) = A219585(n, k). - Andrew Howroyd, Dec 17 2018

A335125 Number of anti-run permutations of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 0, 2, 0, 1, 0, 6, 2, 0, 0, 6, 0, 0, 1, 24, 0, 12, 0, 2, 0, 0, 0, 36, 2, 0, 30, 0, 0, 10, 0, 120, 0, 0, 1, 84, 0, 0, 0, 24, 0, 3, 0, 0, 38, 0, 0, 240, 2, 18, 0, 0, 0, 246, 0, 6, 0, 0, 0, 96, 0, 0, 24, 720, 0, 0, 0, 0, 0, 14, 0, 660, 0, 0, 74, 0, 1, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Jul 01 2020

Keywords

Comments

A multiset whose multiplicities are the prime indices of n (such as row n of A305936) is 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}.
An anti-run is a sequence with no adjacent equal parts.

Examples

			The a(n) permutations for n = 2, 4, 42, 8, 30, 18:
  (1)  (12)  (1212131)  (123)  (121213)  (12123)
       (21)  (1213121)  (132)  (121231)  (12132)
             (1312121)  (213)  (121312)  (12312)
                        (231)  (121321)  (12321)
                        (312)  (123121)  (13212)
                        (321)  (131212)  (21213)
                               (132121)  (21231)
                               (212131)  (21312)
                               (213121)  (21321)
                               (312121)  (23121)
                                         (31212)
                                         (32121)
		

Crossrefs

Positions of zeros are A335126.
Positions of nonzeros are A335127.
The version for the prime indices themselves is A335452.
Anti-run compositions are A003242.
Anti-runs are ranked by A333489.
Separable partitions are ranked by A335433.
Separable factorizations are A335434.
Inseparable partitions are ranked by A335448.
Patterns contiguously matched by compositions are A335457.
Strict permutations of prime indices are A335489.

Programs

  • Mathematica
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Table[Length[Select[Permutations[nrmptn[n]],!MatchQ[#,{_,x_,x_,_}]&]],{n,100}]

A318285 Number of non-isomorphic multiset partitions of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 2, 2, 3, 4, 5, 3, 7, 7, 7, 9, 11, 12, 16, 5, 15, 17, 22, 16, 29, 19, 30, 16, 21, 30, 23, 29, 42, 52, 56, 7, 47, 45, 57, 43, 77, 67, 77, 31, 101, 98, 135, 47, 85, 97, 176, 29, 66, 64, 118, 77, 231, 69, 97, 57, 181, 139, 297, 137, 385, 195, 166, 11, 162, 171, 490, 118
Offset: 1

Views

Author

Gus Wiseman, Aug 23 2018

Keywords

Examples

			Non-isomorphic representatives of the a(12) = 9 multiset partitions of {1,1,2,3}:
  {{1,1,2,3}}
  {{1},{1,2,3}}
  {{2},{1,1,3}}
  {{1,1},{2,3}}
  {{1,2},{1,3}}
  {{1},{1},{2,3}}
  {{1},{2},{1,3}}
  {{2},{3},{1,1}}
  {{1},{1},{2},{3}}
		

Crossrefs

Programs

  • PARI
    \\ See links in A339645 for combinatorial species functions.
    sig(n)={my(f=factor(n), sig=vector(primepi(vecmax(f[,1])))); for(i=1, #f~, sig[primepi(f[i,1])]=f[i,2]); sig}
    C(sig)={my(n=sum(i=1, #sig, i*sig[i]), A=Vec(symGroupSeries(n)-1), B=O(x*x^n), c=prod(i=1, #sig, if(sig[i], sApplyCI(A[sig[i]], sig[i], A[i], i), 1))); polcoef(OgfSeries(sCartProd(c*x^n + B, sExp(x*Ser(A) + B))), n)}
    a(n)={if(n==1, 1, C(sig(n)))} \\ Andrew Howroyd, Jan 17 2023

Formula

a(n) = A317791(A181821(n)).

Extensions

Terms a(31) and beyond from Andrew Howroyd, Jan 17 2023

A095684 Triangle read by rows. There are 2^(m-1) rows of length m, for m = 1, 2, 3, ... The rows are in lexicographic order. The rows have the property that the first entry is 1, the second distinct entry (reading from left to right) is 2, the third distinct entry is 3, etc.

Original entry on oeis.org

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

Views

Author

N. J. A. Sloane, Jun 25 2004

Keywords

Comments

Row k is the unique multiset that covers an initial interval of positive integers and has multiplicities equal to the parts of the k-th composition in standard order (graded reverse-lexicographic, A066099). This composition is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. For example, the 13th composition is (1,2,1), so row 13 is {1,2,2,3}. - Gus Wiseman, Apr 26 2020

Examples

			1, 11, 12, 111, 112, 122, 123, 1111, 1112, 1122, 1123, 1222, 1223, 1233, ...
The 8 strings of length 4 are 1111, 1112, 1122, 1123, 1222, 1223, 1233, 1234.
From _Gus Wiseman_, Apr 26 2020: (Start)
The triangle read by columns begins:
  1:{1}  2:{1,1}  4:{1,1,1}   8:{1,1,1,1}  16:{1,1,1,1,1}
         3:{1,2}  5:{1,1,2}   9:{1,1,1,2}  17:{1,1,1,1,2}
                  6:{1,2,2}  10:{1,1,2,2}  18:{1,1,1,2,2}
                  7:{1,2,3}  11:{1,1,2,3}  19:{1,1,1,2,3}
                             12:{1,2,2,2}  20:{1,1,2,2,2}
                             13:{1,2,2,3}  21:{1,1,2,2,3}
                             14:{1,2,3,3}  22:{1,1,2,3,3}
                             15:{1,2,3,4}  23:{1,1,2,3,4}
                                           24:{1,2,2,2,2}
                                           25:{1,2,2,2,3}
                                           26:{1,2,2,3,3}
                                           27:{1,2,2,3,4}
                                           28:{1,2,3,3,3}
                                           29:{1,2,3,3,4}
                                           30:{1,2,3,4,4}
                                           31:{1,2,3,4,5}
(End)
		

Crossrefs

See A096299 for another version.
The number of distinct parts in row n is A000120(n), also the maximum part.
Row sums are A029931.
Heinz numbers of rows are A057335.
Row lengths are A070939.
Row products are A284001.
The version for prime indices is A305936.
There are A333942(n) multiset partitions of row n.
Multisets of compositions are counted by A034691.
Combinatory separations of normal multisets are A269134.
All of the following pertain to compositions in standard order (A066099):
- Necklaces are A065609.
- Strict compositions are A233564.
- Constant compositions are A272919.
- Lyndon words are A275692.
- Dealings are counted by A333939.
- Distinct parts are counted by A334028.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    ptnToNorm[y_]:=Join@@Table[ConstantArray[i,y[[i]]],{i,Length[y]}];
    Table[ptnToNorm[stc[n]],{n,15}] (* Gus Wiseman, Apr 26 2020 *)

A323307 Number of ways to fill a matrix with the parts of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 2, 4, 2, 6, 3, 12, 18, 12, 2, 36, 4, 10, 20, 72, 2, 60, 4, 40, 60, 24, 3, 120, 80, 14, 360, 120, 4, 240, 2, 240, 42, 32, 70, 720, 6, 27, 112, 480, 2, 210, 4, 84, 420, 40, 4, 1440, 280, 280, 108, 224, 5, 1260, 224, 420, 180, 22, 2, 840, 6, 72, 1680, 2880
Offset: 1

Views

Author

Gus Wiseman, Jan 13 2019

Keywords

Comments

This multiset (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}.

Examples

			The a(22) = 24 matrices:
  [111112] [111121] [111211] [112111] [121111] [211111]
.
  [111] [111] [111] [112] [121] [211]
  [112] [121] [211] [111] [111] [111]
.
  [11] [11] [11] [11] [12] [21]
  [11] [11] [12] [21] [11] [11]
  [12] [21] [11] [11] [11] [11]
.
  [1] [1] [1] [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] [1]
  [2] [1] [1] [1] [1] [1]
		

Crossrefs

Programs

  • Mathematica
    primeMS[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]]}]];
    ptnmats[n_]:=Union@@Permutations/@Select[Union@@(Tuples[Permutations/@#]&/@Map[primeMS,facs[n],{2}]),SameQ@@Length/@#&];
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Array[Length[ptnmats[Times@@Prime/@nrmptn[#]]]&,30]

Formula

a(n) = A318762(n) * A000005(A056239(n)).

A382773 Number of ways to permute a multiset whose multiplicities are the prime indices of n so that the run-lengths are all different.

Original entry on oeis.org

1, 1, 1, 0, 1, 2, 1, 0, 0, 2, 1, 0, 1, 2, 2, 0, 1, 0, 1, 0, 4, 4, 1, 0, 4, 4, 0, 0, 1, 6, 1, 0, 4, 6, 4, 0, 1, 6, 4, 0, 1, 6, 1, 0, 0, 8, 1, 0, 4, 0, 6, 0, 1, 0, 6, 0, 6, 8, 1, 0, 1, 10, 0, 0, 8, 6, 1, 0, 8, 6, 1, 0, 1, 10, 0, 0, 6, 6, 1, 0, 0, 12, 1, 0, 16
Offset: 1

Views

Author

Gus Wiseman, Apr 09 2025

Keywords

Comments

This described multiset (row n of A305936, Heinz number A181821) is generally not the same as the multiset of prime indices of n (A112798). 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}.

Examples

			The a(n) partitions for n = 6, 21, 30, 46:
  (1,1,2)  (1,1,1,1,2,2)  (1,1,1,2,2,3)  (1,1,1,1,1,1,1,1,1,2)
  (2,1,1)  (1,1,1,2,2,1)  (1,1,1,3,2,2)  (1,1,1,1,1,1,1,2,1,1)
           (1,2,2,1,1,1)  (2,2,1,1,1,3)  (1,1,1,1,1,1,2,1,1,1)
           (2,2,1,1,1,1)  (2,2,3,1,1,1)  (1,1,1,1,1,2,1,1,1,1)
                          (3,1,1,1,2,2)  (1,1,1,1,2,1,1,1,1,1)
                          (3,2,2,1,1,1)  (1,1,1,2,1,1,1,1,1,1)
                                         (1,1,2,1,1,1,1,1,1,1)
                                         (2,1,1,1,1,1,1,1,1,1)
		

Crossrefs

Positions of 1 are A008578.
For anti-run permutations we have A335125.
For just prime indices we have A382771, firsts A382772, equal A382857.
These permutations for factorials are counted by A382774, equal A335407.
For equal instead of distinct run-lengths we have A382858.
Positions of 0 are A382912, complement A382913.
A044813 lists numbers whose binary expansion has distinct run-lengths, equal A140690.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A098859 counts partitions with distinct multiplicities, ordered A242882.
A239455 counts Look-and-Say partitions, ranks A351294, conjugate A381432.
A329738 counts compositions with equal run-lengths, ranks A353744.
A329739 counts compositions with distinct run-lengths, ranks A351596.
A351293 counts non-Look-and-Say partitions, ranks A351295, conjugate A381433.

Programs

  • Mathematica
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Table[Length[Select[Permutations[nrmptn[n]],UnsameQ@@Length/@Split[#]&]],{n,100}]

Formula

a(n) = A382771(A181821(n)) = A382771(A304660(n)).

A284001 a(n) = A005361(A283477(n)).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 4, 6, 1, 2, 4, 6, 8, 12, 18, 24, 1, 2, 4, 6, 8, 12, 18, 24, 16, 24, 36, 48, 54, 72, 96, 120, 1, 2, 4, 6, 8, 12, 18, 24, 16, 24, 36, 48, 54, 72, 96, 120, 32, 48, 72, 96, 108, 144, 192, 240, 162, 216, 288, 360, 384, 480, 600, 720, 1, 2, 4, 6, 8, 12, 18, 24, 16, 24, 36, 48, 54, 72, 96, 120, 32, 48, 72, 96, 108, 144, 192, 240, 162, 216, 288, 360, 384, 480
Offset: 0

Views

Author

Antti Karttunen, Mar 18 2017

Keywords

Comments

a(n) is the product of elements of the multiset that covers an initial interval of positive integers with multiplicities equal to the parts of the n-th composition in standard order (graded reverse-lexicographic, A066099). This composition is obtained by taking the set of positions of 1's in the reversed binary expansion of n, prepending 0, taking first differences, and reversing again. For example, the 13th composition is (1,2,1) giving the multiset {1,2,2,3} with product 12, so a(13) = 12. - Gus Wiseman, Apr 26 2020

Crossrefs

Row products of A095684.
All of the following pertain to compositions in standard order (A066099):
- Length is A000120.
- Weighted sum is A029931.
- Necklaces are A065609.
- Sum is A070939.
- Strict compositions are A233564.
- Constant compositions are A272919.
- Lyndon words are A275692.
- Distinct parts are counted by A334028.

Programs

  • Mathematica
    Table[Times @@ FactorInteger[#][[All, -1]] &[Times @@ Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e == 1 :> {Times @@ Prime@ Range@ PrimePi@ p, e}] &[Times @@ Prime@ Flatten@ Position[#, 1] &@ Reverse@ IntegerDigits[n, 2]]], {n, 0, 93}] (* Michael De Vlieger, Mar 18 2017 *)
  • PARI
    A005361(n) = factorback(factor(n)[, 2]); \\ From A005361
    A034386(n) = prod(i=1, primepi(n), prime(i));
    A108951(n) = { my(f=factor(n)); prod(i=1, #f~, A034386(f[i, 1])^f[i, 2]) };  \\ From A108951
    A019565(n) = {my(j,v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ From A019565
    A283477(n) = A108951(A019565(n));
    A284001(n) = A005361(A283477(n));
    
  • Scheme
    (define (A284001 n) (A005361 (A283477 n)))

Formula

a(n) = A005361(A283477(n)).
a(n) = A003963(A057335(n)). - Gus Wiseman, Apr 26 2020
a(n) = A284005(A053645(n)) for n > 0 with a(0) = 1. - Mikhail Kurkov, Jun 05 2021 [verification needed]

A382858 Number of ways to permute a multiset whose multiplicities are the prime indices of n so that the run-lengths are all equal.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 6, 4, 0, 1, 6, 1, 0, 1, 24, 1, 12, 1, 2, 1, 0, 1, 36, 4, 0, 36, 0, 1, 10, 1, 120, 0, 0, 1, 84, 1, 0, 0, 24, 1, 3, 1, 0, 38, 0, 1, 240, 6, 18, 0, 0, 1, 246, 0, 6, 0, 0, 1, 96, 1, 0, 30, 720, 1, 0, 1, 0, 0, 14, 1, 660, 1, 0, 74, 0, 1, 0, 1
Offset: 1

Views

Author

Gus Wiseman, Apr 09 2025

Keywords

Comments

This described multiset (row n of A305936, Heinz number A181821) is generally not the same as the multiset of prime indices of n (A112798). 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}.

Examples

			The a(9) = 4 permutations are:
  (1,1,2,2)
  (1,2,1,2)
  (2,1,2,1)
  (2,2,1,1)
		

Crossrefs

The anti-run case is A335125.
These permutations for factorials are counted by A335407, distinct A382774.
For distinct instead of equal run-lengths we have A382773.
For prime indices we have A382857 (firsts A382878), distinct A382771 (firsts A382772).
Positions of 0 are A382914, signature restriction of A382915.
A003963 gives product of prime indices.
A140690 lists numbers whose binary expansion has equal run-lengths, distinct A044813.
A047966 counts partitions with equal multiplicities, distinct A098859.
A056239 adds up prime indices, row sums of A112798.
A304442 counts partitions with equal run-sums, ranks A353833.
A329738 counts compositions with equal run-lengths, ranks A353744.
A329739 counts compositions with distinct run-lengths, ranks A351596, complement A351291.
A382913 ranks Look-and-Say partitions by signature, complement A382912.

Programs

  • Mathematica
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Table[Length[Select[Permutations[nrmptn[n]],SameQ@@Length/@Split[#]&]],{n,100}]

Formula

a(n) = A382857(A181821(n)) = A382857(A304660(n)).

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
Previous Showing 11-20 of 53 results. Next