OFFSET
1,2
COMMENTS
Compare to: [x^n] exp( n^4 * x ) = n^3 * [x^(n-1)] exp( n^4 * x ) for n>=1.
It is remarkable that this sequence should consist entirely of integers.
LINKS
Paul D. Hanna, Table of n, a(n) for n = 1..200
FORMULA
O.g.f. equals the logarithm of the e.g.f. of A300596.
From Seiichi Manyama, Apr 13 2026: (Start)
G.f. A(x) satisfies [x^n] exp(n^4*A(x)) * (1 - n^4*x) = 0, for n >= 1.
a(n) = n^(4*n-5) - (1/n) * Sum_{k=1..n-1} k * c_n(k) * e_n(n-k),
where c_n(k) = a(k) - n^(4*k-4)/k for 1 <= k <= n-1,
and e_n(0) = 1, e_n(k) = (n^4/k) * Sum_{j=1..k} j * c_n(j) * e_n(k-j) for 1 <= k <= n-1. (End)
EXAMPLE
O.g.f.: A(x) = x + 8*x^2 + 2187*x^3 + 2351104*x^4 + 6153518125*x^5 + 31779658925496*x^6 + 287364845865893467*x^7 + 4200677982722915635200*x^8 + ...
where
exp(A(x)) = 1 + x + 17*x^2/2! + 13171*x^3/3! + 56479849*x^4/4! + 738706542221*x^5/5! + 22885801082965201*x^6/6! + 1448479282286023114807*x^7/7! + ... + A300596(n)*x^n/n! + ...
such that: [x^n] exp( n^4 * A(x) ) = n^4 * [x^(n-1)] exp( n^4 * A(x) ).
PROG
(PARI) {a(n) = my(A=[1]); for(i=1, n+1, A=concat(A, 0); V=Vec(Ser(A)^((#A-1)^4)); A[#A] = ((#A-1)^4*V[#A-1] - V[#A])/(#A-1)^4 ); polcoeff( log(Ser(A)), n)}
for(n=1, 20, print1(a(n), ", "))
(Ruby)
def A300597(n)
a = [0]
(1..n).each{|i|
c = [0] + (1..i - 1).map{|k| a[k] - i ** (4 * k - 4) / k.to_r}
e = [1]
(1..i - 1).each{|k| e << i ** 4 / k.to_r * (1..k).inject(0){|s, j| s + j * c[j] * e[k - j]}}
a << i ** (4 * i - 5) - (1..i - 1).inject(0){|s, k| s + k * c[k] * e[i - k]}.to_i / i
}
a[1..-1]
end
p A300597(20) # Seiichi Manyama, Apr 13 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Mar 09 2018
STATUS
approved