OFFSET
1,1
COMMENTS
a(16) <= 1846273228800. - David A. Corneth, Dec 11 2018
EXAMPLE
264 is in the sequence since its sum of semi-unitary divisors is susigma(264) = 528 = 2 * 264.
MATHEMATICA
f[p_, e_] := (p^Floor[(e+1)/2] - 1)/(p-1) + p^e; susigma[n_] := If[n==1, 1, Times @@ (f @@@ FactorInteger[n])]; aQ[n_] := susigma[n]==2n; Select[Range[10000], aQ]
PROG
(PARI) ssu(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); } \\ A322485
isok(n) = ssu(n) == 2*n; \\ Michel Marcus, Dec 14 2018
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Dec 11 2018
STATUS
approved