login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Variation on A029834: a discrete version of the Mangoldt function. If n is prime then floor(log(prime(n))) else 0.
3

%I #16 Jun 01 2023 16:19:41

%S 0,1,1,0,2,0,2,0,0,0,3,0,3,0,0,0,4,0,4,0,0,0,4,0,0,0,0,0,4,0,4,0,0,0,

%T 0,0,5,0,0,0,5,0,5,0,0,0,5,0,0,0,0,0,5,0,0,0,0,0,5,0,5,0,0,0,0,0,5,0,

%U 0,0,5,0,5,0,0,0,0,0,5,0,0,0,6,0,0,0,0,0,6,0,0,0,0,0,0,0,6,0,0,0,6,0,6,0,0

%N Variation on A029834: a discrete version of the Mangoldt function. If n is prime then floor(log(prime(n))) else 0.

%H Antti Karttunen, <a href="/A062590/b062590.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = delta(tau(n), 2) * floor(log(prime(n))) = A010051(n) * A029835(n), where delta is the Kronecker delta function and tau is the number of divisors function. - _Alonso del Arte_, Sep 11 2013

%e a(5) = 2 because the fifth prime is 11, the logarithm of which is 2.397895...

%e a(6) = 0 because 6 is not prime.

%e a(7) = 2 because the seventh prime is 17, the logarithm of which is 2.833213344...

%t Table[Boole[PrimeQ[n]] Floor[Log[Prime[n]]], {n, 105}] (* _Alonso del Arte_, Sep 07 2013 *)

%o (PARI) v=[]; for(n=1,150,v=concat(v, if(isprime(n),floor(log(prime(n))),))); v

%Y Cf. A029834.

%K nonn

%O 1,5

%A _Jason Earls_, Jul 03 2001