A236485 Primes p such that gpf(lpf(m)-1) = gpf(gpf(m)-1), where m = 2^p-1 is Mersenne composite.
11, 29, 53
Offset: 1
Links
- GIMPS, Exponent Status
Programs
-
PARI
isok(p) = isprime(p) && (mc = (2^p-1)) && (mcpf = (factor(2^p-1))[, 1]) && (length(mcpf) > 1) && (gpf = vecmax(mcpf)) && (lpf = vecmin(mcpf)) && (vecmax(factor(gpf-1)[, 1]) == vecmax(factor(lpf-1)[, 1])); \\ Michel Marcus, Jan 27 2014
Comments