PiBBPBench.py

#!/usr/bin/python
import math
import time

def n16modk (n, k):
    r = 1
    t = 1

    while t <= n:
        t *= 2

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

        t = t / 2

        if t < 1:
            break;

        r = (r * r) % k

    return r

def d16Sj (d, j):
    d16Sj = 0.0
    k = 0

    for k in range(0, d):
        d16Sj += float(n16modk(d - k, 8 * k + j) / 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
d = 0

delta = 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)
delta = time.time() - delta
print(delta)