(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(1, len(pi_digits)):
sk = str(k)
if sk == pi_digits[3*k:3*k+len(sk)]:
print(k, end=", ")
|