OFFSET
1,1
COMMENTS
Here we are considering powerful numbers defined in A001694.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..5000 from Donovan Johnson)
Eric Weisstein's World of Mathematics, Powerful Number.
EXAMPLE
72 = 2^3*3^2 is powerful and 71 and 73 are twin primes, so 72 is in the sequence.
MATHEMATICA
pwfQ[n_] := n == 1 || Min[Transpose[FactorInteger@n][[2]]] > 1; lst={}; Do[If[PrimeQ[n-1] && PrimeQ[n+1] && pwfQ[n], AppendTo[lst, n]], {n, 4, 2*10^6, 4}]; lst
PROG
(PARI) isA113839(n)=isprime(n-1)&&isprime(n+1)&&vecmax(factor(n)[, 2])>1 \\ Charles R Greathouse IV, Jun 29 2011
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Giovanni Resta, Jan 24 2006
STATUS
approved