login
A140452
2^(a(n))-1 contains an overpseudoprime divisor.
1
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,1
COMMENTS
If p is a prime then p is in the sequence iff 2^p-1 is a composite number.
LINKS
PROG
(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
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Shevelev, Jun 26 2008
EXTENSIONS
More terms from Michel Marcus, Dec 09 2018
STATUS
approved