|
%I #26 Jan 16 2023 06:12:26
%S 1,2,3,5,6,8,9,11,13,14,20,21,26,29,30,33,35,36,39,41,43,49,54,55,62,
%T 68,69,75,76,79,81,89,90,105,110,113,117,119,126,134,141,146,154,162,
%U 174,176,178,179,186,191,195,207,209,215,216,222,225,227,230,231,234,237
%N Numbers k such that (k+1)*prime(k) + k is a prime.
%H Robert Israel, <a href="/A277124/b277124.txt">Table of n, a(n) for n = 1..10000</a>
%p p:= 1: count:= 0: Res:= NULL:
%p for n from 1 while count < 100 do
%p p:= nextprime(p);
%p if isprime(n*p+n+p) then
%p count:= count+1; Res:= Res, n
%p fi
%p od:
%p Res; # _Robert Israel_, Oct 25 2018
%t npQ[n_]:=Module[{p=Prime[n]},PrimeQ[n*p+n+p]]; Select[Range[250],npQ] (* _Harvey P. Dale_, Mar 27 2022 *)
%o (PARI) is(n) = isprime((n + 1) * (prime(n) + 1) - 1); \\ _Altug Alkan_, Oct 01 2016
%Y Cf. A000040, A086172.
%K nonn
%O 1,2
%A _Alex Ratushnyak_, Sep 30 2016
|
