I fixed previous errors, but now it says 2^41-1 is a prime :( Anyway, I send the new version -- A.Z.
Attached file: strongprime-old.lisp