[go: up one dir, main page]

login
Revision History for A318153 (Bold, blue-underlined text is an addition; faded, red-underlined text is a deletion.)

Showing all changes.
Number of antichain covers of the free pure symmetric multifunction (with empty expressions allowed) with e-number n.
(history; published version)
#9 by Alois P. Heinz at Tue Sep 11 21:50:48 EDT 2018
STATUS

proposed

approved

#8 by Gus Wiseman at Tue Sep 11 18:35:29 EDT 2018
STATUS

editing

proposed

#7 by Gus Wiseman at Tue Sep 11 02:19:26 EDT 2018
NAME

Number of antichain covers of the orderless expression free pure symmetric multifunction (with empty expressions allowed) with e-number n.

COMMENTS

If n = 1 let e(n) be the leaf symbol "o". Given a positive integer n > 1 we construct a unique orderless expression free pure symmetric multifunction (with empty expressions allowed) e(n) (as can be represented in functional programming languages such as Mathematica) 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). The a(n) is the number of ways to partition e(n) into disjoint subexpressions such that all leaves are covered by exactly one of them.

STATUS

approved

editing

Discussion
Tue Sep 11
18:35
Gus Wiseman: Removed "Mathematica" from the name.
#6 by Susanna Cuyler at Mon Aug 20 07:55:05 EDT 2018
STATUS

proposed

approved

#5 by Gus Wiseman at Mon Aug 20 03:17:48 EDT 2018
STATUS

editing

proposed

#4 by Gus Wiseman at Mon Aug 20 03:17:33 EDT 2018
NAME

Number of antichain covers of the orderless Mathematica expression with e-number n.

COMMENTS

If n = 1 let e(n) be the leaf symbol "o". Given a positive integer n > 1 we construct a unique orderless Mathematica expression e(n) (as can be represented in functional programming languages such as Mathematica) 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). The a(n) is the number of ways to partition e(n) into disjoint subexpressions such that all leaves are covered by exactly one of them.

#3 by Gus Wiseman at Sun Aug 19 10:12:51 EDT 2018
#2 by Gus Wiseman at Sun Aug 19 10:05:56 EDT 2018
NAME

allocated for Gus WisemanNumber of antichain covers of the orderless Mathematica expression with e-number n.

DATA

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, 5, 4, 7, 7, 8, 4, 5, 10, 6, 3, 5, 8, 8, 9, 5, 6, 11, 7, 4, 6, 9, 9, 5, 10, 6, 7, 12, 8, 5, 7, 10, 10, 6, 11, 7, 8, 13, 3, 9, 6, 8, 11, 11, 7, 12, 8, 9, 14, 4, 10, 7, 9, 12, 12, 3, 8

OFFSET

1,2

COMMENTS

If n = 1 let e(n) be the leaf symbol "o". Given a positive integer n > 1 we construct a unique orderless Mathematica expression 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). The a(n) is the number of ways to partition e(n) into disjoint subexpressions such that all leaves are covered by exactly one of them.

FORMULA

If n = rad(x)^(Product_i prime(y_i)^z_i) where rad = A007916 then a(n) = 1 + a(x) * Product_i a(y_i)^z_i.

EXAMPLE

441 is the e-number of o[o,o][o] which has antichain covers {o[o,o][o]}, {o[o,o], o}, {o, o, o, o}}, corresponding to the leaf-colorings 1[1,1][1], 1[1,1][2], 1[2,3][4], so a(441) = 3.

MATHEMATICA

nn=20000;

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]]}, 1+a[radPi[n^(1/g)]]*Product[a[PrimePi[pr[[1]]]]^pr[[2]], {pr, If[g==1, {}, FactorInteger[g]]}]]];

Array[a, 100]

KEYWORD

allocated

nonn

AUTHOR

Gus Wiseman, Aug 19 2018

STATUS

approved

editing

#1 by Gus Wiseman at Sun Aug 19 10:05:56 EDT 2018
NAME

allocated for Gus Wiseman

KEYWORD

allocated

STATUS

approved