primebits.frink

Download or view primebits.frink in plain text format


// Program to test for prime numbers following a pattern
// mentioned by John Olsen.
//
//  This tests n-bit primes of the form:
//  1101 1111 1111 1111 1111 1111
//  1111 1111 0111 1111 1111 1111
//  1111 1111 1101 1111 1111 1111
//  1111 1111 1111 1011 1111 1111
//  1111 1111 1111 1111 1101 1111
//  1111 1111 1111 1111 1110 1111
//  1111 1111 1111 1111 1111 1101

// where bits are successively zeroed, and counts the number
// of primes found for each bit-length.

for bits = 3 to 1000
{
   numPrimes = 0
   base = 2^bits - 1
   for clearBit = 1 to bits-2
   {
      num = base - 2^clearBit
      p = isPrime[num]
      if p == true
         numPrimes = numPrimes + 1

     // Comment in the line below to show the actual numbers found.
     // println[(num -> binary) + "\t$num\t$p"]
   }

   percent = format[numPrimes/(bits-2), 1., 12]
   tested = bits-2
   println["$bits\t$numPrimes\t$percent"]
}

Download or view primebits.frink in plain text format

This is a program written in the programming language Frink.
For more information, view the Frink Documentation or see More Sample Frink Programs.

Alan Eliasen was born 20145 days, 13 hours, 41 minutes ago.