OFFSET
1,1
COMMENTS
Since k^7 - 1 = (k-1)*(k^6 + k^5 + k^4 + k^3 + k^2 + k + 1) and k^7 + 1 = (k+1)*(k^6 - k^5 + k^4 - k^3 + k^2 - k + 1) (and since there is no term less than 3, so k-1 must have at least one prime factor), this sequence lists the numbers k such that k-1, k+1, k^6 + k^5 + k^4 + k^3 + k^2 + k + 1, and k^6 - k^5 + k^4 - k^3 + k^2 - k + 1 are all prime. - Jon E. Schoenfield, Dec 14 2016
MATHEMATICA
Select[Range[100000], PrimeOmega[#^7 - 1] == PrimeOmega[#^7 + 1]== 2 &]
PROG
(Magma) IsSemiprime:=func<n | &+[d[2]: d in Factorization(n)] eq 2>; [n: n in [4..10000] | IsSemiprime(n^7-1)and IsSemiprime(n^7+1)]
CROSSREFS
KEYWORD
nonn
AUTHOR
Vincenzo Librandi, Dec 09 2016
EXTENSIONS
More terms from Jon E. Schoenfield, Dec 14 2016
STATUS
approved