OFFSET
1,2
LINKS
PROG
(PARI) A328617(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, 2] = q*f[k, 1] + lift(1/Mod(m, f[k, 1])))); factorback(f); };
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Antti Karttunen, Oct 23 2019
STATUS
approved