A174651 Numbers n such that Sum_{k=1..n} k*Mobius(k) = n.
1, 88, 125, 5246
Offset: 1
Keywords
Programs
-
Mathematica
muSums = Accumulate@Table[MoebiusMu[k] k, {k, 10^7}]; A174651 = {}; For[i = 1, i <= Length@muSums, i++, If[muSums[[i]] == i, AppendTo[A174651, i]]]; A174651 smQ[x_]:=x==Total[Table[n*MoebiusMu[n],{n,x}]]; Select[Range[6000],smQ] (* Harvey P. Dale, Oct 05 2019 *)
Comments