The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A155200 G.f.: A(x) = exp( Sum_{n>=1} 2^(n^2) * x^n/n ), a power series in x with integer coefficients. 49
 1, 2, 10, 188, 16774, 6745436, 11466849412, 80444398636280, 2306003967992402758, 268654794629082985019564, 126765597346260977505891041836, 241678070948246232010898235031930952, 1858395916567787793818891330877931472153500, 57560683587056536617649234722821582390470430186648 (list; graph; refs; listen; history; text; internal format)
 OFFSET 0,2 COMMENTS More generally, it appears that for m integer, exp( Sum_{n >= 1} m^(n^2) * x^n/n ) is a power series in x with integer coefficients. This is correct: if b(n) = m^(n^2) then by the little Fermat theorem the Gauss congruences hold: b(n*p^k) == b(n*p^(k-1)) ( mod p^k ) for all prime p and positive integers n and k. Then apply Stanley, Ch. 5, Ex. 5.2(a). - Peter Bala, Dec 25 2019 Conjecture: highest exponent of 2 dividing a(n) = A000120(n) = number of 1's in binary expansion of n, so that a(n)/2^A000120(n) is odd for n >= 0. - Paul D. Hanna, Sep 01 2009 REFERENCES R. P. Stanley. Enumerative combinatorics, Vol. 2. Volume 62 of Cambridge Studies in Advanced Mathematics, Cambridge University Press, Cambridge, 1999. LINKS Muniru A Asiru, Table of n, a(n) for n = 0..56 Sawian Jaidee, Patrick Moss, Tom Ward, Time-changes preserving zeta functions, arXiv:1809.09199 [math.DS], 2018. FORMULA Equals column 0 of triangle A155810. G.f. satisfies: 2*A(x)*A(4x) + 8*x*A(x)*A'(4x) - A'(x)*A(4x) = 0. - Paul D. Hanna, Feb 24 2009 From Paul D. Hanna, Mar 11 2009: (Start) The differential equation implies recurrence: n*a(n) = 2*a(n-1) + sum(k = 1, n - 1, 4^k*a(k)*[2*(k+1)*a(n-1-k) - (n-k)*a(n-k)] for n > 0, with a(0) = 1. G.f. A(x) generates A156631: A156631(n) = [x^n] A(x)^(2^n) for n >= 0, where the g.f. of A156631 = Sum_{n >= 0} [Sum_{k >= 1} (2^n*2^k*x)^k/k]^n/n!. (End) a(n) = (1/n)*Sum_{k = 1..n} 2^(k^2)*a(n-k), a(0) = 1. - Vladeta Jovovic, Feb 04 2009 Euler transform of A159034. - Vladeta Jovovic, Apr 02 2009 a(n) = B_n( 0!*2^(1^2), 1!*2^(2^2), 2!*2^(3^2), ..., (n-1)!*2^(n^2) ) / n!, where B_n() is the complete Bell polynomial. - Max Alekseyev, Oct 10 2014 a(n) ~ 2^(n^2) / n. - Vaclav Kotesovec, Oct 09 2019 EXAMPLE G.f.: A(x) = 1 + 2*x + 10*x^2 + 188*x^3 + 16774*x^4 + 6745436*x^5 +... log(A(x)) = 2*x + 2^4*x^2/2 + 2^9*x^3/3 + 2^16*x^4/4 + 2^25*x^5/5 +... MAPLE seq(coeff(series(exp(add(2^(k^2)*x^k/k, k=1..n)), x, n+1), x, n), n = 0 .. 15); # Muniru A Asiru, Dec 19 2018 MATHEMATICA nmax = 14; Exp[Sum[2^(n^2) x^n/n, {n, 1, nmax}]] + O[x]^nmax // CoefficientList[#, x]& (* Jean-François Alcover, Feb 14 2019 *) PROG (PARI) {a(n)=polcoeff(exp(sum(m=1, n+1, 2^(m^2)*x^m/m)+x*O(x^n)), n)} (PARI) {a(n)=if(n==0, 1, (1/n)*(2*a(n-1) + sum(k=1, n-1, 4^k*a(k)*(2*(k+1)*a(n-1-k) - (n-k)*a(n-k)))))} \\ Paul D. Hanna, Mar 11 2009 (PARI) {a(n)=if(n==0, 1, (1/n)*sum(k=1, n, 2^(k^2)*a(n-k)))} \\ Paul D. Hanna, Sep 01 2009 CROSSREFS Cf. A155201, A155202, A155810 (triangle), variants: A155203, A155207. Cf. A000120, A156631, A159034. Sequence in context: A037267 A177399 A194971 * A264563 A156510 A246532 Adjacent sequences: A155197 A155198 A155199 * A155201 A155202 A155203 KEYWORD nonn AUTHOR Paul D. Hanna, Feb 04 2009 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified June 15 11:51 EDT 2024. Contains 373407 sequences. (Running on oeis4.)