OFFSET
1,1
COMMENTS
Prime terms are in A058383.
See A266276(n) = the smallest numbers k such that phi(k) = n * phi(k-1) for n >=1: 2, 3, 7, 1261, 11242771, ...
Number of terms < 10^k: 1, 7, 17, 29, 41, 86, 205, 446, 1001, 2295, ..., . - Robert G. Wilson v, Jan 24 2016
All terms are == +-1 (mod 6) but mostly 1 (> 95%). - Robert G. Wilson v, Jan 24 2016
LINKS
Robert G. Wilson v, Table of n, a(n) for n = 1..2467 (first 405 terms from G. C. Greubel)
FORMULA
a(n) = A067143(n) + 1.
EXAMPLE
19 is in the sequence because phi(19) = 18 = 3*phi(18) = 3*6.
MATHEMATICA
Select[Range[5000], EulerPhi[ # ]==3*EulerPhi[ #-1]&] (* G. C. Greubel, Dec 26 2015 *)
PROG
(Magma) [n: n in [2..2*10^5] | EulerPhi(n) eq 3*EulerPhi(n-1)]; // Vincenzo Librandi, Dec 26 2015
(PARI) isok(n) = eulerphi(n) == 3*eulerphi(n-1); \\ Michel Marcus, Dec 27 2015
(PARI) lista(nn) = for(n=1, nn, if(eulerphi(n) == 3*eulerphi(n-1), print1(n, ", "))); \\ Altug Alkan, Jan 24 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Dec 26 2015
STATUS
approved