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 A236168 #9 Nov 17 2014 14:41:02 %S A236168 2,3,11,23,29,31,41,59,71,113,151,163,191,239,241,269,359,431,433,499, %T A236168 503,521,541,563,661,683,701,751,773,829,883,983,1039,1259,1483,1499, %U A236168 1511,1549,1571,1609,1693,1721,1759,1913 %N A236168 Primes p such that p^3 - p - 1 is prime. %C A236168 Primes in A126421. %H A236168 Harvey P. Dale, <a href="/A236168/b236168.txt">Table of n, a(n) for n = 1..1000</a> %e A236168 269 is prime and 269^3 - 269 - 1 is also prime. So, 269 is a member of this sequence. %t A236168 Select[Prime[Range[300]],PrimeQ[#^3-#-1]&] (* _Harvey P. Dale_, Nov 17 2014 *) %o A236168 (Python) %o A236168 import sympy %o A236168 from sympy import isprime %o A236168 {print(p) for p in range(10**4) if isprime(p) and isprime(p**3-p-1)} %o A236168 (PARI) %o A236168 s=[]; forprime(p=2, 2000, if(isprime(p^3-p-1), s=concat(s, p))); s \\ _Colin Barker_, Jan 20 2014 %Y A236168 Cf. A126421. %K A236168 nonn %O A236168 1,1 %A A236168 _Derek Orr_, Jan 19 2014