OFFSET
1,4
COMMENTS
Prime shifted analog of A160595.
LINKS
MAPLE
f:= proc(n) local F, x, p, t;
F:= ifactors(n)[2];
x:= mul(nextprime(t[1])^t[2], t=F);
p:= numtheory:-phi(x);
p/igcd(x-1, p)
end proc:
map(f, [$1..100]); # Robert Israel, Dec 28 2020
MATHEMATICA
a[n_] := Module[{x, p, e, phi}, x = Product[{p, e} = pe; NextPrime[p]^e, {pe, FactorInteger[n]}]; phi = EulerPhi[x]; phi/GCD[x-1, phi]];
Array[a, 100] (* Jean-François Alcover, Jan 04 2022 *)
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Dec 28 2020
STATUS
approved