A366246 The number of infinitary divisors of n that are "Fermi-Dirac primes" (A050376) and terms of A366242.
0, 1, 1, 0, 1, 2, 1, 1, 0, 2, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 2, 2, 1, 2, 0, 2, 1, 1, 1, 3, 1, 2, 2, 2, 2, 0, 1, 2, 2, 2, 1, 3, 1, 1, 1, 2, 1, 2, 0, 1, 2, 1, 1, 2, 2, 2, 2, 2, 1, 2, 1, 2, 1, 1, 2, 3, 1, 1, 2, 3, 1, 1, 1, 2, 1, 1, 2, 3, 1, 2, 1, 2, 1, 2, 2, 2, 2
Offset: 1
Links
- Amiram Eldar, Table of n, a(n) for n = 1..10000
Programs
-
Mathematica
s[0] = 0; s[n_] := s[n] = s[Floor[n/4]] + If[OddQ[Mod[n, 4]], 1, 0]; f[p_, e_] := s[e]; a[1] = 0; a[n_] := Plus @@ f @@@ FactorInteger[n]; Array[a, 100]
-
PARI
s(e) = if(e>3, s(e\4)) + e%2 \\ after Charles R Greathouse IV at A139351 a(n) = vecsum(apply(s, factor(n)[, 2]));
Comments