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.

A331964 Number of semi-lone-child-avoiding rooted identity trees with n vertices.

Original entry on oeis.org

1, 1, 0, 1, 0, 1, 1, 2, 2, 4, 6, 10, 16, 27, 44, 74, 123, 209, 353, 602, 1026, 1760, 3019, 5203, 8977, 15538, 26930, 46792, 81415, 141939, 247795, 433307, 758672, 1330219, 2335086, 4104064, 7220937, 12718694, 22424283, 39574443, 69903759, 123584852, 218668323
Offset: 1

Views

Author

Gus Wiseman, Feb 04 2020

Keywords

Comments

A rooted tree is semi-lone-child-avoiding if there are no vertices with exactly one child unless that child is an endpoint/leaf. It is an identity tree if the branches of any given vertex are all distinct.

Examples

			The a(9) = 2 through a(12) = 10 semi-lone-child-avoiding rooted identity trees:
  ((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(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(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(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(o))((o)(o(o))))
                                                       (o((o)((o)(o(o)))))
		

Crossrefs

The non-semi version is A000007.
Matula-Goebel numbers of these trees are A331963.
Rooted identity trees are A004111.
Semi-lone-child-avoiding rooted trees are A331934.

Programs

  • Mathematica
    ssei[n_]:=Switch[n,1,{{}},2,{{{}}},_,Join@@Function[c,Select[Union[Sort/@Tuples[ssei/@c]],UnsameQ@@#&]]/@Rest[IntegerPartitions[n-1]]];
    Table[Length[ssei[n]],{n,15}]
  • PARI
    WeighT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v,n,(-1)^(n-1)/n))))-1,-#v)}
    seq(n)={my(v=[1,1]); for(n=2, n-1, v=concat(v, WeighT(v)[n] - v[n])); v} \\ Andrew Howroyd, Feb 09 2020

Extensions

Terms a(36) and beyond from Andrew Howroyd, Feb 09 2020