OFFSET
0,2
COMMENTS
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..200
EXAMPLE
O.g.f.: A(x) = 1 + 3*x + 69*x^2 + 5005*x^3 + 806148*x^4 + 239220375*x^5 +...
such that the coefficients satisfy:
1 = 1*(1-x)^3 + 3*x*(1-2^3*x+x^2)^3 + 69*x^2*(1-3^3*x+3^3*x^2-x^3)^3 + 5005*x^3*(1-4^3*x+6^3*x^2-4^3*x^3+x^4)^3 + 806148*x^4*(1-5^3*x+10^3*x^2-10^3*x^3+5^3*x^4-x^5)^3 +...
Note that the cube-root of the o.g.f., A(x)^(1/3), is an integer series:
A(x)^(1/3) = 1 + x + 22*x^2 + 1624*x^3 + 264962*x^4 + 79136637*x^5 + 38671111558*x^6 + 28642761340956*x^7 + 30413158977739302*x^8 +...+ A245658(n)*x^n +...
CONJECTURE: given
G(x,m,j) = Sum_{n>=0} a(n) * (m*x)^n * [ Sum_{k=0..n+1} C(n+1, k)^3 * (j*x)^k ]^3
then G(x,m,j)^(1/3) is an integer series in x whenever m, j, are integers.
OBSERVATIONS.
The terms of this sequence are congruent to 2 modulo 3 at positions:
[12, 36, 93, 108, 111, 114, 120, 174, 255, 279, ...].
The terms of A245658 are congruent to 2 modulo 3 at positions:
[4, 12, 31, 36, 37, 38, 40, 58, 85, 93, ...].
PROG
(PARI) {a(n)=if(n==0, 1, -polcoeff(sum(m=0, n-1, a(m)*x^m*sum(k=0, m+1, binomial(m+1, k)^3 * (-x)^k +x*O(x^n) )^3 ), n))}
for(n=0, 20, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Sep 09 2014
STATUS
approved