login
A variant of the cuban primes: primes p = (x^3 - y^3)/(x - y) where x = y + 2.
(Formerly M4910 N2105)
7

%I M4910 N2105 #60 Dec 22 2023 08:34:14

%S 13,109,193,433,769,1201,1453,2029,3469,3889,4801,10093,12289,13873,

%T 18253,20173,21169,22189,28813,37633,43201,47629,60493,63949,65713,

%U 69313,73009,76801,84673,106033,108301,112909,115249,129793,139969,142573,147853,169933

%N A variant of the cuban primes: primes p = (x^3 - y^3)/(x - y) where x = y + 2.

%C Primes p such that p = 1 + 3*m^2 for some integer m (A111051). - _Michael Somos_, Sep 15 2005

%D A. J. C. Cunningham, Binomial Factorisations, Vols. 1-9, Hodgson, London, 1923-1929; see Vol. 1, pp. 245-259.

%D N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).

%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

%H Vincenzo Librandi, <a href="/A002648/b002648.txt">Table of n, a(n) for n = 1..5000</a>

%H A. J. C. Cunningham, <a href="/A001912/a001912.pdf">Binomial Factorisations</a>, Vols. 1-9, Hodgson, London, 1923-1929. [Annotated scans of a few pages from Volumes 1 and 2]

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

%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Cuban_prime">Cuban prime</a>.

%F a(n) = 3*A111051(n)^2 + 1. - _Paul F. Marrero Romero_, Nov 03 2023

%e 193 is a term since 193 = (9^3 - 7^3)/(9 - 7) is a prime.

%t Select[Table[3n^2+1,{n,0,700}],PrimeQ] (* _Vincenzo Librandi_, Dec 02 2011 *)

%o (PARI) {a(n)= local(m, c); if(n<1, 0, c=0; m=1; while( c<n, m++; if( isprime(m)&issquare((m-1)/3), c++)); m)} /* _Michael Somos_, Sep 15 2005 */

%o (Magma) [a: n in [0..400] | IsPrime(a) where a is 3*n^2+1]; // _Vincenzo Librandi_, Dec 02 2011

%Y Cf. A002407, A111051 (values of m).

%Y A subsequence of A007645.

%K nonn

%O 1,1

%A _N. J. A. Sloane_

%E Entry revised by _N. J. A. Sloane_, Jan 29 2013