|
%I #8 Jan 06 2021 09:07:27
%S -1,0,0,2,0,7,0,6,6,13,0,13,0,19,22,10,0,19,0,23,32,31,0,-3,20,37,24,
%T 33,0,21,0,6,52,49,58,-36,0,55,62,-9,0,29,0,53,52,67,0,-87,42,53,82,
%U 63,0,-29,94,-15,92,85,0,-219,0,91,74,-22,112,45,0,83,112,45,0,-257,0,109,82,93,136,53,0,-165,42,121
%N Dirichlet inverse of sequence b(n) = 1-A318833(n).
%H Antti Karttunen, <a href="/A340367/b340367.txt">Table of n, a(n) for n = 1..16384</a>
%F a(1) = -1, for n > 1, a(n) = Sum_{d|n, d<n} (1-A318833(n/d)) * a(d).
%o (PARI)
%o up_to = 65537;
%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 (correctly!)
%o A023900(n) = factorback(apply(p -> 1-p, factor(n)[, 1]));
%o A318833(n) = (n+A023900(n));
%o v340367 = DirInverseCorrect(vector(up_to, n, 1-A318833(n)));
%o A340367(n) = v340367[n];
%o (PARI)
%o \\ Or as:
%o A340367(n) = if(1==n, -1, sumdiv(n, d, if(d<n, (1-A318833(n/d))*A340367(d), 0)));
%Y Cf. A000010, A023900, A318833, A340198.
%Y Cf. also A340140, A340197 for similar definitions.
%K sign
%O 1,4
%A _Antti Karttunen_, Jan 05 2021
|
