login
Number of primes with prime subscripts dividing n.
7

%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