A383016 Numbers k such that A382883(k) = 1.
1, 4, 6, 8, 9, 10, 14, 15, 21, 22, 25, 26, 27, 32, 33, 34, 35, 38, 39, 46, 49, 51, 55, 57, 58, 62, 65, 69, 74, 77, 82, 85, 86, 87, 91, 93, 94, 95, 106, 111, 115, 118, 119, 121, 122, 123, 125, 128, 129, 133, 134, 141, 142, 143, 145, 146, 155, 158, 159, 161, 166
Offset: 1
Keywords
Programs
-
Mathematica
V[n_, e_] := If[e == 1, 1, IntegerExponent[n, e]]; f[n_] := f[n] = -DivisorSum[n, V[n, #] * f[#] &, # < n &]; f[1] = 1; Select[Range[200], f[#] == 1 &] (* Amiram Eldar, Apr 29 2025 *)
-
SageMath
def A383016List(upto): return [n for n in srange(1, upto) if A382883(n) == 1]
Comments