OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
For n=4, the divisors of 4 are [1, 2, 4] and summing them as increasing powers yields: 1^1+2^2+4^3 = 69.
For n=12, the divisors of 12 are [1, 2, 3, 4, 6, 12] and summing them as increasing powers yields: 1^1+2^2+3^3+4^4+6^5+12^6 = 2994048.
MAPLE
f:= proc(n) local D, k;
D:=sort(convert(numtheory:-divisors(n), list));
add(D[k]^k, k=1..nops(D))
end proc:
map(f, [$1..100]); # Robert Israel, Sep 11 2020
MATHEMATICA
Total[Divisors[#]^Range[DivisorSigma[0, #]]]&/@Range[40] (* Harvey P. Dale, Aug 16 2011 *)
PROG
(PARI) a(n) = my(d = divisors(n)); sum(k=1, #d, d[k]^k); \\ Michel Marcus, Jan 01 2016
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Jason Earls, Sep 21 2010
EXTENSIONS
Name made precise by Peter Munn, Sep 19 2024
STATUS
approved