A322485 The sum of the semi-unitary divisors of n.
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
Offset: 1
Examples
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.
References
- J. Chidambaraswamy, Sum functions of unitary and semi-unitary divisors, J. Indian Math. Soc., Vol. 31 (1967), pp. 117-126.
Links
- Amiram Eldar, Table of n, a(n) for n = 1..10000
- Pentti Haukkanen, Basic properties of the bi-unitary convolution and the semi-unitary convolution, Indian J. Math, Vol. 40 (1998), pp. 305-315.
- D. Suryanarayana and V. Siva Rama Prasad, Sum functions of k-ary and semi-k-ary divisors, Journal of the Australian Mathematical Society, Vol. 15, No. 2 (1973), pp. 148-162.
- Laszlo Tóth, Sum functions of certain generalized divisors, Ann. Univ. Sci. Budap. Rolando Eötvös, Sect. Math., Vol. 41 (1998), pp. 165-180.
Programs
-
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]
-
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
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.
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
Comments