OFFSET
1,1
COMMENTS
Subsequence of A006881.
LINKS
Peter J. C. Moses, Table of n, a(n) for n = 1..5000
Since 591 = 3*197 and numbers 591 + 3 - 1 = 593, 591 + 197 - 1 = 787 are both primes, then 591 is in the sequence.
FORMULA
A226770(a(n)-1) = 2.
MATHEMATICA
Select[Range[10000], (Last[#2]=={1, 1}&&And@@PrimeQ[#1+First[#2]-1]&)[#1, Transpose[FactorInteger[#1]]]&] (* Peter J. C. Moses, Jul 03 2013 *)
spQ[n_]:=Module[{fi=Transpose[FactorInteger[n]]}, fi[[2]]=={1, 1}&&AllTrue[ n-1+fi[[1]], PrimeQ]]; Select[Range[3000], spQ] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Sep 24 2015 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Shevelev, Jul 02 2013
STATUS
approved