OFFSET
1,1
EXAMPLE
Powers of 2 observable in A070556[this sequence] = {1, 2, 4, 8, 16, 64, 128, 256, 512, 4096, 8192, 32768, 65536, 262144, 524288, ...}. For F(m), Fermat prime:phi[F(m)]=2^m, cototient[2^m]=2^(m-1); if n=113: phi[113]=112, cototient[112]=112-48=64, so 113 is in this sequence.
MATHEMATICA
Do[s= EulerPhi[n]-EulerPhi[EulerPhi[n]]; If[IntegerQ[Log[2, s]]&&PrimeQ[n], Print[n]], {n, 1, 10000000}]
PROG
(PARI) ispow2(n)=n==1<<valuation(n, 2);
forprime(p=2, 4e9, if(ispow2(p-1-eulerphi(p-1)), print1(p", "))) \\ Charles R Greathouse IV, May 17 2011
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Labos Elemer, May 08 2002
EXTENSIONS
a(20)-a(27) from Donovan Johnson, Feb 06 2010
a(28)-a(29) from Charles R Greathouse IV, May 17 2011
STATUS
approved