cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A296564 Decimal expansion of lim_{k->infinity} (1/k)*Sum_{i=1..k} A293630(i).

Original entry on oeis.org

1, 2, 7, 5, 2, 6, 1, 8, 4, 2, 0, 9, 1, 1, 7, 2, 1, 3, 5, 9, 2, 8, 4, 7, 7, 2, 0, 4, 7, 8, 0, 1, 5, 1, 5, 1, 4, 9, 3, 4, 7, 6, 0, 0, 3, 7, 1, 0, 7, 4, 9, 0, 7, 5, 4, 2, 7, 6, 0, 2, 6, 3, 7, 6, 4, 9, 3, 5, 5, 3, 7, 1, 6, 7, 4, 1, 8, 5, 8, 7, 6, 2, 1, 9, 0, 0, 4
Offset: 1

Views

Author

Iain Fox, Dec 15 2017

Keywords

Comments

From Jon E. Schoenfield, Dec 23 2017: (Start)
Starting with the sequence S_0 = {1,2} and extending it one pass at a time as described at A293630 (obtaining S_1 = {1,2,1,1}, S_2 = {1,2,1,1,1,2,1}, etc.), let n_j be the number of terms in S_j; then for j=0,1,2,..., n_j = 2, 4, 7, 13, 37, 73, 145, 289, 865, 1729, 3457, 10369, 20737, 41473, 82945, 248833, 497665, ... (see A291481).
In the algorithm implemented in the PARI program, the variable "build" specifies the number of passes during which the terms of S_j are actually built and stored. The algorithm then uses the terms of S_build to compute the number (n_j) of terms in S_j and their total value (t_j) for each j in build+1..build+n_build. For build=0,1,2,..., the number of decimal digits to which the final ratio t_j/n_j at j = build + n_build matches the actual limit 1.275261842... is 2, 3, 4, 7, 15, 29, 54, 105, 306, 608, 1213, 3629, 7253, 14501, 28995, 86974, 173941, ...
Thus, for example, using build=7, the number of 1s and 2s in the last sequence actually stored, i.e., S_7, is 289, but the number of terms n_j and their total value t_j are computed for every j up through j = build+n_build = 7 + n_7 = 7 + 289 = 296 (both n_296 and t_296 are 104-digit numbers) and the final ratio t_296/n_296 matches the actual limit to 105 decimal digits. (End)
From Iain Fox, Dec 23 2017: (Start)
This is the average value of A293630 on the interval n = 1..infinity.
Is this number transcendental? (End)

Examples

			Equals 1.2752618420911721359284772047801515149347600371...
After generating k steps of A293630:
  k = 0:        [1, 2];                  1.500000000000...
  k = 1:        [1, 2, 1, 1];            1.250000000000...
  k = 2:        [1, 2, 1, 1, 1, 2, 1];   1.285714285714...
  k = 3:        [1, 2, 1, 1, 1, 2, ...]; 1.307692307692...
  k = 4:        [1, 2, 1, 1, 1, 2, ...]; 1.270270270270...
  k = 5:        [1, 2, 1, 1, 1, 2, ...]; 1.273972602739...
  k = 6:        [1, 2, 1, 1, 1, 2, ...]; 1.275862068965...
  ...
  k = infinity: [1, 2, 1, 1, 1, 2, ...]; 1.275261842091...
		

Crossrefs

Cf. A293630.

Programs

  • PARI
    gen(build) = {
    my(S = [1, 2], n = 2, t = 3, L, nPrev, E);
    print(S);
    print(1.0*t/n);
    for(j = 1, build, L = S[#S]; n = n*(1+L)-L; t = t*(1+L)-L^2; nPrev = #S; for(r = 1, L, for(i = 1, nPrev-1, S = concat(S, S[i]))); print(S); print(1.0*t/n));
    E = S;
    for(j = build + 1, build + #E, L = E[#E+1-(j-build)]; n = n*(1+L)-L; t = t*(1+L)-L^2; print(1.0*t/n));
    } \\ (gradually increase build to get more precise answers) Iain Fox, Dec 23 2017 with help of Jon E. Schoenfield