OFFSET
1,1
COMMENTS
For p>3, p == 1 mod 6 and p(p+2) == 0 mod 3, hence, except for the first term, a(n) = subsequence of A014641 Odd octagonal numbers: (2n+1)(6n+1).
EXAMPLE
MATHEMATICA
# (# + 2)/3 & /@ Select[Prime@Range@140, PrimeQ[# + 4] &]
PROG
(PARI) {forprime (p=3, 10^3, isprime(p+4)&print1(p*(p+2)/3, ", "))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Jan 23 2011
EXTENSIONS
More terms from Michel Marcus, Mar 04 2014
STATUS
approved