OFFSET
1,2
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = n if n is in A289509, otherwise a(n) = 1.
MAPLE
g:= proc(n) uses numtheory; igcd(op(map(pi, factorset(n))))=1 end proc:
seq(`if`(g(n), n, 1), n=1..100); # Robert Israel, Sep 19 2019
MATHEMATICA
Table[Max[Select[Divisors[n], #==1||GCD@@PrimePi/@First/@FactorInteger[#]==1&]], {n, 100}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Sep 17 2019
STATUS
approved