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 April 19 14:10 EDT 2024. Contains 371792 sequences. (Running on oeis4.)