OFFSET
1,1
COMMENTS
If k is an even term greater than 2 of A051487 then 2k is another term.
This sequence lists the initial term k_0 of each infinite subsequence that is solution of the equation phi(k) = phi(k - phi(k)).
About 2: one could argue that 2 is primitive since it is not the double of any previous term of A051487, but as 2^k is not solution for n>1, 2 is not primitive.
Each k_0 is of the form k_0 = 6*m with m odd.
If p > 3 is a Sophie Germain prime, then every m = 2*3*p^q, q >=2 is a term because phi(m) = phi(m-phi(m)) = 2*(p-1)*p^(q-1); the first terms that are not of this form are 6, 2310, 5874, ... (see examples).
REFERENCES
Richard K. Guy, Unsolved Problems in Number Theory, 3rd Edition, Springer, 2004, Section B42, p. 150.
EXAMPLE
a(1) = 6 because every k = 3*2^m, m >= 1 satisfies phi(k) = phi(k-phi(k)) = 2^m, and k_0 = 6 is the smallest term of this subsequence of A051487.
a(2) = 150 because every k = 3*5^2*2^m, m >= 1 satisfies phi(k) = phi(k-phi(k)) = 5*2^(m+2) and k_0 = 150 is the smallest term of this subsequence of A051487.
a(3) = 726 because every k = 3*11^2*2^m, m >= 1 satisfies phi(k) = phi(k-phi(k)) = 5*11*2^(m+1) and k_0 = 726 is the smallest term of this subsequence of A051487.
a(5) = 2310 because every k = 3*5*7*11*2^m, m >= 1 satisfies phi(k) = phi(k-phi(k)) = 3*5*2^(m+4) and k_0 = 2310 is the smallest term of this subsequence of A051487.
MAPLE
with(numtheory):
for q from 0 to 13800 do
m := 6*(2*q+1);
if phi(m) = phi(m-phi(m)) then print(m); else fi; od:
PROG
(PARI) isdouble(n, list)= {my(v = Vecrev(list)); for(k=1, #v, if (n == 2*v[k], return(1)); ); }
lista(nn) = {my(list = List(), listp = List()); for (n=3, nn, if (eulerphi(n) == eulerphi(n - eulerphi(n)), if (!isdouble(n, list), listput(listp, n)); listput(list, n); ); ); Vec(listp); } \\ Michel Marcus, Aug 06 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Bernard Schott, Aug 06 2021
STATUS
approved