a(n) = Sum_{d|n} d^rad(d).

Terms

a(0) =1a(1) =5a(2) =28a(3) =21a(4) =3126a(5) =46688a(6) =823544a(7) =85a(8) =757a(9) =10000003130a(10) =285311670612a(11) =3032688a(15) =341a(17) =34059641