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

A318812 Number of total multiset partitions of the multiset of prime indices of n. Number of total factorizations of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 1, 6, 1, 3, 1, 3, 1, 1, 1, 11, 1, 1, 2, 3, 1, 4, 1, 20, 1, 1, 1, 15, 1, 1, 1, 11, 1, 4, 1, 3, 3, 1, 1, 51, 1, 3, 1, 3, 1, 11, 1, 11, 1, 1, 1, 21, 1, 1, 3, 90, 1, 4, 1, 3, 1, 4, 1, 80, 1, 1, 3, 3, 1, 4, 1, 51, 6, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Sep 04 2018

Keywords

Comments

A total multiset partition of m is either m itself or a total multiset partition of a multiset partition of m that is neither minimal nor maximal.
a(n) depends only on the prime signature of n. - Andrew Howroyd, Dec 30 2019

Examples

			The a(24) = 11 total multiset partitions:
  {1,1,1,2}
  {{1},{1,1,2}}
  {{2},{1,1,1}}
  {{1,1},{1,2}}
  {{1},{1},{1,2}}
  {{1},{2},{1,1}}
  {{{1}},{{1},{1,2}}}
  {{{1}},{{2},{1,1}}}
  {{{2}},{{1},{1,1}}}
  {{{1,2}},{{1},{1}}}
  {{{1,1}},{{1},{2}}}
The a(24) = 11 total factorizations:
  24,
  (2*12), (3*8), (4*6),
  (2*2*6), (2*3*4),
  ((2)*(2*6)), ((6)*(2*2)), ((2)*(3*4)), ((3)*(2*4)), ((4)*(2*3)).
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    totfac[n_]:=1+Sum[totfac[Times@@Prime/@f],{f,Select[facs[n],1
    				
  • PARI
    MultEulerT(u)={my(v=vector(#u)); v[1]=1; for(k=2, #u, forstep(j=#v\k*k, k, -k, my(i=j, e=0); while(i%k==0, i/=k; e++; v[j]+=binomial(e+u[k]-1, e)*v[i]))); v}
    seq(n)={my(v=vector(n, i, isprime(i)), u=vector(n), m=logint(n,2)+1); for(r=1, m, u += v*sum(j=r, m, (-1)^(j-r)*binomial(j-1, r-1)); v=MultEulerT(v)); u[1]=1; u} \\ Andrew Howroyd, Dec 30 2019

Formula

a(product of n distinct primes) = A005121(n).
a(prime^n) = A318813(n).

A318813 Number of balanced reduced multisystems with n atoms all equal to 1.

Original entry on oeis.org

1, 1, 2, 6, 20, 90, 468, 2910, 20644, 165874, 1484344, 14653890, 158136988, 1852077284, 23394406084, 317018563806, 4587391330992, 70598570456104, 1151382852200680, 19835976878704628, 359963038816096924, 6863033015330999110, 137156667020252478684, 2867083618970831936826
Offset: 1

Views

Author

Gus Wiseman, Sep 04 2018

Keywords

Comments

For n > 1, also the number of balanced reduced multisystems whose atoms are an integer partition of n with at least one part > 1. 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. - Gus Wiseman, Dec 31 2019

Examples

			The a(5) = 20 balanced reduced multisystems (with n written in place of 1^n):
  5  (14)  (23)  (113)      (122)      (1112)
                 ((1)(13))  ((1)(22))  ((1)(112))
                 ((3)(11))  ((2)(12))  ((2)(111))
                                       ((11)(12))
                                       ((1)(1)(12))
                                       ((1)(2)(11))
                                       (((1))((1)(12)))
                                       (((1))((2)(11)))
                                       (((2))((1)(11)))
                                       (((12))((1)(1)))
                                       (((11))((1)(2)))
		

Crossrefs

Programs

  • Mathematica
    normize[m_]:=m/.Rule@@@Table[{Union[m][[i]],i},{i,Length[Union[m]]}];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    totfact[n_]:=totfact[n]=1+Sum[totfact[Times@@Prime/@normize[f]],{f,Select[facs[n],1
    				
  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
    seq(n)={my(v=vector(n, i, i==1), u=vector(n)); for(r=1, #v, u += v*sum(j=r, #v, (-1)^(j-r)*binomial(j-1, r-1)); v=EulerT(v)); u} \\ Andrew Howroyd, Dec 30 2019

Formula

a(n > 1) = A330679(n)/2. - Gus Wiseman, Dec 31 2019

Extensions

Terms a(14) and beyond from Andrew Howroyd, Dec 30 2019
Terminology corrected by Gus Wiseman, Dec 31 2019

A292505 Number of complete orderless tree-factorizations of n >= 2.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 1, 5, 1, 3, 1, 3, 1, 1, 1, 9, 1, 1, 2, 3, 1, 4, 1, 12, 1, 1, 1, 12, 1, 1, 1, 9, 1, 4, 1, 3, 3, 1, 1, 29, 1, 3, 1, 3, 1, 9, 1, 9, 1, 1, 1, 17, 1, 1, 3, 33, 1, 4, 1, 3, 1, 4, 1, 44, 1, 1, 3, 3, 1, 4, 1, 29, 5, 1, 1, 17, 1
Offset: 2

Views

Author

Gus Wiseman, Sep 17 2017

Keywords

Comments

An orderless tree-factorization (see A292504 for definition) is complete if all leaves are prime numbers. This sequence first differs from A281119 at a(64)=33.
a(n) depends only on the prime signature of n. - Andrew Howroyd, Nov 18 2018

Examples

			The a(60)=17 complete orderless tree-factorizations are: (2(2(35))), (2(3(25))), (2(5(23))), (2(235)), (3(2(25))), (3(5(22))), (3(225)), (5(2(23))), (5(3(22))), (5(223)), ((22)(35)), ((23)(25)), (22(35)), (23(25)), (25(23)), (35(22)), (2235).
		

Crossrefs

Programs

  • Mathematica
    postfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[postfacs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    oltfacs[n_]:=If[n<=1,{{}},Prepend[Union@@Function[q,Sort/@Tuples[oltfacs/@q]]/@DeleteCases[postfacs[n],{n}],n]];
    Table[Length[Select[oltfacs[n],FreeQ[#,_Integer?(!PrimeQ[#]&)]&]],{n,2,100}]
  • PARI
    seq(n)={my(v=vector(n), w=vector(n)); v[1]=1; for(k=2, n, w[k]=v[k]+isprime(k); forstep(j=n\k*k, k, -k, my(i=j, e=0); while(i%k==0, i/=k; e++; v[j]+=binomial(e+w[k]-1, e)*v[i]))); w[2..n]} \\ Andrew Howroyd, Nov 18 2018

Formula

a(p^n) = A000669(n) for prime p. - Andrew Howroyd, Nov 18 2018

A319122 Number of phylogenetic plane trees on n labels.

Original entry on oeis.org

1, 3, 25, 387, 8521, 241683, 8383705, 343826787, 16273985641, 873119718963, 52360707915385, 3470858539699587, 252000934472119561, 19888355652445884243, 1695252683833578455065, 155208762048402360698787, 15190477481877333732410281, 1582657042668691276257233523
Offset: 1

Views

Author

Gus Wiseman, Sep 11 2018

Keywords

Comments

A phylogenetic plane tree on n labels is either the set of labels itself or a finite sequence of at least two phylogenetic plane trees, one on each block of an ordered set partition of the labels.

Examples

			The a(2) = 3 phylogenetic plane trees are {1,2}, ({1},{2}), ({2},{1}).
		

Crossrefs

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    t[n_]:=t[n]=1+Sum[Times@@t/@f,{f,Join@@Permutations/@Select[sps[Range[n]],Length[#]>1&]}];
    Array[t,8]

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

A277130 Number of planar branching factorizations of n.

Original entry on oeis.org

0, 1, 1, 2, 1, 3, 1, 6, 2, 3, 1, 14, 1, 3, 3, 24, 1, 14, 1, 14, 3, 3, 1, 78, 2, 3, 6, 14, 1, 25, 1, 112, 3, 3, 3, 110, 1, 3, 3, 78, 1, 25, 1, 14, 14, 3, 1, 464, 2, 14, 3, 14, 1, 78, 3, 78, 3, 3, 1, 206, 1, 3, 14, 568, 3, 25, 1, 14, 3, 25, 1, 850, 1, 3, 14, 14
Offset: 1

Views

Author

Michel Marcus, Oct 01 2016

Keywords

Comments

A planar branching factorization of n is either the number n itself or a sequence of at least two planar branching factorizations, one of each factor in an ordered factorization of n. - Gus Wiseman, Sep 11 2018

Examples

			From _Gus Wiseman_, Sep 11 2018: (Start)
The a(12) = 14 planar branching factorizations:
  12,
  (2*6), (3*4), (4*3), (6*2), (2*2*3), (2*3*2), (3*2*2),
  (2*(2*3)), (2*(3*2)), (3*(2*2)), ((2*2)*3), ((2*3)*2), ((3*2)*2).
(End)
		

Crossrefs

Programs

  • C
    #include 
    #include 
    #include 
    #define MAX 10000
    /* Number of planar branching factorizations of n. */
    #define lu unsigned long
    lu nbr[MAX]; /* number of branching */
    lu a, b, d, e; /* temporary variables */
    lu n; lu m, p; // factors of n
    lu x; // square root of n
    void main(unsigned argc, char *argv[])
    {
      memset(nbr, 0, MAX*sizeof(lu));
      for (b=0, n=1; nDaniel Mondot, May 19 2017 */
  • Mathematica
    ordfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#1,d]&)/@ordfacs[n/d],{d,Rest[Divisors[n]]}]]
    otfs[n_]:=Prepend[Join@@Table[Tuples[otfs/@f],{f,Select[ordfacs[n],Length[#]>1&]}],n];
    Table[Length[otfs[n]],{n,20}] (* Gus Wiseman, Sep 11 2018 *)

Formula

a(prime^n) = A118376(n). a(product of n distinct primes) = A319122(n). - Gus Wiseman, Sep 11 2018

Extensions

Terms a(65) and beyond from Daniel Mondot, May 19 2017

A295281 Number of complete strict tree-factorizations of n > 1.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Nov 19 2017

Keywords

Comments

A strict tree-factorization (see A295279 for definition) is complete if its leaves are all prime numbers.
From Andrew Howroyd, Nov 18 2018: (Start)
a(n) depends only on the prime signature of n.
This sequence is very similar but not identical to the number of complete orderless identity tree-factorizations of n. The first difference is at n=900 (square of three primes). Here a(n) = 191 whereas the other sequence would have 197. (End)

Examples

			The a(72) = 6 complete strict tree-factorizations are: 2*3*(2*(2*3)), 2*(2*3*(2*3)), 2*(2*(3*(2*3))), 2*(3*(2*(2*3))), 3*(2*(2*(2*3))), (2*3)*(2*(2*3)).
		

Crossrefs

Programs

  • Mathematica
    postfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[postfacs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    sftc[n_]:=Prepend[Join@@Function[fac,Tuples[sftc/@fac]]/@Select[postfacs[n],And[Length[#]>1,UnsameQ@@#]&],n];
    Table[Length[Select[sftc[n],FreeQ[#,_Integer?(!PrimeQ[#]&)]&]],{n,2,100}]
  • PARI
    seq(n)={my(v=vector(n), w=vector(n)); v[1]=1; for(k=2, n, w[k]=v[k]+isprime(k); forstep(j=n\k*k, k, -k, v[j]+=w[k]*v[j/k])); w[2..n]} \\ Andrew Howroyd, Nov 18 2018

Formula

a(product of n distinct primes) = A000311(n).
Positions of zeros are proper prime powers A025475. Positions of nonzero entries are A085971.

A318848 Number of complete tree-partitions of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 1, 1, 2, 3, 5, 4, 12, 9, 12, 17, 34, 29, 44, 26, 92, 90, 277, 68, 171, 93, 806, 144, 197, 309, 581, 269, 2500, 428, 7578, 236, 631, 1025, 869, 954, 24198, 3463, 2402, 712, 75370, 1957, 243800, 1040, 3200, 11705, 776494, 1612, 4349, 2358, 8862, 3993, 2545777
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 sequence of tree-partitions, one of each part of a multiset partition of m with at least two parts. A tree-partition is complete if the leaves are all multisets of length 1.

Examples

			The a(12) = 17 complete tree-partitions of {1,1,2,3} with the leaves (x) replaced with just x:
  (1(1(23)))
  (1(2(13)))
  (1(3(12)))
  (2(1(13)))
  (2(3(11)))
  (3(1(12)))
  (3(2(11)))
  ((11)(23))
  ((12)(13))
  (1(123))
  (2(113))
  (3(112))
  (11(23))
  (12(13))
  (13(12))
  (23(11))
  (1123)
		

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}]]]]];
    allmsptrees[m_]:=Prepend[Join@@Table[Tuples[allmsptrees/@p],{p,Select[mps[m],Length[#]>1&]}],m];
    Table[Length[Select[allmsptrees[nrmptn[n]],FreeQ[#,{?AtomQ,_}]&]],{n,20}]

Formula

a(n) = A281119(A181821(n)).
a(prime(n)) = A196545(n)
a(2^n) = A000311(n).

Extensions

More terms from Jinyuan Wang, Jun 26 2020

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

Original entry on oeis.org

1, 1, 2, 2, 4, 6, 12, 8, 28, 20, 32, 38, 112, 76, 116, 58, 352, 236, 1296, 176, 540, 288, 4448, 374, 612, 1144, 1812, 824, 16640, 1316, 59968, 612, 2336, 4528, 3208, 2924, 231168, 18320, 10632, 2168, 856960, 7132, 3334400, 3776, 11684, 74080, 12679424, 4919, 19192
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 sequence of tree-partitions, one of each part of a multiset partition of m with at least two parts.

Examples

			The a(6) = 6 tree-partitions of {1,1,2}:
  (112)
  ((1)(12))
  ((2)(11))
  ((1)(1)(2))
  ((1)((1)(2)))
  ((2)((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}]]]]];
    allmsptrees[m_]:=Prepend[Join@@Table[Tuples[allmsptrees/@p],{p,Select[mps[m],Length[#]>1&]}],m];
    Table[Length[allmsptrees[nrmptn[n]]],{n,20}]

Formula

a(n) = A281118(A181821(n)).
a(prime(n)) = A289501(n).
a(2^n) = A005804(n).

Extensions

More terms from Jinyuan Wang, Jun 26 2020

A319118 Number of multimin tree-factorizations of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 6, 2, 2, 1, 8, 1, 2, 2, 24, 1, 6, 1, 8, 2, 2, 1, 42, 2, 2, 6, 8, 1, 8, 1, 112, 2, 2, 2, 38, 1, 2, 2, 42, 1, 8, 1, 8, 8, 2, 1, 244, 2, 6, 2, 8, 1, 24, 2, 42, 2, 2, 1, 58, 1, 2, 8, 568, 2, 8, 1, 8, 2, 8, 1, 268, 1, 2, 6, 8, 2, 8, 1, 244, 24
Offset: 1

Views

Author

Gus Wiseman, Sep 10 2018

Keywords

Comments

A multimin factorization of n is an ordered factorization of n into factors greater than 1 such that the sequence of minimal primes dividing each factor is weakly increasing. A multimin tree-factorization of n is either the number n itself or a sequence of multimin tree-factorizations, one of each factor in a multimin factorization of n with at least two factors.

Examples

			The a(12) = 8 multimin tree-factorizations:
  12,
  (2*6), (4*3), (6*2), (2*2*3),
  (2*(2*3)), ((2*2)*3), ((2*3)*2).
Or as series-reduced plane trees of multisets:
  112,
  (1,12), (11,2), (12,1), (1,1,2),
  (1,(1,2)), ((1,1),2), ((1,2),1).
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#1,d]&)/@Select[facs[n/d],Min@@#1>=d&],{d,Rest[Divisors[n]]}]];
    mmftrees[n_]:=Prepend[Join@@(Tuples[mmftrees/@#]&/@Select[Join@@Permutations/@Select[facs[n],Length[#]>1&],OrderedQ[FactorInteger[#][[1,1]]&/@#]&]),n];
    Table[Length[mmftrees[n]],{n,100}]

Formula

a(prime^n) = A118376(n).
a(product of n distinct primes) = A005804(n).
Showing 1-10 of 15 results. Next