login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A070806
Primes p such that cototient(totient(p)) = A070556(p) is a power of 2.
4
3, 5, 7, 13, 17, 29, 97, 113, 193, 257, 449, 509, 769, 7937, 12289, 65537, 114689, 520193, 786433, 7340033, 8388593, 33292289, 33550337, 469762049, 2130706433, 3221225473, 8588886017, 137438691329, 206158430209
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
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