A328845 The first Fibonacci based variant of arithmetic derivative: a(p) = A000045(p) for prime p, a(u*v) = a(u)*v + u*a(v), with a(0) = a(1) = 0.

0, 0, 1, 2, 4, 5, 7, 13, 12, 12, 15, 89, 20, 233, 33, 25, 32, 1597, 33, 4181, 40, 53, 189, 28657, 52, 50, 479, 54, 80, 514229, 65, 1346269, 80, 289, 3211, 100, 84, 24157817, 8381, 725, 100, 165580141, 127, 433494437, 400, 105, 57337, 2971215073, 128, 182, 125, 4825, 984, 53316291173, 135, 500, 188, 12581, 1028487, 956722026041, 160

Offset: 0

Author

Antti Karttunen, Oct 28 2019

Keywords

• nonn

Links

• Antti Karttunen, Table of n, a(n) for n = 0..1001

Crossrefs

Cf. A000040, A000045, A007895, A030426, A113175, A328847, A374201.

Cf. A374046 (indices of even terms), A374047 (of odd terms), A374122 (of multiples of 3), A374202 (2-adic valuation), A374203 (3-adic valuation), A374205 (5-adic valuation), A374125 [a(n) mod 360].

Cf. A374106 [gcd(a(n), A113177(n))], A374035 [gcd(a(n), A328846(n))], A374116 [gcd(a(n), A328768(n))].

For variants of the same formula, see A003415, A258851, A328768, A328769, A328846, A371192.

Programs

Mathematica

A328845[n_] := If[n <= 1, 0, n*Total[MapApply[#2*Fibonacci[#]/# &, FactorInteger[n]]]];
Array[A328845, 100, 0] (* Paolo Xausa, Dec 16 2024 *)

PARI

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

Formula

a(n) = n * Sum e_j * A000045(p_j)/p_j for n = Product p_j^e_j.

a(A000040(n)) = A030426(n).

A007895(a(n)) = A328847(n).