OFFSET
0,1
COMMENTS
The ratio of adjacent terms is 2 except for five terms (if there are exactly five Fermat primes). - T. D. Noe, Jun 21 2012
LINKS
Amiram Eldar, Table of n, a(n) for n = 0..1000 (terms 0..100 from T. D. Noe)
Max Alekseyev, PARI/GP Scripts for Miscellaneous Math Problems (invphi.gp).
FORMULA
Assuming there are only 5 Fermat primes (A019434), a(n) = 2^(n-30)*(2^32-1) for n >= 31.
EXAMPLE
For n = 6, 2^n = 64; the solutions of phi(x) = 64 are {85,128,136,160,170,192,204,240}; the largest is a(6) = 240.
MATHEMATICA
phiinv[ n_, pl_ ] := Module[ {i, p, e, pe, val}, If[ pl=={}, Return[ If[ n==1, {1}, {} ] ] ]; val={}; p=Last[ pl ]; For[ e=0; pe=1, e==0||Mod[ n, (p-1)pe/p ]==0, e++; pe*=p, val=Join[ val, pe*phiinv[ If[ e==0, n, n*p/pe/(p-1) ], Drop[ pl, -1 ] ] ] ]; Sort[ val ] ]; phiinv[ n_ ] := phiinv[ n, Select[ 1+Divisors[ n ], PrimeQ ] ]; Table[ phiinv[ 2^n ][ [ -1 ] ], {n, 0, 30} ] (* phiinv[ n, pl ] = list of x with phi(x)=n and all prime divisors of x in list pl. phiinv[ n ] = list of x with phi(x)=n *)
PROG
(PARI) a(n) = invphiMax(2^n); \\ Amiram Eldar, Nov 11 2024, using Max Alekseyev's invphi.gp
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Nov 30 2000
EXTENSIONS
Edited by Dean Hickerson, Jan 25 2002
STATUS
approved