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 A179549 #11 Dec 26 2018 13:04:01 %S A179549 13,53,127,233,257,293,457,521,541,547,587,857,983,1151,1153,1429, %T A179549 1481,1489,1511,1553,1559,1579,1621,1753,1861,2099,2129,2143,2239, %U A179549 2273,2341,2347,2383,2411,2417,2459,2549,2657,2677,2729,2741,2789,2837,2897,2927 %N A179549 Primes p such that p minus the sum of the square of its digits yields a prime. %C A179549 This sequence differs from A076163 which takes the absolute value of the difference. %H A179549 Harvey P. Dale, <a href="/A179549/b179549.txt">Table of n, a(n) for n = 1..1000</a> %e A179549 a(5)=257 since 257-(2^2+5^2+7^2)=179 is a prime. %t A179549 ssdQ[n_]:=Module[{c=n-Total[IntegerDigits[n]^2]},Positive[c]&&PrimeQ[c]]; Select[Prime[Range[500]],ssdQ] (* _Harvey P. Dale_, Dec 26 2018 *) %o A179549 (PARI) isok(n) = {if (isprime(n), digs = digits(n, 10); isprime(n - sum(i=1, #digs, digs[i]^2));, 0)} \\ _Michel Marcus_, Jul 18 2013 %Y A179549 Cf. A076162, A076163. %K A179549 nonn,base %O A179549 1,1 %A A179549 _Carmine Suriano_, Jul 19 2010 %E A179549 A-number typo and signs in the example corrected - _R. J. Mathar_, Oct 23 2010