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

A300352 Number of strict trees of weight n with distinct leaves.

Original entry on oeis.org

1, 1, 2, 2, 3, 6, 8, 11, 17, 40, 48, 76, 109, 159, 400, 470, 745, 1057, 1576, 2103, 5267, 6022, 9746, 13390, 20099, 26542, 39396, 82074, 101387, 152291, 215676, 308937, 423587, 596511, 799022, 1623311, 1960223, 2947722, 4048704, 5845982, 7794809, 11028888
Offset: 1

Views

Author

Gus Wiseman, Mar 03 2018

Keywords

Comments

A strict tree of weight n > 0 is either a single node of weight n, or a sequence of two or more strict trees with strictly decreasing weights summing to n.

Examples

			The a(8) = 11 strict trees with distinct leaves: 8, (71), ((52)1), ((43)1), (62), ((51)2), (53), ((41)3), (5(21)), (521), (431).

Crossrefs

Cf. A000009, A056239, A063834, A196545, A246867, A273873, A281145, A289501, A294018, A294079, A296150, A299201, A299202, A299203, A300353, A300354, A300355.

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=
Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
str[q_]:=str[q]=If[Length[q]===1,1,Total[Times@@@Map[str,Select[sps[q],And[Length[#]>1,UnsameQ@@Total/@#]&],{2}]]];
Table[Total[str/@Select[IntegerPartitions[n],UnsameQ@@#&]],{n,1,20}]

Formula

a(n) = Sum_{i=1..A000009(n)} A294018(A246867(n,i)).

A300353 Number of strict trees of weight n with odd leaves.

Original entry on oeis.org

1, 1, 0, 1, 1, 2, 2, 4, 7, 14, 24, 46, 92, 186, 368, 750, 1529, 3160, 6510, 13590, 28374, 59780, 125732, 266468, 564188, 1202842, 2560106, 5484304, 11732400, 25229068, 54187918, 116938702, 252039411, 545593378, 1179545874, 2560009400, 5550315640, 12075064432
Offset: 0

Views

Author

Gus Wiseman, Mar 03 2018

Keywords

Comments

This sequence is initially dominated by A300352 but eventually becomes much greater.
A strict tree of weight n > 0 is either a single node of weight n, or a sequence of two or more strict trees with strictly decreasing weights summing to n.

Examples

			The a(8) = 7 strict trees with odd leaves: (71), (53), (((51)1)1), (((31)3)1), (((31)1)3), ((31)31), (((((31)1)1)1)1).

Links

Crossrefs

Cf. A000009, A063834, A078408, A089259, A196545, A279374, A279785, A289501, A294018, A294079, A299203, A300300, A300301, A300351, A300352, A300354, A300355.

Programs

  • Mathematica
    d[n_]:=d[n]=If[EvenQ[n],0,1]+Sum[Times@@d/@y,{y,Select[IntegerPartitions[n],Length[#]>1&&UnsameQ@@#&]}];
Table[d[n],{n,40}]
  • PARI
    seq(n)={my(v=vector(n)); v[1]=1; for(n=2, n, v[n] = polcoef(x/(1-x^2) + prod(k=1, n-1, 1 + v[k]*x^k + O(x*x^n)), n)); concat([1], v)} \\ Andrew Howroyd, Aug 25 2018

Formula

O.g.f: (1 + x/(1-x^2) + Product_{i>0} (1 + a(i)x^i))/2.
a(n) = Sum_{i=1..A000009(n)} A294018(A300351(n,i)).

A300354 Number of enriched p-trees of weight n with distinct leaves.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 8, 8, 13, 17, 54, 56, 98, 125, 195, 500, 606, 921, 1317, 1912, 2635, 6667, 7704, 12142, 16958, 24891, 33388, 47792, 106494, 126475, 195475, 268736, 393179, 523775, 750251, 979518, 2090669, 2457315, 3759380, 5066524, 7420874, 9726501, 13935546
Offset: 0

Views

Author

Gus Wiseman, Mar 03 2018

Keywords

Comments

An enriched p-tree of weight n > 0 is either a single node of weight n, or a sequence of two or more enriched p-trees with weakly decreasing weights summing to n.

Examples

			The a(6) = 8 enriched p-trees with distinct leaves: 6, (42), (51), ((31)2), ((32)1), (3(21)), ((21)3), (321).

Crossrefs

Cf. A000009, A000041, A063834, A196545, A246867, A273873, A281145, A289501, A290261, A294018, A296150, A299201, A299202, A299203, A300352, A300353, A300355.

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
ept[q_]:=ept[q]=If[Length[q]===1,1,Total[Times@@@Map[ept,Join@@Function[sptn,Join@@@Tuples[Permutations/@GatherBy[sptn,Total]]]/@Select[sps[q],Length[#]>1&],{2}]]];
Table[Total[ept/@Select[IntegerPartitions[n],UnsameQ@@#&]],{n,1,30}]

Formula

a(n) = Sum_{i=1..A000009(n)} A299203(A246867(n,i)).

A294079 Strict Moebius function of the multiorder of integer partitions indexed by Heinz numbers.

Original entry on oeis.org

0, 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, 1, 1, 0, -1, -1, -1, -1, 1, -1, -1, -1, 1, 2, 1, 1, 1, -1, 1, 1, 0, 1, -1, 1, 1, -1, -1, -1, -1, -1, 1, -3, 1, -1, 1, 0, -1, 2, 1, 1, -1, 1, 1, 2, 1, -1, 1, 1, -1, 2, 1, 1, 0, -1, 1, -2, -1, -1, -1, -1, 1, -3, -1
Offset: 1

Views

Author

Gus Wiseman, Feb 07 2018

Keywords

Comments

By convention a(1) = 0.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Crossrefs

Cf. A000041, A000720, A056239, A063834, A196545, A273873, A289501, A294018, A294019, A296150, A299201, A299202, A299203.

Programs

  • Mathematica
    nn=120;
ptns=Table[If[n===1,{},Join@@Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]],{n,nn}];
tris=Join@@Map[Tuples[IntegerPartitions/@#]&,ptns];
qmu[y_]:=qmu[y]=If[Length[y]===1,1,-Sum[Times@@qmu/@t,{t,Select[tris,And[Length[#]>1,Sort[Join@@#,Greater]===y,UnsameQ@@#]&]}]];
qmu/@ptns

Formula

mu(y) = Sum_{g(t)=y} (-1)^d(t), where the sum is over all strict trees (A273873) whose multiset of leaves is the integer partition y, and d(t) is the number of non-leaf nodes in t.

A294080 Same-tree Moebius function of the multiorder of integer partitions indexed by Heinz numbers.

Original entry on oeis.org

0, 1, 1, -1, 1, 0, 1, -1, -1, 0, 1, 2, 1, 0, 0, -2, 1, 0, 1, 0, 0, 0, 1, 0, -1, 0, -1, 0, 1, 0, 1, -1, 0, 0, 0, 3, 1, 0, 0, 2, 1, 0, 1, 0, 0, 0, 1, -3, -1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 2, -1, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, -2, 0, 1, -4, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 8
Offset: 1

Views

Author

Gus Wiseman, Feb 07 2018

Keywords

Comments

By convention a(1) = 0.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Links

Crossrefs

Cf. A000041, A056239, A063834, A196545, A273873, A281145, A289501, A294018, A294019, A296150, A299202, A299203.

Programs

  • Mathematica
    nn=120;
ptns=Table[If[n===1,{},Join@@Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]],{n,nn}];
tris=Join@@Map[Tuples[IntegerPartitions/@#]&,ptns];
rmu[y_]:=rmu[y]=If[Length[y]===1,1,-Sum[Times@@rmu/@t,{t,Select[tris,And[Length[#]>1,Sort[Join@@#,Greater]===y,SameQ@@Total/@#]&]}]];
rmu/@ptns
  • PARI
    A056239(n) = { my(f); if(1==n, 0, f=factor(n); sum(i=1, #f~, f[i,2] * primepi(f[i,1]))); }
muifbalancedfactorization(v) = if(!#v, 1, my(pw=A056239(v[1]), m=-1); for(i=1,#v,if(A056239(v[i])!=pw,return(0), m *= A294080(v[i]))); (m));
A294080aux(n, m, facs) = if(1==n, muifbalancedfactorization(Vec(facs)), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A294080aux(n/d, m, newfacs))); (s));
A294080(n) = if(1==n,0,if(isprime(n),1,A294080aux(n, n-1, List([]))));
\\ A memoized implementation:
map294080 = Map();
A294080(n) = if(1==n,0,if(isprime(n),1,if(mapisdefined(map294080,n), mapget(map294080,n), my(v=A294080aux(n, n-1, List([]))); mapput(map294080,n,v); (v)))); \\ Antti Karttunen, Sep 22 2018

Formula

mu(y) = Sum_{g(t)=y} (-1)^d(t), where the sum is over all same-trees (A281145, A294019) whose multiset of leaves is the integer partition y, and d(t) is the number of non-leaf nodes in t.

A294019 Number of same-trees whose leaves are the parts of the integer partition with Heinz number n.

Original entry on oeis.org

0, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 2, 1, 0, 0, 2, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 3, 1, 0, 0, 2, 1, 0, 1, 0, 0, 0, 1, 3, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 2, 3, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 2, 0, 1, 4, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 8
Offset: 1

Views

Author

Gus Wiseman, Feb 07 2018

Keywords

Comments

By convention a(1) = 0.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			The a(108) = 8 same-trees: ((22)(2(11))), ((22)((11)2)), ((2(11))(22)), (((11)2)(22)), (222(11)), (22(11)2), (2(11)22), ((11)222).
From _Antti Karttunen_, Sep 22 2018: (Start)
For 12 = prime(1)^2 * prime(2)^1, we have the following two cases: 2(11) and (11)2, thus a(12) = 2.
For 36 = prime(1)^2 * prime(2)^2, we have the following cases: (11)22, 2(11)2, 22(11), thus a(36) = 3.
For 144  = prime(1)^4 * prime(2)^2, we have the following 14 cases: (1111)(22), (22)(1111); ((11)(11))(22), (22)((11)(11)); (11)(11)22, (11)2(11)2, (11)22(11), 2(11)2(11), 2(11)(11)2, 22(11)(11); ((11)2)(11(2)), ((11)2)(2(11)), (2(11))((11)2), (2(11))(2(11)), thus a(144) = 14.
For n = 8775 = 3^3 * 5^2 * 13^1 = prime(2)^3 * prime(3)^2 * prime(6)^1, we have the following six cases: (222)(33)6, (222)6(33), (33)(222)6, (33)6(222), 6(222)(33), 6(33)(222), thus a(8775) = 6.
(End)

Links

Crossrefs

Cf. A000005, A000041, A000720, A001222, A006241, A056239, A063834, A196545, A215366, A273873, A281145, A289501, A296150, A299200, A299201, A299202, A299203, A294018, A294080.

Programs

  • Mathematica
    nn=120;
ptns=Table[If[n===1,{},Join@@Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]],{n,nn}];
tris=Join@@Map[Tuples[IntegerPartitions/@#]&,ptns];
qci[y_]:=qci[y]=If[Length[y]===1,1,Sum[Times@@qci/@t,{t,Select[tris,And[Length[#]>1,Sort[Join@@#,Greater]===y,SameQ@@Total/@#]&]}]];
qci/@ptns
  • PARI
    A056239(n) = { my(f); if(1==n, 0, f=factor(n); sum(i=1, #f~, f[i,2] * primepi(f[i,1]))); }
productifbalancedfactorization(v) = if(!#v, 1, my(pw=A056239(v[1]), m=1); for(i=1,#v,if(A056239(v[i])!=pw,return(0), m *= A294019(v[i]))); (m));
A294019aux(n, m, facs) = if(1==n, productifbalancedfactorization(Vec(facs)), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A294019aux(n/d, m, newfacs))); (s));
A294019(n) = if(1==n,0,if(isprime(n),1,A294019aux(n, n-1, List([]))));
\\ A memoized implementation:
map294019 = Map();
A294019(n) = if(1==n,0,if(isprime(n),1,if(mapisdefined(map294019,n), mapget(map294019,n), my(v=A294019aux(n, n-1, List([]))); mapput(map294019,n,v); (v)))); \\ Antti Karttunen, Sep 22 2018

Formula

A281145(n) = Sum_{i=1..A000041(n)} a(A215366(n,i)).
a(p^n) = A006241(n) for any prime p and exponent n >= 1. - Antti Karttunen, Sep 22 2018

A301365 Regular triangle where T(n,k) is the number of strict trees of weight n with k leaves.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 2, 2, 1, 0, 1, 2, 4, 4, 1, 0, 1, 3, 7, 9, 7, 1, 0, 1, 3, 10, 19, 20, 11, 1, 0, 1, 4, 15, 35, 51, 43, 16, 1, 0, 1, 4, 18, 55, 104, 123, 84, 22, 1, 0, 1, 5, 25, 84, 196, 298, 284, 153, 29, 1, 0, 1, 5, 30, 120, 331, 624, 783, 614, 260, 37
Offset: 1

Views

Author

Gus Wiseman, Mar 19 2018

Keywords

Comments

A strict tree of weight n > 0 is either a single node of weight n, or a finite sequence of two or more strict trees with strictly decreasing weights summing to n.

Examples

			Triangle begins:
  1
  1   0
  1   1   0
  1   1   1   0
  1   2   2   1   0
  1   2   4   4   1   0
  1   3   7   9   7   1   0
  1   3  10  19  20  11   1   0
  1   4  15  35  51  43  16   1   0
The T(7,3) = 7 strict trees: ((51)1), ((42)1), ((41)2), ((32)2), (4(21)), ((31)3), (421).

Links

Crossrefs

Row sums are A273873.
Cf. A004111, A008284, A032305, A055277, A063834, A281145, A289501, A294018, A294079, A300352, A300442, A300443, A301342, A301364-A301368.

Programs

  • Mathematica
    strtrees[n_]:=Prepend[Join@@Table[Tuples[strtrees/@ptn],{ptn,Select[IntegerPartitions[n],Length[#]>1&&UnsameQ@@#&]}],n];
Table[Length[Select[strtrees[n],Count[#,_Integer,{-1}]===k&]],{n,12},{k,n}]
  • PARI
    A(n)={my(v=vector(n)); for(n=1, n, v[n] = y + polcoef(prod(k=1, n-1, 1 + v[k]*x^k + O(x*x^n)), n)); vector(n, k, Vecrev(v[k]/y, k))}
my(T=A(10));for(n=1, #T, print(T[n])) \\ Andrew Howroyd, Aug 26 2018
Showing 1-7 of 7 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.