PiBBp_bench.py

#!/usr/bin/python

import math
import time

def n16modk(n, k):

  t = 1;

  while t <= n:
    t*=2;

  r = 1;

  while True:
    if n >= t:
      r = (16*r) % k
      n = n - t

    t = t/2;

    if t < 1:
      break;

    r = (r**2) % k;

  return r




def d16Sj(d, j):
    d16sj = 0.0
    var=0;
    for k in range(0,d):
      var = n16modk(d - k, 8 * k + j)
      d16sj += var/float(8 * k + j)

    return d16sj - math.floor(d16sj)


d16Pi = 0.0;
d16S1t = 0.0;
d16S4t = 0.0;
d16S5t = 0.0;
d16S6t = 0.0;

jegy = 0;

kezdes = time.time()

for d in range(100000000, 100000001):

  d16Pi = 0.0

  d16S1t = d16Sj(d, 1)
  d16S4t = d16Sj(d, 4)
  d16S5t = d16Sj(d, 5)
  d16S6t = d16Sj(d, 6)

  d16Pi = 4.0*d16S1t - 2.0*d16S4t - d16S5t - d16S6t

  d16Pi = d16Pi - math.floor(d16Pi)

  jegy = int(math.floor(16.0*d16Pi))

print (jegy)

print (time.time()-kezdes)