MersenneFermat(lb,ub)={
forprime(p=lb,ub,
if(!(lift(Mod(2,2^p-1)^(2^(p-1)-1))==1),print(p)))
}