login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

Numbers k such that (k-1)*k^2 + 1 and k^2 + (k-1) are both prime.
3

%I #38 Sep 08 2022 08:46:07

%S 2,3,5,6,8,13,21,24,26,28,35,45,48,50,55,76,83,89,93,96,100,101,115,

%T 120,138,140,148,149,181,191,195,203,206,209,215,230,258,259,281,285,

%U 294,301,309,323,330,349,358,373,380,386,393,395,423,428,433,474,495

%N Numbers k such that (k-1)*k^2 + 1 and k^2 + (k-1) are both prime.

%C Numbers k such that (k^3 - k^2 + 1)*(k^2 + k - 1) is semiprime.

%C Intersection of A045546 and A111501.

%C Primes in this sequence: 2, 3, 5, 13, 83, 89, 101, 149, 181, 191, ...

%H Daniel Starodubtsev, <a href="/A239135/b239135.txt">Table of n, a(n) for n = 1..10000</a>

%e 2 is in this sequence because (2-1)*2^2+1=5 and 2^2+(2-1)=5 are both prime.

%t Select[Range[600],PrimeQ[#^2+#-1]&&PrimeQ[#^2(#-1)+1]&] (* _Farideh Firoozbakht_, Mar 17 2014 *)

%o (Magma) k := 1;

%o for n in [1..10000] do

%o if IsPrime(k*(n - 1)*n^2 + 1) and IsPrime(k*n^2 + n - 1) then

%o n;

%o end if;

%o end for;

%Y Cf. A239115.

%K nonn

%O 1,1

%A _Ilya Lopatin_ following a suggestion from _Juri-Stepan Gerasimov_, Mar 15 2014