|
|
A322485
|
|
The sum of the semi-unitary divisors of n.
|
|
5
|
|
|
1, 3, 4, 5, 6, 12, 8, 11, 10, 18, 12, 20, 14, 24, 24, 19, 18, 30, 20, 30, 32, 36, 24, 44, 26, 42, 31, 40, 30, 72, 32, 39, 48, 54, 48, 50, 38, 60, 56, 66, 42, 96, 44, 60, 60, 72, 48, 76, 50, 78, 72, 70, 54, 93, 72, 88, 80, 90, 60, 120, 62, 96, 80, 71, 84, 144
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
A semi-unitary divisor of n is defined as the largest divisor d of n such that the largest divisor of d that is a unitary divisor of n/d is 1 (see A322483).
|
|
REFERENCES
|
J. Chidambaraswamy, Sum functions of unitary and semi-unitary divisors, J. Indian Math. Soc., Vol. 31 (1967), pp. 117-126.
|
|
LINKS
|
|
|
FORMULA
|
Multiplicative with a(p^e) = sigma(p^floor((e-1)/2)) + p^e = (p^floor((e+1)/2) - 1)/(p-1) + p^e.
In particular a(p) = p + 1, a(p^2) = p^2 + 1, a(p^3) = p^3 + p + 1.
a(n) <= A000203(n) with equality if and only if n is squarefree (A005117).
Sum_{k=1..n} a(k) ~ c * n^2, where c = (zeta(2)*zeta(3)/2) * Product_{p prime} (1 - 2/p^3 + 1/p^5) = 0.7004703314... . - Amiram Eldar, Nov 24 2022
|
|
EXAMPLE
|
The semi-unitary divisors of 8 are 1, 2, 8 (4 is not semi-unitary divisor since the largest divisor of 4 that is a unitary divisor of 8/4 = 2 is 2 > 1), and their sum is 11, thus a(8) = 11.
|
|
MATHEMATICA
|
f[p_, e_] := (p^Floor[(e+1)/2] - 1)/(p-1) + p^e; susigma[n_] := If[n==1, 1, Times @@ (f @@@ FactorInteger[n])]; Array[susigma, 100]
|
|
PROG
|
(PARI) a(n) = {my(f = factor(n)); for (k=1, #f~, my(p=f[k, 1], e=f[k, 2]); f[k, 1] = (p^((e+1)\2) - 1)/(p-1) + p^e; f[k, 2] = 1; ); factorback(f); } \\ Michel Marcus, Dec 14 2018
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,mult
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|