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 A071394 #31 Sep 01 2018 17:44:44 %S A071394 4,6,33,335,355,3073,8408,64690,481044,1304693,1304719,3524318, %T A071394 3524654,9559785,9559905,70115803,189963234,189963918,514278263, %U A071394 1394194660,3779856591,10246935974,75370122456,204475052725,204475053325,1505578023783,1505578024917 %N A071394 Numbers n divisible by pi(n) [A057809] with prime pi(n); i.e., largest prime factor of n equals pi(n). %H A071394 Giovanni Resta, <a href="/A071394/b071394.txt">Table of n, a(n) for n = 1..49</a> %F A071394 A000720(a(n)) = A006530(a(n)) = A256394(n). - _Jonathan Sondow_, Apr 15 2015 %e A071394 pi(8408) = 1051 and 8408 = 2*2*2*1051. %t A071394 c = 0; lpf[n_] := If[ PrimeQ[n], c++; n, Transpose[ FactorInteger[n]][[1, -1]]]; Do[ If[ lpf[n] == c, Print[n]], {n, 2, 10^7}] %t A071394 Select[Select[Range[2,10^6],IntegerQ[#/PrimePi[#]]&],PrimeQ[PrimePi[#]]&] (* _Ivan N. Ianakiev_, Apr 15 2015 *) %t A071394 Select[Range[10^6], FactorInteger[#][[-1, 1]] == PrimePi@ # &] (* _Michael De Vlieger_, Jul 30 2017 *) %o A071394 (PARI) isok(n) = isprime(p=primepi(n)) && !(n % p); \\ _Michel Marcus_, Jul 31 2017 %Y A071394 Cf. A000720, A006530, A256394. %K A071394 nonn %O A071394 1,1 %A A071394 _Jason Earls_, Jun 12 2002 %E A071394 Edited and extended by _Robert G. Wilson v_, Jun 13 2002 %E A071394 More terms from _Hans Havermann_, Jul 02 2002 %E A071394 a(26)-a(27) from _Giovanni Resta_, Mar 28 2017