A140452 2^(a(n))-1 contains an overpseudoprime divisor.
11, 22, 23, 25, 28, 29, 33, 35, 36, 37, 39, 41, 43, 44, 45, 46, 47, 48, 50, 51, 52, 53, 55, 56, 57, 58, 59, 60, 63, 64, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 81, 82, 83, 84, 86, 87, 88, 90, 91, 92, 94, 95, 96, 97, 99, 100, 101, 102, 103, 104, 105, 106, 108, 109
Offset: 1
Keywords
Links
- V. Shevelev, Overpseudoprimes, Mersenne Numbers and Wieferich Primes, arXiv:0806.3412 [math.NT], 2008-2012.
Programs
-
PARI
f(n) = my(t); sumdiv(2*n+1, d, eulerphi(d)/(t=znorder(Mod(2, d))))*t-t+1; \\ A137576 isopp(n) = (n>1) && !isprime(n) && (n == f((n-1)/2)); \\ A141232 isok(n) = {fordiv(2^n-1, d, if (isopp(d), return (1));); return (0);} \\ Michel Marcus, Dec 09 2018
Extensions
More terms from Michel Marcus, Dec 09 2018
Comments