%I M1430 N0564 #146 Feb 01 2024 05:04:22
%S 2,5,13,17,29,37,41,53,61,73,89,97,101,109,113,137,149,157,173,181,
%T 193,197,229,233,241,257,269,277,281,293,313,317,337,349,353,373,389,
%U 397,401,409,421,433,449,457,461,509,521,541,557,569,577,593,601,613,617
%N Primes congruent to 1 or 2 modulo 4; or, primes of form x^2 + y^2; or, -1 is a square mod p.
%C Or, primes p such that x^2 - p*y^2 represents -1.
%C Primes which are not Gaussian primes (meaning not congruent to 3 mod 4).
%C Every Fibonacci prime (with the exception of F(4) = 3) is in the sequence. If p = 2n+1 is the prime index of the Fibonacci prime, then F(2n+1) = F(n)^2 + F(n+1)^2 is the unique representation of the prime as sum of two squares. - _Sven Simon_, Nov 30 2003
%C Except for 2, primes of the form x^2 + 4y^2. See A140633. - _T. D. Noe_, May 19 2008
%C Primes p such that for all p > 2, p XOR 2 = p + 2. - _Brad Clardy_, Oct 25 2011
%C Greatest prime divisor of r^2 + 1 for some r. - _Michel Lagneau_, Sep 30 2012
%C Empirical result: a(n), as a set, compose the prime factors of the family of sequences produced by A005408(j)^2 + A005408(j+k)^2 = (2j+1)^2 + (2j+2k+1)^2, for j >= 0, and a given k >= 1 for each sequence, with the addition of the prime factors of k if not already in a(n). - _Richard R. Forberg_, Feb 09 2015
%C Primes such that when r is a primitive root then p-r is also a primitive root. - _Emmanuel Vantieghem_, Aug 13 2015
%C Primes of the form (x^2 + y^2)/2. Note that (x^2 + y^2)/2 = ((x+y)/2)^2 + ((x-y)/2)^2 = a^2 + b^2 with x = a + b and y = a - b. More generally, primes of the form (x^2 + y^2) / A001481(n) for every fixed n > 1. - _Thomas Ordowski_, Jul 03 2016
%C Numbers n such that ((n-2)!!)^2 == -1 (mod n). - _Thomas Ordowski_, Jul 25 2016
%C Primes p such that (p-1)!! == (p-2)!! (mod p). - _Thomas Ordowski_, Jul 28 2016
%C The product of 2 different terms (x^2 + y^2)(z^2 + v^2) = (xz + yv)^2 + (xv - yz)^2 is sum of 2 squares (A000404) because (xv - yz)^2 > 0. If x were equal to yz/v then (x^2 + y^2)/(z^2 + v^2) would be equal to ((yz/v)^2 + y^2)/(z^2 + v^2) = y^2/v^2 which is not possible because (x^2 + y^2) and (z^2 + v^2) are prime numbers. For example, (2^2 + 5^2)(4^2 + 9^2) = (2*4 + 5*9)^2 + (2*9 - 5*4)^2. - _Jerzy R Borysowicz_, Mar 21 2017
%D M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 872.
%D David A. Cox, "Primes of the Form x^2 + n y^2", Wiley, 1989.
%D G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford Univ. Press, 1979, p. 219, th. 251, 252.
%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 Zak Seidov, <a href="/A002313/b002313.txt">Table of n, a(n) for n = 1..10000</a> (first 1000 terms from T. D. Noe)
%H Dario Alpern, <a href="https://www.alpertron.com.ar/FSQUARES.HTM">Online program that calculates sum of squares representation</a>
%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 J. Todd, <a href="https://www.jstor.org/stable/2305526">A problem on arc tangent relations</a>, Amer. Math. Monthly, 56 (1949), 517-528.
%H Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/Fermats4nPlus1Theorem.html">Fermat's 4n+1 Theorem</a>
%H G. Xiao, <a href="https://wims.univ-cotedazur.fr/~wims/en_tool~number~twosquares.en.html">Two squares</a>
%H <a href="/index/Ga#gaussians">Index entries for Gaussian integers and primes</a>
%F a(n) ~ 2n log n. - _Charles R Greathouse IV_, Jul 04 2016
%F a(n) = A002331(n)^2 + A002330(n)^2. See crossrefs. - _Wolfdieter Lang_, Dec 11 2016
%e 13 is in the sequence since it is prime and 13 = 4*3 + 1. Also 13 = 2^2 + 3^2. And -1 is a square (mod 13): -1 + 2*13 = 25 = 5^2. Of course, only the first term is congruent to 2 (mod 4). - _Michael B. Porter_, Jul 04 2016
%p with(numtheory): for n from 1 to 300 do if ithprime(n) mod 4 = 1 or ithprime(n) mod 4 = 2 then printf(`%d,`,ithprime(n)) fi; od:
%p # alternative
%p A002313 := proc(n)
%p option remember ;
%p local a;
%p if n = 1 then
%p 2;
%p elif n = 2 then
%p 5;
%p else
%p for a from procname(n-1)+4 by 4 do
%p if isprime(a) then
%p return a ;
%p end if;
%p end do:
%p end if;
%p end proc:
%p seq(A002313(n),n=1..100) ; # _R. J. Mathar_, Feb 01 2024
%t Select[ Prime@ Range@ 115, Mod[#, 4] != 3 &] (* _Robert G. Wilson v_ *)
%t fQ[n_] := Solve[x^2 + 1 == n*y^2, {x, y}, Integers] == {}; Select[ Prime@ Range@ 115, fQ] (* _Robert G. Wilson v_, Dec 19 2013 *)
%o (PARI) select(p->p%4!=3, primes(1000)) \\ _Charles R Greathouse IV_, Feb 11 2011
%o (Haskell)
%o a002313 n = a002313_list !! (n-1)
%o a002313_list = filter ((`elem` [1,2]) . (`mod` 4)) a000040_list
%o -- _Reinhard Zumkeller_, Feb 04 2014
%o (Magma) [p: p in PrimesUpTo(700) | p mod 4 in {1,2}]; // _Vincenzo Librandi_, Feb 18 2015
%Y Apart from initial term, same as A002144. For values of x and y see A002330 and A002331.
%Y Cf. A033203, A038873, A038874, A045331, A008784, A057129, A084163, A084165, A137351.
%K nonn,easy,nice
%O 1,1
%A _N. J. A. Sloane_
%E More terms from _Henry Bottomley_, Aug 10 2000
%E More terms from _James A. Sellers_, Aug 22 2000
|