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-4 of 4 results.

A204922 Ordered differences of Fibonacci numbers.

Original entry on oeis.org

1, 2, 1, 4, 3, 2, 7, 6, 5, 3, 12, 11, 10, 8, 5, 20, 19, 18, 16, 13, 8, 33, 32, 31, 29, 26, 21, 13, 54, 53, 52, 50, 47, 42, 34, 21, 88, 87, 86, 84, 81, 76, 68, 55, 34, 143, 142, 141, 139, 136, 131, 123, 110, 89, 55, 232, 231, 230, 228, 225, 220, 212, 199, 178
Offset: 1

Views

Author

Clark Kimberling, Jan 21 2012

Keywords

Comments

For a guide to related sequences, see A204892. For numbers not in A204922, see A050939.
From Emanuele Munarini, Mar 29 2012: (Start)
Diagonal elements = Fibonacci numbers F(n+1) (A000045)
First column = Fibonacci numbers - 1 (A000071);
Second column = Fibonacci numbers - 2 (A001911);
Row sums = n*F(n+3) - F(n+2) + 2 (A014286);
Central coefficients = F(2*n+1) - F(n+1) (A096140).
(End)

Examples

			a(1) = s(2) - s(1) = F(3) - F(2) = 2-1 = 1, where F=A000045;
a(2) = s(3) - s(1) = F(4) - F(2) = 3-1 = 2;
a(3) = s(3) - s(2) = F(4) - F(3) = 3-2 = 1;
a(4) = s(4) - s(1) = F(5) - F(2) = 5-1 = 4.
From _Emanuele Munarini_, Mar 29 2012: (Start)
Triangle begins:
   1;
   2,  1;
   4,  3,  2;
   7,  6,  5,  3;
  12, 11, 10,  8,  5;
  20, 19, 18, 16, 13,  8;
  33, 32, 31, 29, 26, 21, 13;
  54, 53, 52, 50, 47, 42, 34, 21;
  88, 87, 86, 84, 81, 76, 68, 55, 34;
  ... (End)
		

Crossrefs

Programs

  • Magma
    /* As triangle */ [[Fibonacci(n+2)-Fibonacci(k+1) : k in [1..n]]: n in [1.. 15]]; // Vincenzo Librandi, Aug 04 2015
    
  • Mathematica
    (See the program at A204924.)
  • Maxima
    create_list(fib(n+3)-fib(k+2),n,0,20,k,0,n); /* Emanuele Munarini, Mar 29 2012 */
    
  • PARI
    {T(n,k) = fibonacci(n+2) - fibonacci(k+1)};
    for(n=1,15, for(k=1,n, print1(T(n,k), ", "))) \\ G. C. Greubel, Feb 03 2019
    
  • Sage
    [[fibonacci(n+2) - fibonacci(k+1) for k in (1..n)] for n in (1..15)] # G. C. Greubel, Feb 03 2019

Formula

From Emanuele Munarini, Mar 29 2012: (Start)
T(n,k) = Fibonacci(n+2) - Fibonacci(k+1).
T(n,k) = Sum_{i=k..n} Fibonacci(i+1). (End)

A005013 a(n) = 3*a(n-2) - a(n-4), a(0)=0, a(1)=1, a(2)=1, a(3)=4. Alternates Fibonacci (A000045) and Lucas (A000032) sequences for even and odd n.

Original entry on oeis.org

0, 1, 1, 4, 3, 11, 8, 29, 21, 76, 55, 199, 144, 521, 377, 1364, 987, 3571, 2584, 9349, 6765, 24476, 17711, 64079, 46368, 167761, 121393, 439204, 317811, 1149851, 832040, 3010349, 2178309, 7881196, 5702887, 20633239, 14930352, 54018521, 39088169, 141422324
Offset: 0

Views

Author

Keywords

Comments

S(n,sqrt(5)), with the Chebyshev polynomials A049310, is an integer sequence in the real quadratic number field Q(sqrt(5)) with basis numbers <1,phi>, phi:=(1+sqrt(5))/2. S(n,sqrt(5)) = A(n) + 2*B(n)*phi, with A(n)= a(n+1)*(-1)^n and B(n)= A147600(n-1), n>=0, with A147600(-1):=0.
a(n) = p(n+1) where p(x) is the unique degree-(n-1) polynomial such that p(k) = Fibonacci(k) for k = 1, ..., n. - Michael Somos, Jan 08 2012
Row sums of A227431. - Richard R. Forberg, Jul 29 2013
This is the sequence of Lehmer numbers u_n(sqrt(R),Q) with the parameters R = 5 and Q = 1. It is a strong divisibility sequence, that is, gcd(a(n), a(m)) = a(gcd(n,m)) for all natural numbers n and m. The sequence satisfies a linear recurrence of order four. - Peter Bala, Apr 18 2014
The sequence of convergents of the 2-periodic continued fraction [0; 1, -5, 1, -5, ...] = 1/(1 - 1/(5 - 1/(1 - 1/(5 - ...)))) = (1/2)*(5 - sqrt(5)) begins [0/1, 1/1, 5/4, 4/3, 15/11, 11/8, 40/29, ...]; the denominators give the present sequence. The sequence of numerators [0, 1, 5, 4, 15, 11, 40, ...] is A203976. Cf. A108412 and A026741. - Peter Bala, May 19 2014
Define a binary operation o on the real numbers by x o y = x*sqrt(1 + y^2) + y*sqrt(1 + x^2). The operation o is commutative and associative with identity 0. We have (1/2)*a(2*n + 1) = 1/2 o 1/2 o ... o 1/2 (2*n + 1 terms) and (1/2)*sqrt(5)* a(2*n) = 1/2 o 1/2 o ... o 1/2 (2*n terms). Cf. A084068 and A049629. - Peter Bala, Mar 23 2018

Examples

			G.f. = x + x^2 + 4*x^3 + 3*x^4 + 11*x^5 + 8*x^6 + 29*x^7 + 21*x^8 + 76*x^9 + ...
a(3) = 4 since p(x) = (x^2 - 3*x + 4) / 2 interpolates p(1) = 1, p(2) = 1, p(3) = 2, and p(4) = 4. - _Michael Somos_, Jan 08 2012
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • GAP
    a:=[0,1,1,4];; for n in [5..40] do a[n]:=3*a[n-2]-a[n-4]; od; a; # Muniru A Asiru, Oct 21 2018
  • Haskell
    a005013 n = a005013_list !! n
    a005013_list = alt a000045_list a000032_list where
       alt (f::fs) (:l:ls) = f : l : alt fs ls
    -- Reinhard Zumkeller, Jan 10 2012
    
  • Magma
    I:=[0,1,1,4]; [n le 4 select I[n]  else 3*Self(n-2) - Self(n-4): n in [1..40]]; // Vincenzo Librandi, Feb 09 2016
    
  • Maple
    with(combinat): A005013 := n-> if n mod 2 = 0 then fibonacci(n) else fibonacci(n+1)+fibonacci(n-1); fi;
    A005013:=z*(z**2+z+1)/((z**2+z-1)*(z**2-z-1)); # Simon Plouffe in his 1992 dissertation
  • Mathematica
    CoefficientList[Series[(x + x^2 + x^3)/(1 - 3x^2 + x^4), {x, 0, 40}], x]
    f[n_] = Product[(1 + 4*Sin[k*Pi/n]^2), {k, 1, Floor[(n - 1)/2]}]; a = Table[f[n], {n, 0, 30}]; Round[a]; FullSimplify[ExpandAll[a]] (* Roger L. Bagula and Gary W. Adamson, Nov 26 2008 *)
    LinearRecurrence[{0, 3, 0, -1}, {0, 1, 1, 4}, 100] (* G. C. Greubel, Feb 08 2016 *)
  • PARI
    {a(n) = if( n%2, fibonacci(n+1) + fibonacci(n-1), fibonacci(n))}; /* Michael Somos, Jan 08 2012 */
    
  • PARI
    {a(n) = if( n<0, -a(-n), subst( polinterpolate( vector( n, k, fibonacci(k))), x, n+1))}; /* Michael Somos, Jan 08 2012 */
    

Formula

a(1) = a(2) = 1, a(3) = 4, a(n) = (a(n-1) * a(n-2) - 1) / a(n-3), unless n=3. a(-n) = -a(n).
a(2n) = A001906(n), a(2n+1) = A002878(n). a(n)=F(n+1)+(-1)^(n+1)F(n-1). - Mario Catalani (mario.catalani(AT)unito.it), Sep 20 2002
G.f.: x*(1+x+x^2)/((1-x-x^2)*(1+x-x^2)).
a(n) = Product_{k=1..floor((n-1)/2)} (1 + 4*sin(k*Pi/n)^2). - Roger L. Bagula and Gary W. Adamson, Nov 26 2008
Binomial transform is A096140. - Michael Somos, Apr 13 2012
From Peter Bala, Apr 18 2014: (Start)
a(n) = (alpha^n - beta^n)/(alpha - beta) for n odd, and a(n) = (alpha^n - beta^n)/(alpha^2 - beta^2) for n even, where alpha = (1/2)*(sqrt(5) + 1) and beta = (1/2)*(sqrt(5) - 1). Equivalently, a(n) = U(n-1, sqrt(5)/2) for n odd and a(n) = (1/sqrt(5))*U(n-1, sqrt(5)/2) for n even, where U(n,x) is the Chebyshev polynomial of the second kind. (End)
E.g.f.: (Phi/sqrt(5))*exp(-Phi*x)*(exp(x)-1)*(exp(sqrt(5)*x) - 1/(Phi)^2), where Phi = (1+sqrt(5))/2. - G. C. Greubel, Feb 08 2016
a(n) = (5^floor((n-1)/2)/2^(n-1))*Sum_{k=0..n-1} binomial(n-1,k)/5^floor(k/2). - Tony Foster III, Oct 21 2018
a(n) = hypergeom([(1 - n)/2, (n + 1) mod 2 - n/2], [1 - n], -4) for n >= 2. - Peter Luschny, Sep 03 2019

Extensions

Additional comments from Michael Somos, Jun 01 2000

A128619 Triangle T(n, k) = A127647(n,k) * A128174(n,k), read by rows.

Original entry on oeis.org

1, 0, 1, 2, 0, 2, 0, 3, 0, 3, 5, 0, 5, 0, 5, 0, 8, 0, 8, 0, 8, 13, 0, 13, 0, 13, 0, 13, 0, 21, 0, 21, 0, 21, 0, 21, 34, 0, 34, 0, 34, 0, 34, 0, 34, 0, 55, 0, 55, 0, 55, 0, 55, 0, 55
Offset: 1

Views

Author

Gary W. Adamson, Mar 14 2007

Keywords

Comments

This triangle is different from A128618, which is equal to A128174 * A127647.

Examples

			First few rows of the triangle are:
   1;
   0,  1;
   2,  0,  2;
   0,  3,  0,  3;
   5,  0,  5,  0,  5;
   0,  8,  0,  8,  0,  8;
  13,  0, 13,  0, 13,  0, 13;
   0, 21,  0, 21,  0, 21,  0, 21,
  ...
		

Crossrefs

Programs

  • Magma
    [((n+k+1) mod 2)*Fibonacci(n): k in [1..n], n in [1..15]]; // G. C. Greubel, Mar 17 2024
    
  • Mathematica
    Table[Fibonacci[n]*Mod[n+k+1,2], {n,15}, {k,n}]//Flatten (* G. C. Greubel, Mar 16 2024 *)
  • SageMath
    flatten([[((n+k+1)%2)*fibonacci(n) for k in range(1,n+1)] for n in range(1,16)]) # G. C. Greubel, Mar 17 2024

Formula

T(n, k) = A127647 * A128174, an infinite lower triangular matrix. In odd rows, n terms of F(n), 0, F(n),...; in the n-th row. In even rows, n terms of 0, F(n), 0,...; in the n-th row.
Sum_{k=1..n} T(n, k) = A128620(n-1).
From G. C. Greubel, Mar 16 2024: (Start)
T(n, k) = Fibonacci(n)*(1 + (-1)^(n+k))/2.
Sum_{k=1..n} (-1)^(k-1)*T(n, k) = (-1)^n*A128620(n-1).
Sum_{k=1..floor((n+1)/2)} T(n-k+1, k) = (1/2)*(1-(-1)^n)*A096140(floor((n + 1)/2)).
Sum_{k=1..floor((n+1)/2)} (-1)^(k-1)*T(n-k+1, k) = (1/2)*(1 - (-1)^n)*( Fibonacci(n-1) + (-1)^floor((n-1)/2) * Fibonacci(floor((n-3)/2)) ). (End)

A199512 Triangle T(n,k) = Fibonacci(n+k+1), related to A000045 (Fibonacci numbers).

Original entry on oeis.org

1, 1, 2, 2, 3, 5, 3, 5, 8, 13, 5, 8, 13, 21, 34, 8, 13, 21, 34, 55, 89, 13, 21, 34, 55, 89, 144, 233, 21, 34, 55, 89, 144, 233, 377, 610, 34, 55, 89, 144, 233, 377, 610, 987, 1597, 55, 89, 144, 233, 377, 610, 987, 1597, 2584, 4181
Offset: 0

Views

Author

Philippe Deléham, Nov 07 2011

Keywords

Examples

			Triangle begins :
1
1, 2
2, 3, 5
3, 5, 8, 13
5, 8, 13, 21, 34
8, 13, 21, 34, 55, 89
		

Crossrefs

Rows sums : A096140, Diagonal sums : A128620.

Programs

  • PARI
    T(n, k) = fibonacci(n+k+1);
    tabl(nn) = for (n=0, nn, for (k=0, n, print1(T(n, k), ", ")); print); \\ Michel Marcus, Aug 01 2017

Formula

T(n,k) = T(n,k-1) + T(n-1,k-1) = T(n-1,k-1) + T(n-1,k).
T(n,0) = Fibonacci(n+1) = A000045(n+1).

Extensions

More terms from Michel Marcus, Aug 01 2017
Showing 1-4 of 4 results.