OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[p_, e_] := (-1)^DigitCount[e, 2, 1]; imu[1] = 1; imu[n_] := Times @@ f @@@ FactorInteger[n]; Position[Accumulate[Array[imu, 1000]], 0] // Flatten
PROG
(PARI) imu(n) = vecprod(apply(x -> (-1)^hammingweight(x), factor(n)[, 2]));
lista(n) = {my(s = 0); for(k = 1, n, s+ = imu(k); if(s == 0, print1(k, ", "))); }
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Dec 23 2023
STATUS
approved