A291691 Primes p such that gpf(lpf(2^p - 1) - 1) = p.
2, 3, 5, 7, 11, 13, 23, 29, 37, 43, 47, 53, 73, 79, 83, 97, 113, 131, 151, 173, 179, 181, 191, 197, 211, 223, 233, 239, 251, 263, 277, 281, 283, 307, 317, 337, 353, 359, 367, 383, 397, 419, 431, 439, 443, 457, 461, 463, 467, 487, 491, 499
Offset: 1
Keywords
Examples
We have gpf(lpf(2^11 - 1) - 1) = gpf(23 - 1) = 11, so 11 is a term.
Links
- Charles R Greathouse IV, Table of n, a(n) for n = 1..119
Programs
-
Mathematica
lpf[n_] := FactorInteger[n][[1, 1]]; gpf[n_] := FactorInteger[n][[-1, 1]]; Select[ Prime@ Range@ 45, gpf[lpf[2^# - 1] - 1] == # &] (* Giovanni Resta, Aug 30 2017 *)
-
PARI
listp(nn) = forprime(p=2, nn, if (vecmax(factor(vecmin(factor(2^p-1)[,1])-1)[,1]) == p, print1(p, ", "));); \\ Michel Marcus, Aug 30 2017
Extensions
a(17)-a(26) from Michel Marcus, Aug 30 2017
a(27)-a(34) from Giovanni Resta, Aug 30 2017
a(35)-a(52) from Charles R Greathouse IV, Aug 30 2017
Comments