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.

A331038 Residues of the Lucas-Lehmer primality test for M(127) = 2^127 - 1.

Original entry on oeis.org

3, 7, 47, 2207, 4870847, 23725150497407, 562882766124611619513723647, 9932388036497706472820043948129789713, 102423269049837077051675109560558766898, 7949236499829405891753012242872011683, 119093374737774941856311333667076322210
Offset: 0

Views

Author

Sergio Pimentel, Jan 08 2020

Keywords

Comments

Since a(125) = 0, 2^127 - 1 = 170141183460469231731687303715884105727 is prime. This calculation was carried out by hand by Edouard Lucas. It took him 19 years from 1857 to 1876. The method works with a(0) = 3 since M(127) == 3 (mod 4). It also works with a(0) = 4 or a(0) = 10.

Crossrefs

Programs

  • Mathematica
    NestList[Mod[#^2-2,2^127-1]&, 3,10] (* Stefano Spezia, Mar 28 2025 *)

Formula

a(n) = (a(n-1)^2 - 2) mod (2^127-1) with a(0) = 3; a(125) is the final term.