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 A116886 #16 Jan 20 2015 02:21:44 %S A116886 3,17,103,137,277,313,677,743,1117,1627,2003,2143,3407,3677,4483,5087, %T A116886 5903,7177,7333,8087,8093,8147,8537,8573,9293,9473,10177,10477,11173, %U A116886 13807,14897,15107,16657,19753,21563,22307,24113,26113,26417,26633 %N A116886 Primes p that remain prime through at least 2 iterations of the function f(p) = p^2 + 4. %C A116886 Numbers p with the property that p, q = p^2 + 4, and r = q^2 + 4 are all prime. - _Zak Seidov_, Sep 08 2009 %C A116886 a(n) = sqrt(A165218(n) - 4). - _Zak Seidov_, Sep 08 2009 %H A116886 Zak Seidov, <a href="/A116886/b116886.txt">Table of n, a(n) for n = 1..1000</a> %e A116886 17 is prime, 17^2 + 4 = 293 is prime and 293^2 + 4 = 85853 is prime. %t A116886 Select[Prime[Range[2*7! ]],PrimeQ[ #^2+4]&&PrimeQ[(#^2+4)^2+4]&] (* _Vladimir Joseph Stephan Orlovsky_, Feb 26 2010 *) %t A116886 fQ[n_]:=AllTrue[Rest[NestList[#^2+4&,n,2]],PrimeQ]; Select[Prime[ Range[ 3000]],fQ] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Sep 21 2014 *) %o A116886 (PARI) is(n)=my(q);isprime(p) && isprime(q=p^2+4) && isprime(q^2+4) \\ _Charles R Greathouse IV_, Nov 06 2013 %Y A116886 Cf. A062324, A116887, A116888, A116889, A045637, A062324, A165218. %K A116886 nonn %O A116886 1,1 %A A116886 _Giovanni Resta_, Feb 27 2006 %E A116886 Edited by _N. J. A. Sloane_, Sep 18 2009 at the suggestion of _R. J. Mathar_