OFFSET
1,2
LINKS
PROG
(PARI) A328619(n) = { my(f = factor(n), m, q); for(k=1, #f~, q = (f[k, 2]\f[k, 1]); m = (f[k, 2]%f[k, 1]); if(m&&(f[k, 1]!=2), f[k, 2] = q*f[k, 1] + lift(Mod(m, f[k, 1])/2))); factorback(f); };
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Antti Karttunen, Oct 23 2019
STATUS
approved