login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A099429 A Jacobsthal-Lucas convolution. 2
0, 0, 2, 3, 12, 25, 66, 147, 344, 765, 1710, 3751, 8196, 17745, 38234, 81915, 174768, 371365, 786438, 1660239, 3495260, 7340025, 15379122, 32156323, 67108872, 139810125, 290805086, 603979767, 1252698804, 2594876065, 5368709130, 11095332171, 22906492256 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

COMMENTS

Inverse binomial transform is (-1)^n * a(n). - Michael Somos, Jun 02 2014

If we concatenate the lexicographically ordered bit strings of length n, then a(n) is the number of times 11 appears as a substring, if overlapping substrings are not considered as being separate. - John M. Campbell, Jan 18 2019

LINKS

G. C. Greubel, Table of n, a(n) for n = 0..1000

Index entries for linear recurrences with constant coefficients, signature (2,3,-4,-4).

FORMULA

G.f.: x^2*(2-x)/(1-x-2*x^2)^2. [Typo corrected by Colin Barker, Jun 16 2012]

a(n) = Sum_{k=0..n} J(n-k)*(2^(k-1) -(-1)^k +0^k/2).

a(n) = Sum_{k=0..n+1} J(n-k)*binomial(n-k+1, k)*binomial(1, (k+1)/2)*(1-(-1)^k)/2.

a(n) = A036289(n)/6 +(-1)^n*n/3. - R. J. Mathar, Sep 21 2012

a(-n) = (-2)^(-n-1) * A193449(n) for all n in Z. - Michael Somos, Jun 02 2014

EXAMPLE

G.f. = 2*x^2 + 3*x^3 + 12*x^4 + 25*x^5 + 66*x^6 + 147*x^7 + 344*x^8 + ...

If we concatenate the lexicographically ordered bit strings of length 4, we obtain the expression 0000000100100011010001010110011110001001101010111100110111101111, and we see that the substring 11 appears a total of a(4) = 12 times, with overlapping substrings not being considered as being separate. - John M. Campbell, Jan 18 2019

MATHEMATICA

CoefficientList[Series[x^2*(2-x)/(1-x-2x^2)^2, {x, 0, 32}], x] (* Michael De Vlieger, Jan 18 2019 *)

PROG

(PARI) {a(n) = if( n>=0, polcoeff( x^2*(2-x)/((1+x)*(1-2*x))^2 + x*O(x^n), n), polcoeff( x*(1-2*x)/((1+x)*(2-x))^2 + x*O(x^-n), -n) )}; /* Michael Somos, Jun 02 2014 */

(MAGMA) m:=40; R<x>:=PowerSeriesRing(Integers(), m); [0, 0] cat Coefficients(R!( x^2*(2-x)/(1-x-2*x^2)^2 )); // G. C. Greubel, Feb 25 2019

(Sage) (x^2*(2-x)/(1-x-2*x^2)^2).series(x, 40).coefficients(x, sparse=False) # G. C. Greubel, Feb 25 2019

CROSSREFS

Cf. A001045, A014551, A193449.

Sequence in context: A130337 A293697 A084728 * A049601 A148060 A148061

Adjacent sequences:  A099426 A099427 A099428 * A099430 A099431 A099432

KEYWORD

easy,nonn

AUTHOR

Paul Barry, Oct 15 2004

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified July 30 04:45 EDT 2021. Contains 346348 sequences. (Running on oeis4.)