The Buenos Aires constant is 2.92005097731613…
What’s so special about this number? Let’s see when we use it to initialize the following Python script.
s = 2.920050977316134
for _ in range(10):
i = int(s)
print(i)
s = i*(1 + s - i)
What does this print?
2, 3, 5, 7, 11, 13, 17, 19, 23, 29
If we started with the exact Buenos Aires constant and carried out our operations exactly, the procedure above would generate all prime numbers. In actual IEEE 754 arithmetic, it breaks down around Douglas Adams’ favorite number.
The Buenos Aires constant is defined by the infinite sum
As you can tell, the primes are baked into the definition of λ, so the series can’t generate new primes.
I used mpmath to calculate λ to 100 decimal places:
2.920050977316134712092562917112019468002727899321426719772682533107733772127766124190178112317583742
This required carrying out the series defining λ for the first 56 primes. When I carried out the iteration above also to 100 decimal places, it failed on the 55th prime. So basically I got about as many primes out of the computation as I put into it.
Reference: Beyond Pi and e: a Collection of Constants. James Grime, Kevin Knudson, Pamela Pierce, Ellen Veomett, Glen Whitney. Math Horizons, Vol. 29, No. 1 (September 2021), pp. 8–12
If you write the loop as
for _ in range(1, 11):
It avoids the confusion of using i for two things. Or simpler
for _ in range(10):
This reminds me of the compression scheme where you replace text with the offset into pi that contains the ascii encoding of the text…