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.

Showing 1-10 of 15 results. Next

A192232 Constant term of the reduction of n-th Fibonacci polynomial by x^2 -> x+1. (See Comments.)

Original entry on oeis.org

1, 0, 2, 1, 6, 7, 22, 36, 89, 168, 377, 756, 1630, 3353, 7110, 14783, 31130, 65016, 136513, 285648, 599041, 1254456, 2629418, 5508097, 11542854, 24183271, 50674318, 106173180, 222470009, 466131960, 976694489, 2046447180, 4287928678, 8984443769, 18825088134
Offset: 1

Views

Author

Clark Kimberling, Jun 26 2011

Keywords

Comments

Polynomial reduction: an introduction
...
We begin with an example. Suppose that p(x) is a polynomial, so that p(x)=(x^2)t(x)+r(x) for some polynomials t(x) and r(x), where r(x) has degree 0 or 1. Replace x^2 by x+1 to get (x+1)t(x)+r(x), which is (x^2)u(x)+v(x) for some u(x) and v(x), where v(x) has degree 0 or 1. Continuing in this manner results in a fixed polynomial w(x) of degree 0 or 1. If p(x)=x^n, then w(x)=x*F(n)+F(n-1), where F=A000045, the sequence of Fibonacci numbers.
In order to generalize, write d(g) for the degree of an arbitrary polynomial g(x), and suppose that p, q, s are polynomials satisfying d(s)s in this manner until reaching w such that d(w)s.
The coefficients of (reduction of p by q->s) comprise a vector of length d(q)-1, so that a sequence p(n,x) of polynomials begets a sequence of vectors, such as (F(n), F(n-1)) in the above example. We are interested in the component sequences (e.g., F(n-1) and F(n)) for various choices of p(n,x).
Following are examples of reduction by x^2->x+1:
n-th Fibonacci p(x) -> A192232+x*A112576
n-th cyclotomic p(x) -> A192233+x*A051258
n-th 1st-kind Chebyshev p(x) -> A192234+x*A071101
n-th 2nd-kind Chebyshev p(x) -> A192235+x*A192236
x(x+1)(x+2)...(x+n-1) -> A192238+x*A192239
(x+1)^n -> A001519+x*A001906
(x^2+x+1)^n -> A154626+x*A087635
(x+2)^n -> A020876+x*A030191
(x+3)^n -> A192240+x*A099453
...
Suppose that b=(b(0), b(1),...) is a sequence, and let p(n,x)=b(0)+b(1)x+b(2)x^2+...+b(n)x^n. We define (reduction of sequence b by q->s) to be the vector given by (reduction of p(n,x) by q->s), with components in the order of powers, from 0 up to d(q)-1. For k=0,1,...,d(q)-1, we then have the "k-sequence of (reduction of sequence b by q->s)". Continuing the example, if b is the sequence given by b(k)=1 if k=n and b(k)=0 otherwise, then the 0-sequence of (reduction of b by x^2->x+1) is (F(n-1)), and the 1-sequence is (F(n)).
...
For selected sequences b, here are the 0-sequences and 1-sequences of (reduction of b by x^2->x+1):
b=A000045, Fibonacci sequence (1,1,2,3,5,8,...) yields
0-sequence A166536 and 1-sequence A064831.
b=(1,A000045)=(1,1,1,2,3,5,8,...) yields
0-sequence A166516 and 1-sequence A001654.
b=A000027, natural number sequence (1,2,3,4,...) yields
0-sequence A190062 and 1-sequence A122491.
b=A000032, Lucas sequence (1,3,4,7,11,...) yields
0-sequence A192243 and 1-sequence A192068.
b=A000217, triangular sequence (1,3,6,10,...) yields
0-sequence A192244 and 1-sequence A192245.
b=A000290, squares sequence (1,4,9,16,...) yields
0-sequence A192254 and 1-sequence A192255.
More examples: A192245-A192257.
...
More comments:
(1) If s(n,x)=(reduction of x^n by q->s) and
p(x)=p(0)x^n+p(1)x^(n-1)+...+p(n)x^0, then
(reduction of p by q->s)=p(0)s(n,x)+p(1)s(n-1,x)
+...+p(n-1)s(1,x)+p(n)s(0,x). See A192744.
(2) For any polynomial p(x), let P(x)=(reduction of p(x)
by q->s). Then P(r)=p(r) for each zero r of
q(x)-s(x). In particular, if q(x)=x^2 and s(x)=x+1,
then P(r)=p(r) if r=(1+sqrt(5))/2 (golden ratio) or
r=(1-sqrt(5))/2.

Examples

			The first four Fibonacci polynomials and their reductions by x^2->x+1 are shown here:
F1(x)=1 -> 1 + 0x
F2(x)=x -> 0 + 1x
F3(x)=x^2+1 -> 2+1x
F4(x)=x^3+2x -> 1+4x
F5(x)=x^4+3x^2+1 -> (x+1)^2+3(x+1)+1 -> 6+6x.
From these, read A192232=(1,0,1,1,6,...) and A112576=(0,1,1,4,6,...).
		

Crossrefs

Programs

  • Mathematica
    q[x_] := x + 1;
    reductionRules = {x^y_?EvenQ -> q[x]^(y/2),  x^y_?OddQ -> x q[x]^((y - 1)/2)};
    t = Table[FixedPoint[Expand[#1 /. reductionRules] &, Fibonacci[n, x]], {n, 1, 40}];
    Table[Coefficient[Part[t, n], x, 0], {n, 1, 40}]
      (* A192232 *)
    Table[Coefficient[Part[t, n], x, 1], {n, 1, 40}]
    (* A112576 *)
    (* Peter J. C. Moses, Jun 25 2011 *)
    LinearRecurrence[{1, 3, -1, -1}, {1, 0, 2, 1}, 60] (* Vladimir Joseph Stephan Orlovsky, Feb 08 2012 *)
  • PARI
    Vec((1-x-x^2)/(1-x-3*x^2+x^3+x^4)+O(x^99)) \\ Charles R Greathouse IV, Jan 08 2013

Formula

Empirical G.f.: -x*(x^2+x-1)/(x^4+x^3-3*x^2-x+1). - Colin Barker, Sep 11 2012
The above formula is correct. - Charles R Greathouse IV, Jan 08 2013
a(n) = A265752(A206296(n)). - Antti Karttunen, Dec 15 2015
a(n) = A112576(n) -A112576(n-1) -A112576(n-2). - R. J. Mathar, Dec 16 2015

Extensions

Example corrected by Clark Kimberling, Dec 18 2017

A005248 Bisection of Lucas numbers: a(n) = L(2*n) = A000032(2*n).

Original entry on oeis.org

2, 3, 7, 18, 47, 123, 322, 843, 2207, 5778, 15127, 39603, 103682, 271443, 710647, 1860498, 4870847, 12752043, 33385282, 87403803, 228826127, 599074578, 1568397607, 4106118243, 10749957122, 28143753123, 73681302247, 192900153618, 505019158607, 1322157322203
Offset: 0

Views

Author

Keywords

Comments

Drop initial 2; then iterates of A050411 do not diverge for these starting values. - David W. Wilson
All nonnegative integer solutions of Pell equation a(n)^2 - 5*b(n)^2 = +4 together with b(n)=A001906(n), n>=0. - Wolfdieter Lang, Aug 31 2004
a(n+1) = B^(n)AB(1), n>=0, with compositions of Wythoff's complementary A(n):=A000201(n) and B(n)=A001950(n) sequences. See the W. Lang link under A135817 for the Wythoff representation of numbers (with A as 1 and B as 0 and the argument 1 omitted). E.g., 3=`10`, 7=`010`, 18=`0010`, 47=`00010`, ..., in Wythoff code. a(0) = 2 = B(1) in Wythoff code.
Output of Tesler's formula (as well as that of Lu and Wu) for the number of perfect matchings of an m X n Möbius band where m and n are both even specializes to this sequence for m=2. - Sarah-Marie Belcastro, Jul 04 2009
Numbers having two 1's in their base-phi representation. - Robert G. Wilson v, Sep 13 2010
Pisano period lengths: 1, 3, 4, 3, 2, 12, 8, 6, 12, 6, 5, 12, 14, 24, 4, 12, 18, 12, 9, 6, ... - R. J. Mathar, Aug 10 2012
From Wolfdieter Lang, Feb 18 2013: (Start)
a(n) is also one half of the total number of round trips, each of length 2*n, on the graph P_4 (o-o-o-o) (the simple path with 4 points (vertices) and 3 lines (or edges)). See the array and triangle A198632 for the general case of the graph P_N (there N is n and the length is l=2*k).
O.g.f. for w(4,l) (with zeros for odd l): y*(d/dy)S(4,y)/S(4,y) with y=1/x and Chebyshev S-polynomials (coefficients A049310). See also A198632 for a rewritten form. One half of this o.g.f. for x -> sqrt(x) produces the g.f. (2-3x)/(1-3x+x^2) given below. (End)
Solutions (x, y) = (a(n), a(n+1)) satisfying x^2 + y^2 = 3xy - 5. - Michel Lagneau, Feb 01 2014
Except for the first term, positive values of x (or y) satisfying x^2 - 7xy + y^2 + 45 = 0. - Colin Barker, Feb 16 2014
Except for the first term, positive values of x (or y) satisfying x^2 - 18xy + y^2 + 320 = 0. - Colin Barker, Feb 16 2014
a(n) are the numbers such that a(n)^2-2 are Lucas numbers. - Michel Lagneau, Jul 22 2014
All sequences of this form, b(n+1) = 3*b(n) - b(n-1), regardless of initial values, which includes this sequence, yield this sequence as follows: a(n) = (b(j+n) + b(j-n))/b(j), for any j, except where b(j) = 0. Also note formula below relating this a(n) to all sequences of the form G(n+1) = G(n) + G(n-1). - Richard R. Forberg, Nov 18 2014
A non-simple continued fraction expansion for F(2n*(k+1))/F(2nk) k>=1 is a(n) + (-1)/(a(n) + (-1)/(a(n) + ... + (-1)/a(n))) where a(n) appears exactly k times (F(n) denotes the n-th Fibonacci number). E.g., F(16)/F(12) equals 7 + (-1)/(7 + (-1)/7). Furthermore, these a(n) are exactly the positive integers k such that the non-simple infinite continued fraction k + (-1)/(k + (-1)/(k + (-1)/(k + ...))) belongs to Q(sqrt(5)). Compare to Benoit Cloitre and Thomas Baruchel's comments at A002878. - Greg Dresden, Aug 13 2019
For n >= 1, a(n) is the number of cyclic up-down words of length 2*n over an alphabet of size 3. - Sela Fried, Apr 08 2025

Examples

			G.f. = 2 + 3*x + 7*x^2 + 18*x^3 + 47*x^4 + 123*x^5 + 322*x^6 + 843*x^7 + ... - _Michael Somos_, Aug 11 2009
		

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • Richard P. Stanley, Enumerative combinatorics, Vol. 2. Volume 62 of Cambridge Studies in Advanced Mathematics. Cambridge University Press, Cambridge, 1999.

Crossrefs

Cf. A000032, A002878 (odd-indexed Lucas numbers), A001906 (Chebyshev S(n-1, 3)), a(n) = sqrt(4+5*A001906(n)^2), A228842.
a(n) = A005592(n)+1 = A004146(n)+2 = A065034(n)-1.
First differences of A002878. Pairwise sums of A001519. First row of array A103997.
Cf. A153415, A201157. Also Lucas(k*n): A000032 (k = 1), A014448 (k = 3), A056854 (k = 4), A001946 (k = 5), A087215 (k = 6), A087281 (k = 7), A087265 (k = 8), A087287 (k = 9), A065705 (k = 10), A089772 (k = 11), A089775 (k = 12).

Programs

  • Haskell
    a005248 n = a005248_list !! n
    a005248_list = zipWith (+) (tail a001519_list) a001519_list
    -- Reinhard Zumkeller, Jan 11 2012
  • Magma
    [Lucas(2*n) : n in [0..100]]; // Vincenzo Librandi, Apr 14 2011
    
  • Maple
    a:= n-> (<<2|3>>. <<3|1>, <-1|0>>^n)[1$2]: seq(a(n), n=0..30); # Alois P. Heinz, Jul 31 2008
    with(combinat): seq(5*fibonacci(n)^2+2*(-1)^n, n= 0..26);
  • Mathematica
    a[0] = 2; a[1] = 3; a[n_] := 3a[n - 1] - a[n - 2]; Table[ a[n], {n, 0, 27}] (* Robert G. Wilson v, Jan 30 2004 *)
    Fibonacci[1 + 2n] + 1/2 (-Fibonacci[2n] + LucasL[2n]) (* Tesler. Sarah-Marie Belcastro, Jul 04 2009 *)
    LinearRecurrence[{3, -1}, {2, 3}, 50] (* Sture Sjöstedt, Nov 27 2011 *)
    LucasL[Range[0,60,2]] (* Harvey P. Dale, Sep 30 2014 *)
  • PARI
    {a(n) = fibonacci(2*n + 1) + fibonacci(2*n - 1)}; /* Michael Somos, Jun 23 2002 */
    
  • PARI
    {a(n) = 2 * subst( poltchebi(n), x, 3/2)}; /* Michael Somos, Jun 28 2003 */
    
  • Sage
    [lucas_number2(n,3,1) for n in range(37)] # Zerinvary Lajos, Jun 25 2008
    

Formula

a(n) = Fibonacci(2*n-1) + Fibonacci(2*n+1).
G.f.: (2-3*x)/(1-3*x+x^2). - Simon Plouffe in his 1992 dissertation.
a(n) = S(n, 3) - S(n-2, 3) = 2*T(n, 3/2) with S(n-1, 3) = A001906(n) and S(-2, x) = -1. U(n, x)=S(n, 2*x) and T(n, x) are Chebyshev's U- and T-polynomials.
a(n) = a(k)*a(n - k) - a(n - 2k) for all k, i.e., a(n) = 2*a(n) - a(n) = 3*a(n - 1) - a(n - 2) = 7*a(n - 2) - a(n - 4) = 18*a(n - 3) - a(n - 6) = 47*a(n - 4) - a(n - 8) etc., a(2n) = a(n)^2 - 2. - Henry Bottomley, May 08 2001
a(n) = A060924(n-1, 0) = 3*A001906(n) - 2*A001906(n-1), n >= 1. - Wolfdieter Lang, Apr 26 2001
a(n) ~ phi^(2*n) where phi=(1+sqrt(5))/2. - Joe Keane (jgk(AT)jgk.org), May 15 2002
a(0)=2, a(1)=3, a(n) = 3*a(n-1) - a(n-2) = a(-n). - Michael Somos, Jun 28 2003
a(n) = phi^(2*n) + phi^(-2*n) where phi=(sqrt(5)+1)/2, the golden ratio. E.g., a(4)=47 because phi^(8) + phi^(-8) = 47. - Dennis P. Walsh, Jul 24 2003
With interpolated zeros, trace(A^n)/4, where A is the adjacency matrix of path graph P_4. Binomial transform is then A049680. - Paul Barry, Apr 24 2004
a(n) = (floor((3+sqrt(5))^n) + 1)/2^n. - Lekraj Beedassy, Oct 22 2004
a(n) = ((3-sqrt(5))^n + (3+sqrt(5))^n)/2^n (Note: substituting the number 1 for 3 in the last equation gives A000204, substituting 5 for 3 gives A020876). - Creighton Dement, Apr 19 2005
a(n) = (1/(n+1/2))*Sum_{k=0..n} B(2k)*L(2n+1-2k)*binomial(2n+1, 2k) where B(2k) is the (2k)-th Bernoulli number. - Benoit Cloitre, Nov 02 2005
a(n) = term (1,1) in the 1 X 2 matrix [2,3] . [3,1; -1,0]^n. - Alois P. Heinz, Jul 31 2008
a(n) = 2*cosh(2*n*psi), where psi=log((1+sqrt(5))/2). - Al Hakanson, Mar 21 2009
From Sarah-Marie Belcastro, Jul 04 2009: (Start)
a(n) - (a(n) - F(2n))/2 - F(2n+1) = 0. (Tesler)
Product_{r=1..n} (1 + 4*(sin((4r-1)*Pi/(4n)))^2). (Lu/Wu) (End)
a(n) = Fibonacci(2n+6) mod Fibonacci(2n+2), n > 1. - Gary Detlefs, Nov 22 2010
a(n) = 5*Fibonacci(n)^2 + 2*(-1)^n. - Gary Detlefs, Nov 22 2010
a(n) = A033888(n)/A001906(n), n > 0. - Gary Detlefs, Dec 26 2010
a(n) = 2^(2*n) * Sum_{k=1..2} (cos(k*Pi/5))^(2*n). - L. Edson Jeffery, Jan 21 2012
From Peter Bala, Jan 04 2013: (Start)
Let F(x) = Product_{n>=0} (1 + x^(4*n+1))/(1 + x^(4*n+3)). Let alpha = 1/2*(3 - sqrt(5)). This sequence gives the simple continued fraction expansion of 1 + F(alpha) = 2.31829 56058 81914 31334 ... = 2 + 1/(3 + 1/(7 + 1/(18 + ...))).
Also F(-alpha) = 0.64985 97768 07374 32950 has the continued fraction representation 1 - 1/(3 - 1/(7 - 1/(18 - ...))) and the simple continued fraction expansion 1/(1 + 1/((3-2) + 1/(1 + 1/((7-2) + 1/(1 + 1/((18-2) + 1/(1 + ...))))))).
F(alpha)*F(-alpha) has the simple continued fraction expansion 1/(1 + 1/((3^2-4) + 1/(1 + 1/((7^2-4) + 1/(1 + 1/((18^2-4) + 1/(1 + ...))))))).
Added Oct 13 2019: 1/2 + 1/2*F(alpha)/F(-alpha) = 1.5142923542... has the simple continued fraction expansion 1 + 1/((3 - 2) + 1/(1 + 1/((18 - 2) + 1/(1 + 1/(123 - 2) + 1/(1 + ...))))). (End)
G.f.: (W(0)+6)/(5*x), where W(k) = 5*x*k + x - 6 + 6*x*(5*k-9)/W(k+1) (continued fraction). - Sergei N. Gladkovskii, Aug 19 2013
Sum_{n >= 1} 1/( a(n) - 5/a(n) ) = 1. Compare with A001906, A002878 and A023039. - Peter Bala, Nov 29 2013
0 = a(n) * a(n+2) - a(n+1)^2 - 5 for all n in Z. - Michael Somos, Aug 24 2014
a(n) = (G(j+2n) + G(j-2n))/G(j), for n >= 0 and any j, positive or negative, except where G(j) = 0, and for any sequence of the form G(n+1) = G(n) + G(n-1) with any initial values for G(0), G(1), including non-integer values. G(n) includes Lucas, Fibonacci. Compare with A081067 for odd number offsets from j. - Richard R. Forberg, Nov 16 2014
a(n) = [x^n] ( (1 + 3*x + sqrt(1 + 6*x + 5*x^2))/2 )^n for n >= 1. - Peter Bala, Jun 23 2015
From J. M. Bergot, Oct 28 2015: (Start)
For n>0, a(n) = F(n-1) * L(n) + F(2*n+1) - (-1)^n with F(k) = A000045(k).
For n>1, a(n) = F(n+1) * L(n) + F(2*n-1) - (-1)^n.
For n>2, a(n) = 5*F(2*n-3) + 2*L(n-3) * L(n) + 8*(-1)^n. (End)
For n>1, a(n) = L(n-2)*L(n+2) -7*(-1)^n. - J. M. Bergot, Feb 10 2016
a(n) = 6*F(n-1)*L(n-1) - F(2*n-6) with F(n)=A000045(n) and L(n)=A000032(n). - J. M. Bergot, Apr 21 2017
a(n) = F(2*n) + 2*F(n-1)*L(n) with F(n)=A000045(n) and L(n)=A000032(n). - J. M. Bergot, May 01 2017
E.g.f.: exp(4*x/(1+sqrt(5))^2) + exp((1/4)*(1+sqrt(5))^2*x). - Stefano Spezia, Aug 13 2019
From Peter Bala, Oct 14 2019: (Start)
a(n) = F(2*n+2) - F(2*n-2) = A001906(n+1) - A001906(n-1).
a(n) = trace(M^n), where M is the 2 X 2 matrix [0, 1; 1, 1]^2 = [1, 1; 1, 2].
Consequently the Gauss congruences hold: a(n*p^k) = a(n*p^(k-1)) ( mod p^k ) for all prime p and positive integers n and k. See Zarelua and also Stanley (Ch. 5, Ex. 5.2(a) and its solution).
Sum_{n >= 1} (-1)^(n+1)/( a(n) + 1/a(n) ) = 1/5.
Sum_{n >= 1} (-1)^(n+1)/( a(n) + 3/(a(n) + 2/(a(n))) ) = 1/6.
Sum_{n >= 1} (-1)^(n+1)/( a(n) + 9/(a(n) + 4/(a(n) + 1/(a(n)))) ) = 1/9.
x*exp(Sum_{n >= 1} a(n)*x^/n) = x + 3*x^2 + 8*x^3 + 21*x^4 + ... is the o.g.f. for A001906. (End)
a(n) = n + 2 + Sum_{k=1..n-1} k*a(n-k). - Yu Xiao, May 30 2020
Sum_{n>=1} 1/a(n) = A153415. - Amiram Eldar, Nov 11 2020
Sum_{n>=0} 1/(a(n) + 3) = (2*sqrt(5) + 1)/10 (André-Jeannin, 1991). - Amiram Eldar, Jan 23 2022
a(n) = 2*cosh(2*n*arccsch(2)) = 2*cosh(2*n*asinh(1/2)). - Peter Luschny, May 25 2022
a(n) = (5/2)*(Sum_{k=-n..n} binomial(2*n, n+5*k)) - (1/2)*4^n. - Greg Dresden, Jan 05 2023
a(n) = Sum_{k>=0} Lucas(2*n*k)/(Lucas(2*n)^(k+1)). - Diego Rattaggi, Jan 12 2025

Extensions

Additional comments from Michael Somos, Jun 23 2001

A039717 Row sums of convolution triangle A030523.

Original entry on oeis.org

1, 4, 15, 55, 200, 725, 2625, 9500, 34375, 124375, 450000, 1628125, 5890625, 21312500, 77109375, 278984375, 1009375000, 3651953125, 13212890625, 47804687500, 172958984375, 625771484375, 2264062500000, 8191455078125
Offset: 1

Views

Author

Keywords

Comments

Number of (s(0), s(1), ..., s(2n)) such that 0 < s(i) < 10 and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n, s(0) = 3, s(2n) = 5.
With offset 0 = INVERT transform of A001792: (1, 3, 8, 20, 48, 112, ...). - Gary W. Adamson, Oct 26 2010
From Tom Copeland, Nov 09 2014: (Start)
The array belongs to a family of arrays associated to the Catalan A000108 (t=1), and Riordan, or Motzkin sums A005043 (t=0), with the o.g.f. (1-sqrt(1-4x/(1+(1-t)x)))/2 and inverse x*(1-x)/(1 + (t-1)*x*(1-x)). See A091867 for more info on this family. Here t = -4 (mod signs in the results).
Let C(x) = (1 - sqrt(1-4x))/2, an o.g.f. for the Catalan numbers A000108, with inverse Cinv(x) = x*(1-x) and P(x,t) = x/(1+t*x) with inverse P(x,-t).
O.g.f.: G(x) = x*(1-x)/(1 - 5x*(1-x)) = P(Cinv(x),-5).
Inverse O.g.f.: Ginv(x) = (1 - sqrt(1 - 4*x/(1+5x)))/2 = C(P(x,5)) (signed A026378). Cf. A030528. (End)
p-INVERT of (2^n), where p(s) = 1 - s - s^2; see A289780. - Clark Kimberling, Aug 10 2017

Crossrefs

Programs

  • Mathematica
    CoefficientList[Series[(1 - x) / (1 - 5 x + 5 x^2), {x, 0, 40}], x] (* Vincenzo Librandi, Nov 10 2014 *)
  • PARI
    Vec(x*(1-x)/(1-5*x+5*x^2) + O(x^40)) \\ Altug Alkan, Nov 20 2015

Formula

G.f.: x*(1-x)/(1-5*x+5*x^2) = g1(3, x)/(1-g1(3, x)), g1(3, x) := x*(1-x)/(1-2*x)^2 (g.f. first column of A030523).
From Paul Barry, Apr 16 2004: (Start)
Binomial transform of Fibonacci(2n+2).
a(n) = (sqrt(5)/2 + 5/2)^n*(3*sqrt(5)/10 + 1/2) - (5/2 - sqrt(5)/2)^n*(3*sqrt(5)/10 - 1/2). (End)
a(n) = (1/5)*Sum_{r=1..9} sin(3*r*Pi/10)*sin(r*Pi/2)*(2*cos(r*Pi/10))^(2*n).
a(n) = 5*a(n-1) - 5*a(n-2).
a(n) = Sum_{k=0..n} Sum_{i=0..n} binomial(n, i)*binomial(k+i+1, 2k+1). - Paul Barry, Jun 22 2004
From Johannes W. Meijer, Jul 01 2010: (Start)
Limit_{k->oo} a(n+k)/a(k) = (A020876(n) + A093131(n)*sqrt(5))/2.
Limit_{n->oo} A020876(n)/A093131(n) = sqrt(5).
(End)
From Benito van der Zander, Nov 19 2015: (Start)
Limit_{k->oo} a(k+1)/a(k) = 1 + phi^2 = (5 + sqrt(5)) / 2.
a(n) = a(n-1) * 3 + A081567(n-2) (not proved).
(End)
E.g.f.: exp(x*5/2) * (cosh(x*sqrt(5)/2) + (3/sqrt(5))*sinh(x*sqrt(5)/2)). - Fabian Pereyra, Oct 29 2024

A093129 Binomial transform of Fibonacci(2n-1) (A001519).

Original entry on oeis.org

1, 2, 5, 15, 50, 175, 625, 2250, 8125, 29375, 106250, 384375, 1390625, 5031250, 18203125, 65859375, 238281250, 862109375, 3119140625, 11285156250, 40830078125, 147724609375, 534472656250, 1933740234375, 6996337890625
Offset: 0

Views

Author

Paul Barry, Mar 23 2004

Keywords

Crossrefs

Programs

  • GAP
    a:=[1,2];; for n in [3..30] do a[n]:=5*(a[n-1]-a[n-2]); od; a; # G. C. Greubel, Dec 27 2019
  • Magma
    I:=[1,2]; [n le 2 select I[n] else 5*(Self(n-1) - Self(n-2)): n in [1..30]]; // G. C. Greubel, Dec 27 2019
    
  • Maple
    a:= n-> (<<0|1>, <-5|5>>^n. <<1,2>>)[1,1]:
    seq(a(n), n=0..30);  # Alois P. Heinz, Aug 29 2015
  • Mathematica
    LinearRecurrence[{5, -5}, {1, 2}, 25] (* Jean-François Alcover, May 11 2019 *)
    Table[If[EvenQ[n], 5^(n/2)*Fibonacci[n-1], 5^((n-1)/2)*LucasL[n-1]], {n,0,30}] (* G. C. Greubel, Dec 27 2019 *)
  • PARI
    my(x='x+O('x^30)); Vec((1-3*x)/(1-5*x+5*x^2)) \\ G. C. Greubel, Dec 27 2019
    
  • Sage
    [lucas_number2(n,5,5) for n in range(-1,25)] # Zerinvary Lajos, Jul 08 2008
    

Formula

G.f.: (1-3*x)/(1-5*x+5*x^2).
a(n) = (5-sqrt(5))*((5+sqrt(5))/2)^n/10 + (5+sqrt(5))*((5-sqrt(5))/2)^n/10.
a(n) = A093123(n)/2^n.
a(n) = A020876(n-1). - R. J. Mathar, Sep 05 2008
a(n) = A030191(n) - 3*A030191(n-1). - R. J. Mathar, Jun 29 2012
a(2*n) = 5^n*Fibonacci(2*n-1), a(2*n+1) = 5^n*Lucas(2*n). - G. C. Greubel, Dec 27 2019
E.g.f.: (1/10)*exp((1/2)*(5-sqrt(5))*x)*(5 + sqrt(5) + (5 - sqrt(5))*exp(sqrt(5)*x)). - Stefano Spezia, Dec 28 2019

A093131 Binomial transform of Fibonacci(2n).

Original entry on oeis.org

0, 1, 5, 20, 75, 275, 1000, 3625, 13125, 47500, 171875, 621875, 2250000, 8140625, 29453125, 106562500, 385546875, 1394921875, 5046875000, 18259765625, 66064453125, 239023437500, 864794921875, 3128857421875, 11320312500000, 40957275390625, 148184814453125
Offset: 0

Views

Author

Paul Barry, Mar 23 2004

Keywords

Comments

Second binomial transform of Fibonacci(n). - Paul Barry, Apr 22 2005

Crossrefs

Programs

  • GAP
    a:=[0,1];; for n in [3..30] do a[n]:=5*(a[n-1]-a[n-2]); od; a; # G. C. Greubel, Dec 27 2019
  • Magma
    I:=[0,1]; [n le 2 select I[n] else 5*(Self(n-1) - Self(n-2)): n in [1..30]]; // G. C. Greubel, Dec 27 2019
    
  • Maple
    seq(coeff(series(x/(1-5*x+5*x^2), x, n+1), x, n), n = 0..30); # G. C. Greubel, Dec 27 2019
  • Mathematica
    CoefficientList[Series[x/(1-5x+5x^2), {x,0,30}], x] (* Michael De Vlieger, Dec 22 2019 *)
    Table[If[EvenQ[n], 5^(n/2)*Fibonacci[n], 5^((n-1)/2)*LucasL[n]], {n,0,30}] (* G. C. Greubel, Dec 27 2019 *)
    LinearRecurrence[{5,-5},{0,1},30] (* Harvey P. Dale, Mar 21 2023 *)
  • PARI
    my(x='x+O('x^30)); concat([0], Vec(x/(1-5*x+5*x^2))) \\ G. C. Greubel, Dec 27 2019
    
  • Sage
    def A093131_list(prec):
        P. = PowerSeriesRing(ZZ, prec)
        return P( x/(1-5*x+5*x^2) ).list()
    A093131_list(30) # G. C. Greubel, Dec 27 2019
    

Formula

G.f.: x/(1 - 5*x + 5*x^2).
a(n) = 5*a(n-1) - 5*a(n-2).
a(n) = (((5 + sqrt(5))/2)^n - ((5 - sqrt(5))/2)^n)/sqrt(5).
a(n) = A093130(n)/2^n.
a(n) = Sum_{k=0..n} Sum_{j=0..n} C(n, j)*C(j, k)*Fibonacci(j-k). - Paul Barry, Feb 15 2005
a(n) = Sum_{k=0..n} C(n, k)*2^k*Fibonacci(n-k) = Sum_{k=0..n} C(n, k)*2^(n-k) * Fibonacci(k). - Paul Barry, Apr 22 2005
a(n) = A030191(n-1), n > 0. - R. J. Mathar, Sep 05 2008
E.g.f.: 2*exp(5*x/2)*sinh(sqrt(5)*x/2)/sqrt(5). - Ilya Gutkovskiy, Aug 11 2017
From Kai Wang, Dec 22 2019: (Start)
a(n) = Sum_{i=0..n-1; j=0..n-1; i+2*j=n-1} 5^i*((i+j)!/(i!*j!)).
a(n*k)/a(k) = Sum_{i=0..n-1; j=0..n-1; i+2*j=n-1} (-1)^(j*(k-1))*b(k)^i*((i+j)!/(i!*j!)).
a((2*m+1)*k)/a(k) = Sum_{i=0..m-1} (-1)^(i*k)*A020876((2*m-2*i)*k) + 5^(m*k).
a(2*m*k)/a(k) = Sum_{i=0..m-1} (-1)^(i*k)*A020876((2*m-2*i-1)*k}.
a(m+r)*a(n+s) - a(m+s)*a(n+r) = -5^(n+s)*a(m-n)*a(r-s).
a(m+r)*a(n+s) + a(m+s)*a(n+r) = (2*A020876(m+n+r+s) - 5^(n+s)*A020876(m-n)*A020876(r-s))/5.
A020876(m+r)*A020876(n+s) - A020876(m+s)*A020876(n+r) = 5^(n+s+1)*a(m-n)*a(r-s).
A020876(m+r)*A020876(n+s) - 5*a(m+s)*a(n+r) = 5^(n+s)*A020876(m-n)*A020876(r-s).
A020876(m+r)*A020876(n+s) + 5*a(m+s)*a(n+r) = 2*A020876(m+n+r+s) + 5^(n+s+1)*a(m-n)*a(r-s).
a(n)^2 - a(n+1)*a(n-1) = 5^(n-1).
a(n)^2 - a(n+r)*a(n-r) = 5^(n-r)*a(r)^2.
a(m)*a(n+1) - a(m+1)*a(n) = 5^n*a(m-n).
a(m-n) = (a(m)*A020876(n) - A020876(m)*a(n))/(2*5^n).
a(m+n) = (a(m)*A020876(n) + A020876(m)*a(n))/2.
A020876(n)^2 - A020876(n+r)*A020876(n-r) = -5^(n-r+1)*a(r)^2.
A020876(m)*A020876(n+1) - A020876(m+1)*A020876(n) = -5^(n+1)*a(m-n).
A020876(m+n) - 5^(n)*A020876(m-n) = 5*a(m)*a(n).
A020876(m-n) = (A020876(m)*A020876(n) - 5*a(m)*a(n))/(2*5^n).
A020876(m+n) = (A020876(m)*A020876(n) + 5*a(m)*a(n))/2. (End)
a(2*n) = 5^n*Fibonacci(2*n), a(2*n+1) = 5^n*Lucas(2*n+1). - G. C. Greubel, Dec 27 2019
a(n) = Sum_{k=0..n} (-1)^(k+1)*binomial(2*n, n+k)*(k|5), where (k|5) is the Legendre symbol. - Greg Dresden, Oct 14 2022

A109106 a(n) = (1/sqrt(5))*((sqrt(5) + 1)*((15 + 5*sqrt(5))/2)^(n-1) + (sqrt(5) - 1)*((15 - 5*sqrt(5))/2)^(n-1)).

Original entry on oeis.org

2, 20, 250, 3250, 42500, 556250, 7281250, 95312500, 1247656250, 16332031250, 213789062500, 2798535156250, 36633300781250, 479536132812500, 6277209472656250, 82169738769531250, 1075615844726562500
Offset: 1

Views

Author

Emeric Deutsch, Jun 19 2005

Keywords

Comments

Kekulé numbers for certain benzenoids.

References

  • S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (p. 215, K{T_m}).

Crossrefs

Cf. A179135. - Johannes W. Meijer, Jul 01 2010

Programs

  • Maple
    a:=n->(1/sqrt(5))*((sqrt(5)+1)*((15+5*sqrt(5))/2)^(n-1)+(sqrt(5)-1)*((15-5*sqrt(5))/2)^(n-1)): seq(expand(a(n)),n=1..19);

Formula

G.f.: 2z(1-5z)/(1 - 15z + 25z^2).
From Johannes W. Meijer, Jul 01 2010: (Start)
a(n) = A178381(4*n+2).
Lim_{k->infinity} a(n+k)/a(k) = (A020876(2*n) + 5*A039717(2*n-2)*sqrt(5))/2.
Lim_{n->infinity} A020876(2*n)/(5*A039717(2*n-2)) = sqrt(5).
(End)
a(n) = 2*5^(n-1)*Fibonacci(2*n-1). - Ehren Metcalfe, Apr 21 2018

A094825 Number of (s(0), s(1), ..., s(2n)) such that 0 < s(i) < 10 and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n, s(0) = 1, s(2n) = 7.

Original entry on oeis.org

1, 7, 35, 153, 624, 2444, 9333, 35055, 130207, 479941, 1759616, 6427032, 23412105, 85121783, 309062619, 1121050449, 4063463728, 14721293860, 53313308477, 193023319071, 698715633111, 2528895064637, 9152032060800, 33118656195888
Offset: 3

Views

Author

Herbert Kociemba, Jun 15 2004

Keywords

Crossrefs

Cf. A094865 (partial sums).

Programs

  • Mathematica
    LinearRecurrence[{8,-21,20,-5},{1,7,35,153},30] (* Harvey P. Dale, Jan 16 2015 *)

Formula

a(n) = (1/5)*Sum_{r=1..9} sin(r*Pi/10)*sin(7*r*Pi/10)*(2*cos(r*Pi/10))^(2n).
a(n) = 8*a(n-1) - 21*a(n-2) + 20*a(n-3) - 5*a(n-4).
G.f.: x^3*(1-x)/( (1-3*x+x^2)*(1-5*x+5*x^2) ).
a(n) = -A001906(n)/2 + A020876(n)/10. - R. J. Mathar, Jun 24 2011

A057282 Coefficient triangle of polynomials (falling powers) related to Fibonacci convolutions. Companion triangle to A057281.

Original entry on oeis.org

2, 5, 17, 15, 120, 225, 50, 700, 3050, 4080, 175, 3775, 28625, 89225, 94440, 625, 19225, 223175, 1208975, 3006000, 2666880, 2250, 93500, 1537100, 12689800, 54824650, 115299900, 89016480, 8125, 438250, 9670750, 112454500, 737744125
Offset: 1

Views

Author

Wolfdieter Lang, Sep 13 2000

Keywords

Comments

The row polynomials are q(k,x) := sum(a(k,m)*x^(k-m),m=0..k), k=0,1,2,..
The k-th convolution of F0(n) := A000045(n+1), n >= 0, (Fibonacci numbers starting with F0(0)=1) with itself is Fk(n) := A037027(n+k,k) = (p(k-1,n)*(n+1)*F0(n+1) + q(k-1,n)*(n+2)*F0(n))/(k!*5^k), k=1,2,..., where the companion polynomials p(k,n) := sum(b(k,m)*n^(k-m),m=0..k) are the row polynomials of triangle b(k,m)= A057281(k,m).
a(k,0)= A020876(k), k >= 0.

Examples

			k=2: F2(n)=((5*n^2+21*n+16)*F(n+2)+(5*n^2+27*n+34)*F(n+1))/50, F(n) := A000045(n); see A001628.
2; 5,17; 15,120,225; 50,700,3050,4080; 175,3775,28625,89225,94440; ...
		

Crossrefs

A087424 a(n) = S(4*n,4)/S(n,4) where S(n,m) = Sum_{k=0..n} binomial(n,k)*Fibonacci(m*k).

Original entry on oeis.org

567, 239841, 114082668, 55125843489, 26697877691247, 12934267027240356, 6266540498895923463, 3036106030479071781249, 1470978970343729016987852, 712682440446248640284336721, 345291321126117622870522555983, 167292036479044881831300837903684, 81052212349412217472309893818152407
Offset: 1

Views

Author

Benoit Cloitre, Oct 22 2003

Keywords

Crossrefs

Cf. A020876.

Programs

  • Magma
    [27^n*Fibonacci(8*n)/Fibonacci(2*n): n in [1..15]]; // Vincenzo Librandi, Aug 04 2018
  • Mathematica
    Table[(27^n Fibonacci[8 n] / Fibonacci[2 n]), {n, 15}] (* Vincenzo Librandi, Aug 04 2018 *)
    LinearRecurrence[{567,-40824,413343,-531441},{567,239841,114082668,55125843489},20] (* Harvey P. Dale, Jun 23 2020 *)

Formula

a(n) = (243+108*sqrt(5))^n+(243-108*sqrt(5))^n+((81+27*sqrt(5))/2)^n+((81-27*sqrt(5))/2)^n.
G.f.: -81*x*(26244*x^3-15309*x^2+1008*x-7) / ((729*x^2-486*x+1)*(729*x^2-81*x+1)). - Colin Barker, Dec 01 2012
a(n)/3^(3*n) = L(2*n)*L(4*n) = L(2*n) + L(6*n), L=A000032. - Ehren Metcalfe, Apr 21 2018
a(n) = 27^n*F(8*n)/F(2*n), F=A000045. - Ehren Metcalfe, Aug 03 2018

A087425 a(n) = S(5*n,5)/S(n,5) where S(n,m) = Sum_{k=0..n} binomial(n,k)*Fibonacci(m*k).

Original entry on oeis.org

23105, 459119455, 9758296035305, 208416652653910655, 4452963734477926435505, 95143212432467064852443605, 2032859482921447476046969568705, 43434715031065603778576465510557055
Offset: 1

Views

Author

Benoit Cloitre, Oct 22 2003

Keywords

Crossrefs

Cf. A020876.
Cf. A087423 (m=3), A087424 (m=4).

Programs

  • Maple
    S:=proc(n,m) add(binomial(n,k)*combinat:-fibonacci(m*k),k=0..n) end: m:=5: seq(S(m*n,m)/S(n,m),n=1..16); # Georg Fischer, Jul 07 2021

Formula

a(n) = 121^n+{(21367+9555*sqrt(5))/2}^n+{(21367-9555*sqrt(5))/2}^n+{(1617+715*sqrt(5))/2}^n+{(1617-715*sqrt(5))/2}^n.
a(n) = (x_1)^n+(x_2)^n+(x_3)^n+(x_4)^n+(x_5)^n where (x_i) (1<=i<=5) are the roots of X^5-23105*X^4+37360785*X^3-4520654985*X^2+40931916905*X-25937424601.
Showing 1-10 of 15 results. Next