OFFSET
0,5
COMMENTS
a(p) = 1 for primes p since we have a(p) = (p')^(p') = 1^1 = 1.
MAPLE
a:= n-> (t-> t^t)(n*add(i[2]/i[1], i=ifactors(n)[2])):
seq(a(n), n=0..23); # Alois P. Heinz, Oct 20 2021
MATHEMATICA
Array[#^# &@ If[# < 2, 0, # Total[#2/#1 & @@@ FactorInteger[#]]] &, 19, 2] (* Michael De Vlieger, Oct 18 2021 *)
PROG
(PARI) ad(n) = vecsum([n/f[1]*f[2]|f<-factor(n+!n)~]); \\ A003415
a(n) = my(d=ad(n)); d^d; \\ Michel Marcus, Oct 19 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Oct 18 2021
STATUS
approved