OFFSET
1,3
COMMENTS
Not multiplicative: a(3)*a(5) != a(15), for example.
FORMULA
G.f.: Sum_{k>0} ((2*k-1) * x)^(2*k-1) / (1 - x^(2*k-1)).
a(2^n) = 1.
MATHEMATICA
a[n_] := DivisorSum[n, #^# &, OddQ[#] &]; Array[a, 20] (* Amiram Eldar, Jul 26 2023 *)
PROG
(PARI) a(n) = sumdiv(n, d, (d%2==1)*d^d);
(Python)
from sympy import divisors
def A363590(n): return sum(d**d for d in divisors(n>>(~n & n-1).bit_length(), generator=True)) # Chai Wah Wu, Jul 09 2023
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Seiichi Manyama, Jul 08 2023
STATUS
approved