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 26 results. Next

A001906 F(2n) = bisection of Fibonacci sequence: a(n) = 3*a(n-1) - a(n-2).

Original entry on oeis.org

0, 1, 3, 8, 21, 55, 144, 377, 987, 2584, 6765, 17711, 46368, 121393, 317811, 832040, 2178309, 5702887, 14930352, 39088169, 102334155, 267914296, 701408733, 1836311903, 4807526976, 12586269025, 32951280099, 86267571272, 225851433717, 591286729879, 1548008755920
Offset: 0

Views

Author

N. J. A. Sloane

Keywords

Comments

Apart from initial term, same as A088305.
Second column of array A102310 and of A028412.
Numbers k such that 5*k^2 + 4 is a square. - Gregory V. Richardson, Oct 13 2002
Apart from initial terms, also Pisot sequences E(3,8), P(3,8), T(3,8). See A008776 for definitions of Pisot sequences.
Binomial transform of A000045. - Paul Barry, Apr 11 2003
Number of walks of length 2n+1 in the path graph P_4 from one end to the other one. Example: a(2)=3 because in the path ABCD we have ABABCD, ABCBCD and ABCDCD. - Emeric Deutsch, Apr 02 2004
Simplest example of a second-order recurrence with the sixth term a square.
Number of (s(0), s(1), ..., s(2n)) such that 0 < s(i) < 5 and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n, s(0) = 1, s(2n) = 3. - Lekraj Beedassy, Jun 11 2004
a(n) (for n > 0) is the smallest positive integer that cannot be created by summing at most n values chosen among the previous terms (with repeats allowed). - Andrew Weimholt, Jul 20 2004
All nonnegative integer solutions of Pell equation b(n)^2 - 5*a(n)^2 = +4 together with b(n) = A005248(n), n >= 0. - Wolfdieter Lang, Aug 31 2004
a(n+1) is a Chebyshev transform of 3^n (A000244), where the sequence with g.f. G(x) is sent to the sequence with g.f. (1/(1+x^2))G(x/(1+x^2)). - Paul Barry, Oct 25 2004
a(n) is the number of distinct products of matrices A, B, C, in (A+B+C)^n where commutator [A,B] = 0 but C does not commute with A or B. - Paul D. Hanna and Max Alekseyev, Feb 01 2006
Number of binary words with exactly k-1 strictly increasing runs. Example: a(3)=F(6)=8 because we have 0|0,1|0,1|1,0|01,01|0,1|01,01|1 and 01|01. Column sums of A119900. - Emeric Deutsch, Jul 23 2006
See Table 1 on page 411 of Lukovits and Janezic paper. - Parthasarathy Nambi, Aug 22 2006
Inverse: With phi = (sqrt(5) + 1)/2, log_phi((sqrt(5) a(n) + sqrt(5 a(n)^2 + 4))/2) = n. - David W. Cantrell (DWCantrell(AT)sigmaxi.net), Feb 19 2007
[1,3,8,21,55,144,...] is the Hankel transform of [1,1,4,17,75,339,1558,...](see A026378). - Philippe Deléham, Apr 13 2007
The Diophantine equation a(n) = m has a solution (for m >= 1) if and only if floor(arcsinh(sqrt(5)*m/2)/log(phi)) <> floor(arccosh(sqrt(5)*m/2)/log(phi)) where phi is the golden ratio. An equivalent condition is A130259(m) = A130260(m). - Hieronymus Fischer, May 25 2007
a(n+1) = AB^(n)(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., 1=`1`, 3=`10`, 8=`100`, 21=`1000`, ..., in Wythoff code.
Equals row sums of triangles A140069, A140736 and A140737. - Gary W. Adamson, May 25 2008
a(n) is also the number of idempotent order-preserving partial transformations (of an n-element chain) of width n (width(alpha) = max(Im(alpha))). Equivalently, it is the number of idempotent order-preserving full transformations (of an n-element chain). - Abdullahi Umar, Sep 08 2008
a(n) is the number of ways that a string of 0,1 and 2 of size (n-1) can be arranged with no 12-pairs. - Udita Katugampola, Sep 24 2008
Starting with offset 1 = row sums of triangle A175011. - Gary W. Adamson, Apr 03 2010
As a fraction: 1/71 = 0.01408450... or 1/9701 = 0.0001030821.... - Mark Dols, May 18 2010
Sum of the products of the elements in the compositions of n (example for n=3: the compositions are 1+1+1, 1+2, 2+1, and 3; a(3) = 1*1*1 + 1*2 + 2*1 + 3 = 8). - Dylon Hamilton, Jun 20 2010, Geoffrey Critzer, Joerg Arndt, Dec 06 2010
a(n) relates to regular polygons with even numbers of edges such that Product_{k=1..(n-2)/2} (1 + 4*cos^2 k*Pi/n) = even-indexed Fibonacci numbers with a(n) relating to the 2*n-gons. The constants as products = roots to even-indexed rows of triangle A152063. For example: a(5) = 55 satisfies the product formula relating to the 10-gon. - Gary W. Adamson, Aug 15 2010
Alternatively, product of roots to x^4 - 12x^3 + 51x^2 - 90x + 55, (10th row of triangle A152063) = (4.618...)*(3.618...)*(2.381...)*(1.381...) = 55. - Gary W. Adamson, Aug 15 2010
a(n) is the number of generalized compositions of n when there are i different types of i, (i=1,2,...). - Milan Janjic, Aug 26 2010
Starting with "1" = row sums of triangle A180339, and eigensequence of triangle A137710. - Gary W. Adamson, Aug 28 2010
a(2) = 3 is the only prime.
Number of nonisomorphic graded posets with 0 and uniform hasse graph of rank n > 0, with exactly 2 elements of each rank level above 0. (Uniform used in the sense of Retakh, Serconek, and Wilson. Graded used in Stanley's sense that every maximal chain has the same length n.) - David Nacin, Feb 13 2012
Pisano period lengths: 1, 3, 4, 3, 10, 12, 8, 6, 12, 30, 5, 12, 14, 24, 20, 12, 18, 12, 9, 30, ... - R. J. Mathar, Aug 10 2012
Solutions (x, y) = (a(n), a(n+1)) satisfying x^2 + y^2 = 3xy + 1. - Michel Lagneau, Feb 01 2014
For n >= 1, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,2}. - Milan Janjic, Jan 25 2015
With a(0) = 0, for n > 1, a(n) is the smallest number not already in the sequence such that a(n)^2 - a(n-1)^2 is a Fibonacci number. - Derek Orr, Jun 08 2015
Let T be the tree generated by these rules: 0 is in T, and if p is in T, then p + 1 is in T and x*p is in T and y*p is in T. The n-th generation of T consists of A001906(n) polynomials, for n >= 0. - Clark Kimberling, Nov 24 2015
For n > 0, a(n) = exactly the maximum area of a quadrilateral with sides in order of lengths F(n), F(n), L(n), and L(n) with L(n)=A000032(n). - J. M. Bergot, Jan 20 2016
a(n) = twice the area of a triangle with vertices at (L(n+1), L(n+2)), (F(n+1), F(n+1)), and (L(n+2), L(n+1)), with L(n)=A000032(n). - J. M. Bergot, Apr 20 2016
Except for the initial 0, this is the p-INVERT of (1,1,1,1,1,...) for p(S) = 1 - S - S^2; see A291000. - Clark Kimberling, Aug 24 2017
a(n+1) is the number of spanning trees of the graph T_n, where T_n is a sequence of n triangles, where adjacent triangles share an edge. - Kevin Long, May 07 2018
a(n) is the number of ways to partition [n] such that each block is a run of consecutive numbers, and each block has a fixed point, e.g., for n=3, 12|3 with 1 and 3 as fixed points is valid, but 13|2 is not valid as 1 and 3 do not form a run. Consequently, a(n) also counts the spanning trees of the graph given by taking a path with n vertices and adding another vertex adjacent to all of them. - Kevin Long, May 11 2018
From Wolfdieter Lang, May 31 2018: (Start)
The preceding comment can be paraphrased as follows. a(n) is the row sum of the array A305309 for n >= 1. The array A305309(n, k) gives the sum of the products of the block lengths of the set partition of [n] := {1, 2, ..., n} with A048996(n, k) blocks of consecutive numbers, corresponding to the compositions obtained from the k-th partition of n in Abramowitz-Stegun order. See the comments and examples at A305309.
{a(n)} also gives the infinite sequence of nonnegative numbers k for which k * ||k*phi|| < 1/sqrt(5), where the irrational number phi = A001622 (golden section), and ||x|| is the absolute value of the difference between x and the nearest integer. See, e.g., the Havil reference, pp. 171-172. (End)
This Chebyshev sequence a(n) = S(n-1, 3) (see a formula below) is related to the bisection of Fibonacci sequences {F(a,b;n)}{n>=0} with input F(a,b;0) = a and F(a,b;1) = b, by F(a,b;2*k) = (a+b)*S(k-1, 3) - a*S(k-2, 3) and F(a,b;2*k+1) = b*S(k, 3) + (a-b)*S(k-1, 3), for k >= 0, and S(-2, 3) = -1. Proof via the o.g.f.s GFeven(a,b,t) = (a - t*(2*a-b))/(1 - 3*t + t^2) and GFodd(a,b,t) = (b + t*(a-b))/(1 - 3*t + t^2). The special case a = 0, b = 1 gives back F(2*k) = S(k-1, 3) = a(k). - _Wolfdieter Lang, Jun 07 2019
a(n) is the number of tilings of two n X 1 rectangles joined orthogonally at a common end-square (so to have 2n-1 squares in a right-angle V shape) with only 1 X 1 and 2 X 1 tiles. This is a consequence of F(2n) = F(n+1)*F(n) + F(n)*F(n-1). - Nathaniel Gregg, Oct 10 2021
These are the denominators of the upper convergents to the golden ratio, tau; they are also the numerators of the lower convergents (viz. 1/1 < 3/2 < 8/5 < 21/13 < ... < tau < ... 13/8 < 5/3 < 2/1). - Clark Kimberling, Jan 02 2022
For n > 1, a(n) is the smallest Fibonacci number of unit equilateral triangle tiles needed to make an isosceles trapezoid of height F(n) triangles. - Kiran Ananthpur Bacche, Sep 01 2024

Examples

			G.f. = x + 3*x^2 + 8*x^3 + 21*x^4 + 55*x^5 + 144*x^6 + 377*x^7 + 987*x^8 + ...
a(3) = 8 because there are exactly 8 idempotent order-preserving full transformations on a 3-element chain, namely: (1,2,3)->(1,1,1),(1,2,3)->(2,2,2),(1,2,3)->(3,3,3),(1,2,3)->(1,1,3),(1,2,3)->(2,2,3),(1,2,3)->(1,2,2),(1,2,3)->(1,3,3),(1,2,3)->(1,2,3)-mappings are coordinate-wise. - _Abdullahi Umar_, Sep 08 2008

References

  • Mohammad K. Azarian, The Generating Function for the Fibonacci Sequence, Missouri Journal of Mathematical Sciences, Vol. 2, No. 2, Spring 1990, pp. 78-79. Zentralblatt MATH, Zbl 1097.11516.
  • Mohammad K. Azarian, A Generalization of the Climbing Stairs Problem II, Missouri Journal of Mathematical Sciences, Vol. 16, No. 1, Winter 2004, pp. 12-17.
  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 2,5,6,14,33,55.
  • R. J. Douglas, Tournaments that admit exactly one Hamiltonian cycle, Proc. London Math. Soc., 21 (1970), 716-730.
  • G. Everest, A. van der Poorten, I. Shparlinski and T. Ward, Recurrence Sequences, Amer. Math. Soc., 2003; see esp. p. 255.
  • A. Gerardin, Reply to Query 4389, L'Intermédiaire des Mathématiciens, 22 (1915), 23.
  • Julian Havil, The Irrationals, Princeton University Press, Princeton and Oxford, 2012, pp. 171-172.
  • Howie, J. M. Combinatorial and probabilistic results in transformation semigroups. Words, languages and combinatorics, II (Kyoto, 1992), 200--206, World Sci. Publ., River Edge, NJ, (1994).
  • Laradji, A. and Umar, A. Combinatorial results for semigroups of order-preserving full transformations. Semigroup Forum 72 (2006), 51-62.
  • I. Lukovits, A. Graovac, E. Kalman, G. Kaptay, P. Nagy, S. Nikolic, J. Sytchev and N. Trinajstich, "Nanotubes: Number of Kekulé Structures and Aromaticity", J. Chem. Inf. Comput. Sci, vol. 43 (2003), pp. 609-614. See Equation 6 on page 611.
  • T. Mansour, M. Shattuck, A statistic on n-color compositions and related sequences, Proc. Indian Acad. Sci. (Math. Sci.) Vol. 124, No. 2, May 2014, pp. 127-140.
  • H. Mathieu, Query 3932, L'Intermédiaire des Mathématiciens, 18 (1911), 222. - N. J. A. Sloane, Mar 08 2022
  • I. Niven and H. S. Zuckerman, An Introduction to the Theory of Numbers. 2nd ed., Wiley, NY, 1966, p. 101.
  • Paulo Ribenboim, Primes in Lucas sequences (Chap 4), in 'My Numbers, My Friends', Springer-Verlag 2000 NY, page 27.
  • 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).
  • R. Stanley, Enumerative combinatorics, Vol. 1, Cambridge University Press, Cambridge, 1997, pp. 96-100.

Links

Crossrefs

Fibonacci A000045 = union of this sequence and A001519.
Cf. A000041, A000032, A001622, A048996, A052529, A055991, A078812, A305309, A058038 (pairwise products), A153386.
Inverse sequences A130259 and A130260.
Cf. A033888, A033890, A249450, A344212.

Programs

  • Haskell
    a001906 n = a001906_list !! n
a001906_list =
   0 : 1 : zipWith (-) (map (* 3) $ tail a001906_list) a001906_list
-- Reinhard Zumkeller, Oct 03 2011
  • Magma
    [Fibonacci(2*n): n in [0..30]]; // Vincenzo Librandi, Sep 10 2014
  • Maple
    with(combstruct): SeqSeqSeqL := [T, {T=Sequence(S, card > 0), S=Sequence(U, card > 1), U=Sequence(Z, card >0)}, unlabeled]: seq(count(SeqSeqSeqL, size=n+1), n=0..28); # Zerinvary Lajos, Apr 04 2009
H := (n, a, b) -> hypergeom([a - n/2, b - n/2], [1 - n], -4):
a := n -> `if`(n = 0, 0, H(2*n, 1, 1/2)):
seq(simplify(a(n)), n=0..30); # Peter Luschny, Sep 03 2019
A001906 := proc(n)
    combinat[fibonacci](2*n) ;
end proc:
seq(A001906(n),n=0..20) ; # R. J. Mathar, Jan 11 2024
  • Mathematica
    f[n_] := Fibonacci[2n]; Array[f, 28, 0] (* or *)
LinearRecurrence[{3, -1}, {0, 1}, 28] (* Robert G. Wilson v, Jul 13 2011 *)
Take[Fibonacci[Range[0,60]],{1,-1,2}] (* Harvey P. Dale, May 23 2012 *)
Table[ ChebyshevU[n-1, 3/2], {n, 0, 30}] (* Jean-François Alcover, Jan 25 2013, after Michael Somos *)
CoefficientList[Series[(x)/(1 - 3x + x^2), {x, 0, 30}], x] (* Vincenzo Librandi, Sep 10 2014 *)
  • Maxima
    makelist(fib(2*n),n,0,30); /* Martin Ettl, Oct 21 2012 */
  • MuPAD
    numlib::fibonacci(2*n) $ n = 0..35; // Zerinvary Lajos, May 09 2008
  • PARI
    {a(n) = fibonacci(2*n)}; /* Michael Somos, Dec 06 2002 */
  • PARI
    {a(n) = subst( poltchebi(n+1)*4 - poltchebi(n)*6, x, 3/2)/5}; /* Michael Somos, Dec 06 2002 */
  • PARI
    {a(n) = polchebyshev( n-1, 2, 3/2)}; /* Michael Somos Jun 18 2011 */
  • PARI
    Vec(x/(1-3*x+x^2)+O(x^99)) \\ Charles R Greathouse IV, Oct 24 2012
  • Python
    def a(n, adict={0:0, 1:1}):
    if n in adict:
        return adict[n]
    adict[n]=3*a(n-1) - a(n-2)
    return adict[n] # David Nacin, Mar 04 2012
  • Sage
    [lucas_number1(n,3,1) for n in range(27)] # Zerinvary Lajos, Jun 25 2008
  • Sage
    [fibonacci(2*n) for n in range(0, 28)] # Zerinvary Lajos, May 15 2009

Formula

G.f.: x / (1 - 3*x + x^2). - Simon Plouffe in his 1992 dissertation
a(n) = 3*a(n-1) - a(n-2) = A000045(2*n).
a(n) = -a(-n).
a(n) = A060921(n-1, 0), n >= 1.
a(n) = sqrt((A005248(n)^2 - 4)/5).
a(n) = A007598(n) - A007598(n-2), n > 1.
a(n) = (ap^n - am^n)/(ap-am), with ap := (3+sqrt(5))/2, am := (3-sqrt(5))/2.
Invert transform of natural numbers: a(n) = Sum_{k=1..n} k*a(n-k), a(0) = 1. - Vladeta Jovovic, Apr 27 2001
a(n) = S(n-1, 3) with S(n, x) = U(n, x/2) Chebyshev's polynomials of the 2nd kind, see A049310.
a(n) = Sum_{k=0..n} binomial(n, k)*F(k). - Benoit Cloitre, Sep 03 2002
Limit_{n->infinity} a(n)/a(n-1) = 1 + phi = (3 + sqrt(5))/2. This sequence includes all of the elements of A033888 combined with A033890.
a(0)=0, a(1)=1, a(2)=3, a(n)*a(n-2) + 1 = a(n-1)^2. - Benoit Cloitre, Dec 06 2002
a(n) = n + Sum_{k=0..n-1} Sum_{i=0..k} a(i) = n + A054452(n). - Benoit Cloitre, Jan 26 2003
a(n) = Sum_{k=1..n} binomial(n+k-1, n-k). - Vladeta Jovovic, Mar 23 2003
E.g.f.: (2/sqrt(5))*exp(3*x/2)*sinh(sqrt(5)*x/2). - Paul Barry, Apr 11 2003
Second diagonal of array defined by T(i, 1) = T(1, j) = 1, T(i, j) = Max(T(i-1, j) + T(i-1, j-1); T(i-1, j-1) + T(i, j-1)). - Benoit Cloitre, Aug 05 2003
a(n) = F(n)*L(n) = A000045(n)*A000032(n). - Lekraj Beedassy, Nov 17 2003
F(2n+2) = 1, 3, 8, ... is the binomial transform of F(n+2). - Paul Barry, Apr 24 2004
Partial sums of A001519(n). - Lekraj Beedassy, Jun 11 2004
a(n) = Sum_{i=0..n-1} binomial(2*n-1-i, i)*5^(n-i-1)*(-1)^i. - Mario Catalani (mario.catalani(AT)unito.it), Jul 23 2004
a(n) = Sum_{k=0..n} binomial(n+k, n-k-1) = Sum_{k=0..n} binomial(n+k, 2k+1).
a(n+1) = Sum_{k=0..floor(n/2)} binomial(n-k, k)*(-1)^k*3^(n-2*k). - Paul Barry, Oct 25 2004
a(n) = (n*L(n) - F(n))/5 = Sum_{k=0..n-1} (-1)^n*L(2*n-2*k-1).
The i-th term of the sequence is the entry (1, 2) in the i-th power of the 2 X 2 matrix M = ((1, 1), (1, 2)). - Simone Severini, Oct 15 2005
Computation suggests that this sequence is the Hankel transform of A005807. The Hankel transform of {a(n)} is Det[{{a(1), ..., a(n)}, {a(2), ..., a(n+1)}, ..., {a(n), ..., a(2n-1)}}]. - John W. Layman, Jul 21 2000
a(n+1) = (A005248(n+1) - A001519(n))/2. - Creighton Dement, Aug 15 2004
a(n+1) = Sum_{i=0..n} Sum_{j=0..n} binomial(n-i, j)*binomial(n-j, i). - N. J. A. Sloane, Feb 20 2005
a(n) = (2/sqrt(5))*sinh(2*n*psi), where psi:=log(phi) and phi=(1+sqrt(5))/2. - Hieronymus Fischer, Apr 24 2007
a(n) = ((phi+1)^n - A001519(n))/phi with phi=(1+sqrt(5))/2. - Reinhard Zumkeller, Nov 22 2007
Row sums of triangle A135871. - Gary W. Adamson, Dec 02 2007
a(n)^2 = Sum_{k=1..n} a(2*k-1). This is a property of any sequence S(n) such that S(n) = B*S(n-1) - S(n-2) with S(0) = 0 and S(1) = 1 including {0,1,2,3,...} where B = 2. - Kenneth J Ramsey, Mar 23 2008
a(n) = 1/sqrt(5)*(phi^(2*n+2) - phi^(-2*n-2)), where phi = (1+sqrt(5))/2, the golden ratio. - Udita Katugampola (SIU), Sep 24 2008
If p[i] = i and if A is Hessenberg matrix of order n defined by: A[i,j] = p[j-i+1], (i<=j), A[i,j] = -1, (i = j+1), and A[i,j] = 0 otherwise. Then, for n >= 1, a(n) = det(A). - Milan Janjic, May 02 2010
If p[i] = Stirling2(i,2) and if A is the Hessenberg matrix of order n defined by: A[i,j] = p[j-i+1], (i<=j), A[i,j] = -1, (i = j+1), and A[i,j] = 0 otherwise. Then, for n >= 1, a(n-1) = det(A). - Milan Janjic, May 08 2010
a(n) = F(2*n+10) mod F(2*n+5).
a(n) = 1 + a(n-1) + Sum_{i=1..n-1} a(i), with a(0)=0. - Gary W. Adamson, Feb 19 2011
a(n) is equal to the permanent of the (n-1) X (n-1) Hessenberg matrix with 3's along the main diagonal, i's along the superdiagonal and the subdiagonal (i is the imaginary unit), and 0's everywhere else. - John M. Campbell, Jun 09 2011
a(n), n > 1 is equal to the determinant of an (n-x) X (n-1) tridiagonal matrix with 3's in the main diagonal, 1's in the super and subdiagonals, and the rest 0's. - Gary W. Adamson, Jun 27 2011
a(n) = b such that Integral_{x=0..Pi/2} sin(n*x)/(3/2-cos(x)) dx = c + b*log(3). - Francesco Daddi, Aug 01 2011
a(n+1) = Sum_{k=0..n} A101950(n,k)*2^k. - Philippe Deléham, Feb 10 2012
G.f.: A(x) = x/(1-3*x+x^2) = G(0)/sqrt(5); where G(k)= 1 -(a^k)/(1 - b*x/(b*x - 2*(a^k)/G(k+1))), a = (7-3*sqrt(5))/2, b = 3+sqrt(5), if |x|<(3-sqrt(5))/2 = 0.3819660...; (continued fraction 3 kind, 3-step ). - Sergei N. Gladkovskii, Jun 25 2012
a(n) = 2^n*b(n;1/2) = -b(n;-1), where b(n;d), n=0,1,...,d, denote the delta-Fibonacci numbers defined in comments to A000045 (see also Witula's et al. papers). - Roman Witula, Jul 12 2012
Product_{n>=1} (1 + 1/a(n)) = 1 + sqrt(5). - Peter Bala, Dec 23 2012
Product_{n>=2} (1 - 1/a(n)) = (1/6)*(1 + sqrt(5)). - Peter Bala, Dec 23 2012
G.f.: x/(1-2*x) + x^2/(1-2*x)/(Q(0)-x) where Q(k) = 1 - x/(x*k+1)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Feb 23 2013
G.f.: G(0)/2 - 1, where G(k) = 1 + 1/( 1 - x/(x + (1-x)^2/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jul 16 2013
G.f.: x*G(0)/(2-3*x), where G(k) = 1 + 1/( 1 - x*(5*k-9)/(x*(5*k-4) - 6/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jul 17 2013
Sum_{n>=1} 1/(a(n) + 1/a(n)) = 1. Compare with A001519, A049660 and A049670. - Peter Bala, Nov 29 2013
a(n) = U(n-1,3/2) where U(n-1,x) is Chebyshev polynomial of the second kind. - Milan Janjic, Jan 25 2015
The o.g.f. A(x) satisfies A(x) + A(-x) + 6*A(x)*A(-x) = 0. The o.g.f. for A004187 equals -A(sqrt(x))*A(-sqrt(x)). - Peter Bala, Apr 02 2015
For n > 1, a(n) = (3*F(n+1)^2 + 2*F(n-2)*F(n+1) - F(n-2)^2)/4. - J. M. Bergot, Feb 16 2016
For n > 3, a(n) = floor(MA) - 4 for n even and floor(MA) + 5 for n odd. MA is the maximum area of a quadrilateral with lengths of sides in order L(n), L(n), F(n-3), F(n+3), with L(n)=A000032(n). The ratio of the longer diagonal to the shorter approaches 5/3. - J. M. Bergot, Feb 16 2016
a(n+1) = Sum_{j=0..n} Sum_{k=0..j} binomial(n-j,k)*binomial(j,k)*2^(j-k). - Tony Foster III, Sep 18 2017
a(n) = Sum_{k=0..n-1} Sum_{i=0..n-1} C(k+i,k-i). - Wesley Ivan Hurt, Sep 21 2017
a(n) = Sum_{k=1..A000041(n)} A305309(n, k), n >= 1. Also row sums of triangle A078812.- Wolfdieter Lang, May 31 2018
a(n) = H(2*n, 1, 1/2) for n > 0 where H(n, a, b) -> hypergeom([a - n/2, b - n/2], [1 - n], -4). - Peter Luschny, Sep 03 2019
Sum_{n>=1} 1/a(n) = A153386. - Amiram Eldar, Oct 04 2020
a(n) = A249450(n) + 2. - Leo Tavares, Oct 10 2021
a(n) = -2/(sqrt(5)*tan(2*arctan(phi^(2*n)))), where phi = A001622 is the golden ratio. - Diego Rattaggi, Nov 21 2021
a(n) = sinh(2*n*arcsinh(1/2))/sqrt(5/4). - Peter Luschny, May 21 2022
From Amiram Eldar, Dec 02 2024: (Start)
Product_{n>=1} (1 - (-1)^n/a(n)) = 1 + 1/sqrt(5) (A344212).
Product_{n>=2} (1 + (-1)^n/a(n)) = (5/6) * (1 + 1/sqrt(5)). (End)
a(n) = Sum_{k>=0} Fibonacci(2*n*k)/(Lucas(2*n)^(k+1)). - Diego Rattaggi, Jan 12 2025
Sum_{n>=0} a(n)/3^n = 3. - Diego Rattaggi, Jan 20 2025

A008586 Multiples of 4.

Original entry on oeis.org

0, 4, 8, 12, 16, 20, 24, 28, 32, 36, 40, 44, 48, 52, 56, 60, 64, 68, 72, 76, 80, 84, 88, 92, 96, 100, 104, 108, 112, 116, 120, 124, 128, 132, 136, 140, 144, 148, 152, 156, 160, 164, 168, 172, 176, 180, 184, 188, 192, 196, 200, 204, 208, 212, 216, 220, 224, 228
Offset: 0

Views

Author

N. J. A. Sloane

Keywords

Comments

Apart from initial term(s), dimension of the space of weight 2n cusp forms for Gamma_0( 14 ).
A000466(n), a(n) and A053755(n) are Pythagorean triples. - Zak Seidov, Jan 16 2007
If X is an n-set and Y and Z disjoint 2-subsets of X then a(n-3) is equal to the number of 3-subsets of X intersecting both Y and Z. - Milan Janjic, Aug 26 2007
Number of n-permutations (n>=1) of 5 objects u, v, z, x, y with repetition allowed, containing n-1 u's. Example: if n=1 then n-1 = zero (0) u, a(1)=4 because we have v, z, x, y. If n=2 then n-1 = one (1) u, a(2)=8 because we have vu, zu, xu, yu, uv, uz, ux, uy. A038231 formatted as a triangular array: diagonal: 4, 8, 12, 16, 20, 24, 28, 32, ... - Zerinvary Lajos, Aug 06 2008
For n > 0: numbers having more even than odd divisors: A048272(a(n)) < 0. - Reinhard Zumkeller, Jan 21 2012
A214546(a(n)) < 0 for n > 0. - Reinhard Zumkeller, Jul 20 2012
A090418(a(n)) = 0 for n > 0. - Reinhard Zumkeller, Aug 06 2012
Terms are the differences of consecutive centered square numbers (A001844). - Mihir Mathur, Apr 02 2013
a(n)*Pi = nonnegative zeros of the cycloid generated by a circle of radius 2 rolling along the positive x-axis from zero. - Wesley Ivan Hurt, Jul 01 2013
Apart from the initial term, number of vertices of minimal path on an n-dimensional cubic lattice (n>1) of side length 2, until a self-avoiding walk gets stuck. A004767 + 1. - Matthew Lehman, Dec 23 2013
The number of orbits of Aut(Z^7) as function of the infinity norm n of the representative lattice point of the orbit, when the cardinality of the orbit is equal to 2688. - Philippe A.J.G. Chevalier, Dec 29 2015
First differences of A001844. - Robert Price, May 13 2016
Numbers k such that Fibonacci(k) is a multiple of 3 (A033888). - Bruno Berselli, Oct 17 2017

Links

Crossrefs

Cf. A000290, A000466, A001844, A004767, A008574, A030308, A033888, A035008, A038231, A048272, A053755, A090418, A214546.
Number of orbits of Aut(Z^7) as function of the infinity norm A000579, A154286, A102860, A002412, A045943, A115067, A008585, A005843, A001477, A000217.

Programs

Formula

a(n) = A008574(n), n>0. - R. J. Mathar, Oct 28 2008
a(n) = Sum_{k>=0} A030308(n,k)*2^(k+2). - Philippe Deléham, Oct 17 2011
a(n+1) = A000290(n+2) - A000290(n). - Philippe Deléham, Mar 31 2013
G.f.: 4*x/(1-x)^2. - David Wilding, Jun 21 2014
E.g.f.: 4*x*exp(x). - Stefano Spezia, May 18 2021

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

N. J. A. Sloane

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.

Links

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

A004187 a(n) = 7*a(n-1) - a(n-2) with a(0) = 0, a(1) = 1.

Original entry on oeis.org

0, 1, 7, 48, 329, 2255, 15456, 105937, 726103, 4976784, 34111385, 233802911, 1602508992, 10983760033, 75283811239, 516002918640, 3536736619241, 24241153416047, 166151337293088, 1138818207635569, 7805576116155895, 53500214605455696, 366695926122033977
Offset: 0

Views

Author

N. J. A. Sloane, R. K. Guy

Keywords

Comments

Define the sequence T(a_0,a_1) by a_{n+2} is the greatest integer such that a_{n+2}/a_{n+1}= 0 . A004187 (with initial 0 omitted) is T(1,7).
This is a divisibility sequence.
For n>=2, a(n) equals the permanent of the (n-1) X (n-1) tridiagonal matrix with 7's along the main diagonal, and i's along the superdiagonal and the subdiagonal (i is the imaginary unit). - John M. Campbell, Jul 08 2011
a(n) and b(n) := A056854(n) are the proper and improper nonnegative solutions of the Pell equation b(n)^2 - 5*(3*a(n))^2 = +4. see the cross-reference to A056854 below. - Wolfdieter Lang, Jun 26 2013
For n>=1, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,2,3,4,5,6}. - Milan Janjic, Jan 25 2015
The digital root is A253298, which shares its digital root with A253368. - Peter M. Chema, Jul 04 2016
Lim_{n->oo} a(n+1)/a(n) = 2 + 3*phi = 1+ A090550 = 6.854101... - Wolfdieter Lang, Nov 16 2023

Examples

			a(2) = 7*a(1) - a(0) = 7*7 - 1 = 48. - _Michael B. Porter_, Jul 04 2016

Links

Crossrefs

Cf. A000027, A001906, A001353, A004254, A001109, A049685, A033888. a(n)=sqrt((A056854(n)^2 - 4)/45).
Second column of array A028412.

Programs

  • Magma
    [Fibonacci(4*n)/3 : n in [0..30]]; // Vincenzo Librandi, Jun 07 2011
  • Magma
    /* By definition: */ [n le 2 select n-1 else 7*Self(n-1)-Self(n-2): n in [1..23]]; // Bruno Berselli, Dec 24 2012
  • Maple
    seq(combinat:-fibonacci(4*n)/3, n = 0 .. 30); # Robert Israel, Jan 26 2015
  • Mathematica
    LinearRecurrence[{7,-1},{0,1},30] (* Harvey P. Dale, Jul 13 2011 *)
CoefficientList[Series[x/(1 - 7*x + x^2), {x, 0, 50}], x] (* Vincenzo Librandi, Dec 23 2012 *)
  • Maxima
    a[0]:0$ a[1]:1$ a[n]:=7*a[n-1] - a[n-2]$ A004187(n):=a[n]$ makelist(A004187(n),n,0,30); /* Martin Ettl, Nov 11 2012 */
  • MuPAD
    numlib::fibonacci(4*n)/3 $ n = 0..25; // Zerinvary Lajos, May 09 2008
  • PARI
    a(n)=fibonacci(4*n)/3 \\ Charles R Greathouse IV, Mar 09 2012
  • PARI
    concat(0, Vec(x/(1-7*x+x^2) + O(x^99))) \\ Altug Alkan, Jul 03 2016
  • Sage
    [lucas_number1(n,7,1) for n in range(27)] # Zerinvary Lajos, Jun 25 2008
  • Sage
    [fibonacci(4*n)/3 for n in range(0, 21)] # Zerinvary Lajos, May 15 2009

Formula

G.f.: x/(1-7*x+x^2).
a(n) = F(4*n)/3 = A033888(n)/3, where F=A000045 (the Fibonacci sequence).
a(n) = S(2*n-1, sqrt(9))/sqrt(9) = S(n-1, 7); S(n, x) := U(n, x/2), Chebyshev polynomials of the 2nd kind, A049310.
a(n) = Sum_{i = 0..n-1} C(2*n-1-i, i)*5^(n-i-1). - Mario Catalani (mario.catalani(AT)unito.it), Jul 23 2004
[A049685(n-1), a(n)] = [1,5; 1,6]^n * [1,0]. - Gary W. Adamson, Mar 21 2008
a(n) = A167816(4*n). - Reinhard Zumkeller, Nov 13 2009
a(n) = (((7+sqrt(45))/2)^n-((7-sqrt(45))/2)^n)/sqrt(45). - Noureddine Chair, Aug 31 2011
a(n+1) = Sum_{k = 0..n} A101950(n,k)*6^k. - Philippe Deléham, Feb 10 2012
a(n) = (A081072(n)/3)-1. - Martin Ettl, Nov 11 2012
From Peter Bala, Dec 23 2012: (Start)
Product {n >= 1} (1 + 1/a(n)) = (1/5)*(5 + 3*sqrt(5)).
Product {n >= 2} (1 - 1/a(n)) = (1/14)*(5 + 3*sqrt(5)). (End)
From Peter Bala, Apr 02 2015: (Start)
Sum_{n >= 1} a(n)*x^(2*n) = -A(x)*A(-x), where A(x) = Sum_{n >= 1} Fibonacci(2*n)* x^n.
1 + 5*Sum_{n >= 1} a(n)*x^(2*n) = F(x)*F(-x) = G(x)*G(-x), where F(x) = 1 + A(x) and G(x) = 1 + 5*A(x).
1 + Sum_{n >= 1} a(n)*x^(2*n) = H(x)*H(-x) = I(x)*I(-x), where H(x) = 1 + Sum_{n >= 1} Fibonacci(2*n + 3)*x^n and I(x) = 1 + x + x*Sum_{n >= 1} Fibonacci(2*n - 1)*x^n. (End)
E.g.f.: 2*exp(7*x/2)*sinh(3*sqrt(5)*x/2)/(3*sqrt(5)). - Ilya Gutkovskiy, Jul 03 2016
a(n) = Sum_{k = 0..n-1} (-1)^(n+k+1)*9^k*binomial(n+k, 2*k+1). - Peter Bala, Jul 17 2023
a(n) = Sum_{k = 0..floor(n/2)} (-1)^k*7^(n-2*k)*binomial(n-k, k). - Greg Dresden, Aug 03 2024
From Peter Bala, Jul 22 2025: (Start)
The following products telescope:
Product {n >= 2} (1 + (-1)^n/a(n)) = (3/14)*(3 + sqrt(5)).
Product {n >= 1} (1 - (-1)^n/a(n)) = (1/3)*(3 + sqrt(5)).
Product_{n >= 1} (a(2*n) + 1)/(a(2*n) - 1) = (3/5)*sqrt(5). (End)

Extensions

Entry improved by comments from Michael Somos and Wolfdieter Lang, Aug 02 2000

A033890 a(n) = Fibonacci(4*n + 2).

Original entry on oeis.org

1, 8, 55, 377, 2584, 17711, 121393, 832040, 5702887, 39088169, 267914296, 1836311903, 12586269025, 86267571272, 591286729879, 4052739537881, 27777890035288, 190392490709135, 1304969544928657, 8944394323791464, 61305790721611591, 420196140727489673
Offset: 0

Views

Author

N. J. A. Sloane

Keywords

Comments

(x,y) = (a(n), a(n+1)) are solutions of (x+y)^2/(1+xy)=9, the other solutions are in A033888. - Floor van Lamoen, Dec 10 2001
This sequence consists of the odd-indexed terms of A001906 (whose terms are the values of x such that 5*x^2 + 4 is a square). The even-indexed terms of A001906 are in A033888. Limit_{n->infinity} a(n)/a(n-1) = phi^4 = (7 + 3*sqrt(5))/2. - Gregory V. Richardson, Oct 13 2002
General recurrence is a(n) = (a(1)-1)*a(n-1) - a(n-2), a(1) >= 4, lim_{n->infinity} a(n) = x*(k*x+1)^n, k = a(1) - 3, x = (1 + sqrt((a(1)+1)/(a(1)-3)))/2. Examples in OEIS: a(1)=4 gives A002878. a(1)=5 gives A001834. a(1)=6 gives A030221. a(1)=7 gives A002315. a(1)=8 gives A033890. a(1)=9 gives A057080. a(1)=10 gives A057081. - Ctibor O. Zizka, Sep 02 2008
Indices of square numbers which are also 12-gonal. - Sture Sjöstedt, Jun 01 2009
For positive n, a(n) equals the permanent of the (2n) X (2n) tridiagonal matrix with 3's along the main diagonal, and i's along the superdiagonal and the subdiagonal (i is the imaginary unit). - John M. Campbell, Jul 08 2011
If we let b(0) = 0 and, for n >= 1, b(n) = A033890(n-1), then the sequence b(n) will be F(4n-2) and the first difference is L(4n) or A056854. F(4n-2) is also the ratio of golden spiral length (rounded to the nearest integer) after n rotations. L(4n) is also the pitch length ratio. See illustration in links. - Kival Ngaokrajang, Nov 03 2013
The aerated sequence (b(n))n>=1 = [1, 0, 8, 0, 55, 0, 377, 0, ...] is a fourth-order linear divisibility sequence; that is, if n | m then b(n) | b(m). It is the case P1 = 0, P2 = -5, Q = -1 of the 3-parameter family of divisibility sequences found by Williams and Guy. See A100047. - Peter Bala, Mar 22 2015
Solutions y of Pell equation x^2 - 5*y^2 = 4; corresponding x values are in A342710 (see A342709). - Bernard Schott, Mar 19 2021

Links

Crossrefs

Cf. A000045, A001906, A100047.
Cf. A342709, A342710.

Programs

  • Magma
    [Fibonacci(4*n +2): n in [0..100]]; // Vincenzo Librandi, Apr 17 2011
  • Maple
    A033890 := proc(n)
    option remember;
    if n <= 1 then
        op(n+1,[1,8]);
    else
        7*procname(n-1)-procname(n-2) ;
    end if;
end proc: # R. J. Mathar, Apr 30 2017
  • Mathematica
    Table[Fibonacci[4n + 2], {n, 0, 14}] (* Vladimir Joseph Stephan Orlovsky, Jul 21 2008 *)
LinearRecurrence[{7, -1}, {1, 8}, 50] (* G. C. Greubel, Jul 13 2017 *)
a[n_] := (GoldenRatio^(2 (1 + 2 n)) - GoldenRatio^(-2 (1 + 2 n)))/Sqrt[5]
Table[a[n] // FullSimplify, {n, 0, 21}] (* Gerry Martens, Aug 20 2025 *)
  • PARI
    a(n)=fibonacci(4*n+2);

Formula

G.f.: (1+x)/(1-7*x+x^2).
a(n) = 7*a(n-1) - a(n-2), n > 1; a(0)=1, a(1)=8.
a(n) = S(n,7) + S(n-1,7) = S(2*n,sqrt(9) = 3), where S(n,x) = U(n,x/2) are Chebyshev's polynomials of the 2nd kind. Cf. A049310. S(n,7) = A004187(n+1), S(n,3) = A001906(n+1).
a(n) = ((7+3*sqrt(5))^n - (7-3*sqrt(5))^n + 2*((7+3*sqrt(5))^(n-1) - ((7-3*sqrt(5))^(n-1)))) / (3*(2^n)*sqrt(5)). - Gregory V. Richardson, Oct 13 2002
Let q(n, x) = Sum_{i=0..n} x^(n-i)*binomial(2*n-i, i); then a(n) = (-1)^n*q(n, -9). - Benoit Cloitre, Nov 10 2002
a(n) = L(n,-7)*(-1)^n, where L is defined as in A108299; see also A049685 for L(n,+7). - Reinhard Zumkeller, Jun 01 2005
Define f(x,s) = s*x + sqrt((s^2-1)*x^2+1); f(0,s)=0. a(n) = f(a(n-1),7/2) + f(a(n-2),7/2). - Marcos Carreira, Dec 27 2006
a(n+1) = 8*a(n) - 8*a(n-1) + a(n-2); a(1)=1, a(2)=8, a(3)=55. - Sture Sjöstedt, May 27 2009
a(n) = A167816(4*n+2). - Reinhard Zumkeller, Nov 13 2009
a(n)=b such that (-1)^n*Integral_{0..Pi/2} (cos((2*n+1)*x))/(3/2-sin(x)) dx = c + b*log(3). - Francesco Daddi, Aug 01 2011
a(n) = A000045(A016825(n)). - Michel Marcus, Mar 22 2015
a(n) = A001906(2*n+1). - R. J. Mathar, Apr 30 2017
E.g.f.: exp(7*x/2)*(5*cosh(3*sqrt(5)*x/2) + 3*sqrt(5)*sinh(3*sqrt(5)*x/2))/5. - Stefano Spezia, Apr 14 2025
From Peter Bala, Jun 08 2025: (Start)
Sum_{n >= 1} (-1)^(n+1)/(a(n) - 1/a(n)) = 1/9 [telescoping series: 3/(a(n) - 1/a(n)) = 1/Fibonacci(4*n+4) + 1/Fibonacci(4*n)].
Product_{n >= 1} (a(n) + 3)/(a(n) - 3) = 5/2 [telescoping product:
(a(n) + 3)/(a(n) - 3) = b(n)/b(n-1), where b(n) = (Lucas(4*n+4) - 3)/(Lucas(4*n+4) + 3)].
Product_{n >= 1} (a(n) + 1)/(a(n) - 1) = sqrt(9/5) [telescoping product:
(a(n) + 1)/(a(n) - 1) = c(n)/c(n-1) for n >= 1, where c(n) = Fibonacci(2*n+2)/Lucas(2*n+2)]. (End)
From Gerry Martens, Aug 20 2025: (Start)
a(n) = ((3 + sqrt(5))^(1 + 2*n) - (3 - sqrt(5))^(1 + 2*n)) / (2^(1 + 2*n)*sqrt(5)).
a(n) = Sum_{k=0..2*n} binomial(2*n + k + 1, 2*k + 1). (End)

A058038 a(n) = Fibonacci(2*n)*Fibonacci(2*n+2).

Original entry on oeis.org

0, 3, 24, 168, 1155, 7920, 54288, 372099, 2550408, 17480760, 119814915, 821223648, 5628750624, 38580030723, 264431464440, 1812440220360, 12422650078083, 85146110326224, 583600122205488, 4000054745112195, 27416783093579880, 187917426909946968
Offset: 0

Views

Author

N. J. A. Sloane, Jun 09 2002

Keywords

Comments

Partial sums of A033888, i.e., a(n) = Sum_{k=0..n} Fibonacci(4*k). - Vladeta Jovovic, Jun 09 2002
From Paul Weisenhorn, May 17 2009: (Start)
a(n) is the solution of the 2 equations a(n)+1=A^2 and 5*a(n)+1=B^2
which are equivalent to the Pell equation (10*a(n)+3)^2-5*(A*B)^2=4.
(End)
Numbers a(n) such as a(n)+1 and 5*a(n)+1 are perfect squares. - Sture Sjöstedt, Nov 03 2011

Examples

			G.f. = 3*x + 24*x^2 + 168*x^3 + 1155*x^4 + 7920*x^5 + 54288*x^6 + ... - _Michael Somos_, Jan 23 2025

References

  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 29.
  • H. J. H. Tuenter, Fibonacci summation identities arising from Catalan's identity, Fib. Q., 60:4 (2022), 312-319.

Links

Crossrefs

Cf. A033888, A004187, A094874.
Bisection of A059929, A064831 and A080097.
Related to sum of fibonacci(kn) over n; cf. A000071, A099919, A027941, A138134, A053606.

Programs

  • Magma
    [Fibonacci(2*n)*Fibonacci(2*n+2): n in [0..30]]; // Vincenzo Librandi, Apr 18 2011
  • Maple
    fs4:=n->sum(fibonacci(4*k),k=0..n):seq(fs4(n),n=0..21); # Gary Detlefs, Dec 07 2010
  • Mathematica
    Table[Fibonacci[2 n]*Fibonacci[2 n + 2], {n, 0, 100}] (* Vladimir Joseph Stephan Orlovsky, Jul 01 2011 *)
Accumulate[Fibonacci[4*Range[0,30]]] (* or *) LinearRecurrence[{8,-8,1},{0,3,24},30] (* Harvey P. Dale, Jul 25 2013 *)
  • PARI
    a(n)=fibonacci(2*n)*fibonacci(2*n+2) \\ Charles R Greathouse IV, Jul 02 2013

Formula

a(n) = -3/5 + (1/5*sqrt(5)+3/5)*(2*1/(7+3*sqrt(5)))^n/(7+3*sqrt(5)) + (1/5*sqrt(5)-3/5)*(-2*1/(-7+3*sqrt(5)))^n/(-7+3*sqrt(5)). Recurrence: a(n) = 8*a(n-1) - 8*a(n-2) + a(n-3). G.f.: 3*x/(1-7*x+x^2)/(1-x). - Vladeta Jovovic, Jun 09 2002
a(n) = A081068(n) - 1.
a(n) is the next integer from ((3+sqrt(5))*((7+3*sqrt(5))/2)^(n-1)-6)/10. - Paul Weisenhorn, May 17 2009
a(n) = 7*a(n-1) - a(n-2) + 3, n>1. - Gary Detlefs, Dec 07 2010
a(n) = sum_{k=0..n} Fibonacci(4k). - Gary Detlefs, Dec 07 2010
a(n) = (Lucas(4n+2)-3)/5, where Lucas(n)= A000032(n). - Gary Detlefs, Dec 07 2010
a(n) = (1/5)*(Fibonacci(4n+4) - Fibonacci(4n)-3). - Gary Detlefs, Dec 08 2010
a(n) = 3*A092521(n). - R. J. Mathar, Nov 03 2011
a(0)=0, a(1)=3, a(2)=24, a(n) = 8*a(n-1) - 8*a(n-2) + a(n-3). - Harvey P. Dale, Jul 25 2013
a(n) = A001906(n)*A001906(n+1). - R. J. Mathar, Jul 09 2019
Sum_{n>=1} 1/a(n) = 2/(3 + sqrt(5)) = A094874 - 1. - Amiram Eldar, Oct 05 2020
a(n) = a(-1-n) for all n in Z. - Michael Somos, Jan 23 2025

A103134 a(n) = Fibonacci(6n+4).

Original entry on oeis.org

3, 55, 987, 17711, 317811, 5702887, 102334155, 1836311903, 32951280099, 591286729879, 10610209857723, 190392490709135, 3416454622906707, 61305790721611591, 1100087778366101931, 19740274219868223167, 354224848179261915075, 6356306993006846248183
Offset: 0

Views

Author

Creighton Dement, Jan 24 2005

Keywords

Comments

Gives those numbers which are Fibonacci numbers in A103135.
Generally, for any sequence where a(0)= Fibonacci(p), a(1) = F(p+q) and Lucas(q)*a(1) +- a(0) = F(p+2q), then a(n) = L(q)*a(n-1) +- a(n-2) generates the following Fibonacci sequence: a(n) = F(q(n)+p). So for this sequence, a(n) = 18*a(n-1) - a(n-2) = F(6n+4): q=6, because 18 is the 6th Lucas number (L(0) = 2, L(1)=1); F(4)=3, F(10)=55 and F(16)=987 (F(0)=0 and F(1)=1). See Lucas sequence A000032. This is a special case where a(0) and a(1) are increasing Fibonacci numbers and Lucas(m)*a(1) +- a(0) is another Fibonacci. - Bob Selcoe, Jul 08 2013
a(n) = x + y where x and y are solutions to x^2 = 5*y^2 - 1. (See related sequences with formula below.) - Richard R. Forberg, Sep 05 2013

Links

Crossrefs

Subsequence of A033887.
Cf. A000032, A000045, A001906, A001519, A015448, A014445, A033888, A033889, A033890, A033891, A049310, A049660, A102312, A099100, A134490, A134491, A134492, A134493, A134494, A134495, A103134, A134497, A134498, A134499, A134500, A134501, A134502, A134503, A134504.
Cf. A103135.

Programs

Formula

G.f.: (x+3)/(x^2-18*x+1).
a(n) = 18*a(n-1) - a(n-2) for n>1; a(0)=3, a(1)=55. - Philippe Deléham, Nov 17 2008
a(n) = A007805(n) + A075796(n), as follows from comment above. - Richard R. Forberg, Sep 05 2013
a(n) = ((15-7*sqrt(5)+(9+4*sqrt(5))^(2*n)*(15+7*sqrt(5))))/(10*(9+4*sqrt(5))^n). - Colin Barker, Jan 24 2016
a(n) = S(3*n+1, 3) = 3*S(n,18) + S(n-1,18), with the Chebyshev S polynomials (A049310), S(-1, x) = 0, and S(n, 18) = A049660(n+1). - Wolfdieter Lang, May 08 2023

Extensions

Edited by N. J. A. Sloane, Aug 10 2010

A215466 Expansion of x*(1-4*x+x^2) / ( (x^2-7*x+1)*(x^2-3*x+1) ).

Original entry on oeis.org

0, 1, 6, 38, 252, 1705, 11628, 79547, 544824, 3733234, 25585230, 175356611, 1201893336, 8237850373, 56462937882, 387002396990, 2652553009008, 18180866487757, 124613506702404, 854113665498719, 5854182112700460
Offset: 0

Views

Author

R. J. Mathar, Aug 11 2012

Keywords

Comments

From Peter Bala, Aug 05 2019: (Start)
Let U(n;P,Q), where P and Q are integer parameters, denote the Lucas sequence of the first kind. Then, excluding the case P = -1, the sequence ( U(n;P,1) + U(2*n;P,1) )/(P + 1) is a fourth-order linear divisibility sequence with o.g.f. x*(1 - 2*(P - 1)*x + x^2)/((1 - P*x + x^2)*(1 - (P^2 - 2)*x + x^2)). This is the case P = 3. See A000027 (P = 2), A165998 (P = -2) and A238536 (P = -3).
More generally, the sequence U(n;P,1) + U(2*n;P,1) + ... + U(k*n;P,1) is a linear divisibility sequence of order 2*k. As an example, see A273625 (P = 3, k = 3 and then sequence normalized with initial term 1). (End)

Links

Crossrefs

Cf. A000032, A000045, A001906, A033888, A165998, A215465, A238536, A273625.

Programs

  • Magma
    I:=[0,1,6,38]; [n le 4 select I[n] else 10*Self(n-1)-23*Self(n-2)+10*Self(n-3)-Self(n-4): n in [1..30]]; // Vincenzo Librandi, Dec 23 2012
  • Magma
    /* By definition: */ m:=20; R:=PowerSeriesRing(Integers(), m); [0] cat Coefficients(R!((1-4*x+x^2)/((x^2-7*x+1)*(x^2-3*x+1)))); // Bruno Berselli, Dec 24 2012
  • Maple
    A215466 := proc(n)
    if type(n,'even') then
        A000032(n)*combinat[fibonacci](3*n)/4 ;
    else
        combinat[fibonacci](n)*A000032(3*n)/4 ;
    end if;
end proc:
  • Mathematica
    CoefficientList[Series[x*(1 - 4*x + x^2)/((x^2 - 7*x + 1)*(x^2 - 3*x + 1)), {x, 0, 40}], x] (* Vincenzo Librandi, Dec 23 2012 *)
LinearRecurrence[{10,-23,10,-1},{0,1,6,38},30] (* Harvey P. Dale, Nov 02 2015 *)
  • PARI
    a(n)=([0,1,0,0; 0,0,1,0; 0,0,0,1; -1,10,-23,10]^n*[0;1;6;38])[1,1] \\ Charles R Greathouse IV, Nov 13 2015
  • PARI
    {a(n) = my(w=quadgen(5)^(2*n)); imag(w^2+w)/4}; /* Michael Somos, Dec 29 2022 */

Formula

a(n) = L(n)*F(3n)/4 if n even, = F(n)*L(3n)/4 if n odd, where L=A000032, F=A000045.
a(n) = 3*A004187(n)/4 + A001906(n)/4.
a(n) = 10*a(n-1) - 23*a(n-2) + 10*a(n-3) - a(n-4), a(0)=0, a(1)=1, a(2)=6, a(3)=38. - Harvey P. Dale, Nov 02 2015
a(n) = (1/4)*(Fibonacci(2*n) + Fibonacci(4*n)) = (1/4)*(A001906(n) + A033888(n)). - Peter Bala, Aug 05 2019
E.g.f.: exp(5*x/2)*(cosh(x)+exp(x)*cosh(sqrt(5)*x))*sinh(sqrt(5)*x/2)/sqrt(5). - Stefano Spezia, Aug 17 2019
a(n) = -a(-n) for all n in Z. - Michael Somos, Dec 29 2022

A134504 a(n) = Fibonacci(7n + 6).

Original entry on oeis.org

8, 233, 6765, 196418, 5702887, 165580141, 4807526976, 139583862445, 4052739537881, 117669030460994, 3416454622906707, 99194853094755497, 2880067194370816120, 83621143489848422977, 2427893228399975082453
Offset: 0

Views

Author

Artur Jasinski, Oct 28 2007

Keywords

Links

Crossrefs

Cf. A000045, A001906, A001519, A033887, A015448, A014445, A033888, A033889, A033890, A033891, A102312, A099100, A134490, A134491, A134492, A134493, A134494, A134495, A103134, A134497, A134498, A134499, A134500, A134501, A134502, A134503, A134504.

Programs

Formula

G.f.: (-8-x) / (-1 + 29*x + x^2). - R. J. Mathar, Jul 04 2011
a(n) = A000045(A017053(n)). - Michel Marcus, Nov 08 2013
a(n) = 29*a(n-1) + a(n-2). - Wesley Ivan Hurt, Mar 15 2023

Extensions

Offset changed from 1 to 0 by Vincenzo Librandi, Apr 17 2011

A221125 T(n,k)=Number of nXk arrays of occupancy after each element stays put or moves to some horizontal or antidiagonal neighbor.

Original entry on oeis.org

1, 3, 1, 8, 21, 1, 21, 200, 144, 1, 55, 1950, 4872, 987, 1, 144, 19005, 179300, 118472, 6765, 1, 377, 185164, 6405243, 16436223, 2880520, 46368, 1, 987, 1803949, 229006550, 2152770692, 1506094707, 70036168, 317811, 1, 2584, 17574675
Offset: 1

Views

Author

R. H. Hardin Jan 02 2013

Keywords

Comments

Table starts
.1........3...........8.............21..............55.............144
.1.......21.........200...........1950...........19005..........185164
.1......144........4872.........179300.........6405243.......229006550
.1......987......118472.......16436223......2152770692....282454178707
.1.....6765.....2880520.....1506094707....723234356389.348270153038959
.1....46368....70036168...138001217828.242959402596630
.1...317811..1702839048.12644783009532
.1..2178309.41402332360
.1.14930352
.1

Examples

			Some solutions for n=3 k=4
..1..1..2..0....2..1..1..0....2..0..2..0....2..1..1..0....1..2..0..0
..0..2..2..0....0..1..2..1....1..1..1..1....0..1..2..0....0..2..2..0
..1..1..0..2....0..0..2..2....0..3..1..0....0..2..3..0....0..1..4..0

Links

Crossrefs

Column 2 is A033888
Row 1 is A001906
Showing 1-10 of 26 results. Next

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.