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.

Showing 1-2 of 2 results.

A382307 Position of start of first run of alternating bit values in the base-2 representation of Pi, or -1 if no such run exists.

Original entry on oeis.org

1, 2, 2, 19, 19, 19, 19, 19, 1195, 1697, 1890, 1890, 1890, 1890, 15081, 63795, 206825, 206825, 206825, 470577, 470577, 557265, 557265, 557265, 557265, 557265, 447666572, 447666572, 699793337, 699793337, 2049646803, 2250772991
Offset: 1

Views

Author

James S. DeArmon, Mar 21 2025

Keywords

Comments

In base-2, Pi is: 11.00100100001111110110101010001... For this sequence, the integer part of Pi is ignored, and the first fractional bit is numbered one.
The bit substring may begin with either 0 or 1.
a(27) > 4*10^6. - Michael S. Branicky, Mar 23 2025
Conjecture: no term is -1 (would follow from a proof that Pi is normal in base 2). - Sean A. Irvine, Mar 30 2025
a(33) > 4*10^9. - Pontus von Brömssen, Apr 06 2025

Examples

			The first alternating bit run is "0", at position 1, so a(1) = 1. The second and third alternating bit runs are "01" and "010", starting at position 2, so a(2) and a(3) are both 2.
a(4)-a(8) = 19 since the binary digits of Pi are "10101010" starting at position 19.
		

Crossrefs

Programs

  • Python
    def binary_not(binary_string):
        return ''.join('1' if bit == '0' else '0' for bit in binary_string)
    # !pip install gmpy2   # may be necessary
    import gmpy2
    gmpy2.get_context().precision = 12000000
    pi = gmpy2.const_pi()
    # Convert Pi to binary representation
    binary_pi = gmpy2.digits(pi, 2)[0] # zero-th element is the string of bits
    binary_pi = binary_pi[2:] # remove leading "11" left of decimal point
    outVec = []
    strSearch0 = "" # this substring starts with "0"
    for lenRun in range(1,30):
      strSearch0 += "0" if lenRun%2==1 else "1"
      strSearch1 = binary_not(strSearch0)
      l0 = binary_pi.find(strSearch0)+1 # incr origin-0 result
      l1 = binary_pi.find(strSearch1)+1
      outVec.append(min(l0,l1))
    print(outVec)

Extensions

a(26)-a(32) from Pontus von Brömssen, Apr 06 2025

A382667 Position of the first instance of prime(n), in base 2, in the binary representation of Pi after the binary point.

Original entry on oeis.org

3, 11, 16, 11, 16, 15, 25, 60, 91, 14, 11, 126, 58, 393, 207, 18, 14, 13, 6, 180, 141, 169, 58, 243, 47, 326, 168, 475, 15, 291, 451, 108, 64, 87, 327, 421, 358, 41, 356, 468, 343, 16, 618, 107, 80, 179, 57, 206, 291, 325, 361, 205, 427, 12, 95, 108, 436, 6, 996
Offset: 1

Views

Author

James S. DeArmon, Apr 02 2025

Keywords

Comments

Positions are numbered starting from 1 for the first bit after the binary point in Pi.

Examples

			For n=19, the bits of Pi and their numbering, after the binary point, begin
          1 2 3 4 5 6 7 8 9 ...
   1 1 .  0 0 1 0 0 1 0 0 0 0 1 1 1 1 ...
                    \-----------/
                    prime(19) = 67
prime(19) = 1000011_2 begins at position a(19) = 6.
prime(58) = 271 = 100001111_2 also starts at 6 => a(58) = 6.
		

Crossrefs

Programs

  • Mathematica
    p=Drop[RealDigits[Pi,2,1010][[1]],2](* increase for n>73 *);a[n_]:=First[SequencePosition[p,IntegerDigits[Prime[n],2]][[1]]] (* James C. McMahon, Apr 26 2025 *)
  • Python
    import gmpy2
    from sympy import isprime
    gmpy2.get_context().precision = 12000000
    gmpy2.get_context().round = gmpy2.RoundDown
    pi = gmpy2.const_pi()
    binary_pi = gmpy2.digits(pi, 2)[0][2:] # Get binary digits and remove "11"
    print([binary_pi.find(bin(cand)[2:])+1 for cand in range(2, 700) if isprime(cand)])

Formula

a(n) = A178707(A000040(n)). - Pontus von Brömssen, Apr 12 2025
Showing 1-2 of 2 results.