%I #12 Apr 16 2022 09:40:42
%S 2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
%T 0,1,0,0,0,0,0,0,0,0,0,0,0,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,
%U 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0,0,0,0,0,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2
%N Sum of A353350 and its Dirichlet inverse.
%C a(720) = -1 is the first negative term.
%H Antti Karttunen, <a href="/A353349/b353349.txt">Table of n, a(n) for n = 1..65537</a>
%H <a href="/index/Pri#prime_indices">Index entries for sequences computed from indices in prime factorization</a>
%F a(n) = A353348(n) + A353350(n).
%F a(1) = 2, and for n > 1, a(n) = -Sum_{d|n, 1<d<n} A353348(d) * A353350(n/d).
%F a(p) = 0 for all primes p.
%F a(n) = a(A003961(n)) = a(A348717(n)), for all n >= 1.
%t f[p_, e_] := e*2^(PrimePi[p] - 1); s[1] = 1; s[n_] := Boole@Divisible[Plus @@ f @@@ FactorInteger[n], 3]; sinv[1] = 1; sinv[n_] := -DivisorSum[n, sinv[#]*s[n/#] &, # < n &]; a[n_] := s[n] + sinv[n]; Array[a, 100] (* _Amiram Eldar_, Apr 15 2022 *)
%o (PARI)
%o up_to = 16384;
%o DirInverseCorrect(v) = { my(u=vector(#v)); u[1] = (1/v[1]); for(n=2, #v, u[n] = (-u[1]*sumdiv(n, d, if(d<n, v[n/d]*u[d], 0)))); (u) }; \\ Compute the Dirichlet inverse of the sequence given in input vector v
%o A048675(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; };
%o A353350(n) = (0==(A048675(n)%3));
%o v353348 = DirInverseCorrect(vector(up_to,n,A353350(n)));
%o A353348(n) = v353348[n];
%o A353349(n) = (A353348(n)+A353350(n));
%Y Cf. A003961, A048675, A348717, A353348, A353350.
%K sign
%O 1,1
%A _Antti Karttunen_, Apr 15 2022