This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A185086 #37 Jun 20 2021 02:47:12 %S A185086 5,13,29,41,53,61,73,89,109,113,137,149,157,173,193,229,233,269,281, %T A185086 293,313,317,349,353,373,389,397,409,433,449,461,509,521,557,569,593, %U A185086 601,613,617,653,673,701,733,761,773,797,809,853,857,877,929,937,941,953 %N A185086 Fouvry-Iwaniec primes: Primes of the form k^2 + p^2 where p is a prime. %C A185086 Sequence is infinite, see Fouvry & Iwaniec. %C A185086 Its intersection with A028916 is A262340, by the uniqueness part of Fermat's two-squares theorem. - _Jonathan Sondow_, Oct 05 2015 %C A185086 Named after the French mathematician Étienne Fouvry (b. 1953) and the Polish-American mathematician Henryk Iwaniec (b. 1947). - _Amiram Eldar_, Jun 20 2021 %H A185086 Charles R Greathouse IV, <a href="/A185086/b185086.txt">Table of n, a(n) for n = 1..10000</a> %H A185086 Étienne Fouvry and Henryk Iwaniec, <a href="http://matwbn.icm.edu.pl/ksiazki/aa/aa79/aa7935.pdf">Gaussian primes</a>, Acta Arithmetica, Vol. 79, No. 3 (1997), pp. 249-287. %H A185086 Lasse Grimmelt, <a href="https://arxiv.org/abs/1809.10008">Vinogradov's Theorem with Fouvry-Iwaniec Primes</a>, arXiv:1809.10008 [math.NT], 2018. %H A185086 Art of Problem Solving, <a href="http://www.artofproblemsolving.com/wiki/index.php/Fermat's_Two_Squares_Theorem">Fermat's Two Squares Theorem</a>. %t A185086 nn = 1000; Union[Reap[Do[n = k^2 + p^2; If[n <= nn && PrimeQ[n], Sow[n]], {k, Sqrt[nn]}, {p, Prime[Range[PrimePi[Sqrt[nn]]]]}]][[2, 1]]] %o A185086 (PARI) is(n)=forprime(p=2,sqrtint(n),if(issquare(n-p^2),return(isprime(n))));0 %o A185086 (PARI) list(lim)=my(v=List(),N,t);forprime(p=2,sqrt(lim), N=p^2; for(n=1,sqrt(lim-N), if(ispseudoprime(t=N+n^2), listput(v,t)))); v=vecsort(Vec(v),,8); v %o A185086 (Haskell) %o A185086 a185086 n = a185086_list !! (n-1) %o A185086 a185086_list = filter (\p -> any ((== 1) . a010052) $ %o A185086 map (p -) $ takeWhile (<= p) a001248_list) a000040_list %o A185086 -- _Reinhard Zumkeller_, Mar 17 2013 %Y A185086 Subsequence of A002144 and hence of A002313. %Y A185086 The positive terms of A240130 form a subsequence. %Y A185086 Cf. A010052, A001248, A000040, A028916, A262340. %K A185086 nonn,nice %O A185086 1,1 %A A185086 _Charles R Greathouse IV_, Feb 18 2011