OFFSET
1,1
COMMENTS
Subsequence of A243818: Primes p for which p^i - 4 is prime for i = 1, 3 and 5.
LINKS
Abhiram R Devesh, Table of n, a(n) for n = 1..100
EXAMPLE
Prime p=971 is in this sequence because p-4 = 967 (prime), p^3-4 = 915498607 (prime), p^5-4 = 863169625893847 (prime), and p^7-4 = 813831713247384370687 (prime).
PROG
(Python)
import sympy.ntheory as snt
n=2
while n>1:
....n1=n-4
....n2=((n**3)-4)
....n3=((n**5)-4)
....n4=((n**7)-4)
....##Check if n1 , n2, n3 and n4 are also primes.
....if snt.isprime(n1)== True and snt.isprime(n2)== True and snt.isprime(n3)== True and snt.isprime(n4)== True:
........print(n, n1, n2, n3, n4)
....n=snt.nextprime(n)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Abhiram R Devesh, Jun 12 2014
STATUS
approved