A249800 - OEIS

%I #24 Feb 25 2023 15:18:13

%S 3,2,3,2,5,2,3,11,3,2,5,2,3,2,3,5,5,3,11,2,5,2,5,2,3,2,3,3,7,5,11,2,5,

%T 2,5,2,3,5,3,5,17,2,3,7,3,2,5,3,3,2,5,2,13,2,5,5,3,3,11,2,5,5,5,2,7,2,

%U 3,5,3,2,7,5,3,2,7,2,5,3,3,2,5,113,5,3,11

%N a(n) is the smallest prime q such that n(q+1)+1 is prime, that is, the smallest prime q such that n = (p-1)/(q+1) with p prime; or a(n) = -1 if no such q exists.

%C Variation on Schinzel's Hypothesis.

%H Paolo P. Lava, <a href="/A249800/b249800.txt">Table of n, a(n) for n = 1..1000</a>

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/SchinzelsHypothesis.html">Schinzel's Hypothesis</a>.

%e For n=1 the minimum primes p and q are 5 and 3: (p-1)/(q+1) = (5-1)/(3+1) = 4/4 = 1. Therefore a(1)=3.

%e For n=2 the minimum primes p and q are 7 and 2: (p-1)/(q+1) = (7-1)/(2+1) = 6/3 = 2. Therefore a(2)=2.

%p with(numtheory): P:=proc(q) local k,n;

%p for n from 1 to q do for k from 1 to q do

%p if isprime(n*(ithprime(k)+1)+1) then print(ithprime(k)); break; fi;

%p od; od; end: P(10^5);

%t a249800[n_Integer] := Module[{q}, q = 2; While[CompositeQ[n (q + 1) + 1], q = NextPrime[q]]; q]; a249800/@Range[120] (* _Michael De Vlieger_, Nov 19 2014 *)

%o (PARI) a(n) = my(q=2); while(! isprime(n*(q+1)+1), q = nextprime(q+1)); q; \\ _Michel Marcus_, Nov 07 2014

%Y Cf. A060324, A062251, A064632, A249801, A249802, A249803.

%K nonn,easy

%O 1,1

%A _Paolo P. Lava_, Nov 06 2014