OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..80
EXAMPLE
PROG
(PARI) a034448(n) = {my(f, i, p, e); f=factor(n); for(i=1, #f~, p=f[i, 1]; e=f[i, 2]; f[i, 1]=p^e+1; f[i, 2]=1); factorback(f)};
a188999(n) = {my(f, i, p, e); f=factor(n); for(i=1, #f~, p=f[i, 1]; e=f[i, 2]; f[i, 1]=if(e%2, (p^(e+1)-1)/(p-1), (p^(e+1)-1)/(p-1)-p^(e/2)); f[i, 2]=1); factorback(f)};
isok(n) = (a034448(a188999(n))%n) == 0;
CROSSREFS
KEYWORD
nonn
AUTHOR
Tomohiro Yamada, Jan 16 2024
STATUS
approved