|
PROG
|
(PARI) /* By definition: */
{a(n)=local(A=1); A=sum(m=0, n, x^m*sum(k=0, m, binomial(m, k)^2*x^(5*k)) +x*O(x^n)); polcoeff(A, n)}
for(n=0, 50, print1(a(n), ", "))
(PARI) /* From closed formula: */
{a(n)=local(A=1); A= 1/sqrt((1 - x + x^6)^2 - 4*x^6 +x*O(x^n)); polcoeff(A, n)}
for(n=0, 50, print1(a(n), ", "))
(PARI) /* From a series identity: */
{a(n)=local(A=1+x); A=sum(m=0, n, (2*m)!/(m!)^2 * x^(6*m) / (1 - x + x^6 +x*O(x^n))^(2*m+1)); polcoeff(A, n)}
for(n=0, 40, print1(a(n), ", "))
(PARI) /* From a binomial series identity: */
{a(n)=local(A=1+x); A=sum(m=0, n, x^m*(1-x^5)^(2*m+1)*sum(k=0, n, binomial(m+k, k)^2*x^(5*k)) +x*O(x^n)); polcoeff(A, n)}
for(n=0, 50, print1(a(n), ", "))
(PARI) /* From a binomial series identity: */
{a(n)=local(A=1+x); A=sum(m=0, n\6, x^(6*m)*sum(k=0, n-5*m, binomial(m+k, k)^2*x^k) +x*O(x^n)); polcoeff(A, n)}
for(n=0, 50, print1(a(n), ", "))
(PARI) /* From a binomial series identity: */
{a(n)=local(A=1+x); A=sum(m=0, n\6, x^(6*m) * sum(k=0, m, binomial(m, k)^2*x^k) / (1-x +x*O(x^n))^(2*m+1) ); polcoeff(A, n)}
for(n=0, 40, print1(a(n), ", "))
(PARI) /* From exponential formula: */
{a(n)=local(A=1); A=exp(sum(m=1, n, sum(k=0, m, binomial(2*m, 2*k)*x^(5*k)) * x^m/m) +x*O(x^n)); polcoeff(A, n)}
for(n=0, 50, print1(a(n), ", "))
(PARI) /* From formula for a(n): */
{a(n)=sum(k=0, n\5, binomial(n-5*k, k)^2)}
for(n=0, 50, print1(a(n), ", "))
|