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 A162873 #28 Apr 07 2021 04:47:30
%S A162873 17,19,149,151,197,199,269,271,521,523,809,811,881,883,1049,1051,1061,
%T A162873 1063,1277,1279,1949,1951,1997,1999,2141,2143,2549,2551,2789,2791,
%U A162873 2969,2971,3257,3259,3329,3331,3581,3583,3821,3823,4049,4051,4157,4159,4229,4231
%N A162873 List of pairs (p,r) of twin primes p and r=p+2 such that p-1, p+1 and r+1 are all divisible by squares > 1.
%H A162873 Amiram Eldar, <a href="/A162873/b162873.txt">Table of n, a(n) for n = 1..10000</a>
%e A162873 (197,199) is a pair because all of 196, 198, and 200 are divisible by squares.
%t A162873 dsQ[n_] := Length[Select[FactorInteger[n][[All, 2]], # > 1 &]] > 0; Select[ Partition[ Prime[ Range[ 250]],2,1],#[[2]]-#[[1]]==2&&AllTrue[ #[[1]]+ {-1,1,3},dsQ]&]//Flatten (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Nov 28 2018 *)
%Y A162873 Cf. A089189, A089194, A162870, A162872
%K A162873 nonn,tabf
%O A162873 1,1
%A A162873 _Vladimir Joseph Stephan Orlovsky_, Jul 15 2009
%E A162873 Corrected and previous (incorrect) Mathematica program replaced by _Harvey P. Dale_, Nov 28 2018