The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A104865 Length of sections with the same initial and final digits in the decimal expansion of Pi. 7
10, 22, 19, 11, 18, 21, 5, 11, 2, 7, 20, 14, 6, 44, 4, 34, 17, 10, 6, 15, 8, 12, 10, 4, 11, 13, 21, 24, 16, 5, 11, 17, 19, 39, 33, 17, 4, 8, 7, 3, 20, 10, 6, 4, 21, 20, 11, 12, 3, 5, 4, 5, 27, 2, 3, 21, 7, 22, 13, 7, 6, 8, 4, 4, 8, 2, 2, 8, 4, 4, 11, 3, 9, 28, 7, 49, 30, 3, 5, 8, 24, 5, 11, 3 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
Start with decimal expansion of Pi: s0=3,1,4,1,5,9,2,6,5,3,5,8,9,7,9,3,2,3,8,4,6,2,6,4,3,3,8,3,2,7,9,5,0,2,8,8,4,1,9,7,1,6,9,3,9,9,3,7,5,1,0,5,8,2,0,9,7,4,9,4,4,5,...
Divide s0 to the sections with the same final digits: s={3,1,4,1,5,9,2,6,5,3}, {5,8,9,7,9,3,2,3,8,4,6,2,6,4,3,3,8,3,2,7,9,5}, {0,2,8,8,4,1,9,7,1,6,9,3,9,9,3,7,5,1,0}, {5,8,2,0,9,7,4,9,4,4,5},...
Then A104865(n) = length of [s(n)]: 10,22,19,11,...
A104866(n) = final digits of [s(n)]: 3,5,0,5,...
A104867 gives records in A104865: 10,22,44,49,57,70,...
A104871 gives positions of records A104867 in A104865: 1,2,14,76,146,311;
A104868(n) = round value of the running average of A104865(n)= A104868(n) = round(sum[A104865(i),{i,1,n}]/n)= 10,16,17,16,16,17,15,15,13,... A104868(n) => 11 as n => Infinity.
A104869 gives integer values in sum[A104865(i),{i,1,n}]/n: 10,16,17,16,15,14,14,13,11,11,11,11,...
A104870 gives values of n such that sum[A104865(i),{i,1,n}]/n is integer: 1,2,3,5,14,30,32,63,174,216,312,324,...
LINKS
MATHEMATICA
pi = RealDigits[Pi, 10, 1000][[1]]; lst = {}; f := Block[{k = Position[pi, pi[[1]]][[2, 1]]}, pi = Drop[pi, k]; k]; Do[lst = {lst, f}, {n, 80}]; Flatten[lst] (* Robert G. Wilson v, Mar 29 2005 *)
PROG
(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(2*10**5))[:-1] # alternative to above
pi_digits = pi_digits.replace(".", "")
def aupton(terms):
alst, idx = [], 0
for n in range(1, terms+1):
digit, idx, c = pi_digits[idx], idx+1, 2
while pi_digits[idx] != digit:
idx += 1; c += 1
assert idx < len(pi_digits), "increase precision"
alst.append(c); idx += 1
return alst
print(aupton(84)) # Michael S. Branicky, Oct 22 2021
CROSSREFS
Sequence in context: A080454 A226789 A110367 * A063555 A228010 A303745
KEYWORD
nonn,base
AUTHOR
Zak Seidov, Mar 29 2005
EXTENSIONS
More terms from Robert G. Wilson v, Mar 29 2005
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified June 16 07:19 EDT 2024. Contains 373423 sequences. (Running on oeis4.)