cp's OEIS Frontend

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.

A063999 Numbers k such that the number of primes <= k is phi(phi(k)).

This page as a plain text file.
%I A063999 #19 Aug 19 2019 04:01:57
%S A063999 2,37,55,91,95,136,178,226,507,723,1166,1168,1284,1323,2658,3082,7564,
%T A063999 8166,11578,12014,12018,12026,15340,40220,50544,55490,56810,69732,
%U A063999 89210,100114,100126,111306,131532,226716,226722,232938,519312,683240,849534
%N A063999 Numbers k such that the number of primes <= k is phi(phi(k)).
%H A063999 Harry J. Smith, <a href="/A063999/b063999.txt">Table of n, a(n) for n = 1..44</a>
%e A063999 37 belongs to the sequence because number of primes <= 37 is equal to 12 (counting 2,3,5,7,11,13,17,19,23,29,31,37), while phi(37)=36 (counting 1,2,...,36) and phi(36)=12 (counting 1,5,7,11,13,17,19,23,25,29,31,35).
%p A063999 with(numtheory): p:=proc(n) if pi(n)=phi(phi(n)) then n else fi end: seq(p(n), n=1..900000); #  _Emeric Deutsch_, Feb 24 2005
%o A063999 (PARI) pi(n) = s=0; for(x=1,n, if(isprime(x),s++)); s; for(n=1,10^6, if(pi(n)==eulerphi(eulerphi(n)),print(n)))
%o A063999 (PARI) { default(primelimit, 2500000); n=0; for (m=1, 10^9, if (primepi(m)==eulerphi(eulerphi(m)), write("b063999.txt", n++, " ", m); if (n==44, break)) ) } \\ _Harry J. Smith_, Sep 05 2009
%Y A063999 Cf. A000720, A010554.
%K A063999 nonn
%O A063999 1,1
%A A063999 _Jason Earls_, Sep 06 2001
%E A063999 More terms from _Emeric Deutsch_, Feb 24 2005