OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..250 (terms below 3*10^19)
Eric Weisstein's World of Mathematics, Twin Primes.
Wikipedia, Powerful number: Generalization.
EXAMPLE
139968 = 2^6 * 3^7 is a term since it is 5-full and 139967 and 139969 are twin primes.
MATHEMATICA
Select[6*Range[2*10^5], PrimeQ[# - 1] && PrimeQ[# + 1] && Min[FactorInteger[#][[;; , 2]]] > 4 &]
PROG
(PARI) is(n) = isprime(n-1) && isprime(n+1) && vecmin(factor(n)[, 2]) > 4;
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Feb 23 2023
STATUS
approved