OFFSET
1,1
COMMENTS
LINKS
Freimut Marschner, Table of n, a(n) for n = 1..100000
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]):
seq(ad(3*n+1), n=1..100); # Robert Israel, Aug 25 2014
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
Freimut Marschner, Aug 06 2014
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