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 A085989 #15 Feb 16 2025 08:32:50
%S A085989 8,13,18,20,25,29,32,34,40,41,45,50,52,53,58,61,65,68,72,73,74,80,85,
%T A085989 89,90,97,98,100,104,106,109,113,116,117,125,128,130,136,137,145,146,
%U A085989 148,149,153,157,160,162,164,169,170,173,178,180,181,185,193,194,200
%N A085989 Numbers that can be expressed as a sum of two squares, each >=2.
%H A085989 Robert Israel, <a href="/A085989/b085989.txt">Table of n, a(n) for n = 1..10000</a>
%H A085989 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/Landau-RamanujanConstant.html">Landau-Ramanujan Constant</a>
%e A085989 8 = 2^2 + 2^2, 13 = 2^2 + 3^2, ...
%p A085989 filter:= n -> ormap(t -> subs(t,x)>1 and subs(t,y) > 1, [isolve(x^2 + y^2 = n)]):
%p A085989 select(filter, [$0..200]); # _Robert Israel_, Mar 06 2017
%t A085989 r[n_] := Reduce[2 <= x <= y && n == x^2 + y^2, {x, y}, Integers]; Select[Range[200], r[#] =!= False &] (* _Jean-François Alcover_, Oct 29 2012 *)
%Y A085989 Cf. A001481.
%K A085989 nonn
%O A085989 1,1
%A A085989 _Eric W. Weisstein_, Jul 06 2003
%E A085989 Offset changed by _Robert Israel_, Mar 06 2017