|
%I #17 Oct 05 2019 10:21:00
%S 1,88,125,5246
%N Numbers n such that Sum_{k=1..n} k*Mobius(k) = n.
%C No more terms < 10^7.
%t 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
%t smQ[x_]:=x==Total[Table[n*MoebiusMu[n],{n,x}]]; Select[Range[6000],smQ] (* _Harvey P. Dale_, Oct 05 2019 *)
%Y Cf. A008683, A002321, A055615, A068340.
%K nonn
%O 1,2
%A _Carl Najafi_, Dec 22 2012
|
