OFFSET
1,1
COMMENTS
These numbers m appear to satisfy cototient(m) > totient(m) or 2*phi(m) < m; they seem to be the missing terms mentioned in A067800. - Labos Elemer, May 08 2003
All elements in this sequence must have 2*phi(m) < m, but not the reverse. See A118700. - Franklin T. Adams-Watters, May 21 2006
The numbers of terms that do not exceed 10^k, for k = 3, 4, ..., are 11, 108, 1139, 11036, 111796, ... . Apparently, the asymptotic density of this sequence exists and equals 0.011... . - Amiram Eldar, Nov 21 2024
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
Max Alekseyev, PARI/GP Scripts for Miscellaneous Math Problems (invphi.gp).
MAPLE
N:= 10^4: # to get all terms <= N
PhiE:= map(numtheory:-phi, [seq(i, i=2..N, 2)]):
A:= NULL:
for n from 1 to N by 2 do
t:= numtheory:-phi(n);
if 2*t < n and member(t, PhiE[1..(n-1)/2]) then A:= A, n fi;
od:
A; # Robert Israel, Jan 06 2017
PROG
(PARI) is(m) = m%2 && #select(k -> !(k%2) && k < m, invphi(eulerphi(m))) > 0; \\ Amiram Eldar, Nov 21 2024, using Max Alekseyev's invphi.gp
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
STATUS
approved