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 A245097 #15 Jan 24 2022 04:34:30 %S A245097 2,67,22,4,27,8,5,0,5,7,1,13,2,0,62,9,7,2,6,0,8,11,2,0,0,0,9,5,1,5,4, %T A245097 0,12,0,4,13,2,0,12,12,14,33,7,0,7,18,3,0,3,0,3,8,2,0,2,0,2,16,6,51,6, %U A245097 0,50,0,11,6,4,0,4,20,7,3,8,0,3,0,32,47,2,0 %N A245097 Number of primes greater than n and not greater than the smallest pseudoprime to base n. %H A245097 Reinhard Zumkeller, <a href="/A245097/b245097.txt">Table of n, a(n) for n = 1..10000</a> %F A245097 a(n) = Sum {k=n+1..A007535(n)} A010051(k). %F A245097 a(n) = A000720(A007535(n)) - A000720(n). %o A245097 (Haskell) %o A245097 a245097 n = sum $ map a010051' [n + 1 .. a007535 n] %o A245097 (PARI) f(n) = forcomposite(m=n+1, , if(Mod(n, m)^(m-1)==1, return(m))); \\ A007535 %o A245097 a(n) = primepi(f(n)) - primepi(n); \\ _Michel Marcus_, Jan 24 2022 %Y A245097 Cf. A000720, A007535, A010051. %K A245097 nonn %O A245097 1,1 %A A245097 _Reinhard Zumkeller_, Nov 12 2014