login
Dirichlet inverse of A358680, where A358680 is the characteristic function of the numbers with even arithmetic derivative (A003415).
16

%I #25 Jan 17 2023 16:31:12

%S 1,0,0,-1,0,0,0,-1,-1,0,0,-1,0,0,-1,0,0,0,0,-1,-1,0,0,-1,-1,0,0,-1,0,

%T 0,0,1,-1,0,-1,1,0,0,-1,-1,0,0,0,-1,0,0,0,1,-1,0,-1,-1,0,0,-1,-1,-1,0,

%U 0,1,0,0,0,1,-1,0,0,-1,-1,0,0,1,0,0,0,-1,-1,0,0,1,0,0,0,1,-1,0,-1,-1,0,0,-1,-1,-1,0,-1,3,0,0,0,1,0,0,0,-1,0

%N Dirichlet inverse of A358680, where A358680 is the characteristic function of the numbers with even arithmetic derivative (A003415).

%C As A358680 is not multiplicative, neither is this sequence.

%C For all odd numbers n with an odd number of prime factors (with mult.), a(n) = 0. Proof: In the convolution formula, when n is any term of A067019, either the divisor (n/d) or d (but not both) has an odd number of prime factors. As A358680 is zero for all A067019 (see A235991), it is easy to show by induction that also a(n) is zero for all such numbers.

%C For all numbers n of the form 4m+2, a(n) = 0. Proof: In the convolution formula, when n is any term of A016825, the other divisor of the pair {(n/d), d} is always odd, and the other is always even (particularly of the form 4u+2). As A358680 is zero for all A016825 (see A235991), it is easy to show by induction that also a(n) is zero for all such numbers.

%C Therefore, nonzero values (including any odd values, see A359783) occur only on a subset of A235992, and A359781(n) <= A358680(n).

%C The above proof can be made simpler by realizing that A235992 is a multiplicative semigroup and therefore for any n that is in its complement A235991, at least the other of the divisors d or n/d must be in A235991. Compare also to a proof given in A353348. - _Antti Karttunen_, Jan 17 2023

%H Antti Karttunen, <a href="/A359780/b359780.txt">Table of n, a(n) for n = 1..65537</a>

%F a(1) = 1, and for n > 1, a(n) = -Sum_{d|n, d<n} A358680(n/d) * a(d).

%o (PARI)

%o A358680(n) = if(n<=1, 1, my(f=factor(n)); 0==((n*sum(i=1, #f~, f[i, 2]/f[i, 1]))%2));

%o memoA359780 = Map();

%o A359780(n) = if(1==n,1,my(v); if(mapisdefined(memoA359780,n,&v), v, v = -sumdiv(n,d,if(d<n,A358680(n/d)*A359780(d),0)); mapput(memoA359780,n,v); (v)));

%Y Cf. A003415, A016825, A067019, A235992, A358680, A359781 (parity of terms), A359782 (positions of even terms), A359783 (of odd terms), A359784.

%Y Cf. also A359763 [= a(A003961(n))], A359793, A359823 and A353348.

%K sign

%O 1,96

%A _Antti Karttunen_, Jan 13 2023