|
%I #15 Nov 03 2022 16:35:52
%S 1,4,9,4,25,36,49,16,27,100,121,9,169,196,225,8,289,108,361,50,441,
%T 484,529,144,125,676,27,49,841,900,961,64,1089,1156,1225,108,1369,
%U 1444,1521,400,1681,1764,1849,121,675,2116,2209,144,343,500,2601,338,2809,36
%N Denominators of arithmetic derivative of 1/n: -A003415(n)/n^2.
%H Alois P. Heinz, <a href="/A068238/b068238.txt">Table of n, a(n) for n = 1..10000</a>
%p d:= n-> n*add(i[2]/i[1], i=ifactors(n)[2]):
%p a:= n-> denom(-d(n)/n^2):
%p seq(a(n), n=1..80); # _Alois P. Heinz_, Jun 07 2015
%t d[n_] := If[n < 2, 0, n Sum[f[[2]]/f[[1]], {f, FactorInteger[n]}]];
%t a[n_] := Denominator[-d[n]/n^2];
%t Array[a, 80] (* _Jean-François Alcover_, Mar 12 2019 *)
%o (Python)
%o from fractions import Fraction
%o from sympy import factorint
%o def A068238(n): return Fraction(sum((Fraction(e,p) for p,e in factorint(n).items())),n).denominator # _Chai Wah Wu_, Nov 03 2022
%Y Cf. A003415, A068237 (numerators).
%K nonn,frac
%O 1,2
%A _Reinhard Zumkeller_, Feb 23 2002
|
