OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
s[n_] := s[n] = Module[{e = FactorInteger[n!][[;; , 2]]}, Sum[DigitCount[e[[k]], 2, 1], {k, 1, Length[e]}]]; Select[Range[2, 300], s[#] > s[# + 1] &]
PROG
(PARI) s(n) = {my(e = factor(n!)[, 2]); sum(k=1, #e, hammingweight(e[k])); }
lista(kmax) = {my(s1 = s(1), s2); for(k = 2, kmax, s2 = s(k); if(s1 > s2, print1(k-1, ", ")); s1 = s2); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jul 16 2024
STATUS
approved