|
%I #9 Feb 09 2020 18:58:24
%S 1,0,1,1,2,3,5,9,16,30,55,105,200,388,754,1483,2923,5807,11575,23190,
%T 46608,94043,190287,386214,785831,1602952,3276845,6712905,13778079,
%U 28330583,58350582,120370731,248676129,514459237,1065696295,2210302177,4589599429,9540623926
%N Number of lone-child-avoiding rooted semi-identity trees with n vertices.
%C Lone-child-avoiding means there are no unary branchings.
%C In a semi-identity tree, the non-leaf branches of any given vertex are distinct.
%H Andrew Howroyd, <a href="/A331966/b331966.txt">Table of n, a(n) for n = 1..1000</a>
%H Gus Wiseman, <a href="https://docs.google.com/document/d/e/2PACX-1vS1zCO9fgAIe5rGiAhTtlrOTuqsmuPos2zkeFPYB80gNzLb44ufqIqksTB4uM9SIpwlvo-oOHhepywy/pub">Sequences counting series-reduced and lone-child-avoiding trees by number of vertices.</a>
%e The a(1) = 1 through a(9) = 16 trees (empty column shown as dot):
%e o . (oo) (ooo) (oooo) (ooooo) (oooooo) (ooooooo) (oooooooo)
%e (o(oo)) (o(ooo)) (o(oooo)) (o(ooooo)) (o(oooooo))
%e (oo(oo)) (oo(ooo)) (oo(oooo)) (oo(ooooo))
%e (ooo(oo)) (ooo(ooo)) (ooo(oooo))
%e (o(o(oo))) (oooo(oo)) (oooo(ooo))
%e ((oo)(ooo)) (ooooo(oo))
%e (o(o(ooo))) ((oo)(oooo))
%e (o(oo(oo))) (o(o(oooo)))
%e (oo(o(oo))) (o(oo)(ooo))
%e (o(oo(ooo)))
%e (o(ooo(oo)))
%e (oo(o(ooo)))
%e (oo(oo(oo)))
%e (ooo(o(oo)))
%e ((oo)(o(oo)))
%e (o(o(o(oo))))
%t ssb[n_]:=If[n==1,{{}},Join@@Function[c,Select[Union[Sort/@Tuples[ssb/@c]],UnsameQ@@DeleteCases[#,{}]&]]/@Rest[IntegerPartitions[n-1]]];
%t Table[Length[ssb[n]],{n,10}]
%o (PARI) WeighT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v,n,(-1)^(n-1)/n))))-1,-#v)}
%o seq(n)={my(v=[0, 0]); for(n=2, n-1, v=concat(v, 1 + vecsum(WeighT(v)) - v[n])); v[1]=1; v} \\ _Andrew Howroyd_, Feb 09 2020
%Y The non-semi case is A000007.
%Y Lone-child-avoiding rooted trees are A001678.
%Y The locally disjoint case is A212804.
%Y Not requiring lone-child-avoidance gives A306200.
%Y Matula-Goebel numbers of these trees are A331965.
%Y The semi-lone-child-avoiding version is A331993.
%Y Cf. A000081, A004111, A291636, A300660, A306202, A316694, A331683, A331686, A331783, A331875, A331964, A331994.
%K nonn
%O 1,5
%A _Gus Wiseman_, Feb 05 2020
%E Terms a(31) and beyond from _Andrew Howroyd_, Feb 09 2020
| 