A348978 Numerator of ratio A332993(n) / sigma(n).
1, 1, 1, 1, 1, 5, 1, 1, 1, 8, 1, 11, 1, 11, 7, 1, 1, 31, 1, 6, 29, 17, 1, 23, 1, 20, 1, 25, 1, 17, 1, 1, 15, 26, 43, 67, 1, 29, 53, 38, 1, 71, 1, 13, 11, 35, 1, 47, 1, 27, 23, 46, 1, 47, 67, 53, 77, 44, 1, 37, 1, 47, 23, 1, 79, 37, 1, 20, 31, 113, 1, 139, 1, 56, 53, 67, 89, 131, 1, 26, 1, 62, 1, 155, 103, 65, 39, 83
Offset: 1
Links
Crossrefs
Programs
-
Mathematica
f[n_] := n/FactorInteger[n][[1, 1]]; g[1] = 1; g[n_] := g[n] = n + g[f[n]]; a[n_] := Numerator[g[n]/DivisorSigma[1, n]]; Array[a, 100] (* Amiram Eldar, Nov 06 2021 *)
-
PARI
A332993(n) = if(1==n,n,n + A332993(n/vecmin(factor(n)[,1]))); A348978(n) = { my(u=A332993(n)); (u/gcd(sigma(n), u)); };
Comments