OFFSET
1,1
COMMENTS
Records are for n's such that A071904(n) = squares of primes.
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
MATHEMATICA
nn=501; With[{ci=Complement[Range[9, nn, 2], Prime[Range[PrimePi[nn]]]]}, FactorInteger[ #][[1, 1]]&/@ci] (* Harvey P. Dale, Nov 30 2012 *)
PROG
(Python)
from sympy import primepi, primefactors
def A162022(n):
if n == 1: return 3
m, k = n, primepi(n) + n + (n>>1)
while m != k:
m, k = k, primepi(k) + n + (k>>1)
return min(primefactors(m)) # Chai Wah Wu, Jul 31 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Jun 25 2009
EXTENSIONS
Corrected example a(4)=5 Francesco Antoni (francesco_antoni(AT)yahoo.com), Aug 04 2010
STATUS
approved