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 A374328 #7 Jul 06 2024 05:29:52 %S A374328 0,0,1,0,0,0,1,0,0,1,0,0,0,2,0,1,0,1,0,0,0,1,0,1,0,0,0,0,0,0,1,0,0,0, %T A374328 0,0,0,1,1,0,0,2,1,1,0,1,0,0,0,0,0,1,0,0,1,0,0,0,1,0,0,0,0,0,1,1,0,0, %U A374328 0,2,2,0,0,1,0,0,0,0,1,0,1,0,0,0,0,1,1 %N A374328 The base-2 logarithm of the maximal exponent in the prime factorization of the numbers whose maximal exponent in their prime factorization is a power of 2. %C A374328 First differs from {A369934(n+1), n>=1} at n = 378. %C A374328 The first occurrence of k = 0, 1, ... is at 1, 3, 14, 224, 57307, ..., which is the position of 2^(2^k) at A369938. %H A374328 Amiram Eldar, <a href="/A374328/b374328.txt">Table of n, a(n) for n = 1..10000</a> %F A374328 a(n) = log_2(A374327(n)). %F A374328 a(n) = log_2(A051903(A369938(n))). %F A374328 Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{k>=0} k * d(k) / Sum_{k>=0} d(k) = 0.35575339438419889972..., where d(k) = 1/zeta(2^k+1) - 1/zeta(2^k) for k>=1, and d(0) = 1/zeta(2). %t A374328 f[n_] := Module[{e = If[n == 1, 0, Max[FactorInteger[n][[;; , 2]]]]}, If[e == 2^IntegerExponent[e, 2], IntegerExponent[e, 2], Nothing]]; Array[f, 150] %o A374328 (PARI) lista(kmax) = {my(e); for(k = 2, kmax, e = vecmax(factor(k)[, 2]); if(e >> valuation(e, 2) == 1, print1(valuation(e, 2), ", ")));} %Y A374328 Cf. A051903, A369934, A369938. %Y A374328 Similar sequences: A374324, A374325, A374326, A374327. %K A374328 nonn,easy %O A374328 1,14 %A A374328 _Amiram Eldar_, Jul 04 2024