%I #10 Jul 31 2022 07:46:25

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

%T 1,2,0,0,1,1,0,1,0,0,3,0,0,1,0,2,2,0,0,3,1,0,1,0,0,2,0,0,2,0,1,1,0,1,

%U 1,1,0,2,0,0,3,0,0,1,0,1,4,0,0,1,2,0,1

%N a(n) is the number of Fermat numbers dividing n, counted with multiplicity.

%C The multiplicity of a divisor d (not necessarily a prime) of n is defined in A169594 (see also the first formula).

%C A000244(n) is the least number k such that a(k) = n.

%C The asymptotic density of occurrences of 0 is 1/2.

%C The asymptotic density of occurrences of 1 is (1/2) * Sum_{k>=0} 1/(2^(2^k)+1) = (1/2) * A051158 = 0.2980315860... .

%H Amiram Eldar, <a href="/A356242/b356242.txt">Table of n, a(n) for n = 1..10000</a>

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/FermatNumber.html">Fermat Number</a>.

%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Fermat_number">Fermat number</a>.

%F a(n) = Sum_{k>=1} v(A000215(k), n), where v(m, n) is the exponent of the largest power of m that divides n.

%F a(A000215(n)) = 1.

%F a(A000244(n)) = a(A000351(n)) = a(A001026(n)) = n.

%F a(A003593(n)) = A112754(n).

%F a(n) >= A356241(n).

%F a(A051179(n)) = n.

%F a(A080307(n)) > 0 and a(A080308(n)) = 0.

%F Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{k>=0} 1/(2^(2^k)) = 0.8164215090... (A007404).

%t f = Table[(2^(2^n) + 1), {n, 0, 5}]; a[n_] := Total[IntegerExponent[n, f]]; Array[a, 100]

%Y Cf. A000244, A000215, A007404, A051158, A051179, A169594, A356241.

%Y Cf. A080307 (positions of nonzeros), A080308 (positions of 0's).

%K nonn

%O 1,9

%A _Amiram Eldar_, Jul 30 2022