Python, https://www.quantamagazine.org/mathematicians-discover-prime-conspiracy-20160313/ test

test.py

# Infinite primes generator and quantum article check:
#   https://www.quantamagazine.org/mathematicians-discover-prime-conspiracy-20160313/

from tests import miller_rabin

def primes_generator():
    i = 2
    while True:
        if miller_rabin(i):
            yield i
        i += 1


if __name__ == '__main__':
    primes = primes_generator()
    nine_ends_pair = 0
    one_ends_pair = 0
    nines = 0
    tmp = None
    primes_length = 0
    for x in primes:
        primes_length += 1
        if primes_length > 10 ** 8: break
        if not tmp:
            tmp = x % 10
            if tmp == 9: nines += 1
        elif tmp:
            a, b = tmp, x % 10
            if tmp == 9: nines += 1
            if a == 9 and b == 9: nine_ends_pair += 1
            if a == 9 and b == 1: one_ends_pair += 1
            tmp = x % 10


    print(nine_ends_pair, one_ends_pair, nines, primes_length)  ## -> Replicates outcomes: 4622916, 7991431, 25000026, 100000001       
                                                                ## For 9 and 1.