miracl.lib库

/*
 *   Program to calculate mersenne primes
 *   using Lucas-Lehmer test - Knuth p.391
 *
 *   Try this only in a 32-bit （or better!） environment
 *
 */

#include 
#include “miracl.h“
#define LIMIT 100000

int main（）
{ /* calculate mersenne primes */
    BOOL compo;
    big LmT;
    int ikqpr;
    miracl *mip=mirsys（50000）;
    L=mirvar（0）;
    m=mirvar（0）;
    T=mirvar（0）;
    gprime（LIMIT）;
    for （k=1;;k++）
    { /* test only prime exponents */
        q=mip->PRIMES[k];
        if （q==0） break;
        expb2（qm）;
        decr（m1m）;    /* m=2^q-1 */

/* try to find a factor. Should perhaps keep trying over a bigger range... */

        compo=FALSE;
        for（i=2;i<16*q;i+=2）
        { /* prime factors （if they exist） are always *
           * of the form i*q+1 and