|
|
A189639
|
|
Numbers n such that n'' = n'+1 where n' and n'' are respectively the first and the second arithmetic derivative of n (A003415).
|
|
3
|
|
|
161, 209, 221, 1935, 4265, 15941, 22217, 24041, 25637, 30377, 38117, 39077, 48617, 49097, 55877, 68441, 73817, 76457, 80357, 88457, 95237, 98117, 99941, 105641, 110057, 115397, 122537, 130217, 131141, 136517, 143237, 147941, 148697, 152357, 154457, 159077
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
The arithmetic derivative of a(n) is a Giuga's number A007850 (solution of n' = n+1).
|
|
LINKS
|
|
|
EXAMPLE
|
161' = 30, 161'' = 30' = 31 ==> 161'' = 161'+1 so 161 is a term.
|
|
PROG
|
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 */
for(n=1, 10^6, d1=A003415(n); d2=A003415(d1); if(d2==d1+1, print1(n, ", "))); /* show terms */
(Haskell)
import Data.List (elemIndices)
a189710 n = a189710_list !! (n-1)
a189710_list = elemIndices 0 $
zipWith (-) (map a003415 a003415_list) (map pred a003415_list)
|
|
CROSSREFS
|
|
|
KEYWORD
|
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|