|
|
A327158
|
|
Unitary multiply-perfect numbers: n divides usigma(n), where usigma is the sum of unitary divisors of n (A034448).
|
|
9
|
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
10^13 < a(6) <= 146361946186458562560000. - Giovanni Resta, Aug 29 2019
|
|
LINKS
|
|
|
PROG
|
(PARI)
A034448(n) = { my(f=factorint(n)); prod(k=1, #f~, 1+(f[k, 1]^f[k, 2])); }; \\ After code in A034448
isA327158(n) = (gcd(n, A034448(n))==n);
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,more
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|