login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A328768
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.
23
0, 0, 1, 2, 4, 6, 7, 30, 12, 12, 17, 210, 20, 2310, 67, 28, 32, 30030, 33, 510510, 44, 104, 431, 9699690, 52, 60, 4633, 54, 148, 223092870, 71, 6469693230, 80, 652, 60077, 192, 84, 200560490130, 1021039, 6956, 108, 7420738134810, 229, 304250263527210, 884, 114, 19399403, 13082761331670030, 128, 420, 145, 90124, 9292, 614889782588491410, 135, 1116, 324
OFFSET
0,4
FORMULA
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).
a(n) = n * Sum e_j * (p_j)#/(p_j^2) for n = Product p_j^e_j with (p_j)# = A034386(p_j).
For all n >= 0, A276150(a(n)) = A328771(n).
PROG
(PARI)
A002110(n) = prod(i=1, n, prime(i));
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]));
(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)));
CROSSREFS
Cf. A042965 (indices of even terms), A016825 (of odd terms), A152822 (antiparity of terms), A373992 (indices of multiples of 3), A374212 (2-adic valuation), A374213 (3-adic valuation), A374123 [a(n) mod 360].
Cf. A374031 [gcd(a(n), A276085(n))], A374116 [gcd(a(n), A328845(n))].
For variants of the same formula, see A003415, A258851, A328769, A328845, A328846, A371192.
Sequence in context: A061937 A065804 A374123 * A275693 A135033 A317647
KEYWORD
nonn
AUTHOR
Antti Karttunen, Oct 28 2019
STATUS
approved

  NODES
orte 1
see 2
Story 1