OFFSET
1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
MATHEMATICA
cubeFreeQ[n_] := Max[FactorInteger[n][[;; , 2]]] < 3; Mod[Select[Range[200], cubeFreeQ], 2]
PROG
(PARI) lista(kmax) = print1(1, ", "); for(k = 2, kmax, if(vecmax(factor(k)[, 2]) < 3, print1(k % 2, ", ")));
(Python)
from sympy import mobius, integer_nthroot
def A375349(n):
def f(x): return n+x-sum(mobius(k)*(x//k**3) for k in range(1, integer_nthroot(x, 3)[0]+1))
m, k = n, f(n)
while m != k:
m, k = k, f(k)
return m&1 # Chai Wah Wu, Aug 13 2024
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Aug 12 2024
STATUS
approved