A359883 a(n) = Sum_{d|n} d^(n-d) * (n/d)^d.

1, 3, 4, 21, 6, 367, 8, 5129, 19693, 106411, 12, 9590989, 14, 105614223, 2439477016, 8590983185, 18, 1788338285659, 20, 44174380774421, 1483406745548512, 584318428289047, 24, 2300697575522919961, 298023223876953151, 2481152876039086107

Offset

1,2

Links

Table of n, a(n) for n=1..26.

Formula

a(n) = [x^n] Sum_{k>0} (n * x / k)^k / (1 - (k * x)^k).

If p is prime, a(p) = 1 + p.

Mathematica

a[n_] := DivisorSum[n, #^(n-#) * (n/#)^# &]; Array[a, 30] (* Amiram Eldar, Aug 09 2023 *)

Prog

(PARI) a(n) = sumdiv(n, d, d^(n-d)*(n/d)^d);

Crossrefs

Cf. A359882.

Sequence in context: A197410 A308689 A339482 * A032830 A225478 A254884

Adjacent sequences: A359880 A359881 A359882 * A359884 A359885 A359886

Keyword

nonn

Author

Seiichi Manyama, Jan 16 2023

Status

approved