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 10 results.

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

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)
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.

Crossrefs

Fibonacci A000045 = union of this sequence and A001519.
Inverse sequences A130259 and A130260.

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

A019952 Decimal expansion of tangent of 54 degrees.

Original entry on oeis.org

1, 3, 7, 6, 3, 8, 1, 9, 2, 0, 4, 7, 1, 1, 7, 3, 5, 3, 8, 2, 0, 7, 2, 0, 9, 5, 8, 1, 9, 1, 0, 8, 8, 7, 6, 7, 9, 5, 2, 5, 8, 9, 9, 3, 3, 6, 0, 0, 8, 1, 5, 8, 6, 6, 3, 3, 6, 5, 6, 7, 5, 7, 6, 5, 6, 1, 9, 0, 9, 5, 1, 9, 3, 7, 6, 7, 1, 7, 2, 9, 8, 5, 0, 6, 5, 9, 5, 2, 9, 9, 3, 1, 1, 0, 0, 7, 0, 1, 9
Offset: 1

Views

Author

Keywords

Comments

Also the decimal expansion of cotangent of 36 degrees. - Mohammad K. Azarian, Jun 30 2013
A quartic number with denominator 5. - Charles R Greathouse IV, Aug 27 2017
Conjecture: Product (2/3) * (8/7) * (12/13) * (18/17) * (22/23) * (32/33) * ... * (a_n/b_n) = sqrt(25 + 10*sqrt(5))/5 = tan(3*Pi/10) = A019952, where a_n even, a_n + b_n = a(n), |a_n - b_n| = 1, n >= 0. - Dimitris Valianatos, Feb 14 2020
Also the limiting value of the distance between the lines F(n)*x + F(n+1)*y = 0 and F(n)*x + F(n+1)*y = F(n+2) (where F(n)=A000045(n) are the Fibonacci numbers and n>0). - Burak Muslu, Apr 03 2021
Decimal expansion of the radius of an inscribed sphere in a rhombic triacontahedron with unit edge length. - Wesley Ivan Hurt, May 11 2021

Examples

			1.376381920471173538207209581910887679525899336...
		

Crossrefs

Cf. A344171 (rhombic triacontahedron surface area).
Cf. A344172 (rhombic triacontahedron volume).
Cf. A344212 (rhombic triacontahedron midradius).

Programs

  • Magma
    SetDefaultRealField(RealField(100)); R:= RealField(); Tan(3*Pi(R)/10); // G. C. Greubel, Nov 22 2018
    
  • Maple
    Digits:=100: evalf(tan(3*Pi/10)); # Wesley Ivan Hurt, Oct 07 2014
  • Mathematica
    RealDigits[Tan[3*Pi/10], 10, 100][[1]] (* Wesley Ivan Hurt, Oct 07 2014 *)
    RealDigits[Tan[54 Degree],10,120][[1]] (* Harvey P. Dale, Jul 16 2016 *)
  • PARI
    tan(3*Pi/10) \\ Charles R Greathouse IV, Aug 27 2017
    
  • Python
    from sympy import sqrt
    [print(i, end=', ') for i in str(sqrt(1+2/sqrt(5)).n(110)) if i!='.'] # Karl V. Keller, Jr., Jun 19 2020
  • Sage
    numerical_approx(tan(3*pi/10), digits=100) # G. C. Greubel, Nov 22 2018
    

Formula

Equals A019863/A019845 = 1/A019934. - R. J. Mathar, Jul 26 2010
The largest positive solution of cos(4*arctan(1/x)) = cos(6*arctan(1/x)). - Thomas Olson, Oct 03 2014
Equals sqrt(25 + 10*sqrt(5))/5. - G. C. Greubel, Nov 22 2018
Equals sqrt(2 + sqrt(5))/5^(1/4). - Burak Muslu, Apr 03 2021
From Wesley Ivan Hurt, May 11 2021: (Start)
Equals phi^2/sqrt(1+phi^2) where phi is the golden ratio.
Equals sqrt(1+2/sqrt(5)). (End)
Equals Product_{k>=1} (1 - (-1)^k/A090772(k)). - Amiram Eldar, Nov 23 2024
Equals 2*A375067. - Hugo Pfoertner, Nov 23 2024

A187798 Decimal expansion of (3-phi)/2, where phi is the golden ratio.

Original entry on oeis.org

6, 9, 0, 9, 8, 3, 0, 0, 5, 6, 2, 5, 0, 5, 2, 5, 7, 5, 8, 9, 7, 7, 0, 6, 5, 8, 2, 8, 1, 7, 1, 8, 0, 9, 4, 1, 1, 3, 9, 8, 4, 5, 4, 1, 0, 0, 9, 7, 1, 1, 8, 5, 6, 8, 9, 3, 2, 2, 7, 5, 6, 8, 8, 6, 4, 7, 3, 6, 9, 7, 6, 8, 5, 9, 0, 5, 4, 8, 7, 7, 5, 1, 4, 6, 3, 9, 6, 3, 9, 7, 9, 0, 5, 3, 0, 4, 4, 3, 1, 2, 5, 7, 6, 2, 2
Offset: 0

Views

Author

Joost Gielen, Aug 30 2013

Keywords

Comments

This is the height h of the isosceles triangle in a regular pentagon inscribed in the unit circle formed from a diagonal as base and two adjacent pentagon sides. h = sqrt(sqrt(3-phi)^2 - (sqrt(2 + phi)/2)^2) = sqrt(10 - 5*phi)/2 = (3 - phi)/2. - Wolfdieter Lang, Jan 07 2018

Examples

			0.6909830056250525758977065828171809411398454100971185689322756886473697685905...
		

Crossrefs

Programs

  • Mathematica
    RealDigits[(3 - GoldenRatio)/2, 10, 111][[1]] (* or *)
    RealDigits[(5 - Sqrt[5])/4, 10, 111][[1]] (* Robert G. Wilson v, Jan 07 2018 *)
  • PARI
    (5-sqrt(5))/4 \\ Charles R Greathouse IV, Aug 31 2013

Formula

Equals (3-phi)/2 = A094874/2 with phi from A001622.
From Amiram Eldar, Nov 28 2024: (Start)
Equals 1/A344212.
Equals Product_{k>=0} (1 - 1/A081011(k)). (End)

Extensions

Extended by Charles R Greathouse IV, Aug 31 2013

A242671 Decimal expansion of k2, a Diophantine approximation constant such that the area of the "critical parallelogram" (in this case a square) is 4*k2.

Original entry on oeis.org

7, 2, 3, 6, 0, 6, 7, 9, 7, 7, 4, 9, 9, 7, 8, 9, 6, 9, 6, 4, 0, 9, 1, 7, 3, 6, 6, 8, 7, 3, 1, 2, 7, 6, 2, 3, 5, 4, 4, 0, 6, 1, 8, 3, 5, 9, 6, 1, 1, 5, 2, 5, 7, 2, 4, 2, 7, 0, 8, 9, 7, 2, 4, 5, 4, 1, 0, 5, 2, 0, 9, 2, 5, 6, 3, 7, 8, 0, 4, 8, 9, 9, 4, 1, 4, 4, 1, 4, 4, 0, 8, 3, 7, 8, 7, 8, 2, 2, 7, 4
Offset: 0

Views

Author

Jean-François Alcover, May 20 2014

Keywords

Comments

Quoting Steven Finch: "The slopes of the 'critical parallelogram' are (1+sqrt(5))/2 [phi] and (1-sqrt(5))/2 [-1/phi]."
Essentially the same as A229780, A134972, A134945, A098317 and A002163. - R. J. Mathar, May 23 2014
Let W_n be the collection of all binary words of length n that do not contain two consecutive 0's. Let r_n be the ratio of the total number of 1's in W_n divided by the total number of letters in W_n. Then lim_{n->oo} r_n = 0.723606... Equivalently, lim_{n->oo} A004798(n)/(n*A000045(n+2)) = 0.723606... - Geoffrey Critzer, Feb 04 2022
The limiting frequency of the digit 0 in the base phi representation of real numbers in the range [0,1], where phi is the golden ratio (A001622) (Rényi, 1957). - Amiram Eldar, Mar 18 2025

Examples

			k2 = 0.723606797749978969640917366873127623544...
		

References

  • Steven R. Finch, Mathematical Constants, Cambridge University Press, 2003, Section 2.23, p. 176.

Crossrefs

Programs

  • Mathematica
    RealDigits[(1+1/Sqrt[5])/2, 10, 100] // First
  • PARI
    (1 + 1/sqrt(5))/2 \\ Stefano Spezia, Dec 07 2024

Formula

Equals (1 + 1/sqrt(5))/2.
Equals 1/A094874. - Michel Marcus, Dec 01 2018
From Amiram Eldar, Feb 11 2022: (Start)
Equals phi/sqrt(5), where phi is the golden ratio (A001622).
Equals lim_{k->oo} Fibonacci(k+1)/Lucas(k). (End)
From Amiram Eldar, Nov 28 2024: (Start)
Equals A344212/2 = A296184/5 = A300074^2 = sqrt(A229780).
Equals Product_{k>=1} (1 - 1/A081007(k)). (End)
Equals 1 - A244847. - Amiram Eldar, Mar 18 2025

A344171 Decimal expansion of 12*sqrt(5).

Original entry on oeis.org

2, 6, 8, 3, 2, 8, 1, 5, 7, 2, 9, 9, 9, 7, 4, 7, 6, 3, 5, 6, 9, 1, 0, 0, 8, 4, 0, 2, 4, 7, 7, 5, 3, 1, 4, 8, 2, 5, 2, 8, 7, 4, 2, 0, 3, 1, 5, 3, 3, 8, 3, 0, 8, 6, 9, 1, 2, 5, 0, 7, 6, 6, 9, 4, 4, 9, 2, 6, 2, 5, 1, 1, 0, 7, 6, 5, 3, 6, 5, 8, 7, 9, 2, 9, 7, 2, 9, 7, 2, 9, 0, 0, 5
Offset: 2

Views

Author

Wesley Ivan Hurt, May 10 2021

Keywords

Comments

Decimal expansion of the surface area of a rhombic triacontahedron with unit edge length.

Examples

			26.83281572999747635691008...
		

Crossrefs

Cf. A344172 (rhombic triacontahedron volume).
Cf. A344212 (rhombic triacontahedron midradius).
Cf. A019952 (rhombic triacontahedron radius of inscribed sphere).

Programs

  • Mathematica
    RealDigits[12 Sqrt[5], 10, 100][[1]] // Flatten

A344172 Decimal expansion of 4*sqrt(5 + 2*sqrt(5)).

Original entry on oeis.org

1, 2, 3, 1, 0, 7, 3, 4, 1, 4, 8, 7, 0, 1, 0, 1, 3, 6, 1, 0, 2, 8, 1, 1, 6, 2, 3, 0, 4, 1, 4, 7, 6, 3, 9, 2, 9, 6, 0, 2, 6, 8, 0, 8, 5, 7, 4, 1, 5, 1, 1, 6, 9, 7, 0, 8, 1, 5, 6, 6, 2, 5, 0, 0, 1, 4, 9, 9, 4, 5, 3, 1, 5, 3, 9, 8, 0, 3, 6, 3, 6, 7, 3, 8, 1, 8, 3, 7, 4, 8, 8, 6, 6
Offset: 2

Views

Author

Wesley Ivan Hurt, May 10 2021

Keywords

Comments

Decimal expansion of the volume of a rhombic triacontahedron with unit edge length.

Examples

			12.310734148701013610281...
		

Crossrefs

Cf. A344171 (rhombic triacontahedron surface area).
Cf. A344212 (rhombic triacontahedron midradius).
Cf. A019952 (rhombic triacontahedron radius of inscribed sphere).

Programs

  • Mathematica
    RealDigits[4 Sqrt[5 + 2 Sqrt[5]], 10, 100][[1]] // Flatten

A127693 Expansion of psi(x^2) + x * psi(x^10) in powers of x where psi() is a Ramanujan theta function.

Original entry on oeis.org

1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0
Offset: 0

Views

Author

Michael Somos, Jan 19 2007

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Examples

			G.f. = 1 + x + x^2 + x^6 + x^11 + x^12 + x^20 + x^30 + x^31 + x^42 + x^56 + x^61 + ...
G.f. = q + q^5 + q^9 + q^25 + q^45 + q^49 + q^81 + q^121 + q^125 + q^169 + q^225 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ (EllipticTheta[ 2, 0, x] + EllipticTheta[ 2, 0, x^5]) / (2 x^(1/4)), {x, 0, n}]; (* Michael Somos, Jul 08 2015 *)
  • PARI
    {a(n) = issquare(4*n + 1) + issquare(20*n + 5)};

Formula

Expansion of f(-x^2, -x^3) * f(x, -x^4) / f(-x^2, -x^2) = f(x^2, -x^3) * f(x, x^4) / f(-x^10, -x^10) where f(,) is Ramanujan's general theta function. - Michael Somos, Jul 30 2012
Euler transform of period 20 sequence [ 1, 0, -1, 0, 0, 1, -1, 0, 1, -1, 1, 0, -1, 1, 0, 0, -1, 0, 1, -1, ...].
a(n) = b(4*n + 1) where b() is multiplicative and b(2^e) = 0^e, b(5^e) = 1, else b(p^e) = (1 + (-1)^e) / 2.
a(9*n + 2) = a(5*n + 1) = a(n), a(5*n + 3) = a(5*n + 4) = a(6*n + 3) = a(6*n + 4) = a(9*n + 5) = a(9*n + 8) = 0.
G.f.: Sum_{k>0} x^(k*(k - 1)) + x^(5*k*(k - 1) + 1).
G.f.: Product_{k>0} (1 - x^(10*k)) * (1 + x^(10*k - 1)) * (1 + x^(10*k-2)) * (1 - x^(10*k - 3)) * (1 + x^(10*k - 4)) * (1 + x^(10*k - 6)) * (1 - x^(10*k - 7)) * (1 + x^(10*k -8)) * (1 + x^(10*k - 9)).
a(2*n) = A010054(n). a(3*n) = A089806(n). a(6*n) = A080995(n).
Sum_{k=1..n} a(k) ~ c * sqrt(n), where c = 1 + 1/sqrt(5) = 1.447213... (A344212). - Amiram Eldar, Dec 29 2023

A214284 Characteristic function of squares or five times squares.

Original entry on oeis.org

1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Michael Somos, Jul 09 2012

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
A195198 is a similar sequence except with three instead of five. - Michael Somos, Oct 22 2017

Examples

			G.f. = 1 + x + x^4 + x^5 + x^9 + x^16 + x^20 + x^25 + x^36 + x^45 + x^49 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ Series[ (EllipticTheta[ 3, 0, q] + EllipticTheta[ 3, 0, q^5]) / 2, {q, 0, n}], {q, 0, n}];
    a[ n_] := If[ n < 0, 0, Boole[ OddQ [ Length @ Divisors @ n] || OddQ [ Length @ Divisors[5 n]]]];
  • PARI
    {a(n) = issquare(n) || issquare(5*n)};
    
  • PARI
    {a(n) = if( n<1, n==0, direuler( p=2, n, if( p==5, 1 + X, 1) / (1 - X^2))[n])};

Formula

Expansion of f(q, q^9) * f(-q^8, -q^12) / f(-q^4, -q^16) in powers of q where f(, ) is Ramanujan's general theta function.
Expansion of f(q^3, q^7) * f(-q^2, -q^3) / f(-q, -q^4) in powers of q where f(, ) is Ramanujan's general theta function.
Euler transform of period 20 sequence [1, -1, 0, 1, 0, 0, 0, -1, 1, -1, 1, -1, 0, 0, 0, 1, 0, -1, 1, -1, ...].
a(n) is multiplicative with a(0) = a(5^e) = 1, a(p^e) = 1 if e is even, 0 otherwise.
G.f.: (theta_3(q) + theta_3(q^5)) / 2 = 1 + (Sum_{k>0} x^(k^2) + x^(5*k^2)).
Dirichlet g.f.: zeta(2*s) * (1 + 5^-s).
a(4*n + 2) = a(4*n + 3) = 0. a(4*n + 1) = A127693(n). a(5*n) = a(n).
Sum_{k=0..n} a(k) ~ c * sqrt(n), where c = 1+1/sqrt(5) = 1.447213... (A344212). - Amiram Eldar, Sep 14 2023

A081005 a(n) = Fibonacci(4n+3) + 1, or Fibonacci(2n+1)*Lucas(2n+2).

Original entry on oeis.org

3, 14, 90, 611, 4182, 28658, 196419, 1346270, 9227466, 63245987, 433494438, 2971215074, 20365011075, 139583862446, 956722026042, 6557470319843, 44945570212854, 308061521170130, 2111485077978051, 14472334024676222, 99194853094755498, 679891637638612259
Offset: 0

Views

Author

R. K. Guy, Mar 01 2003

Keywords

References

  • Hugh C. Williams, Edouard Lucas and Primality Testing, John Wiley and Sons, 1998, p. 75.

Crossrefs

Cf. A000045 (Fibonacci numbers), A000032 (Lucas numbers), A344212.

Programs

  • GAP
    List([0..30], n-> Fibonacci(4*n+3)+1); # G. C. Greubel, Jul 15 2019
  • Magma
    [Fibonacci(4*n+3)+1: n in [0..30]]; // Vincenzo Librandi, Apr 15 2011
    
  • Maple
    with(combinat): for n from 0 to 30 do printf(`%d,`,fibonacci(4*n+3)+1) od: # James Sellers, Mar 03 2003
  • Mathematica
    Fibonacci[4Range[0,30]+3]+1 (* or *) LinearRecurrence[{8,-8,1}, {3,14,90}, 30] (* Harvey P. Dale, Jan 02 2013 *)
  • PARI
    vector(30, n, n--; fibonacci(4*n+3)+1) \\ G. C. Greubel, Jul 15 2019
    
  • Sage
    [fibonacci(4*n+3)+1 for n in (0..30)] # G. C. Greubel, Jul 15 2019
    

Formula

a(n) = 8*a(n-1) - 8*a(n-2) + a(n-3).
G.f.: (3-10*x+2*x^2)/((1-x)*(1-7*x+x^2)). - Colin Barker, Jun 24 2012
Product_{n>=0} (1 + 1/a(n)) = 1 + 1/sqrt(5) = A344212. - Amiram Eldar, Nov 28 2024

Extensions

More terms from James Sellers, Mar 03 2003

A049657 a(n) = (F(8*n+3) - 2)/3, where F = A000045 (the Fibonacci sequence).

Original entry on oeis.org

0, 29, 1393, 65472, 3075821, 144498145, 6788337024, 318907342013, 14981856737617, 703828359326016, 33064951031585165, 1553348870125176769, 72974331944851723008, 3428240252537905804637, 161054317537336721094961
Offset: 0

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    [(Fibonacci(8*n+3) - 2)/3: n in [0..30]]; // G. C. Greubel, Dec 02 2017
  • Mathematica
    LinearRecurrence[{48, -48, 1}, {0, 29, 1393}, 50] (* or *) Table[( Fibonacci[8*n+3] - 2)/3, {n,0,30}] (* G. C. Greubel, Dec 02 2017 *)
  • PARI
    for(n=0,30, print1((fibonacci(8*n+3) - 2)/3, ", ")) \\ G. C. Greubel, Dec 02 2017
    

Formula

G.f.: -x*(29+x) / ( (x-1)*(x^2-47*x+1) ).
Product_{n>=1} (1 - 1/a(n)) = 2*(1+1/sqrt(5))/3 = (2/3) * A344212. - Amiram Eldar, Nov 28 2024

Extensions

Description corrected by and more terms from Michael Somos
Showing 1-10 of 10 results.