|
%I #27 Mar 07 2018 13:24:38
%S 161,209,221,1935,4265,15941,22217,24041,25637,30377,38117,39077,
%T 48617,49097,55877,68441,73817,76457,80357,88457,95237,98117,99941,
%U 105641,110057,115397,122537,130217,131141,136517,143237,147941,148697,152357,154457,159077
%N Numbers n such that n'' = n'+1 where n' and n'' are respectively the first and the second arithmetic derivative of n (A003415).
%C The arithmetic derivative of a(n) is a Giuga's number A007850 (solution of n' = n+1).
%H Reinhard Zumkeller and Donovan Johnson, <a href="/A189639/b189639.txt">Table of n, a(n) for n = 1..1000</a> (first 100 terms from Reinhard Zumkeller)
%e 161' = 30, 161'' = 30' = 31 ==> 161'' = 161'+1 so 161 is a term.
%o (PARI) /* using _Michael B. Porter_'s code from A003415: */
%o A003415(n) = {local(fac); if(n<1, 0, fac=factor(n); sum(i=1, matsize(fac)[1], n*fac[i, 2]/fac[i, 1]))} /* arithmetic derivative */
%o for(n=1,10^6,d1=A003415(n);d2=A003415(d1);if(d2==d1+1,print1(n,", "))); /* show terms */
%o /* _Joerg Arndt_, Apr 25 2011 */
%o (Haskell)
%o import Data.List (elemIndices)
%o a189710 n = a189710_list !! (n-1)
%o a189710_list = elemIndices 0 $
%o zipWith (-) (map a003415 a003415_list) (map pred a003415_list)
%o -- _Reinhard Zumkeller_, May 09 2011
%Y Cf. A003415, A007850, A054377, A132632, A189639.
%K nonn,look
%O 1,1
%A _Giorgio Balzarotti_, Apr 24 2011
|
