/*
*  Riesel's Number Research
*  Sierpinski's Number Research 
*
*  Author	: Rosario Turco
*  Date		: 16/04/2010 
*  Revision	: 1.0           
*/


/*
*   Riesel's numbers Research 
*   Sierpinski's Number Research 
*   R. Turco
*
*/

RSL(start=3, end=11, s=1, n=1000, r=1)=local();{ 
 
 if( start%2 == 0, error("inserire uno start dispari"););
 i = start;

 while(i<end+1,
       RSLR(i,s,n,r);   
       i = i+2;
 );
 print("start : ", start," end : ", end);
}

RSLR(val,s=1,n=1000,r=1)=local(out="",trovato=0); { 

    for(j=s,n,

        if( r == 1,
            a = (val*2^j) - 1;
        );
        if( r == 0,
            a = (val*2^j) + 1;
        );

        if( isprime(a) == 1,
            trovato=1;
            print("Nok at n: ", j);
            j=n+1;
        );

     );

     if(trovato == 0,  
        if( r == 1,  
            out = concat( out, "Probable Riesel's number: k=");
        );

        if( r == 0,  
            out = concat( out, "Probable Sierpinski's number: k=");
        );
        out = concat( out, val);
        out = concat( out, " on ");
        out = concat( out, n);
        out = concat( out, " values ...");

        print(out);

     );


}