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

A382912 Numbers k such that row k of A305936 (a multiset whose multiplicities are the prime indices of k) has no permutation with all distinct run-lengths.

Original entry on oeis.org

4, 8, 9, 12, 16, 18, 20, 24, 27, 28, 32, 36, 40, 44, 45, 48, 50, 52, 54, 56, 60, 63, 64, 68, 72, 75, 76, 80, 81, 84, 88, 90, 92, 96, 98, 99, 100, 104, 108, 112, 116, 117, 120, 124, 125, 126, 128, 132, 135, 136, 140, 144, 148, 150, 152, 153, 156, 160, 162, 164
Offset: 1

Views

Author

Gus Wiseman, Apr 12 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}.
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, sum A056239.

Examples

			The terms, prime indices, and corresponding multisets begin:
   4:       {1,1} {1,2}
   8:     {1,1,1} {1,2,3}
   9:       {2,2} {1,1,2,2}
  12:     {1,1,2} {1,1,2,3}
  16:   {1,1,1,1} {1,2,3,4}
  18:     {1,2,2} {1,1,2,2,3}
  20:     {1,1,3} {1,1,1,2,3}
  24:   {1,1,1,2} {1,1,2,3,4}
  27:     {2,2,2} {1,1,2,2,3,3}
  28:     {1,1,4} {1,1,1,1,2,3}
  32: {1,1,1,1,1} {1,2,3,4,5}
  36:   {1,1,2,2} {1,1,2,2,3,4}
  40:   {1,1,1,3} {1,1,1,2,3,4}
  44:     {1,1,5} {1,1,1,1,1,2,3}
  45:     {2,2,3} {1,1,1,2,2,3,3}
  48: {1,1,1,1,2} {1,1,2,3,4,5}
  50:     {1,3,3} {1,1,1,2,2,2,3}
  52:     {1,1,6} {1,1,1,1,1,1,2,3}
		

Crossrefs

The Look-and-Say partition is ranked by A048767, listed by A381440.
Look-and-Say partitions are counted by A239455, ranks A351294.
Non-Look-and-Say partitions are counted by A351293.
For prime indices instead of signature we have A351295, conjugate A381433.
The complement is A382913.
For equal instead of distinct run-lengths we have A382914, see A382858, A382879, A382915.
A056239 adds up prime indices, row sums of A112798.
A329739 counts compositions with distinct run-lengths, ranks A351596, complement A351291.
A381431 lists the section-sum partition of n, ranks A381436, union A381432.

Programs

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

A382913 Numbers k such that row k of A305936 (a multiset whose multiplicities are the prime indices of k) has a permutation with all distinct run-lengths.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 25, 26, 29, 30, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 46, 47, 49, 51, 53, 55, 57, 58, 59, 61, 62, 65, 66, 67, 69, 70, 71, 73, 74, 77, 78, 79, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 101, 102, 103
Offset: 1

Views

Author

Gus Wiseman, Apr 12 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}.
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, sum A056239.

Examples

			The terms, prime indices, and corresponding multisets begin:
   1:    {} {}
   2:   {1} {1}
   3:   {2} {1,1}
   5:   {3} {1,1,1}
   6: {1,2} {1,1,2}
   7:   {4} {1,1,1,1}
  10: {1,3} {1,1,1,2}
  11:   {5} {1,1,1,1,1}
  13:   {6} {1,1,1,1,1,1}
  14: {1,4} {1,1,1,1,2}
  15: {2,3} {1,1,1,2,2}
  17:   {7} {1,1,1,1,1,1,1}
  19:   {8} {1,1,1,1,1,1,1,1}
  21: {2,4} {1,1,1,1,2,2}
  22: {1,5} {1,1,1,1,1,2}
  23:   {9} {1,1,1,1,1,1,1,1,1}
  25: {3,3} {1,1,1,2,2,2}
  26: {1,6} {1,1,1,1,1,1,2}
		

Crossrefs

Look-and-Say partitions are counted by A239455, ranks A351294.
Non-Look-and-Say partitions are counted by A351293, ranks A351295.
For prime indices instead of signature we have A351294, conjugate A381432.
The Look-and-Say partition of n is listed by A381440, rank A048767.
The complement is A382912.
For equal run-lengths we have the complement of A382914, see A382858, A382879, A382915.
A044813 lists numbers whose binary expansion has distinct run-lengths.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A329739 counts compositions with distinct run-lengths, ranks A351596.
A381431 ranks section-sum partition, listed by A381436.

Programs

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

A386579 Number of permutations of row n of A305936 (a multiset whose multiplicities are the prime indices of n) with k adjacent unequal parts.

Original entry on oeis.org

1, 1, 0, 0, 2, 1, 0, 0, 0, 2, 1, 1, 0, 0, 0, 0, 0, 6, 0, 2, 2, 2, 0, 2, 2, 0, 1, 0, 0, 0, 0, 0, 0, 6, 6, 1, 0, 0, 0, 0, 0, 0, 2, 3, 0, 0, 0, 2, 3, 4, 1, 0, 0, 0, 24, 1, 0, 0, 0, 0, 0, 0, 0, 0, 6, 12, 12, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6, 12, 2, 0, 2, 4, 6, 3, 0
Offset: 2

Views

Author

Gus Wiseman, Aug 04 2025

Keywords

Comments

Row 1 is empty, so offset is 2.
Same as A386578 with rows reversed.
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

			Row n = 21 counts the following permutations:
  .  111122  111221  111212  112121  .
     221111  112211  112112  121121
             122111  121112  121211
             211112  211121
                     211211
                     212111
Triangle begins:
  .
  1
  1  0
  0  2
  1  0  0
  0  2  1
  1  0  0  0
  0  0  6
  0  2  2  2
  0  2  2  0
  1  0  0  0  0
  0  0  6  6
  1  0  0  0  0  0
  0  2  3  0  0
  0  2  3  4  1
  0  0  0 24
  1  0  0  0  0  0  0
  0  0  6 12 12
  1  0  0  0  0  0  0  0
  0  0  6 12  2
  0  2  4  6  3  0
		

Crossrefs

Column k = 0 is A010051.
Row lengths are A056239.
Row sums are A318762.
Column k = last is A335125.
For prime indices we have A374252, reverse A386577.
Reversing all rows gives A386578.
A003242 and A335452 count anti-runs, ranks A333489, patterns A005649.
A025065(n - 2) counts partitions of inseparable type, ranks A335126, sums of A386586.
A124762 gives inseparability of standard compositions, separability A333382.
A305936 is a multiset whose multiplicities are the prime indices of n.
A325534 counts separable multisets, ranks A335433, sums of A386583.
A325535 counts inseparable multisets, ranks A335448, sums of A386584.
A336106 counts partitions of separable type, ranks A335127, sums of A386585.

Programs

  • Mathematica
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    ugt[c_,x_]:=Select[Permutations[c],Function[q,Length[Select[Range[Length[q]-1],q[[#]]!=q[[#+1]]&]]==x]];
    Table[Table[Length[ugt[nrmptn[n],k]],{k,0,Length[nrmptn[n]]-1}],{n,30}]

A386578 Irregular triangle read by rows where T(n,k) is the number of permutations of row n of A305936 (a multiset whose multiplicities are the prime indices of n) with k adjacent equal parts.

Original entry on oeis.org

1, 0, 1, 2, 0, 0, 0, 1, 1, 2, 0, 0, 0, 0, 1, 6, 0, 0, 2, 2, 2, 0, 0, 2, 2, 0, 0, 0, 0, 0, 1, 6, 6, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 3, 2, 0, 1, 4, 3, 2, 0, 24, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 12, 12, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 2, 12, 6, 0, 0, 0, 3, 6, 4, 2, 0
Offset: 2

Views

Author

Gus Wiseman, Aug 04 2025

Keywords

Comments

Row 1 is empty, so offset is 2.
Same as A386579 with rows reversed.
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

			Row n = 21 counts the following permutations:
  .  112121  111212  111221  111122  .
     121121  112112  112211  221111
     121211  121112  122111
             211121  211112
             211211
             212111
Triangle begins
   .
   1
   0  1
   2  0
   0  0  1
   1  2  0
   0  0  0  1
   6  0  0
   2  2  2  0
   0  2  2  0
   0  0  0  0  1
   6  6  0  0
   0  0  0  0  0  1
   0  0  3  2  0
   1  4  3  2  0
  24  0  0  0
   0  0  0  0  0  0  1
  12 12  6  0  0
   0  0  0  0  0  0  0  1
   2 12  6  0  0
   0  3  6  4  2  0
		

Crossrefs

Column k = last is A010051.
Row lengths are A056239.
Initial zeros are counted by A252736 = A001222 - 1.
Row sums are A318762.
Column k = 0 is A335125.
For prime indices we have A386577.
Reversing all rows gives A386579.
A003242 and A335452 count anti-runs, ranks A333489, patterns A005649.
A025065(n - 2) counts partitions of inseparable type, ranks A335126, sums of A386586.
A124762 gives inseparability of standard compositions, separability A333382.
A305936 is a multiset whose multiplicities are the prime indices of n.
A325534 counts separable multisets, ranks A335433, sums of A386583.
A325535 counts inseparable multisets, ranks A335448, sums of A386584.
A336106 counts partitions of separable type, ranks A335127, sums of A386585.

Programs

  • Mathematica
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    aqt[c_,x_]:=Select[Permutations[c],Function[q,Length[Select[Range[Length[q]-1],q[[#]]==q[[#+1]]&]]==x]];
    Table[Table[Length[aqt[nrmptn[n],k]],{k,0,Length[nrmptn[n]]-1}],{n,30}]

A321188 Number of set systems with no singletons whose multiset union is row n of A305936 (a multiset whose multiplicities are the prime indices of n).

Original entry on oeis.org

1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 4, 0, 1, 0, 0, 0, 0, 0, 3, 0, 0, 1, 0, 0, 0, 0, 11, 0, 0, 0, 4, 0, 0, 0, 1
Offset: 1

Views

Author

Gus Wiseman, Oct 29 2018

Keywords

Comments

A set system is a finite set of finite nonempty sets.
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.

Examples

			The a(36) = 4 set systems with no singletons whose multiset union is {1,1,2,2,3,4}:
  {{1,2},{1,2,3,4}}
  {{1,2,3},{1,2,4}}
  {{1,2},{1,3},{2,4}}
  {{1,2},{1,4},{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]]]];
    hyp[m_]:=Select[mps[m],And[And@@UnsameQ@@@#,UnsameQ@@#,Min@@Length/@#>1]&];
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Table[Length[hyp[nrmptn[n]]],{n,30}]

A181821 a(n) = smallest integer with factorization as Product p(i)^e(i) such that Product p(e(i)) = n.

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 16, 30, 36, 24, 32, 60, 64, 48, 72, 210, 128, 180, 256, 120, 144, 96, 512, 420, 216, 192, 900, 240, 1024, 360, 2048, 2310, 288, 384, 432, 1260, 4096, 768, 576, 840, 8192, 720, 16384, 480, 1800, 1536, 32768, 4620, 1296, 1080, 1152, 960, 65536
Offset: 1

Views

Author

Matthew Vandermast, Dec 07 2010

Keywords

Comments

A permutation of A025487. a(n) is the member m of A025487 such that A181819(m) = n. a(n) is also the member of A025487 whose prime signature is conjugate to the prime signature of A108951(n).
If n = Product_i prime(e(i)) with the e(i) weakly decreasing, then a(n) = Product_i prime(i)^e(i). For example, 90 = prime(3) * prime(2) * prime(2) * prime(1), so a(90) = prime(1)^3 * prime(2)^2 * prime(3)^2 * prime(4)^1 = 12600. - Gus Wiseman, Jan 02 2019

Examples

			The canonical factorization of 24 is 2^3*3^1. Therefore, p(e(i)) = prime(3)*prime(1)(i.e., A000040(3)*A000040(1)), which equals 5*2 = 10. Since 24 is the smallest integer for which p(e(i)) = 10, a(10) = 24.
		

Crossrefs

Programs

  • Maple
    a:= n-> (l-> mul(ithprime(i)^l[i], i=1..nops(l)))(sort(map(i->
                 numtheory[pi](i[1])$i[2], ifactors(n)[2]), `>`)):
    seq(a(n), n=1..70);  # Alois P. Heinz, Sep 05 2018
  • Mathematica
    With[{s = Array[If[# == 1, 1, Times @@ Map[Prime@ Last@ # &, FactorInteger@ #]] &, 2^16]}, Array[First@ FirstPosition[s, #] &, LengthWhile[Differences@ Union@ s, # == 1 &]]] (* Michael De Vlieger, Dec 17 2018 *)
    Table[Times@@MapIndexed[Prime[#2[[1]]]^#1&,Reverse[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]],{n,30}] (* Gus Wiseman, Jan 02 2019 *)
  • PARI
    A181821(n) = { my(f=factor(n),p=0,m=1); forstep(i=#f~,1,-1,while(f[i,2], f[i,2]--; m *= (p=nextprime(p+1))^primepi(f[i,1]))); (m); }; \\ Antti Karttunen, Dec 10 2018
    
  • Python
    from math import prod
    from sympy import prime, primepi, factorint
    def A181821(n): return prod(prime(i)**e for i, e in enumerate(sorted(map(primepi,factorint(n,multiple=True)),reverse=True),1)) # Chai Wah Wu, Sep 15 2023

Formula

If A108951(n) = Product p(i)^e(i), then a(n) = Product A002110(e(i)). I.e., a(n) = A108951(A181819(A108951(n))).
a(A181819(n)) = A046523(n). - [See also A124859]. Antti Karttunen, Dec 10 2018
a(n) = A025487(A361808(n)). - Pontus von Brömssen, Mar 25 2023
a(n) = A108951(A122111(n)). - Antti Karttunen, Sep 15 2023

Extensions

Definition corrected by Gus Wiseman, Jan 02 2019

A318283 Sum of elements of the multiset spanning an initial interval of positive integers with multiplicities equal to the prime indices of n in weakly decreasing order.

Original entry on oeis.org

0, 1, 2, 3, 3, 4, 4, 6, 6, 5, 5, 7, 6, 6, 7, 10, 7, 9, 8, 8, 8, 7, 9, 11, 9, 8, 12, 9, 10, 10, 11, 15, 9, 9, 10, 13, 12, 10, 10, 12, 13, 11, 14, 10, 13, 11, 15, 16, 12, 12, 11, 11, 16, 16, 11, 13, 12, 12, 17, 14, 18, 13, 14, 21, 12, 12, 19, 12, 13, 13, 20, 18
Offset: 1

Views

Author

Gus Wiseman, Aug 23 2018

Keywords

Examples

			The multiset spanning an initial interval of positive integers with multiplicities equal to the prime indices of 90 in weakly decreasing order is {1,1,1,2,2,3,3,4}, so a(90) = 1+1+1+2+2+3+3+4 = 17.
		

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}]]]]];
    Total/@Array[nrmptn,100]

Formula

a(n) = A056239(A181821(n)).

A318284 Number of 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, 5, 9, 7, 7, 11, 11, 12, 16, 15, 15, 26, 22, 21, 29, 19, 30, 36, 31, 30, 66, 38, 42, 52, 56, 52, 47, 45, 57, 92, 77, 67, 77, 74, 101, 98, 135, 64, 137, 97, 176, 135, 109, 109, 118, 105, 231, 249, 97, 141, 181, 139, 297, 198, 385, 195, 269
Offset: 1

Views

Author

Gus Wiseman, Aug 23 2018

Keywords

Examples

			The a(12) = 11 multiset partitions of {1,1,2,3}:
  {{1,1,2,3}}
  {{1},{1,2,3}}
  {{2},{1,1,3}}
  {{3},{1,1,2}}
  {{1,1},{2,3}}
  {{1,2},{1,3}}
  {{1},{1},{2,3}}
  {{1},{2},{1,3}}
  {{1},{3},{1,2}}
  {{2},{3},{1,1}}
  {{1},{1},{2},{3}}
		

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}]]]]];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[facs[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(n=vecsum(sig), A=O(x*x^vecmax(sig)), 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]))*permcount(p)); s/n!}
    a(n)={if(n==1, 1, my(s=sig(n)); if(#s==1, numbpart(s[1]), count(sig(n))))} \\ Andrew Howroyd, Dec 10 2018

Formula

a(n) = A001055(A181821(n)).
a(prime(n)^k) = A219727(n,k). - Andrew Howroyd, Dec 10 2018

A335126 A multiset whose multiplicities are the prime indices of n is inseparable.

Original entry on oeis.org

3, 5, 7, 10, 11, 13, 14, 17, 19, 21, 22, 23, 26, 28, 29, 31, 33, 34, 37, 38, 39, 41, 43, 44, 46, 47, 51, 52, 53, 55, 57, 58, 59, 61, 62, 65, 66, 67, 68, 69, 71, 73, 74, 76, 78, 79, 82, 83, 85, 86, 87, 88, 89, 91, 92, 93, 94, 95, 97, 101, 102, 103, 104, 106
Offset: 1

Views

Author

Gus Wiseman, Jul 01 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 of terms together with the corresponding multisets begins:
   3: {1,1}
   5: {1,1,1}
   7: {1,1,1,1}
  10: {1,1,1,2}
  11: {1,1,1,1,1}
  13: {1,1,1,1,1,1}
  14: {1,1,1,1,2}
  17: {1,1,1,1,1,1,1}
  19: {1,1,1,1,1,1,1,1}
  21: {1,1,1,1,2,2}
  22: {1,1,1,1,1,2}
  23: {1,1,1,1,1,1,1,1,1}
  26: {1,1,1,1,1,1,2}
  28: {1,1,1,1,2,3}
  29: {1,1,1,1,1,1,1,1,1,1}
		

Crossrefs

The complement is A335127.
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_,_}]&]=={}&]

A318762 Number of permutations of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 3, 1, 6, 6, 4, 1, 12, 1, 5, 10, 24, 1, 30, 1, 20, 15, 6, 1, 60, 20, 7, 90, 30, 1, 60, 1, 120, 21, 8, 35, 180, 1, 9, 28, 120, 1, 105, 1, 42, 210, 10, 1, 360, 70, 140, 36, 56, 1, 630, 56, 210, 45, 11, 1, 420, 1, 12, 420, 720, 84, 168, 1, 72, 55
Offset: 1

Views

Author

Gus Wiseman, Sep 03 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}.

Examples

			The a(12) = 12 permutations are (1123), (1132), (1213), (1231), (1312), (1321), (2113), (2131), (2311), (3112), (3121), (3211).
		

Crossrefs

Programs

  • Maple
    a:= n-> (l-> add(i, i=l)!/mul(i!, i=l))(map(i->
           numtheory[pi](i[1])$i[2], ifactors(n)[2])):
    seq(a(n), n=1..100);  # Alois P. Heinz, Sep 03 2018
  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Total[primeMS[n]]!/Times@@Factorial/@primeMS[n],{n,100}]
  • PARI
    sig(n)={my(f=factor(n)); concat(vector(#f~, i, vector(f[i, 2], j, primepi(f[i, 1]))))}
    a(n)={if(n==1, 1, my(s=sig(n)); vecsum(s)!/prod(i=1, #s, s[i]!))}  \\ Andrew Howroyd, Dec 17 2018

Formula

If n = Product prime(x_i)^y_i is the prime factorization of n, then a(n) = (Sum x_i * y_i)! / Product (x_i!)^y_i.
a(n) = A008480(A181821(n)).
a(n) = A112624(n) * A124794(n). - Max Alekseyev, Oct 15 2023
Sum_{m in row n of A215366} a(m) = A005651(n).
Sum_{m in row n of A215366} a(m) * A008480(m) = A000670(n).
Sum_{m in row n of A215366} a(m) * A008480(m) / A001222(m)! = A000110(n).
Showing 1-10 of 53 results. Next