login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Numbers n such that the set {2*n+p^2, p any prime} contains exactly one prime.
7

%I #8 May 21 2013 02:29:36

%S 1,4,7,10,16,19,22,25,31,37,40,46,49,52,61,64,70,79,82,85,91,94,109,

%T 112,115,121,124,127,130,136,142,151,154,169,172,175,187,190,196,205,

%U 211,217,220,226,229,235,241,247,250,256,274,277,280,289,292,295,304,316

%N Numbers n such that the set {2*n+p^2, p any prime} contains exactly one prime.

%C The sequence forms a subset of A016777, as explained below:

%C For each prime p<>3 we have p^2 =1 (mod 3), see A024700.

%C (i) For the k where 2*k=2 (mod 3), that is where k=1 (mod 3), this leads to 2*k+p^2=0 (mod 3), so the 2*k+p^2 are divisible by 3 (not prime) unless p=3.

%C The subcase where 2*k+3^2 is prime generates this sequence here; the subcase where it is not generates A138685.

%C (ii) For the k where 2*k=0 (mod 3), that is where k=0 (mod 3), one can select any p^2 =1 (mod 3)

%C to generate a prime 2*k+p^2 = 1 (mod 3), so these k generate many primes (of the form A002476).

%C (iii) For the k where 2*k=1 (mod 3), that is where k=2 (mod 3), one can select any p^2 =1 (mod 3)

%C to generate a prime 2*k+p^2 = 2 (mod 3), so these k generate many primes (of the form A003627).

%C The unique primes associated with each n are in A007528: n=1 associated with A007528(2)=11=2*1+3^2,

%C n=4 associated with A007528(3)=17=2*4+3^2 etc.

%H Vincenzo Librandi, <a href="/A138694/b138694.txt">Table of n, a(n) for n = 1..1000</a>

%F {This sequence here} Union {A138685} = {A016777}.

%e 3 is not in the sequence because {6+2^2, 6+3^3, 6+5^2, 6+7^2,..} = {10, 15, 31, 55,..,127,..,367,..}

%e contains the primes 31, 127, 367,..., generated with p=5,11,19...

%e 4 is in the sequence because {8+2^2, 8+3^3, 8+5^2, 8+7^2,..} = {12, 17, 33, 57,...} contains

%e only one prime (that is, 17), generated with p=3.

%t b = {}; Do[a = {}; Do[If[PrimeQ[2*k + Prime[n]^2], AppendTo[a, k]], {n, 1, 100}]; If[Length[a] < 2, AppendTo[b, a]], {k, 1, 500}]; Union[Flatten[b]]

%Y Cf. A138479, A138685, A138686, A138691, A138692, A138693, A138694, A007528, A138696, A138697, A138698, A138699, A138700.

%K nonn

%O 1,2

%A _Artur Jasinski_, Mar 27 2008

%E Edited by _R. J. Mathar_, May 15 2009