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 A100570 #67 Dec 25 2020 11:02:08 %S A100570 1,2,3,12,17,28,32,72,108,117,297,657 %N A100570 Numbers that are not the sum of a square and a semiprime. %C A100570 No others up to 300000. Computed in collaboration with _Ray Chandler_. It appears that this sequence is finite, that is, that almost every positive integer is the sum of a semiprime and a square number. There are probably no further exceptions after a(12)=657. %C A100570 The statement about the finiteness of this sequence (namely, a(n)<=657) is much stronger than the Goldbach binary conjecture. Indeed, a much weaker conjecture, that this sequence contains no perfect squares >1, already implies the Goldbach conjecture. Cf. comment in A241922. - _Vladimir Shevelev_, May 01 2014 %C A100570 From _Daniel Mikhail_, Nov 23 2020: (Start) %C A100570 There are no new terms in this sequence between 658 and 2^28. %C A100570 Notably, A014090 (numbers that are not the sum of a square and one prime) is a known infinite sequence. (End) %H A100570 Daniel Mikhail, <a href="https://raw.githubusercontent.com/mikhaidn/SemiprimeCalculations/main/Summary%20of%202%5E28%20results">Brief glance of the search for integer solutions to a(n)-k^2 = semiprime</a> %F A100570 An integer is not an element for any integers i, j of the pairwise sum of {A001358(i)} and {A000290(j)}. %e A100570 From _Daniel Mikhail_, Nov 23 2020: (Start) %e A100570 An integer m is in this set if, for any primes, p and q, there does not exist a natural k, such that m-k^2 = p*q. %e A100570 Consider m=12 and all k such that k^2 < 12: k is either 0,1,4, or 9. %e A100570 12 - 0 = 12 = 2*2*2*3 => not semiprime; %e A100570 12 - 1 = 11 => not semiprime; %e A100570 12 - 4 = 8 = 2*2*2 => not semiprime; %e A100570 12 - 9 = 3 => not semiprime. %e A100570 Therefore, 12 is a term. (End) %t A100570 lim = 657; Complement[Range[lim],Select[Flatten[Outer[Plus,Select[Range[lim], PrimeOmega[#] == 2 &],Table[i^2, {i, 0, Sqrt[lim]}]]], # <= lim &]] (* _Robert Price_, Apr 10 2019 *) %Y A100570 Cf. A000290, A001358, A046903, A241922, A014090. %K A100570 nonn,hard,more %O A100570 1,2 %A A100570 _Jonathan Vos Post_, Nov 29 2004