A317653 -id:A317653 - cp's OEIS frontend

A317658 Number of positions in the n-th free pure symmetric multifunction (with empty expressions allowed) with one atom.

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

Offset: 1

Gus Wiseman, Aug 03 2018

nonn

Given a positive integer n > 1 we construct a unique free pure symmetric multifunction e(n) 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)].

Also the number of positions in the orderless Mathematica expression with e-number n.

			The first twenty Mathematica expressions:
   1: o
   2: o[]
   3: o[][]
   4: o[o]
   5: o[][][]
   6: o[o][]
   7: o[][][][]
   8: o[o[]]
   9: o[][o]
  10: o[o][][]
  11: o[][][][][]
  12: o[o[]][]
  13: o[][o][]
  14: o[o][][][]
  15: o[][][][][][]
  16: o[o,o]
  17: o[o[]][][]
  18: o[][o][][]
  19: o[o][][][][]
  20: o[][][][][][][]

Mathematica Reference, Orderless

First differs from A277615 at a(128) = 5, A277615(128) = 6.
Cf. A001003, A052893, A053492, A255906, A277996, A279944, A280000.
Cf. A317652, A317653, A317654, A317655, A317656.

nn=100;
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,x,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[exp[n],{n,1,nn}]

a(rad(x)^(prime(y_1) * ... * prime(y_k))) = a(x) + a(y_1) + ... + a(y_k).
e(2^(2^n)) = o[o,...,o].
e(2^prime(2^prime(2^...))) = o[o[...o[o]]].
e(rad(rad(rad(...)^2)^2)^2) = o[o][o]...[o].

A317776 Number of strict multiset partitions of normal multisets of size n, where a multiset is normal if it spans an initial interval of positive integers.

Original entry on oeis.org

1, 1, 3, 13, 59, 313, 1847, 11977, 84483, 642405, 5228987, 45297249, 415582335, 4021374193, 40895428051, 435721370413, 4850551866619, 56282199807401, 679220819360775, 8508809310177481, 110454586096508563, 1483423600240661781, 20581786429087269819

Offset: 0

Gus Wiseman, Aug 06 2018

nonn

			The a(3) = 13 strict multiset partitions:
  {{1,1,1}}, {{1},{1,1}},
  {{1,2,2}}, {{1},{2,2}}, {{2},{1,2}},
  {{1,1,2}}, {{1},{1,2}}, {{2},{1,1}},
  {{1,2,3}}, {{1},{2,3}}, {{2},{1,3}}, {{3},{1,2}}, {{1},{2},{3}}.

Alois P. Heinz, Table of n, a(n) for n = 0..300

Cf. A001055, A007716, A045778, A255906, A281116, A317449, A317532, A317583, A317653, A317752, A317757, A317775.

Row sums of A327116.

C:= binomial:
b:= proc(n, i, k) option remember; `if`(n=0, 1, `if`(i<1, 0, add(
      b(n-i*j, min(n-i*j, i-1), k)*C(C(k+i-1, i), j), j=0..n/i)))
    end:
a:= n-> add(add(b(n$2, i)*(-1)^(k-i)*C(k, i), i=0..k), k=0..n):
seq(a(n), n=0..23);  # Alois P. Heinz, Sep 16 2019

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]]]];
allnorm[n_Integer]:=Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1];
Table[Length[Select[Join@@mps/@allnorm[n],UnsameQ@@#&]],{n,9}]
(* Second program: *)
c := Binomial;
b[n_, i_, k_] := b[n, i, k] = If[n==0, 1, If[i<1, 0, Sum[b[n - i*j, Min[n - i*j, i-1], k] c[c[k+i-1, i], j], {j, 0, n/i}]]];
a[n_] := Sum[b[n, n, i] (-1)^(k-i) c[k, i], {k, 0, n}, {i, 0, k}];
a /@ Range[0, 23] (* Jean-François Alcover, Dec 17 2020, after Alois P. Heinz *)

a(0), a(8)-a(22) from Alois P. Heinz, Sep 16 2019

A317654 Number of free pure symmetric multifunctions whose leaves are a strongly normal multiset of size n.

Original entry on oeis.org

1, 3, 26, 375, 6696, 159837, 4389226, 144915350, 5377002075, 227624621051, 10632808475596, 550932945236121, 31062550998284221, 1907051034025848314, 126052420069459211076, 8956882232940915920404, 679298518935625486287703, 54868537321267493152151502, 4696952405203792017289469056

Offset: 1

Gus Wiseman, Aug 03 2018

nonn

A multiset is strongly normal if it spans an initial interval of positive integers with weakly decreasing multiplicities. A free pure symmetric multifunction f in EPSM is either (case 1) a positive integer, or (case 2) an expression of the form h[g_1, ..., g_k] where k > 0, h is in EPSM, each of the g_i for i = 1, ..., k is in EPSM, and for i < j we have g_i <= g_j under a canonical total ordering of EPSM, such as the Mathematica ordering of expressions.

			The a(3) = 26 free pure symmetric multifunctions:
1[1[1]], 1[1,1], 1[1][1],
1[1[2]], 1[2[1]], 1[1,2], 2[1[1]], 2[1,1], 1[1][2], 1[2][1], 2[1][1],
1[2[3]], 1[3[2]], 1[2,3], 2[1[3]], 2[3[1]], 2[1,3], 3[1[2]], 3[2[1]], 3[1,2], 1[2][3], 2[1][3], 1[3][2], 3[1][2], 2[3][1], 3[2][1].

Cf. A001003, A052893, A053492, A255906, A277996, A279944, A280000.

Cf. A317652, A317653, A317655, A317656, A317658.

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]]]];
exprUsing[m_]:=exprUsing[m]=If[Length[m]==0,{},If[Length[m]==1,{First[m]},Join@@Cases[Union[Table[PR[m[[s]],m[[Complement[Range[Length[m]],s]]]],{s,Take[Subsets[Range[Length[m]]],{2,-2}]}]],PR[h_,g_]:>Join@@Table[Apply@@@Tuples[{exprUsing[h],Union[Sort/@Tuples[exprUsing/@p]]}],{p,mps[g]}]]]];
got[y_]:=Join@@Table[Table[i,{y[[i]]}],{i,Range[Length[y]]}];
Table[Sum[Length[exprUsing[got[y]]],{y,IntegerPartitions[n]}],{n,6}]

\\ See links in A339645 for combinatorial species functions.
cycleIndexSeries(n)={my(p=O(x)); for(n=1, n, p = x*sv(1) + p*(sExp(p)-1)); p}
StronglyNormalLabelingsSeq(cycleIndexSeries(15)) \\ Andrew Howroyd, Jan 01 2021

Terms a(8) and beyond from Andrew Howroyd, Jan 01 2021

A317652 Number of free pure symmetric multifunctions whose leaves are an integer partition of n.

Original entry on oeis.org

1, 1, 2, 6, 22, 93, 421, 2010, 9926, 50357, 260728, 1372436, 7321982, 39504181, 215168221, 1181540841, 6534058589, 36357935615, 203414689462, 1143589234086, 6457159029573, 36602333187792, 208214459462774, 1188252476400972, 6801133579291811, 39032172166792887

Offset: 0

Gus Wiseman, Aug 03 2018

nonn

A free pure symmetric multifunction f in EPSM is either (case 1) a positive integer, or (case 2) an expression of the form h[g_1, ..., g_k] where k > 0, h is in EPSM, each of the g_i for i = 1, ..., k is in EPSM, and for i < j we have g_i <= g_j under a canonical total ordering of EPSM, such as the Mathematica ordering of expressions.

			The a(4) = 22 free pure symmetric multifunctions:
  1[1[1[1]]]  1[1[2]]  1[3]  2[2]  4
  1[1[1][1]]  1[2[1]]  3[1]
  1[1][1[1]]  2[1[1]]
  1[1[1]][1]  1[1][2]
  1[1][1][1]  1[2][1]
  1[1[1,1]]   2[1][1]
  1[1,1[1]]   1[1,2]
  1[1][1,1]   2[1,1]
  1[1,1][1]
  1[1,1,1]

Andrew Howroyd, Table of n, a(n) for n = 0..200

Cf. A001003, A052893, A053492, A277996, A279944, A280000.

Cf. A317653, A317654, A317655, A317656, A317658.

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]]]];
exprUsing[m_]:=exprUsing[m]=If[Length[m]==0,{{}},If[Length[m]==1,{First[m]},Join@@Cases[Union[Table[PR[m[[s]],m[[Complement[Range[Length[m]],s]]]],{s,Take[Subsets[Range[Length[m]]],{2,-2}]}]],PR[h_,g_]:>Join@@Table[Apply@@@Tuples[{exprUsing[h],Union[Sort/@Tuples[exprUsing/@p]]}],{p,mps[g]}]]]];
Table[Sum[Length[exprUsing[y]],{y,IntegerPartitions[n]}],{n,0,6}]

EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
seq(n)={my(v=[]); for(n=1, n, my(t=EulerT(v)); v=concat(v, 1 + sum(k=1, n-1, v[k]*t[n-k]))); concat([1],v)} \\ Andrew Howroyd, Aug 28 2018

Terms a(12) and beyond from Andrew Howroyd, Aug 28 2018

A317876 Number of free pure symmetric identity multifunctions (with empty expressions allowed) with one atom and n positions.

Original entry on oeis.org

1, 1, 2, 4, 10, 25, 67, 184, 519, 1489, 4342, 12812, 38207, 114934, 348397, 1063050, 3262588, 10064645, 31190985, 97061431, 303165207, 950115502, 2986817742, 9415920424, 29760442192, 94286758293, 299377379027, 952521579944, 3036380284111, 9696325863803

Offset: 1

Gus Wiseman, Aug 09 2018

nonn

A free pure symmetric identity multifunction (with empty expressions allowed) (FOI) 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 FOI, each of the g_i for i = 1, ..., k >= 0 is an FOI, and for i < j we have g_i < g_j under a canonical total ordering such as the Mathematica ordering of expressions. The number of positions in an FOI is the number of brackets [...] plus the number of o's.

Also the number of free orderless identity Mathematica expressions with one atom and n positions.

			The a(5) = 10 FOIs:
  o[o[o]]
  o[o][o]
  o[o[][]]
  o[o,o[]]
  o[][o[]]
  o[][][o]
  o[o[]][]
  o[][o][]
  o[o][][]
  o[][][][]

Andrew Howroyd, Table of n, a(n) for n = 1..200

Cf. A000081, A004111, A052893, A053492, A277996, A280000, A317652, A317653, A317654, A317875.

Cf. A317877, A317878, A317879, A317880, A317881.

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[Union[Sort/@Tuples[allIdExpr/@p]],UnsameQ@@#&]}],{p,IntegerPartitions[g]}]]];
Table[Length[allIdExpr[n]],{n,12}]

WeighT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v,n,(-1)^(n-1)/n))))-1,-#v)}
seq(n)={my(v=[1]); for(n=2, n, my(t=WeighT(v)); v=concat(v, v[n-1] + sum(k=1, n-2, v[k]*t[n-k-1]))); v} \\ Andrew Howroyd, Aug 19 2018

From Ilya Gutkovskiy, Apr 30 2019: (Start)
G.f. A(x) satisfies: A(x) = x * (1 + A(x) * exp(Sum_{k>=1} (-1)^(k+1)*A(x^k)/k)).
G.f.: A(x) = Sum_{n>=1} a(n)*x^n = x * (1 + (Sum_{n>=1} a(n)*x^n) * Product_{n>=1} (1 + x^n)^a(n)). (End)

Terms a(16) and beyond from Andrew Howroyd, Aug 19 2018

A317655 Number of free pure symmetric multifunctions with leaves a multiset whose multiplicities are the integer partition with Heinz number n.

Original entry on oeis.org

0, 1, 1, 2, 3, 8, 10, 15, 50, 35, 37, 96, 144, 160, 299, 184, 589, 840, 2483, 578, 1729, 750, 10746, 1627, 2246, 3578, 9357, 3367, 47420, 6397, 212668, 3155, 9818, 17280, 15666, 18250, 966324, 84232, 54990, 12471, 4439540, 45015

Offset: 1

Gus Wiseman, Aug 03 2018

nonn

The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

A free pure symmetric multifunction f in EPSM is either (case 1) a positive integer, or (case 2) an expression of the form h[g_1, ..., g_k] where k > 0, h is in EPSM, each of the g_i for i = 1, ..., k is in EPSM, and for i < j we have g_i <= g_j under a canonical total ordering of EPSM, such as the Mathematica ordering of expressions.

			The a(6) = 8 free pure symmetric multifunctions:
  1[1[2]]
  1[2[1]]
  2[1[1]]
  1[1][2]
  1[2][1]
  2[1][1]
  1[1,2]
  2[1,1]

Cf. A001003, A052893, A053492, A255906, A277996, A279944, A280000.

Cf. A317652, A317653, A317654, A317656, A317658.

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]]]];
exprUsing[m_]:=exprUsing[m]=If[Length[m]==0,{},If[Length[m]==1,{First[m]},Join@@Cases[Union[Table[PR[m[[s]],m[[Complement[Range[Length[m]],s]]]],{s,Take[Subsets[Range[Length[m]]],{2,-2}]}]],PR[h_,g_]:>Join@@Table[Apply@@@Tuples[{exprUsing[h],Union[Sort/@Tuples[exprUsing/@p]]}],{p,mps[g]}]]]];
got[y_]:=Join@@Table[Table[i,{y[[i]]}],{i,Range[Length[y]]}];
primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[exprUsing[got[Reverse[primeMS[n]]]]],{n,40}]

A317656 Number of free pure symmetric multifunctions whose leaves are the integer partition with Heinz number n.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 3, 1, 2, 1, 8, 1, 2, 2, 10, 1, 8, 1, 8, 2, 2, 1, 35, 1, 2, 3, 8, 1, 15, 1, 37, 2, 2, 2, 50, 1, 2, 2, 35, 1, 15, 1, 8, 8, 2, 1, 160, 1, 8, 2, 8, 1, 35, 2, 35, 2, 2, 1, 96, 1, 2, 8, 144, 2, 15, 1, 8, 2, 15, 1, 299, 1, 2, 8, 8, 2, 15, 1, 160

Offset: 1

Gus Wiseman, Aug 03 2018

nonn

A free pure symmetric multifunction f in EPSM is either (case 1) a positive integer, or (case 2) an expression of the form h[g_1, ..., g_k] where k > 0, h is in EPSM, each of the g_i for i = 1, ..., k is in EPSM, and for i < j we have g_i <= g_j under a canonical total ordering of EPSM, such as the Mathematica ordering of expressions.

			The a(12) = 8 free pure symmetric multifunctions are 1[1[2]], 1[2[1]], 1[1,2], 2[1[1]], 2[1,1], 1[1][2], 1[2][1], 2[1][1].

Cf. A001003, A052893, A053492, A255906, A277996, A279944, A280000.

Cf. A317652, A317653, A317654, A317655, A317658.

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]]]];
exprUsing[m_]:=exprUsing[m]=If[Length[m]==0,{},If[Length[m]==1,{First[m]},Join@@Cases[Union[Table[PR[m[[s]],m[[Complement[Range[Length[m]],s]]]],{s,Take[Subsets[Range[Length[m]]],{2,-2}]}]],PR[h_,g_]:>Join@@Table[Apply@@@Tuples[{exprUsing[h],Union[Sort/@Tuples[exprUsing/@p]]}],{p,mps[g]}]]]];
primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[exprUsing[primeMS[n]]],{n,100}]

A317659 Regular triangle where T(n,k) is the number of distinct free pure symmetric multifunctions (with empty expressions allowed) with one atom, n positions, and k leaves.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 3, 1, 0, 1, 6, 5, 1, 0, 1, 10, 17, 7, 1, 0, 1, 15, 43, 33, 9, 1, 0, 1, 21, 92, 118, 55, 11, 1, 0, 1, 28, 174, 341, 252, 82, 13, 1, 0, 1, 36, 302, 845, 935, 463, 115, 15, 1, 0, 1, 45, 490, 1864, 2921, 2103, 769, 153, 17, 1, 0, 1, 55, 755

Offset: 1

Gus Wiseman, Aug 03 2018

			The T(5,3) = 5 expressions are o[o[o]], o[o,o[]], o[][o,o], o[o][o], o[o,o][].
Triangle begins:
    1
    1    0
    1    1    0
    1    3    1    0
    1    6    5    1    0
    1   10   17    7    1    0
    1   15   43   33    9    1    0
    1   21   92  118   55   11    1    0
    1   28  174  341  252   82   13    1    0
    1   36  302  845  935  463  115   15    1    0
    1   45  490 1864 2921 2103  769  153   17    1    0
    1   55  755 3755 7981 8012 4145 1187  197   19    1    0

Mathematica Reference, Orderless

Cf. A001003, A052893, A053492, A255906, A277996, A279944, A280000.

Cf. A317652, A317653, A317654, A317655, A317656, A317676.

maxUsing[n_]:=If[n==1,{"o"},Join@@Cases[Table[PR[k,n-k-1],{k,n-1}],PR[h_,g_]:>Join@@Table[Apply@@@Tuples[{maxUsing[h],Union[Sort/@Tuples[maxUsing/@p]]}],{p,IntegerPartitions[g]}]]];
Table[Length[Select[maxUsing[n],Length[Position[#,"o"]]==k&]],{n,12},{k,n}]

cp's OEIS Frontend

A317658 Number of positions in the n-th free pure symmetric multifunction (with empty expressions allowed) with one atom.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A317776 Number of strict multiset partitions of normal multisets of size n, where a multiset is normal if it spans an initial interval of positive integers.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Extensions

A317654 Number of free pure symmetric multifunctions whose leaves are a strongly normal multiset of size n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

Extensions

A317652 Number of free pure symmetric multifunctions whose leaves are an integer partition of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A317876 Number of free pure symmetric identity multifunctions (with empty expressions allowed) with one atom and n positions.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A317655 Number of free pure symmetric multifunctions with leaves a multiset whose multiplicities are the integer partition with Heinz number n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A317656 Number of free pure symmetric multifunctions whose leaves are the integer partition with Heinz number n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A317659 Regular triangle where T(n,k) is the number of distinct free pure symmetric multifunctions (with empty expressions allowed) with one atom, n positions, and k leaves.

Views

Author