%I #19 Aug 31 2021 23:18:51
%S 2,4,10,36,56,94,126,224,260,270,300,350,686,716,780,1036,1070,1080,
%T 1156,1174,1210,1394,1416,1434,1440,1460,1524,1550,1576,1616,1654,
%U 1660,1700,1756,1860,1980,2054,2084,2096,2116,2224,2454,2600,2664,2770,2864
%N Positive integers n such that n^6 + 1 is semiprime.
%C n^6+1 can only be prime when n = 1, n^6+1 = 2. This is because the sum of cubes formula gives the polynomial factorization n^6+1 = (n^2+1) * (n^4 - n^2 + 1). Hence n^6+1 can only be semiprime when both (n^2+1) and (n^4 - n^2 + 1) are primes.
%H Robert Price, <a href="/A103854/b103854.txt">Table of n, a(n) for n = 1..1134</a>
%F a(n)^6 + 1 is semiprime. (a(n)^2+1) is prime and (a(n)^4 - a(n)^2 + 1) is prime.
%e n n^6+1 = (n^2+1) * (n^4 - n^2 + 1)
%e 2 65 = 5 * 13
%e 4 4097 = 17 * 241
%e 10 1000001 = 101 * 9901
%e 36 2176782337 = 1297 * 1678321
%e 56 30840979457 = 3137 * 9831361
%e 94 689869781057 = 8837 * 78066061
%e 126 4001504141377 = 15877 * 252031501
%e 224 126324651851777 = 50177 * 2517580801
%t semiprimeQ[n_] := Plus @@ Last /@ FactorInteger[n] == 2; Select[ 2Range@1526, semiprimeQ[ #^6 + 1] &] (* _Robert G. Wilson v_, May 26 2006 *)
%t Select[Range[200000], PrimeQ[#^2 + 1] && PrimeQ[(#^6 + 1)/(#^2 + 1)] &] (* _Robert Price_, Mar 11 2015 *)
%o (PARI) is(n)=my(s=n^2); isprime(s+1) && isprime(s^2-s+1) \\ _Charles R Greathouse IV_, Aug 31 2021
%Y Subsequence of A005574.
%Y Cf. A001358, A001538, A085722, A096173, A186669, A104238, A103854, A105041, A105066, A105078, A105122, A105142, A105237, A104335, A104479, A104494, A104657, A105282.
%K easy,nonn
%O 1,1
%A _Jonathan Vos Post_, Mar 31 2005
%E More terms from _Robert G. Wilson v_, May 26 2006
|