%I #26 Jan 01 2024 09:14:58
%S 0,0,1,0,1,1,0,0,1,1,1,1,0,0,2,0,1,1,0,1,1,1,0,1,1,0,1,0,0,2,1,0,2,1,
%T 1,1,0,0,1,1,1,1,0,1,2,0,0,1,0,1,2,0,0,1,2,0,1,0,1,2,0,1,1,0,1,2,1,1,
%U 1,1,0,1,0,0,2,0,1,1,0,1,1,1,1,1,2,0,1,1,0,2,0,0,2,0,1,1,0,0,2,1,0,2,0,0,2,0,0,1,1,2,1,0,0,1,1,0,1,1,1,2
%N Number of primes with prime subscripts dividing n.
%H G. C. Greubel, <a href="/A279952/b279952.txt">Table of n, a(n) for n = 1..5000</a>
%F G.f.: Sum_{k>=1} x^prime(prime(k))/(1 - x^prime(prime(k))).
%F a(n) = Sum_{d|n} A111406(d-1). - _Ridouane Oudra_, Sep 12 2023
%F Additive with a(p^e) = 1 if primepi(p) is prime, and 0 otherwise. - _Amiram Eldar_, Nov 03 2023
%F Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{n>=1} 1/A006450(n) = 1.04... (see A006450 for a better estimate of this constant). - _Amiram Eldar_, Jan 01 2024
%e a(15) = 2 because 15 has 4 divisors {1,3,5,15} among which 2 divisors {3,5} are primes with prime subscripts.
%p with(numtheory): seq(add(pi(pi(d))-pi(pi(d-1)), d in divisors(n)), n=1..80); # _Ridouane Oudra_, Sep 12 2023
%t Rest[nmax = 120; CoefficientList[Series[Sum[x^Prime[Prime[k]]/(1 - x^Prime[Prime[k]]), {k, 1, nmax}], {x, 0, nmax}], x]]
%t f[p_, e_] := If[PrimeQ[PrimePi[p]], 1, 0]; a[1] = 0; a[n_] := Plus @@ f @@@ FactorInteger[n]; Array[a, 100] (* _Amiram Eldar_, Nov 03 2023 *)
%o (PARI) my(x='x+O('x^120)); concat([0, 0], Vec(sum(k=1, 120, x^prime(prime(k))/(1 - x^prime(prime(k)))))) \\ _Indranil Ghosh_, May 23 2017
%Y Cf. A001221, A006450, A257994.
%Y Cf. A000720, A111406.
%K nonn,easy
%O 1,15
%A _Ilya Gutkovskiy_, Dec 23 2016