OFFSET
0,3
COMMENTS
This sequence is conjectured to consist entirely of integers.
EXAMPLE
G.f.: A(x) = -1 + x + 6*x^2 + 76*x^3 + 74*x^4 + 26058*x^5 - 7794632*x^6 - 14388774920*x^7 + 213380711849598*x^8 + 28487519401892735386*x^9 + ...
The table of coefficients of x^k in A(x)^(n^2), k >= 0, begins:
n=1: [-1, 1, 6, 76, 74, 26058, -7794632, ...];
n=2: [1, -4, -18, -236, 761, -99192, 31519032, ...];
n=3: [-1, 9, 18, 336, -4716, 221688, -72062976, ...];
n=4: [1, -16, 24, -336, 13116, -438576, 130813640, ...];
n=5: [-1, 25, -150, 600, -25800, 855780, -210360800, ...];
n=6: [1, -36, 414, -2316, 46161, -1632600, 316244424, ...];
n=7: [-1, 49, -882, 8036, -97706, 3042606, -459537288, ...]; ...
in which, by definition, the following sums along the columns equal zero:
0 = (6)/1! + (-18)/2! + (18)/3! ;
0 = (76)/1! + (-236)/2! + (336)/3! + (-336)/4! ;
0 = (74)/1! + (761)/2! + (-4716)/3! + (13116)/4! + (-25800)/5! ;
0 = (26058)/1! + (-99192)/2! + (221688)/3! + (-438576)/4! + (855780)/5! + (-1632600)/6! ;
0 = (-7794632)/1! + (31519032)/2! + (-72062976)/3! + (130813640)/4! + (-210360800)/5! + (316244424)/6! + (-459537288)/7! ; ...
One may continue the above pattern to determine all the terms of this sequence.
PROG
(PARI) {a(n) = my(A=[-1, 1, 0]); for(i=1, n, A=concat(A, 0);
A[#A-1] = (-1)^(#A-1) * (#A-2)! * polcoeff( sum(m=0, #A-1, Ser(A)^(m^2) / m!) , #A-2) ); A[n+1]}
for(n=0, 16, print1(a(n), ", "))
CROSSREFS
KEYWORD
sign
AUTHOR
Paul D. Hanna, Dec 29 2021
STATUS
approved