STATUS
proposed
approved
proposed
approved
editing
proposed
Robert Israel, <a href="/A306711/b306711.txt">Table of n, a(n) for n = 1..10000</a>
psi:= k -> mul((t+1)/t, t=numtheory:-factorset(k))*k:
select(t -> igcd(t, psi(t)) <> igcd(t, numtheory:-phi(t)), [$1..1000]); # Robert Israel, Apr 28 2019
approved
editing
reviewed
approved
proposed
reviewed
editing
proposed
Numbers m such that A306695(m) = m is an elementare terms.
proposed
editing
editing
proposed
(PARI) dpsi(n) = n * sumdivmult(n, d, issquarefree(d)/d); \\ A001615
isok(k) = gcd(k, eulerphi(k)) != gcd(k, dpsi(k)); \\ Michel Marcus, Mar 21 2019
The ratio gcd(k, phi(k)) : gcd(k, psi(k)), when reduced, contains at least one prime number that divides a(n).
proposed
editing