OFFSET
1,1
COMMENTS
There are only 52189 exponential abundant numbers less than 50 million, which suggests that these account for approximately 0.1% of all integers.
Includes 36*m for all m coprime to 6 that are not squarefree. - Robert Israel, Feb 19 2019
The asymptotic density of this sequence is Sum_{n>=1} f(A328136(n)) = 0.001043673..., where f(n) = (6/(Pi^2*n))*Product_{prime p|n}(p/(p+1)). - Amiram Eldar, Sep 02 2022
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
Peter Hagis, Jr., Some Results Concerning Exponential Divisors, Internat. J. Math. & Math. Sci., Vol. 11, No. 2, (1988), pp. 343-350.
Eric Weisstein's World of Mathematics, e-Divisor.
EXAMPLE
The third integer that is exceeded by its proper exponential divisor sum is 3600. Hence a(3) = 3600.
MAPLE
filter:= proc(n) local L, m, i, j;
L:= ifactors(n)[2];
m:= nops(L);
mul(add(L[i][1]^j, j=numtheory:-divisors(L[i][2])), i=1..m)>2*n
end proc:
select(filter, [$1..10^5]); # Robert Israel, Feb 19 2019
MATHEMATICA
ExponentialDivisors[1]={1}; ExponentialDivisors[n_]:=Module[{}, {pr, pows}=Transpose@FactorInteger[n]; divpowers=Distribute[Divisors[pows], List]; Sort[Times@@(pr^Transpose[divpowers])]]; properexponentialdivisorsum[k_]:=Plus@@ExponentialDivisors[k]-k; Select[Range[5 10^4], properexponentialdivisorsum[ # ]># &]
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Ant King, Apr 28 2007
STATUS
approved