This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A085731 #49 Feb 23 2024 07:27:47 %S A085731 1,1,1,4,1,1,1,4,3,1,1,4,1,1,1,16,1,3,1,4,1,1,1,4,5,1,27,4,1,1,1,16,1, %T A085731 1,1,12,1,1,1,4,1,1,1,4,3,1,1,16,7,5,1,4,1,27,1,4,1,1,1,4,1,1,3,64,1, %U A085731 1,1,4,1,1,1,12,1,1,5,4,1,1,1,16,27,1,1,4,1,1,1,4,1,3,1,4,1,1,1,16 %N A085731 Greatest common divisor of n and its arithmetic derivative. %C A085731 a(n) = 1 iff n is squarefree (A005117), cf. A068328. %C A085731 This sequence is very probably multiplicative. - _Mitch Harris_, Apr 19 2005 %H A085731 T. D. Noe, <a href="/A085731/b085731.txt">Table of n, a(n) for n = 1..10000</a> %H A085731 Victor Ufnarovski and Bo Ahlander, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL6/Ufnarovski/ufnarovski.html">How to Differentiate a Number</a>, J. Integer Seqs., Vol. 6, 2003, #03.3.4. %F A085731 a(n) = GCD(n, A003415(n)). %F A085731 Multiplicative with a(p^e) = p^e if p divides e; a(p^e) = p^(e-1) otherwise. - _Eric M. Schmidt_, Oct 22 2013 %F A085731 From _Antti Karttunen_, Feb 28 2021: (Start) %F A085731 Thus a(A276086(n)) = A328572(n), by the above formula and the fact that A276086 is a permutation of A048103. %F A085731 a(n) = n / A083346(n) = A190116(n) / A086130(n). (End) %t A085731 d[0] = d[1] = 0; d[n_] := d[n] = n*Total[Apply[#2/#1 &, FactorInteger[n], {1}]]; a[n_] := GCD[n, d[n]]; Table[a[n], {n, 1, 96}] (* _Jean-François Alcover_, Feb 21 2014 *) %t A085731 f[p_, e_] := p^If[Divisible[e, p], e, e - 1]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* _Amiram Eldar_, Oct 31 2023 *) %o A085731 (Haskell) %o A085731 a085731 n = gcd n $ a003415 n -- _Reinhard Zumkeller_, May 10 2011 %o A085731 (PARI) a(n) = {my(f = factor(n)); for (i=1, #f~, if (f[i,2] % f[i,1], f[i,2]--);); factorback(f);} \\ _Michel Marcus_, Feb 14 2016 %Y A085731 Cf. A003415, A005117, A048103, A068328, A083346, A083347, A086130, A129251, A189100, A189036, A189103, A190116, A276086, A327858, A327938, A328572, A340070. %Y A085731 Cf. A072873 (positions of records), A268398 (partial sums). %K A085731 nonn,easy,mult %O A085731 1,4 %A A085731 _Reinhard Zumkeller_, Jul 20 2003