|
%I #16 Oct 28 2019 20:01:01
%S 0,0,1,2,4,6,7,30,12,12,17,210,20,2310,67,28,32,30030,33,510510,44,
%T 104,431,9699690,52,60,4633,54,148,223092870,71,6469693230,80,652,
%U 60077,192,84,200560490130,1021039,6956,108,7420738134810,229,304250263527210,884,114,19399403,13082761331670030,128,420,145,90124,9292,614889782588491410,135,1116,324
%N The first primorial based variant of arithmetic derivative: a(prime(i)) = A002110(i-1), where prime(i) = A000040(i), a(u*v) = a(u)*v + u*a(v), with a(0) = a(1) = 0.
%H Antti Karttunen, <a href="/A328768/b328768.txt">Table of n, a(n) for n = 0..1001</a>
%H <a href="/index/Pri#primorial_numbers">Index entries for sequences related to primorial numbers</a>
%F a(n) = n * Sum e_j * A276085(p_j)/p_j for n = Product p_j^e_j, where for primes p, A276085(p) = A002110(A000720(p)-1).
%F a(n) = n * Sum e_j * (p_j)#/(p_j^2) for n = Product p_j^e_j with (p_j)# = A034386(p_j).
%F For all n >= 0, A276150(a(n)) = A328771(n).
%o (PARI)
%o A002110(n) = prod(i=1,n,prime(i));
%o A328768(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]*A002110(primepi(f[i,1])-1)/f[i, 1]));
%o (PARI) A328768(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]*A002110(primepi(f[i, 1]))/(f[i, 1]^2)));
%Y Cf. A002110, A108951, A276085, A276150, A328771.
%Y Cf. also A003415, A258851, A328769, A328845, A328846.
%K nonn
%O 0,4
%A _Antti Karttunen_, Oct 28 2019
|
