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

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

A305936 Irregular triangle whose n-th row is the multiset spanning an initial interval of positive integers with multiplicities equal to the n-th row of A296150 (the prime indices of n in weakly decreasing order).

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Aug 23 2018

Keywords

Examples

			Row 90 is {1,1,1,2,2,3,3,4} because 90 = prime(3)*prime(2)*prime(2)*prime(1).
Triangle begins:
   1:
   2:  1
   3:  1  1
   4:  1  2
   5:  1  1  1
   6:  1  1  2
   7:  1  1  1  1
   8:  1  2  3
   9:  1  1  2  2
  10:  1  1  1  2
  11:  1  1  1  1  1
  12:  1  1  2  3
  13:  1  1  1  1  1  1
		

Crossrefs

Row lengths are A056239. Number of distinct elements in row n is A001222(n). Number of distinct multiplicities in row n is A001221(n).

Programs

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

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

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

A318287 Number of non-isomorphic 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, 3, 4, 5, 3, 7, 4, 7, 9, 5, 5, 12, 6, 12, 14, 10, 8, 13, 12, 14, 14, 18, 10, 34
Offset: 1

Views

Author

Gus Wiseman, Aug 23 2018

Keywords

Examples

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

Crossrefs

Formula

a(n) = A318357(A181821(n)).

A318362 Number of non-isomorphic set multipartitions (multisets of sets) of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Aug 24 2018

Keywords

Examples

			Non-isomorphic representatives of the a(12) = 5 set multipartitions of {1,1,2,3}:
  {{1},{1,2,3}}
  {{1,2},{1,3}}
  {{1},{1},{2,3}}
  {{1},{2},{1,3}}
  {{1},{1},{2},{3}}
		

Crossrefs

Formula

a(n) = A318369(A181821(n)).

A330664 Number of non-isomorphic balanced reduced multisystems of maximum depth whose degrees (atom multiplicities) are the weakly decreasing prime indices of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 2, 1, 4, 5, 5, 7, 16, 16, 27, 2, 61, 33, 272, 27, 123, 61, 1385, 27, 78, 272, 95, 123, 7936, 362
Offset: 1

Views

Author

Gus Wiseman, Dec 28 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 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}.

Examples

			Non-isomorphic representatives of the a(n) multisystems for n = 2, 3, 6, 9, 10, 12 (commas and outer brackets elided):
  1  11  {1}{12}  {{1}}{{1}{22}}  {{1}}{{1}{12}}  {{1}}{{1}{23}}
         {2}{11}  {{11}}{{2}{2}}  {{11}}{{1}{2}}  {{11}}{{2}{3}}
                  {{1}}{{2}{12}}  {{1}}{{2}{11}}  {{1}}{{2}{13}}
                  {{12}}{{1}{2}}  {{12}}{{1}{1}}  {{12}}{{1}{3}}
                                  {{2}}{{1}{11}}  {{2}}{{1}{13}}
                                                  {{2}}{{3}{11}}
                                                  {{23}}{{1}{1}}
		

Crossrefs

The non-maximal version is A330666.
The case of constant or strict atoms is A000111.
Labeled versions are A330728, A330665 (prime indices), and A330675 (strongly normal).
Non-isomorphic multiset partitions whose degrees are the prime indices of n are A318285.

Formula

For n > 1, a(2^n) = a(prime(n)) = A000111(n - 1).

A318371 Number of non-isomorphic strict set multipartitions (sets of sets) of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Aug 24 2018

Keywords

Examples

			Non-isomorphic representatives of the a(24) = 6 strict set multipartitions of {1,1,2,3,4}:
  {{1},{1,2,3,4}}
  {{1,2},{1,3,4}}
  {{1},{2},{1,3,4}}
  {{1},{1,2},{3,4}}
  {{2},{1,3},{1,4}}
  {{1},{2},{3},{1,4}}
		

Crossrefs

Formula

a(n) = A318370(A181821(n)).

A318560 Number of combinatory separations of a multiset whose multiplicities are the prime indices of n in weakly decreasing order.

Original entry on oeis.org

1, 1, 2, 2, 3, 4, 5, 3, 8, 7, 7, 8, 11, 12, 15, 5, 15, 17, 22, 14, 27, 19, 30, 13, 27, 30, 33, 26, 42, 37, 56, 7, 44, 45, 51, 34, 77, 67, 72, 25
Offset: 1

Views

Author

Gus Wiseman, Aug 28 2018

Keywords

Comments

A multiset is normal if it spans an initial interval of positive integers. The type of a multiset is the unique normal multiset that has the same sequence of multiplicities when its entries are taken in increasing order. For example the type of 335556 is 112223. A (headless) combinatory separation of a multiset m is a multiset of normal multisets {t_1,...,t_k} such that there exist multisets {s_1,...,s_k} with multiset union m and such that s_i has type t_i for each i = 1...k.
The prime indices of n are the n-th row of A296150.

Examples

			The a(18) = 17 combinatory separations of {1,1,2,2,3}:
  {11223}
  {1,1122} {1,1123} {1,1223} {11,112} {12,112} {12,122} {12,123}
  {1,1,112} {1,1,122} {1,1,123} {1,11,11} {1,11,12} {1,12,12}
  {1,1,1,11} {1,1,1,12}
  {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}]]]]];
    normize[m_]:=m/.Rule@@@Table[{Union[m][[i]],i},{i,Length[Union[m]]}];
    Table[Length[Union[Sort/@Map[normize,mps[nrmptn[n]],{2}]]],{n,20}]

A324326 Number of crossing multiset partitions of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 10, 0, 0, 0, 10, 0, 0, 0, 7, 0, 0, 0, 0, 0, 0, 0, 31, 0, 0, 0, 0, 0, 36, 0, 14, 0, 0, 0, 25, 0, 0, 0, 71, 0, 0, 0, 0, 0, 0, 0, 103, 0, 0, 0, 0, 0, 0, 0, 75
Offset: 1

Views

Author

Gus Wiseman, Feb 22 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}.
A multiset partition is crossing if it contains two blocks of the form {{...x...y...},{...z...t...}} with x < z < y < t or z < x < t < y.

Examples

			The a(36) = 10 crossing multiset partitions of {1,1,2,2,3,4}:
  {{1,3},{1,2,2,4}}
  {{2,4},{1,1,2,3}}
  {{1,1,3},{2,2,4}}
  {{1,2,3},{1,2,4}}
  {{1},{1,3},{2,2,4}}
  {{1},{2,4},{1,2,3}}
  {{2},{1,3},{1,2,4}}
  {{2},{1,1,3},{2,4}}
  {{1,2},{1,3},{2,4}}
  {{1},{2},{1,3},{2,4}}
		

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]]}]];
    croXQ[stn_]:=MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x
    				

Formula

a(n) + A324325(n) = A318284(n).
Showing 1-10 of 13 results. Next