login
Numbers n such that n^16+1 and (n+2)^16+1 are both prime.
0

%I #8 Jun 24 2015 14:45:23

%S 74,156,474,476,686,688,774,776,830,832,982,2698,3560,3646,4016,5974,

%T 8542,8732,10730,10858,12164,12902,13022,13154,13554,13838,13840,

%U 20436,20454,20504,21912,24554,24770,25582,30144,31634,31882,32358,32360,32992,33758

%N Numbers n such that n^16+1 and (n+2)^16+1 are both prime.

%e 74 is in the sequence because 74^16+1 = 808551180810136214718004658177 and 76^16+1 = 1238846438084943599707227160577 are both prime.

%t lst={}; Do[p=n^16+1; q=(n+2)^16+1;If[PrimeQ[p] && PrimeQ[q], AppendTo[lst, n]], {n, 0, 34000}];lst

%t Select[Range[34000],AllTrue[{#^16+1,(#+2)^16+1},PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Jun 24 2015 *)

%Y Cf. A006313.

%K nonn,easy

%O 1,1

%A _Michel Lagneau_, Oct 17 2012