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 A319050 #21 Apr 05 2019 17:27:06 %S A319050 7,23,43,47,79,97,151,167,223,241,331,349,359,367,439,523,547,619,691, %T A319050 727,773,823,839,907,1051,1087,1123,1223,1231,1249,1303,1367,1423, %U A319050 1447,1483,1523,1571,1627,1663,1699,1723,1811,1823,1847,1861,1879,1951,1987,2131,2203,2207 %N A319050 Primes p such that neither p + 1 nor p + 2 is squarefree. %H A319050 Seiichi Manyama, <a href="/A319050/b319050.txt">Table of n, a(n) for n = 1..10000</a> %e A319050 8 = 2^3 and 9 = 3^2. So 7 is a term. %e A319050 24 = 2^3*3 and 25 = 5^2. So 23 is a term. %t A319050 Select[Prime[Range[400]],NoneTrue[#+{1,2},SquareFreeQ]&] (* Requires Mathematica version 10 or later *) (* _Harvey P. Dale_, Apr 05 2019 *) %o A319050 (PARI) isok(p) = isprime(p) && !issquarefree(p+1) && !issquarefree(p+2); \\ _Michel Marcus_, Sep 09 2018 %Y A319050 Cf. A000040, A049098, A257108, A318959, A319051. %K A319050 nonn %O A319050 1,1 %A A319050 _Seiichi Manyama_, Sep 08 2018