OFFSET
1,1
COMMENTS
Primes p such that p-2 divides 2^(p-1) - 1. The only member in A006512 is 5. - Robert Israel, Dec 03 2014
N=647089 is the smallest composite number such that (n-1)^2-1 divides 2^(n-1)-1. - M. F. Hasler, Jul 24 2015
LINKS
Robert G. Wilson v, Table of n, a(n) for n = 1..1076 (first 303 terms from Chai Wah Wu)
MAPLE
select(p -> isprime(p) and 2 &^ (p-1) - 1 mod (p-2) = 0, [seq(2*i+1, i=1..10^6)]); # Robert Israel, Dec 03 2014
MATHEMATICA
Select[Prime[Range[2, 81000]], PowerMod[2, #-1, #(#-2)]==1&] (* Harvey P. Dale, Sep 11 2023 *)
PROG
(PARI) lista(nn) = {forprime(p = 3, nn, if (! ((2^(p-1)-1) % (p*(p-2))), print1(p, ", ")); )} \\ Michel Marcus, Dec 02 2013
(Python)
from sympy import prime
from gmpy2 import powmod
A081762_list = [p for p in (prime(n) for n in range(2, 10**5)) if powmod(2, p-1, p*(p-2)) == 1] # Chai Wah Wu, Dec 03 2014
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Apr 09 2003
EXTENSIONS
More terms from Michel Marcus, Dec 02 2013
STATUS
approved