OFFSET
0,3
COMMENTS
Conjectures:
(1) a(2^n) is odd for n>=0.
(2) a(n) is odd iff n is a Fibbinary number: a(A003714(k)) is odd for k>=0.
(3) The number of odd terms between a(2^n) and a(2^(n+1)-1), inclusively, is Fibonacci(n+1), for n>=0.
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..520
FORMULA
G.f. A(x) satisfies:
(1) [x^(n-1)] 1/A(x)^(n^2) = 0 for n>2.
(2) [x^(n-1)] (1/x)*Series_Reversion( x*A(x)^n ) = 0 for n>2.
(3) [x^(n-2)] ( (1/x)*Series_Reversion( x*A(x)^n ) )^(1/n) = 0 for n>3.
a(n) ~ (-1)^n * c * 2^(2*n) * n^(n-5/2) / (exp(n) * d^n * (2-d)^n), where d = -LambertW(-2*exp(-2)) = -A226775 = 0.40637573995995990767695812412483975821... and c = 0.01284812446900190... - Vaclav Kotesovec, Sep 27 2017
EXAMPLE
G.f.: A(x) = 1 + x + 5*x^2 + 34*x^3 + 273*x^4 + 2331*x^5 + 22110*x^6 + 190450*x^7 + 2540975*x^8 - 1071509*x^9 + 1200284739*x^10 - 50263360280*x^11 + 3102388877800*x^12 - 199436004737160*x^13 + 14155468007742978*x^14 - 1088800915851203694*x^15 +...
such that the coefficient of x^n in 1/A(x)^(n^2) equals zero for n>1.
Notice that a(n) seems to be odd only when n is a Fibbinary number (A003714):
[0, 1, 2, 4, 5, 8, 9, 10, 16, 17, 18, 20, 21, 32, 33, 34, 36, 37, 40, 41, 42, 64, 65, 66, 68, 69, 72, 73, 74, 80, 81, 82, 84, 85, 128, 29, 130, 132, 133, ...].
RELATED TABLES.
(1) The table of coefficients in 1/A(x)^(n^2) begins:
n=1: [1, -1, -4, -25, -194, -1603, -15264, -122316, -1897710, ...];
n=2: [1, -4, -10, -56, -427, -3360, -33546, -218880, -5179834, ...];
n=3: [1, -9, 0, -21, -252, -1701, -25992, 2970, -7903413, ...];
n=4: [1, -16, 56, 0, -84, 784, -18656, 384896, -13426530, ...];
n=5: [1, -25, 200, -525, 0, 2695, -38600, 878150, -26292375, ...];
n=6: [1, -36, 486, -3000, 7821, 0, -101322, 1916352, -52357590, ...];
n=7: [1, -49, 980, -10241, 58898, -170079, 0, 4515000, -108626140, ...];
n=8: [1, -64, 1760, -27136, 256048, -1500352, 4979712, 0, -234893352, ...];
n=9: [1, -81, 2916, -61425, 838026, -7720839, 48097152, -184870512, 0, ...]; ...
such that the main diagonal is all zeros after the initial terms.
(2) The table of coefficients in (1/x)*Series_Reversion( x*A(x)^n ) begins:
n=1: [1, -1, -3, -14, -85, -504, -4424, -6796, -878157, ...];
n=2: [1, -2, -3, -8, -40, -60, -2604, 48112, -1747260, ...];
n=3: [1, -3, 0, 2, -15, 189, -3850, 101700, -3340845, ...];
n=4: [1, -4, 6, 0, -35, 396, -7182, 181824, -5817510, ...];
n=5: [1, -5, 15, -30, 0, 714, -13335, 315060, -9679455, ...];
n=6: [1, -6, 27, -104, 315, 0, -19957, 532848, -15864336, ...];
n=7: [1, -7, 42, -238, 1260, -5481, 0, 713796, -25010433, ...];
n=8: [1, -8, 60, -448, 3310, -23352, 136696, 0, -31112163, ...];
n=9: [1, -9, 81, -750, 7065, -66420, 598626, -4474764, 0, ...]; ...
in which the main diagonal is all zeroes after the initial terms.
(3) The table of coefficients in ((1/x)*Series_Reversion( x*A(x)^n ))^(1/n) begins:
n=1: [1, -1, -3, -14, -85, -504, -4424, -6796, -878157, ...];
n=2: [1, -1, -2, -6, -28, -70, -1446, 22302, -855032, ...];
n=3: [1, -1, -1, -1, -7, 49, -1191, 31569, -1051695, ...];
n=4: [1, -1, 0, 1, -6, 78, -1544, 40605, -1328178, ...];
n=5: [1, -1, 1, 0, -9, 117, -2118, 53232, -1699905, ...];
n=6: [1, -1, 2, -4, 0, 141, -2958, 71900, -2216860, ...];
n=7: [1, -1, 3, -11, 37, 0, -3245, 95286, -2941059, ...];
n=8: [1, -1, 4, -21, 118, -581, 0, 99086, -3760182, ...];
n=9: [1, -1, 5, -34, 259, -2002, 13212, 0, -3775221, ...];
n=10: [1, -1, 6, -50, 476, -4788, 47578, -397090, 0, ...]; ...
in which the secondary diagonal is all zeroes after the initial terms.
(4) The table of coefficients in 1/A(x)^n begins:
n=1: [1, -1, -4, -25, -194, -1603, -15264, -122316, ...];
n=2: [1, -2, -7, -42, -322, -2618, -25145, -191580, ...];
n=3: [1, -3, -9, -52, -396, -3168, -30889, -220332, ...];
n=4: [1, -4, -10, -56, -427, -3360, -33546, -218880, ...];
n=5: [1, -5, -10, -55, -425, -3286, -33990, -195585, ...];
n=6: [1, -6, -9, -50, -399, -3024, -32938, -157122, ...];
n=7: [1, -7, -7, -42, -357, -2639, -30968, -108718, ...];
n=8: [1, -8, -4, -32, -306, -2184, -28536, -54368, ...];
n=9: [1, -9, 0, -21, -252, -1701, -25992, 2970, -7903413, ...]; ...
where the main diagonal divided by n begins:
D = [1, -2/2, -9/3, -56/4, -425/5, -3024/6, -30968/7, -54368/8, -7903413/9, ...],
D = [1, -1, -3, -14, -85, -504, -4424, -6796, -878157, 25703710, -1270518018, ...].
Compare D to:
A230218 = [1, 1, -3, 14, -85, 504, -4424, 6796, -878157, -25703710, -1270518018, ...];
the g.f. G(x) of A230218 obeys: [x^n] G(x)^(n^2+n+1) = 0 for n>1.
PROG
(PARI) {a(n) = my(A=[1, 1]); for(i=2, n+1, A=concat(A, 0); A[#A]=Vec(1/Ser(A)^((#A)^2))[#A]/(#A)^2 ); A[n+1]}
for(n=0, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
sign
AUTHOR
Paul D. Hanna, Sep 25 2017
STATUS
approved