has(n)=numdiv(n)==42;
list(lim)={
	my(v=List(), t);
	forprime(p=2, sqrtnint(lim\=1, 41), t=p^41; if(has(t+1), listput(v, t)); if(has(t-1), listput(v, t-1)));
	forprime(p=2, sqrtnint(lim\2, 20), my(p20=p^20); forprime(q=2, lim\p20, if(p==q, next); t=p20*q; if(has(t+1), listput(v, t)); if(has(t-1), listput(v, t-1))));
	forprime(p=2, sqrtnint(lim\2, 13), my(p13=p^13); forprime(q=2, sqrtint(lim\p13), if(p==q, next); t=p13*q^2; if(has(t+1), listput(v, t)); if(has(t-1), listput(v, t-1))));
	forprime(p=2, sqrtnint(lim\32, 6), my(p6=p^6); forprime(q=2, sqrtnint(lim\p6,5), if(p==q, next); t=p6*q^5; if(has(t+1), listput(v, t)); if(has(t-1), listput(v, t-1))));
	forprime(p=2, sqrtnint(lim\12, 6), my(p6=p^6); forprime(q=2, sqrtint(lim\(2*p6)), if(p==q, next); my(p6q2=p6*q^2); forprime(r=2, lim\p6q2, if(p==r || q==r, next); t=p6q2*r; if(has(t+1), listput(v, t)); if(has(t-1), listput(v, t-1)))));
	Set(v);
}