|
PROG
|
(PARI) is(n)=numdiv(n)==34 && numdiv(n+1)==34
(PARI) has(n)=if(n%4==2, ispower(n/2, 16, &n) && isprime(n), bitand(n, 131071)==65536 && isprime(n>>16) && n>65536) \\ check if n is even with 34 divisors
list(lim)=my(v=List(), t); forprime(p=2, sqrtnint(lim\=1, 33), t=p^33; if(has(t+1), listput(v, t)); if(has(t-1), listput(v, t-1))); forprime(p=3, sqrtnint(lim\3, 16), my(p16=p^16); forprime(q=3, lim\p16, if(p==q, next); t=p16*q; if(has(t+1), listput(v, t)); if(has(t-1), listput(v, t-1)))); Set(v)
|