(Python)
from sympy import S
# download https://stuff.mit.edu/afs/sipb/contrib/pi/pi-billion.txt, then
#with open('pi-billion.txt', 'r') as f: pi_digits = f.readline()
pi_digits = str(S.Pi.n(3*10**5+2))[:-2] # alternative to above
pi_digits = pi_digits.replace(".", "")
def afind():
for k in range(6, len(pi_digits)):
sk = str(k)
if sk == pi_digits[k-5:k-5+len(sk)]:
print(k, end=", ")
|