|
PROG
|
(PARI) is(n)=numdiv(n)==38 && numdiv(n+1)==38
(PARI) has(n)=if(n%4==2, ispower(n/2, 18, &n) && isprime(n), bitand(n, 524287)==262144 && isprime(n>>18) && n>262144) \\ check if n is even with 38 divisors
list(lim)=my(v=List(), t); forprime(p=2, sqrtnint(lim\=1, 37), t=p^37; if(has(t+1), listput(v, t)); if(has(t-1), listput(v, t-1))); forprime(p=3, sqrtnint(lim\3, 18), my(p18=p^18); forprime(q=3, lim\p18, if(p==q, next); t=p18*q; if(has(t+1), listput(v, t)); if(has(t-1), listput(v, t-1)))); Set(v)
|