login
Primes of the form 3x^3+x+1.
0

%I #6 Mar 02 2023 15:17:36

%S 5,197,3011,5197,10141,46901,59077,81031,499181,555637,648061,902357,

%T 985597,1119817,1479197,1778197,2491847,2738117,2910997,3183727,

%U 4444747,4562741,4804957,5184121,6145277,7381261,8958097,9146021,11609837

%N Primes of the form 3x^3+x+1.

%H David Broadhurst, <a href="http://physics.open.ac.uk/~dbroadhu/chg_alg.pdf">An algebraic approach to primality proving</a>.

%t Select[Table[3x^3+x+1,{x,200}],PrimeQ] (* _Harvey P. Dale_, Mar 02 2023 *)

%o (PARI) g(n,k) = for(x=1,n,y=k*x^3+x+1;if(isprime(y),print1(y",")))

%K easy,nonn

%O 1,1

%A _Cino Hilliard_, Feb 09 2006