OFFSET
1,1
LINKS
Mauro Fiorentini, Table of n, a(n) for n = 1..154 (all terms for n up to 10^9).
EXAMPLE
phi(629) = 576 = 4*phi(630).
MAPLE
A:= NULL:
y:= numtheory:-phi(1):
for n from 1 to 10^6 do
x:= numtheory:-phi(n+1);
if y = 4*x then A:= A, n fi;
y:= x;
od:
A; # Robert Israel, Apr 15 2015
MATHEMATICA
Select[Range@ 1000000, EulerPhi@ # == 4 EulerPhi[# + 1] &] (* Michael De Vlieger, Apr 13 2015 *)
Position[Partition[EulerPhi[Range[6*10^6]], 2, 1], _?(#[[1]]==4#[[2]]&), {1}, Heads->False]//Flatten (* Harvey P. Dale, Sep 18 2016 *)
PROG
(PARI) s=[]; for(n=1, 1000000, if(eulerphi(n)==4*eulerphi(n+1), s=concat(s, n))); s \\ Colin Barker, Apr 13 2015
(Magma) [n: n in [1..10^7] | EulerPhi(n) eq 4*EulerPhi(n+1)]; // Vincenzo Librandi, Apr 14 2015
(Sage) [n for n in (1..1000000) if euler_phi(n) == 4*euler_phi(n+1)]; # Bruno Berselli, Apr 14 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Mauro Fiorentini, Apr 13 2015
STATUS
approved