OFFSET
1,2
COMMENTS
Squarefree terms are 1, 10, 70, 78, 110, 174, 222, 230, 246, 426, 438, ...
Common terms of this sequence and A033632 are 1, 14406, 544500, 141118050, ...
From Robert Israel, Sep 03 2017: (Start)
If an even number m is in the sequence, then so is 2*m.
Are there any odd terms other than 1? (End)
a(1) = 1 is the only odd term. LHS of equation allows for 1 and 3 but only for k <= 6. RHS allows for 1 and only for k = 1. - Torlach Rush, Jul 28 2023
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
psi:= proc(n) n*mul((1+1/i[1]), i=ifactors(n)[2]) end:
select(psi @ numtheory:-phi = numtheory:-phi @ psi, [$1..1000]); # Robert Israel, Sep 03 2017
MATHEMATICA
f[n_] := n Sum[MoebiusMu[d]^2/d, {d, Divisors@ n}]; Select[Range[10^3], f[EulerPhi@ #] == EulerPhi[f@ #] &] (* Michael De Vlieger, Sep 03 2017 *)
PROG
(PARI) a001615(n) = my(f=factor(n)); prod(i=1, #f~, f[i, 1]^f[i, 2] + f[i, 1]^(f[i, 2]-1));
isok(n) = eulerphi(a001615(n))==a001615(eulerphi(n)); \\ after Charles R Greathouse IV at A001615
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Altug Alkan, Sep 03 2017
STATUS
approved