|
%I #22 Dec 31 2023 06:23:41
%S 0,0,1,0,0,1,1,0,1,0,0,1,0,1,1,0,0,1,0,0,2,0,0,1,0,0,1,1,0,1,1,0,1,0,
%T 1,1,0,0,1,0,0,2,0,0,1,0,0,1,1,0,1,0,0,1,0,1,1,0,0,1,0,1,2,0,0,1,0,0,
%U 1,1,0,1,0,0,1,0,1,1,0,0,1,0,0,2,0,0,1,0,0,1,1,0,2,0,0,1,0,1,1,0
%N Number of distinct Mersenne primes dividing n.
%C a(n) = m first occurs at n = A098918(m). - _Robert Israel_, Feb 03 2020
%H Antti Karttunen, <a href="/A147645/b147645.txt">Table of n, a(n) for n = 1..131072</a> (terms 1..10000 from Robert Israel)
%F From _Antti Karttunen_, May 12 2022: (Start)
%F a(n) = A154402(n) - A353786(n)
%F a(n) = a(2*n) = a(A000265(n)).
%F a(n) <= A331410(n). (End)
%F Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = A173898 = 0.516454... . - _Amiram Eldar_, Dec 31 2023
%e a(21)=2 because 1, 3, 7 and 21 are divisors of 21. Then 21 has two divisors that are Mersenne primes (A000668): 3 and 7.
%p N:= 100: # for a(1)..a(N)
%p V:= Vector(N):
%p for i from 1 do
%p m:= numtheory:-mersenne([i]);
%p if m > N then break fi;
%p for j from m by m to N do
%p V[j]:= V[j]+1
%p od od:
%p convert(V,list); # _Robert Israel_, Feb 03 2020
%o (PARI) A147645(n) = { my(m=3,s=0); while(m<=n, s += (isprime(m)*!(n%m)); m += (m+1)); (s); }; \\ _Antti Karttunen_, May 12 2022
%Y Cf. A000265, A000668, A001221, A080225, A098918, A154402, A173898, A353786.
%Y Coincides with A331410 on A054784.
%K easy,nonn
%O 1,21
%A _Omar E. Pol_, Nov 09 2008
|
