OFFSET
1,2
REFERENCES
Allan C. Wechsler, Posting to Math Fun Mailing List, July 16 2016.
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 6 because sigma(3)/3 = 4/3; prime power factors are {2^2, 3^-1}; 2^2 is larger in magnitude thus p = 2 so 3 * 2 = 6.
a(4) = 28 because sigma(4)/4= 7/4; prime power factors are {2^2, 7^1}; 7 is larger in magnitude thus p = 7 so 4 * 7 = 28. - Michael De Vlieger, Jul 20 2016
MATHEMATICA
Table[n #1^Sign[#2] & @@ Last@ FactorInteger@ First@ MaximalBy[#, Abs@ Log@ # &] &@ Map[#1^#2 & @@ # &, FactorInteger@ #] &[DivisorSigma[1, n]/n], {n, 74}] (* Michael De Vlieger, Jul 20 2016 *)
PROG
(PARI) cmpf(x, y) = f[x, 1]^abs(f[x, 2]) >= f[y, 1]^abs(f[y, 2]);
a(n) = if (n==1, 1, f = factor(sigma(n)/n); vf = vector(#f~, k, k); vsi = vecsort(vf, cmpf, 1); imax = vsi[#f~]; if (f[imax, 2] > 0, n*f[imax, 1], n/f[imax, 1])); \\ Michel Marcus, Jul 20 2016; corrected Jun 13 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Jul 20 2016
EXTENSIONS
More terms from Michel Marcus, Jul 20 2016
STATUS
approved