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.

A001075 a(0) = 1, a(1) = 2, a(n) = 4*a(n-1) - a(n-2).

Original entry on oeis.org

1, 2, 7, 26, 97, 362, 1351, 5042, 18817, 70226, 262087, 978122, 3650401, 13623482, 50843527, 189750626, 708158977, 2642885282, 9863382151, 36810643322, 137379191137, 512706121226, 1913445293767, 7141075053842, 26650854921601, 99462344632562, 371198523608647
Offset: 0

Views

Author

N. J. A. Sloane

Keywords

Comments

Chebyshev's T(n,x) polynomials evaluated at x=2.
x = 2^n - 1 is prime if and only if x divides a(2^(n-2)).
Any k in the sequence is succeeded by 2*k + sqrt{3*(k^2 - 1)}. - Lekraj Beedassy, Jun 28 2002
For all elements x of the sequence, 12*x^2 - 12 is a square. Lim_{n -> infinity} a(n)/a(n-1) = 2 + sqrt(3) = (4 + sqrt(12))/2 which preserves the kinship with the equation "12*x^2 - 12 is a square" where the initial "12" ends up appearing as a square root. - Gregory V. Richardson, Oct 10 2002
This sequence gives the values of x in solutions of the Diophantine equation x^2 - 3*y^2 = 1; the corresponding values of y are in A001353. The solution ratios a(n)/A001353(n) are obtained as convergents of the continued fraction expansion of sqrt(3): either as successive convergents of [2;-4] or as odd convergents of [1;1,2]. - Lekraj Beedassy, Sep 19 2003 [edited by Jon E. Schoenfield, May 04 2014]
a(n) is half the central value in a list of three consecutive integers, the lengths of the sides of a triangle with integer sides and area. - Eugene McDonnell (eemcd(AT)mac.com), Oct 19 2003
a(3+6*k) - 1 and a(3+6*k) + 1 are consecutive odd powerful numbers. See A076445. - T. D. Noe, May 04 2006
The intermediate convergents to 3^(1/2), beginning with 3/2, 12/7, 45/26, 168/97, comprise a strictly increasing sequence; essentially, numerators=A005320, denominators=A001075. - Clark Kimberling, Aug 27 2008
The upper principal convergents to 3^(1/2), beginning with 2/1, 7/4, 26/15, 97/56, comprise a strictly decreasing sequence; numerators=A001075, denominators=A001353. - Clark Kimberling, Aug 27 2008
a(n+1) is the Hankel transform of A000108(n) + A000984(n) = (n+2)*Catalan(n). - Paul Barry, Aug 11 2009
Also, numbers such that floor(a(n)^2/3) is a square: base 3 analog of A031149, A204502, A204514, A204516, A204518, A204520, A004275, A001541. - M. F. Hasler, Jan 15 2012
Pisano period lengths: 1, 2, 2, 4, 3, 2, 8, 4, 6, 6, 10, 4, 12, 8, 6, 8, 18, 6, 5, 12, ... - R. J. Mathar, Aug 10 2012
Except for the first term, positive values of x (or y) satisfying x^2 - 4*x*y + y^2 + 3 = 0. - Colin Barker, Feb 04 2014
Except for the first term, positive values of x (or y) satisfying x^2 - 14*x*y + y^2 + 48 = 0. - Colin Barker, Feb 10 2014
A triangle with row sums generating the sequence can be constructed by taking the production matrix M. Take powers of M, extracting the top rows.
M =
1, 1, 0, 0, 0, 0, ...
2, 0, 3, 0, 0, 0, ...
2, 0, 0, 3, 0, 0, ...
2, 0, 0, 0, 3, 0, ...
2, 0, 0, 0, 0, 3, ...
...
The triangle generated from M is:
1,
1, 1,
3, 1, 3,
11, 3, 3, 9,
41, 11, 9, 9, 27,
...
The left border is A001835 and row sums are (1, 2, 7, 26, 97, ...). - Gary W. Adamson, Jul 25 2016
Even-indexed terms are odd while odd-indexed terms are even. Indeed, a(2*n) = 2*(a(n))^2 - 1 and a(2*n+1) = 2*a(n)*a(n+1) - 2. - Timothy L. Tiffin, Oct 11 2016
For each n, a(0) divides a(n), a(1) divides a(2n+1), a(2) divides a(4*n+2), a(3) divides a(6*n+3), a(4) divides a(8*n+4), a(5) divides a(10n+5), and so on. Thus, a(k) divides a((2*n+1)*k) for each k > 0 and n >= 0. A proof of this can be found in Bhargava-Kedlaya-Ng's first solution to Problem A2 of the 76th Putnam Mathematical Competition. Links to the exam and its solutions can be found below. - Timothy L. Tiffin, Oct 12 2016
From Timothy L. Tiffin, Oct 21 2016: (Start)
If any term a(n) is a prime number, then its index n will be a power of 2. This is a consequence of the results given in the previous two comments. See A277434 for those prime terms.
a(2n) == 1 (mod 6) and a(2*n+1) == 2 (mod 6). Consequently, each odd prime factor of a(n) will be congruent to 1 modulo 6 and, thus, found in A002476.
a(n) == 1 (mod 10) if n == 0 (mod 6), a(n) == 2 (mod 10) if n == {1,-1} (mod 6), a(n) == 7 (mod 10) if n == {2,-2} (mod 6), and a(n) == 6 (mod 10) if n == 3 (mod 6). So, the rightmost digits of a(n) form a repeating cycle of length 6: 1, 2, 7, 6, 7, 2. (End)
a(A298211(n)) = A002350(3*n^2). - A.H.M. Smeets, Jan 25 2018
(2 + sqrt(3))^n = a(n) + A001353(n)*sqrt(3), n >= 0; integers in the quadratic number field Q(sqrt(3)). - Wolfdieter Lang, Feb 16 2018
Yong Hao Ng has shown that for any n, a(n) is coprime with any member of A001834 and with any member of A001835. - René Gy, Feb 26 2018
Positive numbers k such that 3*(k-1)*(k+1) is a square. - Davide Rotondo, Oct 25 2020
a(n)*a(n+1)-1 = a(2*n+1)/2 = A001570(n) divides both a(n)^6+1 and a(n+1)^6+1. In other words, for k = a(2*n+1)/2, (k+1)^6 has divisors congruent to -1 modulo k (cf. A350916). - Max Alekseyev, Jan 23 2022

Examples

			2^6 - 1 = 63 does not divide a(2^4) = 708158977, therefore 63 is composite. 2^5 - 1 = 31 divides a(2^3) = 18817, therefore 31 is prime.
G.f. = 1 + 2*x + 7*x^2 + 26*x^3 + 97*x^4 + 362*x^5 + 1351*x^6 + 5042*x^7 + ...

References

  • Serge Lang, Introduction to Diophantine Approximations, Addison-Wesley, New York, 1966.
  • Eugene McDonnell, "Heron's Rule and Integer-Area Triangles", Vector 12.3 (January 1996) pp. 133-142.
  • 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).
  • P.-F. Teilhet, Reply to Query 2094, L'Intermédiaire des Mathématiciens, 10 (1903), 235-238.

Links

Crossrefs

Bisections are A011943 and A094347.
Cf. A001353, A065918, A071954, A001571, A001834, A003500, A016064, A082840, A026150, A277434.

Programs

  • Haskell
    a001075 n = a001075_list !! n
a001075_list =
   1 : 2 : zipWith (-) (map (4 *) $ tail a001075_list) a001075_list
-- Reinhard Zumkeller, Aug 11 2011
  • Magma
    I:=[1, 2]; [n le 2 select I[n] else 4*Self(n-1) - Self(n-2): n in [1..30]]; // G. C. Greubel, Dec 19 2017
  • Maple
    A001075 := proc(n)
    orthopoly[T](n,2) ;
end proc:
seq(A001075(n),n=0..30) ; # R. J. Mathar, Apr 14 2018
  • Mathematica
    Table[ Ceiling[(1/2)*(2 + Sqrt[3])^n], {n, 0, 24}]
CoefficientList[Series[(1-2*x) / (1-4*x+x^2), {x, 0, 24}], x] (* Jean-François Alcover, Dec 21 2011, after Simon Plouffe *)
LinearRecurrence[{4,-1},{1,2},30] (* Harvey P. Dale, Aug 22 2015 *)
Round@Table[LucasL[2n, Sqrt[2]]/2, {n, 0, 20}] (* Vladimir Reshetnikov, Sep 15 2016 *)
ChebyshevT[Range[0, 20], 2] (* Eric W. Weisstein, May 26 2017 *)
a[ n_] := LucasL[2*n, x]/2 /. x->Sqrt[2]; (* Michael Somos, Sep 05 2022 *)
  • PARI
    {a(n) = subst(poltchebi(abs(n)), x, 2)};
  • PARI
    {a(n) = real((2 + quadgen(12))^abs(n))};
  • PARI
    {a(n) = polsym(1 - 4*x + x^2, abs(n))[1 + abs(n)]/2};
  • PARI
    a(n)=polchebyshev(n,1,2) \\ Charles R Greathouse IV, Nov 07 2016
  • PARI
    my(x='x+O('x^30)); Vec((1-2*x)/(1-4*x+x^2)) \\ G. C. Greubel, Dec 19 2017
  • SageMath
    [lucas_number2(n,4,1)/2 for n in range(0, 25)] # Zerinvary Lajos, May 14 2009
  • SageMath
    def a(n):
    Q = QuadraticField(3, 't')
    u = Q.units()[0]
    return (u^n).lift().coeffs()[0]  # Ralf Stephan, Jun 19 2014

Formula

G.f.: (1 - 2*x)/(1 - 4*x + x^2). - Simon Plouffe in his 1992 dissertation
E.g.f.: exp(2*x)*cosh(sqrt(3)*x).
a(n) = 4*a(n-1) - a(n-2) = a(-n).
a(n) = (S(n, 4) - S(n-2, 4))/2 = T(n, 2), with S(n, x) := U(n, x/2), S(-1, x) := 0, S(-2, x) := -1. U, resp. T, are Chebyshev's polynomials of the second, resp. first, kind. S(n-1, 4) = A001353(n), n >= 0. See A049310 and A053120.
a(n) = A001353(n+2) - 2*A001353(n+1).
a(n) = sqrt(1 + 3*A001353(n)) (cf. Richardson comment, Oct 10 2002).
a(n) = 2^(-n)*Sum_{k>=0} binomial(2*n, 2*k)*3^k = 2^(-n)*Sum_{k>=0} A086645(n, k)*3^k. - Philippe Deléham, Mar 01 2004
a(n) = ((2 + sqrt(3))^n + (2 - sqrt(3))^n)/2; a(n) = ceiling((1/2)*(2 + sqrt(3))^(n)).
a(n) = cosh(n * log(2 + sqrt(3))).
a(n) = Sum_{k=0..floor(n/2)} binomial(n, 2*k)*2^(n-2*k)*3^k. - Paul Barry, May 08 2003
a(n+2) = 2*a(n+1) + 3*Sum_{k>=0} a(n-k)*2^k. - Philippe Deléham, Mar 03 2004
a(n) = 2*a(n-1) + 3*A001353(n-1). - Lekraj Beedassy, Jul 21 2006
a(n) = left term of M^n * [1,0] where M = the 2 X 2 matrix [2,3; 1,2]. Right term = A001353(n). Example: a(4) = 97 since M^4 * [1,0] = [A001075(4), A001353(4)] = [97, 56]. - Gary W. Adamson, Dec 27 2006
Binomial transform of A026150: (1, 1, 4, 10, 28, 76, ...). - Gary W. Adamson, Nov 23 2007
First differences of A001571. - N. J. A. Sloane, Nov 03 2009
Sequence satisfies -3 = f(a(n), a(n+1)) where f(u, v) = u^2 + v^2 - 4*u*v. - Michael Somos, Sep 19 2008
a(n) = Sum_{k=0..n} A201730(n,k)*2^k. - Philippe Deléham, Dec 06 2011
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - x*(3*k - 4)/(x*(3*k - 1) - 2/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 28 2013
a(n) = Sum_{k=0..n} A238731(n,k). - Philippe Deléham, Mar 05 2014
a(n) = (-1)^n*(A125905(n) + 2*A125905(n-1)), n > 0. - Franck Maminirina Ramaharo, Nov 11 2018
a(n) = (tan(Pi/12)^n + tan(5*Pi/12)^n)/2. - Greg Dresden, Oct 01 2020
From Peter Bala, Aug 17 2022: (Start)
a(n) = (1/2)^n * [x^n] ( 4*x + sqrt(1 + 12*x^2) )^n.
The g.f. A(x) satisfies A(2*x) = 1 + x*B'(x)/B(x), where B(x) = 1/sqrt(1 - 8*x + 4*x^2) is the g.f. of A069835.
The Gauss congruences a(n*p^k) == a(n*p^(k-1)) (mod p^k) hold for all primes p >= 3 and positive integers n and k.
Sum_{n >= 1} 1/(a(n) - (3/2)/a(n)) = 1.
Sum_{n >= 1} (-1)^(n+1)/(a(n) + (1/2)/a(n)) = 1/3.
Sum_{n >= 1} 1/(a(n)^2 - 3/2) = 1 - 1/sqrt(3). (End)
a(n) = binomial(2*n, n) + 2*Sum_{k > 0} binomial(2*n, n+2*k)*cos(k*Pi/3). - Greg Dresden, Oct 11 2022
2*a(n) + 2^n = 3*Sum_{k=-n..n} (-1)^k*binomial(2*n, n+6*k). - Greg Dresden, Feb 07 2023

Extensions

More terms from James Sellers, Jul 10 2000
Chebyshev comments from Wolfdieter Lang, Oct 31 2002

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

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