OFFSET
1,4
MAPLE
A166546 := proc(n) option remember; local a; if n = 1 then 1; else for a from procname(n-1)+1 do if isprime(1+numtheory[tau](a)) then return a; end if; end do ; end if; end proc: A166724 := proc(n) numtheory[pi](A166546(n)) ; end proc: A166724 := proc(n) option remember ; numtheory[pi](A166546(n)) ; end proc: A167530 := proc(nmax) ol := 0 ; ct := 1 ; for n from 2 to nmax do if A166724(n) <> ol then printf("%d, ", ct) ; ct := 1 ; else ct := ct+1 ; end if; ol := A166724(n) ; end do ; return ; end proc: A167530(400) ; # R. J. Mathar, Feb 01 2010
CROSSREFS
KEYWORD
nonn
AUTHOR
Giovanni Teofilatto, Nov 05 2009
EXTENSIONS
Definition corrected by Giovanni Teofilatto, Jan 30 2010
More terms from R. J. Mathar, Feb 01 2010
Name clarified by Jason Kimberley, Aug 09 2018
STATUS
approved