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-6 of 6 results.

A316112 Number of leaves in the free pure symmetric multifunction (with empty expressions allowed) with e-number n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Aug 18 2018

Keywords

Comments

If n = 1 let e(n) be the leaf symbol "o". Given a positive integer n > 1 we construct a unique free pure symmetric multifunction e(n) with one atom by expressing n as a power of a number that is not a perfect power to a product of prime numbers: n = rad(x)^(prime(y_1) * ... * prime(y_k)) where rad = A007916. Then e(n) = e(x)[e(y_1), ..., e(y_k)]. For example, e(21025) = o[o[o]][o] because 21025 = rad(rad(1)^prime(rad(1)^prime(1)))^prime(1).

Examples

			e(21025) = o[o[o]][o] has 4 leaves so a(21025) = 4.

Crossrefs

Cf. A007916, A052409, A052410, A109129, A277576, A277996, A300626, A316112, A317056, A317658, A317765, A317994.

Programs

  • Mathematica
    nn=1000;
radQ[n_]:=If[n==1,False,GCD@@FactorInteger[n][[All,2]]==1];
rad[n_]:=rad[n]=If[n==0,1,NestWhile[#+1&,rad[n-1]+1,Not[radQ[#]]&]];
Clear[radPi];Set@@@Array[radPi[rad[#]]==#&,nn];
a[n_]:=If[n==1,1,With[{g=GCD@@FactorInteger[n][[All,2]]},a[radPi[Power[n,1/g]]]+Sum[a[PrimePi[pr[[1]]]]*pr[[2]],{pr,If[g==1,{},FactorInteger[g]]}]]];
Table[a[n],{n,100}]

Formula

a(rad(x)^(prime(y_1) * ... * prime(y_k))) = a(x) + a(y_1) + ... + a(y_k) where rad = A007916.

A317056 Depth of the free pure symmetric multifunction (with empty expressions allowed) with e-number n.

Original entry on oeis.org

0, 1, 2, 1, 3, 2, 4, 2, 2, 3, 5, 3, 3, 4, 6, 1, 4, 4, 5, 7, 2, 5, 5, 6, 3, 8, 2, 3, 6, 6, 7, 3, 4, 9, 3, 2, 4, 7, 7, 8, 4, 5, 10, 4, 3, 5, 8, 8, 4, 9, 5, 6, 11, 5, 4, 6, 9, 9, 5, 10, 6, 7, 12, 2, 6, 5, 7, 10, 10, 6, 11, 7, 8, 13, 3, 7, 6, 8, 11, 11, 2, 7, 12
Offset: 1

Views

Author

Gus Wiseman, Aug 18 2018

Keywords

Comments

If n = 1 let e(n) be the leaf symbol "o". Given a positive integer n > 1 we construct a unique free pure symmetric multifunction e(n) with one atom by expressing n as a power of a number that is not a perfect power to a product of prime numbers: n = rad(x)^(prime(y_1) * ... * prime(y_k)) where rad = A007916. Then e(n) = e(x)[e(y_1), ..., e(y_k)]. For example, e(21025) = o[o[o]][o] because 21025 = rad(rad(1)^prime(rad(1)^prime(1)))^prime(1).

Examples

			e(21025) = o[o[o]][o] has depth 3 so a(21025) = 3.

Crossrefs

Cf. A007916, A052409, A052410, A109082, A277576, A277996, A300626, A316112, A317056, A317658, A317765, A317994.

Programs

  • Mathematica
    nn=1000;
radQ[n_]:=If[n===1,False,GCD@@FactorInteger[n][[All,2]]===1];
rad[n_]:=rad[n]=If[n===0,1,NestWhile[#+1&,rad[n-1]+1,Not[radQ[#]]&]];
Clear[radPi];Set@@@Array[radPi[rad[#]]==#&,nn];
exp[n_]:=If[n===1,"o",With[{g=GCD@@FactorInteger[n][[All,2]]},Apply[exp[radPi[Power[n,1/g]]],exp/@Flatten[Cases[FactorInteger[g],{p_?PrimeQ,k_}:>ConstantArray[PrimePi[p],k]]]]]];
Table[Max@@Length/@Position[exp[n],_],{n,200}]

A317994 Number of inequivalent leaf-colorings of the free pure symmetric multifunction with e-number n.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 1, 4, 2, 2, 2, 1, 4, 2, 2, 2, 2, 1, 2, 4, 2, 2, 2, 2, 2, 1, 2, 5, 4, 2, 2, 2, 2, 2, 1, 2, 5, 4, 2, 2, 2, 2, 2, 2, 1, 2, 5, 4, 2, 2, 2, 2, 2, 2, 1, 5, 2, 5, 4, 2, 2, 2, 2, 2, 2, 1, 5, 2, 5, 4, 2, 2, 4, 2, 2, 2, 2, 1, 5
Offset: 1

Views

Author

Gus Wiseman, Aug 18 2018

Keywords

Comments

If n = 1 let e(n) be the leaf symbol "o". Given a positive integer n > 1 we construct a unique free pure symmetric multifunction (with empty expressions allowed) e(n) with one atom by expressing n as a power of a number that is not a perfect power to a product of prime numbers: n = rad(x)^(prime(y_1) * ... * prime(y_k)) where rad = A007916. Then e(n) = e(x)[e(y_1), ..., e(y_k)]. For example, e(21025) = o[o[o]][o] because 21025 = rad(rad(1)^prime(rad(1)^prime(1)))^prime(1).

Examples

			Inequivalent representatives of the a(441) = 11 colorings of the expression e(441) = o[o,o][o] are the following.
  1[1,1][1]
  1[1,1][2]
  1[1,2][1]
  1[1,2][2]
  1[1,2][3]
  1[2,2][1]
  1[2,2][2]
  1[2,2][3]
  1[2,3][1]
  1[2,3][2]
  1[2,3][4]

Crossrefs

Cf. A007916, A052409, A052410, A277576, A277996, A300626, A316112, A317056, A317658, A317765.

A317765 Number of distinct subexpressions of the free pure symmetric multifunction (with empty expressions allowed) with e-number n.

Original entry on oeis.org

1, 2, 3, 2, 4, 3, 5, 3, 3, 4, 6, 4, 4, 5, 7, 2, 5, 5, 6, 8, 3, 6, 6, 7, 4, 9, 3, 4, 7, 7, 8, 4, 5, 10, 4, 3, 5, 8, 8, 9, 5, 6, 11, 5, 4, 6, 9, 9, 5, 10, 6, 7, 12, 6, 5, 7, 10, 10, 6, 11, 7, 8, 13, 3, 7, 6, 8, 11, 11, 7, 12, 8, 9, 14, 4, 8, 7, 9, 12, 12, 3, 8
Offset: 1

Views

Author

Gus Wiseman, Aug 18 2018

Keywords

Comments

If n = 1 let e(n) be the leaf symbol "o". Given a positive integer n > 1 we construct a unique free pure symmetric multifunction (with empty expressions allowed) e(n) with one atom by expressing n as a power of a number that is not a perfect power to a product of prime numbers: n = rad(x)^(prime(y_1) * ... * prime(y_k)) where rad = A007916. Then e(n) = e(x)[e(y_1), ..., e(y_k)]. For example, e(21025) = o[o[o]][o] because 21025 = rad(rad(1)^prime(rad(1)^prime(1)))^prime(1).

Examples

			The a(12) = 4 subexpressions of o[o[]][] are {o, o[], o[o[]], o[o[]][]}.

Crossrefs

Cf. A007916, A052409, A052410, A277576, A277996, A300626, A316112, A317056, A317658, A317713, A317994.

Programs

  • Mathematica
    nn=1000;
radQ[n_]:=If[n===1,False,GCD@@FactorInteger[n][[All,2]]===1];
rad[n_]:=rad[n]=If[n===0,1,NestWhile[#+1&,rad[n-1]+1,Not[radQ[#]]&]];
Clear[radPi];Set@@@Array[radPi[rad[#]]==#&,nn];
exp[n_]:=If[n===1,"o",With[{g=GCD@@FactorInteger[n][[All,2]]},Apply[exp[radPi[Power[n,1/g]]],exp/@Flatten[Cases[FactorInteger[g],{p_?PrimeQ,k_}:>ConstantArray[PrimePi[p],k]]]]]];
Table[Length[Union[Cases[exp[n],_,{0,Infinity},Heads->True]]],{n,100}]

A304486 Number of inequivalent leaf-colorings of the unlabeled rooted tree with Matula-Goebel number n.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 2, 3, 2, 2, 1, 4, 2, 4, 2, 5, 2, 4, 3, 4, 4, 2, 2, 7, 2, 5, 3, 9, 2, 5, 1, 7, 2, 4, 4, 9, 4, 7, 5, 7, 2, 11, 4, 4, 4, 4, 2, 12, 7, 4, 4, 11, 5, 7, 2, 16, 7, 5, 2, 11, 4, 2, 9, 11, 5, 5, 3, 9, 4, 11
Offset: 1

Views

Author

Gus Wiseman, Aug 17 2018

Keywords

Examples

			Inequivalent representatives of the a(52) = 11 colorings of the tree (oo(o(o))) are the following.
  (11(1(1)))
  (11(1(2)))
  (11(2(1)))
  (11(2(2)))
  (11(2(3)))
  (12(1(1)))
  (12(1(2)))
  (12(1(3)))
  (12(3(1)))
  (12(3(3)))
  (12(3(4)))

Crossrefs

Cf. A000081, A001678, A004111, A007097, A049076, A061773, A061775, A109082, A109129, A300626, A303024, A303431, A317713, A317994.

A304485 Regular triangle where T(n,k) is the number of inequivalent colorings of free pure symmetric multifunctions (with empty expressions allowed) with n positions and k leaves.

Original entry on oeis.org

1, 1, 0, 1, 2, 0, 1, 6, 4, 0, 1, 12, 23, 7, 0, 1, 20, 81, 73, 12, 0, 1, 30, 209, 407, 206, 19, 0, 1, 42, 451, 1566, 1751, 534, 30, 0, 1, 56, 858, 4711, 9593, 6695, 1299, 45, 0, 1, 72, 1494, 11951, 39255, 51111, 23530, 3004, 67, 0, 1, 90, 2430, 26752, 130220, 278570, 245319, 77205, 6664, 97, 0
Offset: 1

Views

Author

Gus Wiseman, Aug 17 2018

Keywords

Comments

A free pure symmetric multifunction (with empty expressions allowed) f in EOME is either (case 1) a positive integer, or (case 2) a possibly empty expression of the form h[g_1, ..., g_k] where k >= 0, h is in EOME, each of the g_i for i = 1, ..., k is in EOME, and for i < j we have g_i <= g_j under a canonical total ordering of EOME, such as the Mathematica ordering of expressions.
T(n,k) is also the number of inequivalent colorings of orderless Mathematica expressions with n positions and k leaves.

Examples

			Inequivalent representatives of the T(5,3) = 23 Mathematica expressions:
  1[][1,1]  1[1,1][]  1[1][1]  1[1[1]]  1[1,1[]]
  1[][1,2]  1[1,2][]  1[1][2]  1[1[2]]  1[1,2[]]
  1[][2,2]  1[2,2][]  1[2][1]  1[2[1]]  1[2,1[]]
  1[][2,3]  1[2,3][]  1[2][2]  1[2[2]]  1[2,2[]]
                      1[2][3]  1[2[3]]  1[2,3[]]
Triangle begins:
    1
    1    0
    1    2    0
    1    6    4    0
    1   12   23    7    0
    1   20   81   73   12    0
    1   30  209  407  206   19    0
    1   42  451 1566 1751  534   30    0

Links

Crossrefs

Row sums are A300626.
Cf. A000612, A007716, A052893, A053492, A277996, A280000, A317676.

Programs

  • PARI
    \\ See links in A339645 for combinatorial species functions.
cycleIndexSeries(n)={my(p=O(x)); for(n=1, n, p = x*sv(1) + x*p*sExp(p)); p}
T(n)={my(v=Vec(InequivalentColoringsSeq(sFuncSubst(cycleIndexSeries(n), i->sv(i)*y^i)))); vector(n, n, Vecrev(v[n]/y, n))}
{ my(A=T(10)); for(n=1, #A, print(A[n])) } \\ Andrew Howroyd, Jan 01 2021

Extensions

Terms a(37) and beyond from Andrew Howroyd, Jan 01 2021
Showing 1-6 of 6 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.