OFFSET
1,1
COMMENTS
EXAMPLE
For k = 3040 = 2^5 * 5 * 19, phi(k) = 2^4 * 4 * 18 = 2^7 * 3^2, tau(k) = 6* 2* 2 = 2^3 * 3, sigma(k) = 2^3 * 3^3 * 5 *7, gcd(k,phi(k)) + gcd(k tau(k)) = 2^5 + 2^3 = 40 and gcd(k,sigma(k)) = 2^3 * 5 = 40.
MATHEMATICA
aQ[n_] := GCD[n, EulerPhi[n]] + GCD[n, DivisorSigma[0, n]] == GCD[n, DivisorSigma[1, n]]; Select[Range[300000], aQ] (* Amiram Eldar, Oct 23 2019 *)
PROG
(Magma) [k: k in [1..300000]| Gcd(k, NumberOfDivisors(k))+Gcd(k, EulerPhi(k)) eq Gcd(k, SumOfDivisors(k))];
(PARI) isok(k) = gcd(k, numdiv(k)) + gcd(k, eulerphi(k)) == gcd(k, sigma(k)); \\ Michel Marcus, Oct 24 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Marius A. Burtea, Oct 18 2019
STATUS
approved