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.
2 = 1^2+1^4, 5 = 2^2+1^4, 17 = 4^2+1^4 = 1^2+2^4, …, 97 = 9^2+2^4 = 4^2+3^4, etc.
mx = 10^12; lst = {}; Do[a = i^2 + j^4; If[ PrimeQ[a], AppendTo[ lst, a]], {i, Sqrt[mx]}, {j, Sqrt[ Sqrt[mx - i^2]]}]; Table[ Length@ Select[lst, # < 10^n &], {n, 12}]
2 = 1^2+1^4, 5 = 2^2+1^4, 17 = 4^2+1^4 = 1^2+2^4, ..., 97 = 9^2+2^4 = 4^2+3^4, etc.
mx = 2^40; lst = {}; Do[a = i^2 + j^4; If[ PrimeQ[a], AppendTo[ lst, a]], {i, Sqrt[mx]}, {j, Sqrt[ Sqrt[mx - i^2]]}]; Table[ Length@ Select[lst, # <2^n &], {n, 40}]
mx = 2^40; lst = {}; Do[a = i^2 + j^4; If[ PrimeQ[a], AppendTo[lst, a]], {i, Sqrt[mx]}, {j, Sqrt[ Sqrt[mx - i^2]]}]; Table[ Length@ Select[ Union@ lst, # < 2^n &], {n, 40}]
a(n)=my(N=2^n,v=List(),t);for(a=1,sqrt(N),forstep(b=a%2+1, sqrtint(sqrtint(N-a^2)), 2, t=a^2+b^4;if(isprime(t),listput(v,t)))); 1+#vecsort(Vec(v),,8) \\ Charles R Greathouse IV, Jun 12 2013
Comments