The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A007645 Generalized cuban primes: primes of the form x^2 + xy + y^2; or primes of the form x^2 + 3*y^2; or primes == 0 or 1 (mod 3).
(Formerly M2637)
86

%I M2637 #106 Aug 21 2021 13:13:49

%S 3,7,13,19,31,37,43,61,67,73,79,97,103,109,127,139,151,157,163,181,

%T 193,199,211,223,229,241,271,277,283,307,313,331,337,349,367,373,379,

%U 397,409,421,433,439,457,463,487,499,523,541,547,571,577,601,607,613

%N Generalized cuban primes: primes of the form x^2 + xy + y^2; or primes of the form x^2 + 3*y^2; or primes == 0 or 1 (mod 3).

%C Also, odd primes p such that -3 is a square mod p. - _N. J. A. Sloane_, Dec 25 2017

%C Equivalently, primes of the form p = (x^3 - y^3)/(x - y). If x=y+1 we get the cuban primes A002407, which is therefore a subsequence.

%C These are not to be confused with the Eisenstein primes, which are the primes in the ring of integers Z[w], where w = (-1+sqrt(-3))/2. The present sequence gives the rational primes which are also Eisenstein primes. - _N. J. A. Sloane_, Feb 06 2008

%C Also primes of the form x^2+3y^2 and, except for 3, x^2+xy+7y^2. See A140633. - _T. D. Noe_, May 19 2008

%C Conjecture: this sequence is Union(A002383,A162471). - _Daniel Tisdale_, Jul 04 2009

%C Primes p such that antiharmonic mean B(p) of the numbers k < p such that gcd(k, p) = 1 is not integer, where B(p) = A053818(p) / A023896(p) = A175505(p) / A175506(p) = (2p - 1) / 3. Primes p such that A175506(p) > 1. Subsequence of A179872. Union a(n) + A179891 = A179872. Example: a(6) = 37 because B(37) = A053818(37) / A023896(37) = A175505(37) / A175506(37) = 16206 / 666 = 73 / 3 (not integer). Cf. A179871, A179872, A179873, A179874, A179875, A179876, A179877, A179878, A179879, A179880, A179882, A179883, A179884, A179885, A179886, A179887, A179890, A179891, A003627, A034934. - _Jaroslav Krizek_, Aug 01 2010

%C Subsequence of Loeschian numbers, cf. A003136 and A024614; A088534(a(n)) > 0. - _Reinhard Zumkeller_, Oct 30 2011

%C Primes such that there exist a unique x, y, with 1 < x <= y < p, x + y == 1 (mod p) and x * y == 1 (mod p). - _Jon Perry_, Feb 02 2014

%C The prime factors of A002061. - _Richard R. Forberg_, Dec 10 2014

%C This sequence gives the primes p which solve s^2 == -3 (mod 4*p) (see Buell, Proposition 4.1., p. 50, for Delta = -3). p = 2 is not a solution. x^2 == -3 (mod 4) has solutions for all odd x. x^2 == -3 (mod p) has for odd primes p, not 3, the solutions of Legendre(-3|p) = +1 which are p == {1, 7} (mod 12). For p = 3 the representative solution is x = 0. Hence the solution of s^2 == -3 (mod 4*p) are the odd primes p = 3 and p == {1, 7} (mod 12) (or the primes p = 0, 1 (mod 3)). - _Wolfdieter Lang_, May 22 2021

%D D. A. Buell, Binary Quadratic Forms. Springer-Verlag, NY, 1989, p. 50.

%D Conway, J. H. and Guy, R. K. The Book of Numbers. New York: Springer-Verlag, pp. 220-223, 1996.

%D David A. Cox, "Primes of the Form x^2 + n y^2", Wiley, 1989, p. 7.

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

%D Wagon, S. "Eisenstein Primes." Section 9.8 in Mathematica in Action. New York: W. H. Freeman, pp. 319-323, 1991.

%H T. D. Noe, <a href="/A007645/b007645.txt">Table of n, a(n) for n = 1..1000</a>

%H U. P. Nair, <a href="http://arXiv.org/abs/math.NT/0408107">Elementary results on the binary quadratic form a^2+ab+b^2</a>, arXiv:math/0408107 [math.NT], 2004.

%H N. J. A. Sloane et al., <a href="https://oeis.org/wiki/Binary_Quadratic_Forms_and_OEIS">Binary Quadratic Forms and OEIS</a> (Index to related sequences, programs, references)

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

%F p == 0 or 1 (mod 3).

%F {3} UNION A002476. - _R. J. Mathar_, Oct 28 2008

%F A007645 UNION A003627 = A000040. - _Juri-Stepan Gerasimov_, Jan 28 2010

%p select(isprime,[3, seq(6*k+1, k=1..1000)]); # _Robert Israel_, Dec 12 2014

%t Join[{3},Select[Prime[Range[150]],Mod[#,3]==1&]] (* _Harvey P. Dale_, Aug 21 2021 *)

%o (PARI) forprime(p=2,1e3,if(p%3<2,print1(p", "))) \\ _Charles R Greathouse IV_, Jun 16 2011

%o (Haskell)

%o a007645 n = a007645_list !! (n-1)

%o a007645_list = filter ((== 1) . a010051) $ tail a003136_list

%o -- _Reinhard Zumkeller_, Jul 11 2013, Oct 30 2011

%Y Apart from initial term, same as A045331.

%Y Cf. A002407 (cuban primes, a subsequence). A002648 and A201477 are also subsequences.

%Y Cf. A001479, A001480 (x and y such that a(n) = x^2 + 3y^2)

%Y Cf. A000040, A003627. - _Juri-Stepan Gerasimov_, Jan 28 2010

%Y Primes in A003136 and A034017.

%K nonn,easy

%O 1,1

%A _N. J. A. Sloane_, _Mira Bernstein_ and _Robert G. Wilson v_

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

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 19 23:42 EDT 2024. Contains 372703 sequences. (Running on oeis4.)