A383103 - OEIS

%I #9 Apr 29 2025 16:52:25

%S 4,8,9,25,27,32,49,121,125,128,169,243,289,343,361,529,841,900,961,

%T 1331,1369,1681,1764,1849,2048,2187,2197,2209,2809,3125,3481,3721,

%U 4356,4489,4900,4913,5041,5329,6084,6241,6859,6889,7921,8192,9409,10201,10404,10609

%N Numbers k such that A008683(k) = 0 and A382883(k) = 1.

%C See the comments in A382883.

%t 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[11000], !SquareFreeQ[#] && f[#] == 1 &] (* _Amiram Eldar_, Apr 29 2025 *)

%o (SageMath)

%o print([n for n in range(1, 10610) if moebius(n) == 0 and A382883(n) == 1])

%Y Cf. A008683, A382883, A053810.

%K nonn

%O 1,1

%A _Peter Luschny_, Apr 17 2025