[go: up one dir, main page]

login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
Search: a277996 -id:a277996
Displaying 1-10 of 32 results found. page 1 2 3 4
     Sort: relevance | references | number | modified | created      Format: long | short | data
A317884 Number of series-reduced achiral free pure multifunctions (with empty expressions allowed) with one atom and n positions. +0
6
1, 1, 1, 2, 4, 8, 14, 26, 47, 87, 160, 295, 540, 997, 1832, 3369, 6197, 11406, 20975, 38594, 70991, 130610, 240275, 442043, 813184, 1496053, 2752251, 5063319, 9314879, 17136632, 31526032, 57998423, 106699160, 196294065, 361120800, 664352454, 1222204958 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,4
COMMENTS
A series-reduced achiral expression (SRAE) is either (case 1) the leaf symbol "o", or (case 2) a possibly empty but not unitary expression of the form h[g, ..., g], where h and g are SRAEs. The number of positions in an SRAE is the number of brackets [...] plus the number of o's.
Also the number of series-reduced achiral Mathematica expressions with one atom and n positions.
LINKS
FORMULA
a(1) = 1; a(n > 1) = a(n-1) + Sum_{0 < k < n-1} a(k) * Sum_{d|(n-k-1), d < n-k-1} a(d).
EXAMPLE
The a(6) = 8 SRAEs:
o[o,o,o,o]
o[o[],o[]]
o[][o,o,o]
o[][][o,o]
o[o,o,o][]
o[][o,o][]
o[o,o][][]
o[][][][][]
MAPLE
a:= proc(n) option remember; `if`(n=1, 1, a(n-1)+add(a(j)*add(
a(d), d=numtheory[divisors](n-j-1) minus {n-j-1}), j=1..n-1))
end:
seq(a(n), n=1..45); # Alois P. Heinz, Sep 05 2018
MATHEMATICA
allAchExprSR[n_] := If[n == 1, {"o"}, Join @@ Cases[Table[PR[k, n - k - 1], {k, n - 1}], PR[h_, g_] :> Join @@ Table[Apply @@@ Tuples[{allAchExprSR[h], Select[Tuples[allAchExprSR /@ p], SameQ @@ # &]}], {p, If[g == 0, {{}}, Join @@ Permutations /@ Rest[IntegerPartitions[g]]]}]]]; Table[Length[allAchExprSR[n]], {n, 12}]
(* Second program: *)
a[n_] := a[n] = If[n == 1, 1, a[n-1] + Sum[a[j]*DivisorSum[
n-j-1, If[# < n-j-1, a[#], 0]&], {j, 1, n-2}]];
Array[a, 45] (* Jean-François Alcover, May 17 2021, after Alois P. Heinz *)
PROG
(PARI) seq(n)={my(p=O(x)); for(n=1, n, p = x + p*x*(1 + sum(k=2, n-2, subst(p + O(x^(n\k+1)), x, x^k)) ) + O(x*x^n)); Vec(p)} \\ Andrew Howroyd, Aug 19 2018
(PARI) seq(n)={my(v=vector(n)); v[1]=1; for(n=2, #v, v[n]=v[n-1] + sum(i=1, n-2, v[i]*sumdiv(n-i-1, d, if(d<n-i-1, v[d], 0)))); v} \\ Andrew Howroyd, Aug 19 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Aug 09 2018
EXTENSIONS
Terms a(13) and beyond from Andrew Howroyd, Aug 19 2018
STATUS
approved
A317879 Number of free pure identity multifunctions (with empty expressions allowed) with one atom and n positions. +0
8
1, 1, 2, 4, 11, 29, 83, 251, 767, 2403, 7652, 24758, 80875, 266803, 887330, 2972108, 10016981, 33942461, 115572864, 395226810, 1356840007, 4674552089, 16156355357, 56003840659, 194651585875, 678220460687, 2368505647624, 8288873657180, 29064904732911 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,3
COMMENTS
A free pure identity multifunction (with empty expressions allowed) (IME) is either (case 1) the leaf symbol "o", or (case 2) a possibly empty expression of the form h[g_1, ..., g_k] where h is an IME, each of the g_i for i = 1, ..., k >= 0 is an IME, and for i != j we have g_i != g_j. The number of positions in an IME is the number of brackets [...] plus the number of o's.
Also the number of identity Mathematica expressions with one atom and n positions.
LINKS
EXAMPLE
The a(5) = 11 IMEs:
o[o[o]]
o[o][o]
o[o[][]]
o[o[],o]
o[o,o[]]
o[][o[]]
o[][][o]
o[o[]][]
o[][o][]
o[o][][]
o[][][][]
MATHEMATICA
allIdExpr[n_]:=If[n==1, {"o"}, Join@@Cases[Table[PR[k, n-k-1], {k, n-1}], PR[h_, g_]:>Join@@Table[Apply@@@Tuples[{allIdExpr[h], Select[Tuples[allIdExpr/@p], UnsameQ@@#&]}], {p, Join@@Permutations/@IntegerPartitions[g]}]]];
Table[Length[allIdExpr[n]], {n, 12}]
PROG
(PARI) seq(n)={my(v=vector(n)); v[1]=1; for(n=2, n, my(p=prod(k=1, n, 1 + sum(i=1, n\k, binomial(v[k], i)*x^(i*k)*y^i) + O(x*x^n))); v[n]=v[n-1]+sum(k=1, n-2, v[n-k-1]*subst(serlaplace(y^0*polcoef(p, k)), y, 1))); v} \\ Andrew Howroyd, Sep 01 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Aug 09 2018
EXTENSIONS
Terms a(13) and beyond from Andrew Howroyd, Sep 01 2018
STATUS
approved
A317056 Depth of the free pure symmetric multifunction (with empty expressions allowed) with e-number n. +0
9
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 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,3
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).
LINKS
EXAMPLE
e(21025) = o[o[o]][o] has depth 3 so a(21025) = 3.
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}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Aug 18 2018
STATUS
approved
A318153 Number of antichain covers of the free pure symmetric multifunction (with empty expressions allowed) with e-number n. +0
4
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 (list; graph; refs; listen; history; text; internal format)
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 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.
LINKS
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]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Aug 19 2018
STATUS
approved
A318152 e-numbers of unlabeled rooted trees. A number n is in the sequence iff n = 2^(prime(y_1) * ... * prime(y_k)) for some k > 0 and y_1, ..., y_k already in the sequence. +0
4
1, 4, 16, 128, 256, 16384, 65536, 268435456, 4294967296, 562949953421312, 9007199254740992, 72057594037927936, 18446744073709551616, 316912650057057350374175801344, 81129638414606681695789005144064, 5192296858534827628530496329220096 (list; graph; refs; listen; history; text; internal format)
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 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 sequence consists of all numbers n such that e(n) contains no empty subexpressions f[] or subexpressions in heads f[x_1, ..., x_k][y_1, ..., y_k] where k,j >= 0.
LINKS
EXAMPLE
The sequence contains 16384 = 2^14 = 2^(prime(1) * prime(4)) because 1 and 4 both already belong to the sequence.
The sequence of unlabeled rooted trees with e-numbers in the sequence begins:
1: o
4: (o)
16: (oo)
128: ((o))
256: (ooo)
16384: (o(o))
65536: (oooo)
. (oo(o))
. (ooooo)
. ((o)(o))
((oo))
(ooo(o))
(oooooo)
(o(o)(o))
(o(oo))
(oooo(o))
(ooooooo)
(oo(o)(o))
MATHEMATICA
baQ[n_]:=Or[n==1, MatchQ[FactorInteger[n], {{2, _?(And@@Cases[FactorInteger[#], {p_, k_}:>baQ[PrimePi[p]]]&)}}]];
Select[2^Range[0, 50], baQ]
CROSSREFS
A subsequence of A000079 and A318151.
KEYWORD
nonn
AUTHOR
Gus Wiseman, Aug 19 2018
STATUS
approved
A318151 e-numbers of unlabeled rooted trees with empty leaves o[] allowed. A number n is in the sequence iff n = 2^(prime(y_1) * ... * prime(y_k)) for some k >= 0 and y_1, ..., y_k already in the sequence. +0
1
1, 2, 4, 8, 16, 64, 128, 256, 512, 4096, 16384, 65536, 262144, 524288, 2097152, 16777216, 134217728, 268435456, 4294967296, 68719476736, 274877906944, 4398046511104, 281474976710656, 562949953421312, 9007199254740992, 18014398509481984, 72057594037927936 (list; graph; refs; listen; history; text; internal format)
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 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 sequence consists of all numbers n such that e(n) contains no subexpressions in heads f[x_1, ..., x_k][y_1, ..., y_k] where k,j >= 0.
LINKS
CROSSREFS
A subsequence of A000079.
KEYWORD
nonn
AUTHOR
Gus Wiseman, Aug 19 2018
STATUS
approved
A318150 e-numbers of free pure functions with one atom. +0
5
1, 4, 36, 128, 2025, 21025, 279936, 4338889, 449482401, 78701569444, 373669453125, 18845583322500, 1347646586640625, 202054211912421649, 6193981883008128893161, 139629322539586311507076, 170147232533595290155627, 355156175404848064835984400 (list; graph; refs; listen; history; text; internal format)
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 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). This sequence consists of all numbers n such that e(n) contains no non-unitary subexpressions f[x_1, ..., x_k] where k != 1.
LINKS
FORMULA
a(1) = 1, and if a and b are in this sequence then so is rad(a)^prime(b). - Charlie Neder, Feb 23 2019
EXAMPLE
The sequence of all free pure functions with one atom together with their e-numbers begins:
1: o
4: o[o]
36: o[o][o]
128: o[o[o]]
2025: o[o][o][o]
21025: o[o[o]][o]
279936: o[o][o[o]]
4338889: o[o][o][o][o]
CROSSREFS
A subsequence of A001597.
KEYWORD
nonn
AUTHOR
Gus Wiseman, Aug 19 2018
EXTENSIONS
More terms from Charlie Neder, Feb 23 2019
STATUS
approved
A318149 e-numbers of free pure symmetric multifunctions with one atom. +0
5
1, 4, 16, 36, 128, 256, 441, 1296, 2025, 16384, 21025, 65536, 77841, 194481, 220900, 279936, 1679616, 1803649, 4100625, 4338889, 268435456, 273571600, 442050625, 449482401, 1801088541, 4294967296, 4334247225, 6059221281 (list; graph; refs; listen; history; text; internal format)
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 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 sequence consists of all numbers n such that e(n) contains no empty subexpressions f[].
LINKS
EXAMPLE
The sequence of free pure symmetric multifunctions with one atom "o", together with their e-numbers begins:
1: o
4: o[o]
16: o[o,o]
36: o[o][o]
128: o[o[o]]
256: o[o,o,o]
441: o[o,o][o]
1296: o[o][o,o]
2025: o[o][o][o]
16384: o[o,o[o]]
21025: o[o[o]][o]
65536: o[o,o,o,o]
77841: o[o,o,o][o]
194481: o[o,o][o,o]
220900: o[o,o][o][o]
279936: o[o][o[o]]
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]]]]]];
Select[Range[nn], FreeQ[exp[#], _[]]&]
PROG
(Python) See Neder link.
CROSSREFS
A subsequence of A001597.
KEYWORD
nonn
AUTHOR
Gus Wiseman, Aug 19 2018
EXTENSIONS
a(16)-a(27) from Charlie Neder, Sep 01 2018
STATUS
approved
A317994 Number of inequivalent leaf-colorings of the free pure symmetric multifunction with e-number n. +0
9
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 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,4
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).
LINKS
EXAMPLE
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
KEYWORD
nonn
AUTHOR
Gus Wiseman, Aug 18 2018
STATUS
approved
A317765 Number of distinct subexpressions of the free pure symmetric multifunction (with empty expressions allowed) with e-number n. +0
8
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 (list; graph; refs; listen; history; text; internal format)
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 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).
LINKS
EXAMPLE
The a(12) = 4 subexpressions of o[o[]][] are {o, o[], o[o[]], o[o[]][]}.
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}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Aug 18 2018
STATUS
approved
page 1 2 3 4

Search completed in 0.020 seconds

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 30 00:57 EDT 2024. Contains 375520 sequences. (Running on oeis4.)