pibbpbench

from math import floor
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
    for k in xrange(0, d + 1):
        d16sj += n16modk(d - k, 8 * k + j) / (float)(8 * k + j)
    return d16sj - floor(d16sj)

jegy = 0
start = time.time()
for d in xrange(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 - floor (d16Pi)
    jegy = int(floor(16.0 * d16Pi))

print jegy
print time.time() - start