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.

A096639 Smallest prime p == 5 mod 8 (A007521) and p > prime(n+2) such that p is a quadratic residue mod the first n odd primes 3, 5, 7, 11, ..., prime(n+1), and p is a quadratic non-residue mod prime(n+2).

This page as a plain text file.
%I A096639 #10 Mar 07 2013 17:12:42
%S A096639 5,13,61,109,421,1621,7309,8941,13381,82021,365509,300301,1336141,
%T A096639 644869,8658589,3462229,6810301,16145221,165163909,43030381,163384621,
%U A096639 249623581,2283397141,1272463669,2055693949
%N A096639 Smallest prime p == 5 mod 8 (A007521) and p > prime(n+2) such that p is a quadratic residue mod the first n odd primes 3, 5, 7, 11, ..., prime(n+1), and p is a quadratic non-residue mod prime(n+2).
%C A096639 Same as smallest prime p == 5 mod 8 with the property that the Legendre symbol (p|q) = 1 for the first n odd primes q = prime(k+1), k = 1, 2, ..., n, and (p|q) = -1 for q = prime(n+2).
%t A096639 f[n_] := Block[{k = 2}, While[ JacobiSymbol[n, Prime[k]] == 1, k++ ]; Prime[k]]; t = Table[0, {50}]; Do[p = Prime[n]; If[Mod[p, 8] == 5, a = f[p]; If[ t[[ PrimePi[a]]] == 0, t[[ PrimePi[a]]] = p; Print[ PrimePi[a], " = ", p]]], {n, 10^9}]; t
%Y A096639 Cf. A096634, A096636.
%K A096639 nonn
%O A096639 0,1
%A A096639 _Robert G. Wilson v_, Jun 24 2004
%E A096639 Better name from _Jonathan Sondow_, Mar 07 2013