cp's OEIS Frontend

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.

A107171 Primes of the form 5x^2 + 13y^2.

This page as a plain text file.
%I A107171 #19 Feb 10 2017 09:48:58
%S A107171 5,13,97,137,193,197,457,593,613,617,733,877,1097,1237,1373,1553,1753,
%T A107171 1877,1913,1997,2273,2293,2333,2377,2477,2593,2917,2953,3037,3373,
%U A107171 3517,3593,3673,3677,3697,3733,3853,4217,4337,4457,4513,4517,4673
%N A107171 Primes of the form 5x^2 + 13y^2.
%C A107171 Discriminant = -260. See A107132 for more information.
%H A107171 Vincenzo Librandi and Ray Chandler, <a href="/A107171/b107171.txt">Table of n, a(n) for n = 1..10000</a> [First 1000 terms from Vincenzo Librandi]
%H A107171 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)
%t A107171 QuadPrimes2[5, 0, 13, 10000] (* see A106856 *)
%o A107171 (PARI) list(lim)=my(v=List(),w,t); for(x=0, sqrtint(lim\5), w=5*x^2; for(y=0, sqrtint((lim-w)\13), if(isprime(t=w+13*y^2), listput(v,t)))); Set(v) \\ _Charles R Greathouse IV_, Feb 10 2017
%K A107171 nonn,easy
%O A107171 1,1
%A A107171 _T. D. Noe_, May 13 2005