This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A375849 #11 Aug 31 2024 18:00:58 %S A375849 1,1,3,3,1,1,7,7,1,1,5,5,11,11,15,15,3,3,1,9,19,19,3,3,23,23,25,25,7, %T A375849 7,31,31,15,15,17,17,35,35,9,9,39,39,41,41,21,21,3,3,47,47,49,49,3,13, %U A375849 53,53,27,27,9,9,57,57,63,63,31,31,31,15,67,67,11,11 %N A375849 The maximum odd exponent in the prime factorization of n!. %C A375849 The sequence of indices of record values, 2, 4, 8, 14, 16, 22, 26, 28, 32, 38, ..., are the odious numbers (A000069) multiplied by 2 (A128309). %H A375849 Amiram Eldar, <a href="/A375849/b375849.txt">Table of n, a(n) for n = 2..10000</a> %H A375849 <a href="/index/Fa#factorial">Index entries for sequences related to factorial numbers</a>. %F A375849 a(n) = A375032(n!). %F A375849 max(a(n), A375850(n)) = A011371(n). %t A375849 a[n_] := Max[Select[FactorInteger[n!][[;; , 2]], OddQ]]; Array[a, 100, 2] %o A375849 (PARI) a(n) = {my(e = select(x -> (x % 2), factor(n!)[, 2])); if(#e > 0, vecmax(e));} %o A375849 (Python) %o A375849 from collections import Counter %o A375849 from sympy import factorint %o A375849 def A375849(n): return max(filter(lambda x: x&1,sum((Counter(factorint(i)) for i in range(2,n+1)),start=Counter()).values())) # _Chai Wah Wu_, Aug 31 2024 %Y A375849 Cf. A000069, A000142, A011371, A128309, A375032, A375850. %K A375849 nonn,easy %O A375849 2,3 %A A375849 _Amiram Eldar_, Aug 31 2024