|
|
A051377
|
|
a(1)=1; for n > 1, a(n) = sum of exponential divisors (or e-divisors) of n.
|
|
63
|
|
|
1, 2, 3, 6, 5, 6, 7, 10, 12, 10, 11, 18, 13, 14, 15, 22, 17, 24, 19, 30, 21, 22, 23, 30, 30, 26, 30, 42, 29, 30, 31, 34, 33, 34, 35, 72, 37, 38, 39, 50, 41, 42, 43, 66, 60, 46, 47, 66, 56, 60, 51, 78, 53, 60, 55, 70, 57, 58, 59, 90, 61, 62, 84, 78, 65, 66, 67, 102, 69, 70, 71
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
The e-divisors (or exponential divisors) of x=Product p(i)^r(i) are all numbers of the form Product p(i)^s(i) where s(i) divides r(i) for all i.
a(n) = n if and only if n is squarefree. - Jon Perry, Nov 13 2012
|
|
LINKS
|
Eric Weisstein's World of Mathematics, e-Divisor
|
|
FORMULA
|
The average order of a(n) is Dn + O(n^e) for any e > 0, due to Fabrykowski & Subbarao, where D is about 0.568. (D >= 0.5 since a(n) >= n.) - Charles R Greathouse IV, Sep 22 2023
|
|
EXAMPLE
|
a(8)=10 because 2 and 2^3 are e-divisors of 8 and 2+2^3=10.
|
|
MAPLE
|
local a, pe, p, e;
a := 1;
for pe in ifactors(n)[2] do
p := pe[1] ;
e := pe[2] ;
add(p^d, d=numtheory[divisors](e)) ;
a := a*% ;
end do:
a ;
end proc:
|
|
MATHEMATICA
|
a[n_] := Times @@ (Sum[ First[#]^d, {d, Divisors[Last[#]]}] & ) /@ FactorInteger[n]; Table[a[n], {n, 1, 71}] (* Jean-François Alcover, Apr 06 2012 *)
|
|
PROG
|
(PARI) a(n)=my(f=factor(n)); prod(i=1, #f[, 1], sumdiv(f[i, 2], d, f[i, 1]^d)) \\ Charles R Greathouse IV, Nov 22 2011
(PARI) ediv(n, f=factor(n))=my(v=List(), D=apply(divisors, f[, 2]~), t=#f~); forvec(u=vector(t, i, [1, #D[i]]), listput(v, prod(j=1, t, f[j, 1]^D[j][u[j]]))); Set(v)
(Haskell)
a051377 n = product $ zipWith sum_e (a027748_row n) (a124010_row n) where
sum_e p e = sum [p ^ d | d <- a027750_row e]
(GAP) A051377:=n->Product(List(Collected(Factors(n)), p -> Sum(DivisorsInt(p[2]), d->p[1]^d))); List([1..10^4], n -> A051377(n)); # Muniru A Asiru, Oct 29 2017
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy,nice,mult
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|