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 A127871 #11 Sep 18 2019 21:04:05 %S A127871 2,16,28,34,82,118,266,296,328,436,778,1126,1238,1280,1486,1496,1612, %T A127871 1952,2102,2192,2312,2414,2578,2690,2770,2834,2872,3100,3410,3550, %U A127871 3620,3752,4012,4016,4240,4264,4450,4772,5084,5458,5732,5798,5864,6704,7208 %N A127871 Numbers n such that n^4+1 and n^4+3 are twin primes. %C A127871 Intersection of A000068 and A125259. %C A127871 Smallest k such that a(k+1) = a(k) + 2 is 364. - _Altug Alkan_, May 15 2018 %H A127871 Robert Israel, <a href="/A127871/b127871.txt">Table of n, a(n) for n = 1..10000</a> %p A127871 select(t -> isprime(t^4+1) and isprime(t^4+3), [seq(i,i=2..10^4,2)]); # _Robert Israel_, May 14 2018 %t A127871 Select[Range[0,200000,2],PrimeQ[ #^4+1]&&PrimeQ[ #^4+3]&] %t A127871 Select[Range[7300],AllTrue[#^4+{1,3},PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Sep 18 2019 *) %o A127871 (PARI) isok(n) = isprime(n^4+1) && isprime(n^4+3); \\ _Michel Marcus_, May 15 2018 %Y A127871 Cf. A000068, A125259. %K A127871 nonn %O A127871 1,1 %A A127871 _Zak Seidov_, Apr 05 2007