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 A182434 #19 Feb 26 2015 17:25:43 %S A182434 0,0,1,1,1,1,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3, %T A182434 3,3,4,4,4,4,4,4,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,6,6, %U A182434 6,6,6,6,7,7,7,7,7,7,8,8,8,8,8,8,8,8,8 %N A182434 Number of primes p < n such that 4*p+1 is also prime. %H A182434 Enrique Pérez Herrero, <a href="/A182434/b182434.txt">Table of n, a(n) for n = 1..10000</a> %F A182434 a(10^n) = A182265(n). %F A182434 a(n) = sum(i=2..n, floor(phi(4*i^2+i)/(4*i^2-4*i))). - _Enrique Pérez Herrero_, May 02 2012. %t A182434 Accumulate[Table[Boole[PrimeQ[n]&&PrimeQ[4n+1]],{n,1,200}]] %t A182434 Accumulate[If[AllTrue[{#,4#+1},PrimeQ],1,0]&/@Range[90]] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Feb 26 2015 *) %Y A182434 Cf. A023212, A182265. %Y A182434 Cf. A005384, A156660, A156874, A092816. %K A182434 nonn %O A182434 1,7 %A A182434 _Enrique Pérez Herrero_, Apr 28 2012