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.

A193423 Primes of the form x^2+y^3 with x, y and y^3 - x^2 >= 0.

This page as a plain text file.
%I A193423 #28 Jan 24 2023 10:07:50
%S A193423 2,31,43,73,89,113,241,337,347,359,379,443,487,521,593,599,733,829,
%T A193423 953,1009,1049,1129,1213,1289,1361,1367,1753,1777,1907,2017,2089,2213,
%U A193423 2297,2341,2393,2521,2689,2753,2953,2969,3221,3391,3571,3631,3797,3833,4051,4133,4159,4177,4217,4457,4721,4937,5237,5813
%N A193423 Primes of the form x^2+y^3 with x, y and y^3 - x^2 >= 0.
%C A193423 Subsequence of A066649.
%H A193423 Charles R Greathouse IV, <a href="/A193423/b193423.txt">Table of n, a(n) for n = 1..10000</a>
%e A193423 2 is in the sequence because 2 is prime, 2=1^2+1^3 and 1^3-1^2=0;
%e A193423 31 is in the sequence because 31 is prime, 31=2^2+3^3 and 3^3-2^2>0;
%e A193423 43 is in the sequence because 43 is prime, 43=4^2+3^3 and 3^3-4^2>0;
%o A193423 (PARI) list(lim)=my(v=List(),t,B); lim\=1; for(b=1,sqrtn(lim+.5,3),B=b^3; for(a=1,min(lim-B,sqrtint(B)),if(isprime(t=B+a^2),listput(v,t)))); vecsort(Vec(v),,8) \\ _Charles R Greathouse IV_, Aug 28 2011
%K A193423 nonn
%O A193423 1,1
%A A193423 _Juri-Stepan Gerasimov_, Aug 27 2011
%E A193423 Corrected by _D. S. McNeil_, Aug 27 2011