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.
%I A097681 #16 Feb 16 2025 08:32:54
%S A097681 1,6,36,216,1296,7776,47376,314496,2612736,28740096,368395776,
%T A097681 4796983296,60300205056,750367328256,10151357239296,164475953381376,
%U A097681 3110937349718016,61410199093641216,1174438559356747776
%N A097681 E.g.f.: (1/(1-x^6))*exp( 6*sum_{i>=0} x^(6*i+1)/(6*i+1) ) for an order-6 linear recurrence with varying coefficients.
%C A097681 Limit_{n->inf} n*n!/a(n) = 6*c = 0.1140186893... where c = 6*exp(psi(1/6)+EulerGamma) = 0.0190031148...(A097671) and EulerGamma is the Euler-Mascheroni constant (A001620) and psi() is the Digamma function (see Mathworld link).
%D A097681 Mohammad K. Azarian, Problem 1218, Pi Mu Epsilon Journal, Vol. 13, No. 2, Spring 2010, p. 116. Solution published in Vol. 13, No. 3, Fall 2010, pp. 183-185.
%D A097681 A. M. Odlyzko, Linear recurrences with varying coefficients, in Handbook of Combinatorics, Vol. 2, R. L. Graham, M. Grotschel and L. Lovasz, eds., Elsevier, Amsterdam, 1995, pp. 1135-1138.
%H A097681 Benoit Cloitre, <a href="/A097679/a097679.pdf">On a generalization of Euler-Gauss formula for the Gamma function</a>, preprint 2004.
%H A097681 Andrew Odlyzko, <a href="http://www.dtc.umn.edu/~odlyzko/doc/asymptotic.enum.pdf">Asymptotic enumeration methods</a>, in Handbook of Combinatorics, vol. 2, 1995, pp. 1063-1229.
%H A097681 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/DigammaFunction.html">Digamma Function</a>.
%F A097681 For n>=6: a(n) = 6*a(n-1) + n!/(n-6)!*a(n-6); for n<6: a(n)=6^n. E.g.f.: 1/(1-x^6)*(1+x)/(1-x)*sqrt((1+x+x^2)/(1-x+x^2))* exp(sqrt(3)*atan(sqrt(3)*x/(1-x^2))).
%e A097681 The sequence {1, 6, 36/2!, 216/3!, 1296/4!, 7776/5!, 47376/6!,...} is generated by a recursion described by Benoit Cloitre's generalized Euler-Gauss formula for the Gamma function (see Cloitre link).
%o A097681 (PARI) {a(n)=n!*polcoeff(1/(1-x^6)*exp(6*sum(i=0,n,x^(6*i+1)/(6*i+1)))+x*O(x^n),n)}
%o A097681 (PARI) a(n)=if(n<0,0,if(n==0,1,6*a(n-1)+if(n<6,0,n!/(n-6)!*a(n-6))))
%Y A097681 Cf. A097671, A097677-A097680, A097682-A097682.
%K A097681 nonn
%O A097681 0,2
%A A097681 _Paul D. Hanna_, Sep 01 2004