login
Difference between largest and smallest positive exponent in prime factorization of n; a(1) = 0 by convention.
5

%I #26 Jan 05 2024 02:49:46

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

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

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

%N Difference between largest and smallest positive exponent in prime factorization of n; a(1) = 0 by convention.

%H Antti Karttunen, <a href="/A062977/b062977.txt">Table of n, a(n) for n = 1..100000</a> (first 4000 terms from Harry J. Smith)

%H <a href="/index/Eu#epf">Index entries for sequences computed from exponents in factorization of n</a>.

%F a(n) = A051903(n) - A051904(n).

%F a(A108951(n)) = A325226(n) = A001222(n) - A071178(n). - _Antti Karttunen_, Nov 17 2019

%F Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = A033150 - 1 = 0.705211... . - _Amiram Eldar_, Jan 05 2024

%e a(24) = 2 since 24 = 2^3*3^1 and max(3,1) - min(3,1) = 3 - 1 = 2;

%e a(25) = 0 since 25 = 5^2 and max(2) - min(2) = 2 - 2 = 0.

%t dlsp[n_]:=Module[{xp=FactorInteger[n][[All,2]]},Max[xp]-Min[xp]]; Join[ {0},Array[ dlsp,120]] (* _Harvey P. Dale_, Jan 28 2021 *)

%o (PARI) { for (n=1, 4000, if (n<2, M=m=0, f=factor(n)~; M=m=f[2, 1]; for (i=2, length(f), M=max(M, f[2, i]); m=min(m, f[2, i]))); write("b062977.txt", n, " ", M - m) ) } \\ _Harry J. Smith_, Aug 14 2009

%o (PARI) A062977(n) = if((1==n),0,n=(factor(n)[, 2]); vecmax(n)-vecmin(n)); \\ _Antti Karttunen_, Nov 17 2019

%Y Cf. A072774 (positions of zeros), A059404 (of nonzeros).

%Y Cf. A001222, A033150, A051903, A051904, A071178, A108951, A325226.

%K nonn,easy

%O 1,24

%A _Henry Bottomley_, Jul 24 2001