OFFSET
1,2
COMMENTS
Also a(n) is the maximum value of the n-th column of A342255. - Jianing Song, Aug 09 2022
LINKS
Jianing Song, Table of n, a(n) for n = 1..10000
FORMULA
a(n) is the largest prime factor of n if n is in A342256, 1 otherwise.
EXAMPLE
gcd(6, Phi_6(a)) = gcd(6, a^2-a+1) = 3 for a == 2 (mod 3), 1 otherwise, so {gcd(6, Phi_6(a)): a in Z} has period 3, hence a(6) = 3.
gcd(12, Phi_12(a)) = gcd(12, a^4-a^2+1) = 1 for all n, so {gcd(12, Phi_12(a)): a in Z} has period 1, hence a(12) = 1.
PROG
(PARI) a(n) = if(n>1, my(L=factor(n), d=omega(n), p=L[d, 1]); if((p-1)%(n/p^L[d, 2])==0, p, 1), 1)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Jianing Song, Mar 07 2021
STATUS
approved