|
|
A245838
|
|
Arithmetic derivative of (3*n + 1), n >= 1, (A016777)'.
|
|
1
|
|
|
4, 1, 7, 1, 32, 1, 13, 10, 32, 1, 19, 1, 68, 1, 25, 14, 56, 16, 31, 1, 192, 1, 59, 1, 80, 1, 43, 22, 140, 20, 49, 1, 140, 1, 55, 1, 240, 28, 61, 22, 128, 1, 101, 26, 212, 1, 73, 34, 152, 1, 113, 1, 432, 1, 85, 26, 176, 95, 91, 1, 284, 28, 143, 1, 252, 1, 103
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Comparing A016777(n) = 3n+1 and its arithmetic derivative a(n) = (A016777(n))' allows us to select the primes in A016777, because (prime(n))' = 1.
|
|
LINKS
|
|
|
FORMULA
|
|
|
EXAMPLE
|
a(2) = (3*2 + 1)' = (7)' = 1, a(5846) = (A016777(5846))' = (3*5846 + 1)' = (17539)' = 1.
|
|
MAPLE
|
ad:= n -> n*add(f[2]/f[1], f=ifactors(n)[2]):
|
|
PROG
|
(PARI) a(n) = my(n=3*n+1); sum(i=1, #f=factor(n)~, n/f[1, i]*f[2, i]); \\ Michel Marcus, Aug 12 2014
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
Edited: in name n>=0 replaced by n>=1. Example corrected. - Wolfdieter Lang, Oct 14 2014
Missing comma in data at a(63) inserted by Andrew Howroyd, Feb 22 2018
|
|
STATUS
|
approved
|
|
|
|