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

A107920 Lucas and Lehmer numbers with parameters (1 +- sqrt(-7))/2.

Original entry on oeis.org

0, 1, 1, -1, -3, -1, 5, 7, -3, -17, -11, 23, 45, -1, -91, -89, 93, 271, 85, -457, -627, 287, 1541, 967, -2115, -4049, 181, 8279, 7917, -8641, -24475, -7193, 41757, 56143, -27371, -139657, -84915, 194399, 364229, -24569, -753027, -703889, 802165, 2209943, 605613, -3814273
Offset: 0

Views

Author

Michael Somos, May 28 2005

Keywords

Comments

The sequences A001607, A077020, A107920, A167433, A169998 are all essentially the same except for signs.
This is an example of a sequence of Lehmer numbers. In this case, the two parameters, alpha and beta, are (1 +- i*sqrt(7))/2. Bilu, Hanrot, Voutier and Mignotte show that all terms of a Lehmer sequence a(n) have a primitive factor for n > 30. Note that for this sequence, a(30) = 24475 = 5*5*11*89 has no primitive factors. - T. D. Noe, Oct 29 2003
Row sums of Riordan array (1/(1+2*x^2), x/(1+2*x^2)). - Paul Barry, Sep 10 2005
Pisano period lengths: 1, 1, 8, 2, 24, 8, 21, 2, 24, 24, 10, 8, 168, 21, 24, 4, 144, 24, 360, 24, ... - R. J. Mathar, Aug 10 2012
This is the Lucas Sequence U_n(P, Q) = U_n(1, 2). V_n(1, 2) = A002249(n). - Raphie Frank, Dec 25 2013
Note that (A002249(n)/2)^2 + 7*(a(n)/2)^2 = 2^n for all n in N. This is a specific case of the Lucas sequence identity (V_n/2)^2 - D*(U_n/2)^2 = Q^n where V_n = (a^n + b^n), U_n = (a^n - b^n)/(a - b), Q = (a*b) = 2 and D = (a - b)^2 = -7; a = (1 + sqrt(-7))/2 and b = (1 - sqrt(-7))/2. - Raphie Frank, Nov 26 2015
For the special case where |a(n)| = 1, true for n if and only if n is in {1, 2, 3, 5, 13} = {A215795(n) + 1} = {A060728(n) - 2}, then, additionally, by the Lucas sequence identity (U_2n = U_n*V_n), we have (a(2n)/2)^2 + 7*(a(n)/2)^2 = 2^n. - Raphie Frank, Nov 26 2015

Examples

			G.f. = x + x^2 - x^3 - 3*x^4 - x^5 + 5*x^6 + 7*x^7 - 3*x^8 - 17*x^9 - 11*x^10 + ...
		

Crossrefs

Programs

  • Magma
    [0] cat [n le 2 select 1 else Self(n-1)-2*Self(n-2): n in [1..45]]; // Vincenzo Librandi, Nov 27 2015
  • Maple
    a:= n-> (Matrix([[1,1],[ -2,0]])^n)[1,2]: seq(a(n), n=0..45); # Alois P. Heinz, Sep 03 2008
  • Mathematica
    LinearRecurrence[{1, -2}, {0, 1}, 50] (* Vincenzo Librandi, Nov 27 2015 *)
    a[ n_] := Im[ ((1 + Sqrt[-7]) / 2)^n // FullSimplify] 2 / Sqrt[7]; (* Michael Somos, Jan 19 2017 *)
    a[n_] := If[n < 2, n, Hypergeometric2F1[1 - n/2, (1 - n)/2, 1 - n, 8]];
    Table[a[n], {n, 0, 45}] (* Peter Luschny, Feb 23 2018 *)
  • PARI
    {a(n) = imag(quadgen(-7)^n)};
    
  • PARI
    my(x='x+O('x^100)); concat(0, Vec(x/(1-x+2*x^2))) \\ Altug Alkan, Dec 04 2015
    
  • Sage
    [lucas_number1(n,1,2) for n in range(0, 46)] # Zerinvary Lajos, Apr 22 2009
    

Formula

G.f.: x / (1 - x + 2*x^2).
a(n) = a(n-1) - 2*a(n-2).
a(n) = -(-1)^n*A001607(n).
From Paul Barry, Sep 10 2005: (Start)
a(n+1) = Sum_{k=0..n} C((n+k)/2, k)*(-2)^((n-k)/2)*(1+(-1)^(n-k))/2.
a(n+1) = Sum_{k=0..floor(n/2)} C(n-k, k)*(-2)^k. (End)
a(n+1) = Sum_{k=0..n} A109466(n,k)*2^(n-k). - Philippe Deléham, Oct 26 2008
a(n) = ((1 - i*sqrt(7))^n - (1 + i*sqrt(7))^n)*i/(2^n*sqrt(7)), where i=sqrt(-1). - Bruno Berselli, Jul 01 2011
(a(2*(A060728(n)) - 4))^2 = (A002249(A060728(n) - 2))^2 = 2^(A060728(n)) - 7 = A227078(n), the Ramanujan-Nagell squares. - Raphie Frank, Dec 25 2013
a(n) = -a(-n) * 2^n for all n in Z. - Michael Somos, Jan 19 2017
G.f.: x / (1 - x / (1 + 2*x / (1 - 2*x))). - Michael Somos, Jan 19 2017
a(n) = S(n-1, 1/sqrt(2))*(sqrt(2))^(n-1), n >= 0, with the Chebyshev S polynomials (coefficients in A049310), and S(-1, x) = 0. - Wolfdieter Lang, Feb 22 2018
a(n) = hypergeom([1-n/2, (1-n)/2], [1-n], 8) for n >= 2. - Peter Luschny, Feb 23 2018

A060728 Numbers n such that Ramanujan's equation x^2 + 7 = 2^n has an integer solution.

Original entry on oeis.org

3, 4, 5, 7, 15
Offset: 1

Views

Author

Lekraj Beedassy, Apr 25 2001

Keywords

Comments

See A038198 for corresponding x. - Lekraj Beedassy, Sep 07 2004
Also numbers such that 2^(n-3)-1 is in A000217, i.e., a triangular number. - M. F. Hasler, Feb 23 2009
With respect to M. F. Hasler's comment above, all terms 2^(n-3) - 1 are known as the Ramanujan-Nagell triangular numbers (A076046). - Raphie Frank, Mar 31 2013
Interestingly enough, all the solutions correspond to noncomposite x, i.e., x = 1 for the first term, and primes 3, 5, 11, 181 for the following terms. - M. F. Hasler, Mar 11 2024

Examples

			The fifth and ultimate solution to Ramanujan's equation is obtained for the 15th power of 2, so that we have x^2 + 7 = 2^15 yielding x = 181.
		

References

  • J.-M. De Koninck, Ces nombres qui nous fascinent, Entry 181, p. 56, Ellipses, Paris 2008.
  • J. Roberts, Lure of the Integers. pp. 90-91, MAA 1992.
  • Ian Stewart & David Tall, Algebraic Number Theory and Fermat's Last Theorem, 3rd Ed. Natick, Massachusetts (2002): 96-98.

Crossrefs

Programs

  • Magma
    [n: n in [0..100] | IsSquare(2^n-7)]; // Vincenzo Librandi, Jan 07 2014
    
  • Mathematica
    ramaNagell[n_] := Reduce[x^2 + 7 == 2^n, x, Integers] =!= False; Select[ Range[100], ramaNagell] (* Jean-François Alcover, Sep 21 2011 *)
  • PARI
    is(n)=issquare(2^n-7) \\ Anders Hellström, Dec 12 2015

Formula

a(n) = log_2(8*A076046(n) + 8) = log_2(A227078(n) + 7)
Empirically, a(n) = Fibonacci(c + 1) + 2 = ceiling[e^((c - 1)/2)] + 2 where {c} is the complete set of positive solutions to {n in N | 2 cos(2*Pi/n) is in Z}; c is in {1,2,3,4,6} (see A217290).

Extensions

Added keyword "full", M. F. Hasler, Feb 23 2009

A227078 The Ramanujan-Nagell squares: A038198(n)^2.

Original entry on oeis.org

1, 9, 25, 121, 32761
Offset: 0

Views

Author

Raphie Frank, Jun 30 2013

Keywords

Comments

a(n) = (2*x - 1)^2 = (sqrt(2)*sqrt(sqrt(6*y^2 - 5) + 1) - 1)^2 = 2^(z + 3) - 7 for x, y, z are the solutions to two Diophantine equations noted by R. K. Guy: 2*x^2*(x^2 - 1) = 3*(y^2 - 1) & x*(x - 1)/2 = 2^z - 1 (see A180445). x = {1, 2, 3, 6, 91} = A180445(n), y = {1, 3, 7, 29, 6761} = A227077(n), and z = {0, 1, 2, 4, 12} = A215795(n).

References

  • J.-M. De Koninck, Ces nombres qui nous fascinent, Entry 181, p. 56, Ellipses, Paris 2008.
  • L. J. Mordell, Diophantine Equations, Academic Press, NY, 1969, p. 205.

Crossrefs

Formula

a(n) + 7 = 2^A060728(n).
(a(n) - 1)/8 = A076046(n).

A180445 All positive solutions, x, for each of two Diophantine equations noted by Richard K. Guy.

Original entry on oeis.org

1, 2, 3, 6, 91
Offset: 1

Views

Author

Jonathan Vos Post, Sep 05 2010

Keywords

Comments

2*(x^2)*((x^2)-1) = 3*((y^2)-1) has only these five positive solutions.
x*(x-1)/2 = (2^z)-1 has only these five positive solutions.
Richard K. Guy notes, as Example 29: "True, but why the coincidence?"
Algebraically, y solutions = {1, 3, 7, 29, 6761} can be derived from x solutions as follows: y = sqrt(((2*x^2 - 1)^2 + 5)/6). From this relationship it becomes clear that the form (((2*x^2 - 1)^2 + 5)/6) can only be an integer square for x is in {1, 2, 3, 6, 91}. Thus, x and y solutions are also unique integer solutions to the following equivalency: (2x^2 - 1)^2 = 6y^2 - 5. From this relationship the following statement naturally follows: ((sqrt(6*y^2 - 5) + 1)/2 - sqrt((sqrt(6*(y^2) - 5) + 1)/2))/2 = (2^z - 1) = {0, 1, 3, 15, 4095} = A076046(n), the Ramanujan-Nagell triangular numbers; z = {0, 1, 2, 4, 12} = (A060728(n) - 3). - Raphie Frank, Jun 26 2013

Crossrefs

Formula

x = sqrt((sqrt(6*(y^2) - 5) + 1)/2) = (sqrt(2^(z + 3) - 7) + 1)/2; y = {1, 3, 7, 29, 6761} and z = (A060728(n) - 3) = A215795(n) = {0, 1, 2, 4, 12}. - Raphie Frank, Jun 23 2013

A328129 Finite cardinalities of equivalence classes of real intervals with respect to the symmetric transitive closure of R(x,y) = "x is an integer multiple of y".

Original entry on oeis.org

1, 2, 4, 7, 12, 24, 41, 58, 75, 92, 109, 214, 319, 424, 529, 634, 1176, 1718, 2260, 2802, 3344, 3886, 4428, 4970, 9411
Offset: 1

Views

Author

Peter Munn, Oct 04 2019

Keywords

Comments

Start with a point marked on an axis at unit distance from the origin. How many distinct points can be marked on the axis, if the offset from the origin of each additional point is a positive integer multiple or submultiple of that of a point already marked, and no less than unit distance nor greater than a given maximum, d? The possible finite answers, which occur only if d < 27/5, are the numbers in this sequence.
If we think of an interval with lower bound 1 as growing in time, equivalence classes of cardinality a(n) appear when the interval reaches [1, r_min(n)] for a rational r_min(n) whose numerator and denominator are 5-smooth, and disappear when the interval reaches [1, r_max(n)) where r_max(n) = (r_min(n+1))^2/r_min(n).
The initially known sequence terms (n <= 20) have the following property: a(n) = 2 * a(n-1) or, for some i < n, a(n) = 2 * a(n-1) - a(i); the former case occurring if r_min(n) equals a prime number, the latter occurring if and only if r_min(n) = r_max(i). Having identified underlying reasons that look like the basis of a proof, the author conjectures this is true for all n > 1.
Infinite equivalence classes occur if the interval is a scaled image of an interval containing (1, 27/5). Intervals that are scaled images of an interval containing [1, 50/9) or (1, 50/9] have no finite equivalence classes; for other intervals, the smallest equivalence class has cardinality 1, 2, 4 or 12 (cf. A215795).
The 5-smooth rationals in the critical semi-open interval [1, 50/9) are therefore partitioned into infinite equivalence classes. I believe these classes can be numbered, ordered by the power(s), k, of 5, for which 5^k*2^j is in the class for at least one integer j. 8/5, 1/1 and 5/1 are in class 0; 25/8 and 125/64 are in class 1; 32/25 and 128/125 are in class -1. If we close the critical interval, the inclusion of 50/9 unites classes 0 and 1 as defined above. If an interval, I, (or any scaled image of I) is a proper superset of the closed interval [1, 50/9], I am confident all 5-smooth rationals in I are in the same equivalence class. - Peter Munn, Jun 29 2022
From Peter Munn, Sep 16 2022: (Start)
The linked figure, "The disappearance of finite classes", shows various interval parts of the interval [1, 27/5] as rectangles. The arrowed lines between the rectangles represent the relation R. It can be seen as summarizing a graph between uncountably many points.
There are 12 rectangles around the circumference annotated as representing points that are on a finite path. These correspond to the points (equivalently real numbers) that are in equivalence classes of cardinality 12 -- the only finite equivalence classes when the interval is large enough to have infinite equivalence classes.
If we animate this figure to show the upper bound of the interval increasing from 27/5, the 12 outer rectangles would shrink. Consider, in particular, the one on the right marked with bounds 27/25 and 10/9. 10/9 derives from (((1 * 5) / 3) * 2) / 3, a progression that can followed around the inner edge of the circuit in the figure. 10/9 stays fixed. 27/25 grows as 1/5 of the upper bound, and the finite classes/paths disappear when the upper bound reaches 50/9 = 5 * 10/9.
(End)

Examples

			The real interval [1.0, 2.0] has equivalence classes including {1.5} and {1.0, 2.0}, so |{1.5}| = 1 and |{1.0, 2.0}| = 2 are in the sequence.
To establish whether 3 is in the sequence, let x be the smallest member (by absolute value) of an equivalence class E on an interval I. E has cardinality greater than 1 only if 2*x is in I. If 3*x is not in I, there are no other multiples or submultiples of {x, 2*x} in I, so E has cardinality less than 3. Otherwise 3*x is in I, 3/2*x is in I, so {x, 2*x, 3*x, 3/2*x} is a subset of E and E has cardinality greater than 3. So 3 is not in the sequence.
Equivalence classes of cardinality 7 take the form {x, 2*x, 3*x, 3/2*x, 4*x, 4/3*x, 8/3*x} for real x <> 0.
Equivalence classes of cardinality 12 take the form {x, 2*x, 3*x, 3/2*x, 4*x, 4/3*x, 8/3*x, 9/2*x, 9/4*x, 9/8*x, 27/8*x, 27/16*x} for real x <> 0.
In the table below the columns are as follows. |E|: cardinality of equivalence class; I_min: smallest interval with lower bound 1 that has such an equivalence class; I_max: largest interval with lower bound 1 that has such an equivalence class; u_max: upper bound of I_max in decimal notation.
  |E|    I_min            I_max                  u_max
    1    [1, 1]           (1, 4)                 4.0
    2    [1, 2]           (1, 9/2)               4.5
    4    [1, 3]           (1, 16/3)              5.3333333333...
    7    [1, 4]           (1, 81/16)             5.0625
   12    [1, 9/2]         (1, 50/9)              5.5555555555...
   24    [1, 5]           (1, 6561/1280)         5.12578125
   41    [1, 81/16]       (1, 531441/102400)     5.189853515625
   58    [1, 6561/1280]   (1, 3^16/(2^16*5^3))   5.2547266845...
  ...
  214    [1, 16/3]        (1, 2^52*5^8/3^43)     5.3592727015...
  ...
   oo    (1, 27/5)        not defined            not defined
		

Crossrefs

Programs

  • PARI
    { my (points = Set([1]), newmax);
      print (1," ",matsize(points)[2]);
      for (n=2,25,
        newmax = vecmin(setminus(setbinop((x,y)->x*y,[2,3,5],points),points));
        points = setunion(setbinop((x,y)->x/y,[newmax],points),points);
    print (n," ",matsize(points)[2]); )} \\ Peter Munn, Jun 29 2022

Formula

a(n) = |S_n| where S_n = S_(n-1) U {min( {2*x, 3*x, 5*x : x in S_(n-1)} \ S_(n-1)) / y : y in S_(n-1)} n >= 2, S_1 = {1.0}.

Extensions

a(21)-a(25) from Peter Munn, Jun 29 2022
Showing 1-5 of 5 results.