cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A100223 G.f. A(x) satisfies: 2^n - 1 = Sum_{k=0..n} [x^k]A(x)^n and also satisfies: (2+z)^n - (1+z)^n + z^n = Sum_{k=0..n} [x^k](A(x)+z*x)^n for all z, where [x^k]A(x)^n denotes the coefficient of x^k in A(x)^n.

This page as a plain text file.
%I A100223 #46 Sep 24 2021 13:40:47
%S A100223 1,0,1,1,0,-2,-3,1,11,15,-13,-77,-86,144,595,495,-1520,-4810,-2485,
%T A100223 15675,39560,6290,-159105,-324805,87075,1592843,2616757,-2136539,
%U A100223 -15726114,-20247800,32296693,152909577,145139491,-417959049,-1460704685,-885536173,4997618808,13658704994,3223741399
%N A100223 G.f. A(x) satisfies: 2^n - 1 = Sum_{k=0..n} [x^k]A(x)^n and also satisfies: (2+z)^n - (1+z)^n + z^n = Sum_{k=0..n} [x^k](A(x)+z*x)^n for all z, where [x^k]A(x)^n denotes the coefficient of x^k in A(x)^n.
%C A100223 The current sequence and A007440, A214649, A343773 form a cluster. The base of the cluster is A343773. We get the current sequence if we shift all the terms in A343773 to the right by 2 and take a(0) = 1. - _Gennady Eremin_, May 15 2021
%H A100223 Gennady Eremin, <a href="/A100223/b100223.txt">Table of n, a(n) for n = 0..800</a>
%H A100223 Gennady Eremin, <a href="https://arxiv.org/abs/2108.10676">Walking in the OEIS: From Motzkin numbers to Fibonacci numbers. The "shadows" of Motzkin numbers</a>, arXiv:2108.10676 [math.CO], 2021.
%F A100223 G.f.: A(x) = (1+x+sqrt(1-2*x+5*x^2))/2.
%F A100223 G.f.: A(x) = x/(series_reversion[x*(1-x)/(1-x-x^2)]).
%F A100223 D-finite with recurrence a(n) = ((2*n-3)*a(n-1) - 5*(n-3)*a(n-2))/n (for n>2), with a(0)=1, a(1)=0, a(2)=1.
%F A100223 Given g.f. A(x), then B(x) = A(x) - 1 + x series reversion is -B(-x). - _Michael Somos_, Sep 07 2005
%F A100223 Given g.f. A(x) and C(x) = reversion of x + x^2, then B(x) = A(x) - 1 + x satisfies B(x) = x + C(x * B(x)). - _Michael Somos_, Sep 07 2005
%F A100223 From _Paul Barry_, Nov 07 2010: (Start)
%F A100223 a(n+1) has Hankel transform (-1)^n*F(n).
%F A100223 a(n+2) has Hankel transform (-1)^comb(n+1,2).
%F A100223 a(n+3) has Hankel transform (-1)^comb(n+1,2)*F(n+2).
%F A100223 Hankel transform of a(n+4) is sum{k=0..n, (-1)^(n-k+1)*(n-k+1)(F(2k+2)-C(1,k)-C(0,k))}. (End)
%F A100223 G.f.: A(x) = x + G(0) = 1 + x^2/G(0) where G(k) = 1 - x + x^2/G(k+1); (continued fraction). - _Sergei N. Gladkovskii_, Dec 25 2011
%F A100223 Lim sup n->infinity |a(n)|^(1/n) = sqrt(5).
%F A100223 From _Gennady Eremin_, May 15 2021: (Start)
%F A100223 a(n+2) = A343773(n), n >= 0.
%F A100223 G.f.: A(x) = 1 + x^2*B(x), where B(x) is the g.f. of A343773.
%F A100223 Lim_{n->infinity} a(n)/A001006(n) = 0. (End)
%e A100223 G.f. = 1 + x^2 + x^3 - 2*x^5 - 3*x^6 + x^7 + 11*x^8 + 15*x^9 - 13*x^10 - 77*x^11 - 86*x^12 + 144*x^13 + ...
%e A100223 From the table of powers of A(x) (A100224), we see that
%e A100223 2^n-1 = Sum of coefficients [x^0] through [x^n] in A(x)^n:
%e A100223 A^1=[1,0],1,1,0,-2,-3,1,11,...
%e A100223 A^2=[1,0,2],2,1,-2,-5,-2,12,...
%e A100223 A^3=[1,0,3,3],3,0,-5,-6,6,...
%e A100223 A^4=[1,0,4,4,6],4,-2,-8,-3,...
%e A100223 A^5=[1,0,5,5,10,10],5,-5,-10,...
%e A100223 A^6=[1,0,6,6,15,18,17],6,-9,...
%e A100223 A^7=[1,0,7,7,21,28,35,28],7,...
%e A100223 A^8=[1,0,8,8,28,40,60,64,46],...
%e A100223 the main diagonal of which is A001610 = [0,2,3,6,10,17,...],
%e A100223 where Sum_{n>=1} A001610(n-1)/n*x^n = log((1-x)/(1-x-x^2)).
%t A100223 CoefficientList[Series[(1+x+Sqrt[1-2*x+5*x^2])/2, {x, 0, 30}], x] (* _Vaclav Kotesovec_, Feb 11 2015 *)
%o A100223 (PARI) {a(n) = if( n<=0, n==0, (2^n - 1 - sum( k=0, n, polcoeff( sum( j=0, min(k, n-1), a(j) * x^j)^n + x * O(x^k), k))) / n)}
%o A100223 (PARI) {a(n) = if( n<=0, n==0, if( n==1, 0, if( n==2, 1, ((2*n-3) * a(n-1) - 5 * (n-3) * a(n-2)) / n)))}
%o A100223 (PARI) {a(n) = if( n<0, 0, polcoeff((1 + x + sqrt(1 - 2*x + 5*x^2 + x * O(x^n))) / 2, n))}
%o A100223 (Python)
%o A100223 A100223 = [1, 0, 1]
%o A100223 for n in range(3, 801):
%o A100223     A100223.append( ((2*n-3)*A100223[-1]
%o A100223       - 5*(n-3)*A100223[-2])//n )  # _Gennady Eremin_, Apr 24 2021
%Y A100223 Cf. A001006, A001610, A007440, A100224, A100225, A214649, A343773.
%K A100223 sign
%O A100223 0,6
%A A100223 _Paul D. Hanna_, Nov 28 2004
%E A100223 Entry revised by _Paul D. Hanna_, Mar 19 2013