For learning

from functools import reduce

def generate(initial, factor, rem, scale=1):
    a = initial
    af = factor
    while True:
        a = (a * af) % rem
        if a % scale == 0: yield a

gen = (1 if (x & 0xFFFF) == (y & 0xFFFF) else 0 for x,y in
       zip(generate(618,16807,2**31-1,4),
           generate(814,48271,2**31-1,8)))
s = reduce(lambda x,y: x+y,(next(gen) for _ in range(5000000)))
print(s)