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.

A325205 a(n) is the number of labeled rooted trees on a set of size n where each node has at most 6 neighbors that are further away from the root than the node itself.

This page as a plain text file.
%I A325205 #9 Dec 22 2024 09:08:09
%S A325205 0,1,2,9,64,625,7776,117649,2097144,43046136,999970020,25936053990,
%T A325205 742947675624,23295384644532,793591829158128,29187143427692250,
%U A325205 1152639088016576160,48646833059722978080,2185150741063924810176,104085328898784937079376,5240461483486301616704160
%N A325205 a(n) is the number of labeled rooted trees on a set of size n where each node has at most 6 neighbors that are further away from the root than the node itself.
%C A325205 A preimage constraint on a function is a set of nonnegative integers such that the size of the inverse image of any element is one of the values in that set. View a labeled rooted tree as an endofunction on the set {1,2,...,n} by sending every non-root node to its neighbor that is closer to the root and sending the root to itself. Thus, a(n) is the number of endofunctions on a set of size n with exactly one cyclic point and such that each preimage has at most 6 entries.
%H A325205 B. Otto, <a href="https://arxiv.org/abs/1903.00542">Coalescence under Preimage Constraints</a>, arXiv:1903.00542 [math.CO], 2019, Corollaries 5.3 and 7.8.
%F A325205 a(n) = (n-1)! * [x^(n-1)] e_6(x)^n, where e_k(x) is the truncated exponential 1 + x + x^2/2! + ... + x^k/k!. The link above yields explicit constants c_k, r_k so that the columns are asymptotically c_6 * n^(-3/2) * r_6^-n.
%o A325205 (Python)
%o A325205 # print first num_entries entries in the sequence
%o A325205 import math, sympy; x=sympy.symbols('x')
%o A325205 k=6; num_entries = 64
%o A325205 P=range(k+1); eP=sum([x**d/math.factorial(d) for d in P]); r = [0,1]; curr_pow = eP
%o A325205 for term in range(1,num_entries-1):
%o A325205     curr_pow=(curr_pow*eP).expand()
%o A325205     r.append(curr_pow.coeff(x**term)*math.factorial(term))
%o A325205 print(r)
%Y A325205 Column k=6 of A325201; see that entry for sequences related to other preimage constraints constructions.
%K A325205 easy,nonn
%O A325205 0,3
%A A325205 _Benjamin Otto_, Apr 11 2019