OFFSET
1,1
COMMENTS
All the terms in this sequence are 1 mod 9.
LINKS
K. D. Bajpai, Table of n, a(n) for n = 1..733
EXAMPLE
a(1) = 524521;
(524521^3 + 2)/3 = 48102471044890921;
(524521^5 + 2)/3 = 13234061480615091039311002201;
(524521^7 + 2)/3 = 3640985160809159281478976663465873196681;
all four are prime.
MATHEMATICA
Select[Prime[Range[10^7]], PrimeQ[(#^3 + 2)/3] && PrimeQ[(#^5 + 2)/3] && PrimeQ[(#^7 + 2)/3] &]
PROG
(PARI) is(n)=n%9==1 && isprime(n) && isprime((n^3+2)/3) && isprime((n^5+2)/3) && isprime((n^7+2)/3) \\ Charles R Greathouse IV, Jan 23 2015
(Magma) [p: p in PrimesInInterval(3, 10000000) | IsPrime((p^3 + 2) div 3) and IsPrime((p^5 + 2) div 3) and IsPrime((p^7 + 2) div 3)]; // Vincenzo Librandi, Mar 27 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
K. D. Bajpai, Jan 23 2015
STATUS
approved