ROSECODE 449
Machin Like Formulae
One of the most famous formulae for computing is
due to John Machin (1706) who used it to compute 100 decimal places.
We can rewrite it as:
It turns out that many integers like 239 have an inverse which can be decomposed as a linear combination of Arctan of inverses.
A few examples:
The irreductible numbers, those which cannot be decomposed this way are called Störmer numbers.
Thus, the inverse of every number which is not a Störmer number can be expressed as a linear combinaison of Arctan of inverse of smaller irreductible integers.
What makes John Machin's formula particularily efficient is the fact that it contains few terms and the numbers involved are not too small: the convergence of the Arctan serie is then quite fast.
Let's define the cost of a formula as the sum of the inverse of the distinct integers without 1.
Numbers like 21 or 31 are not suitable for a computation of as they do not include the number 1 in their decomposition, their cost will be infinity.
What are the 10 most efficient numbers less than 20000?
Answer format: Comma separated list in ascending order of cost.
Example: 41,46,75,17,32,7,68,93,43,57 (for numbers less than 100)
[My timing : < 10 sec]
We can rewrite it as:
It turns out that many integers like 239 have an inverse which can be decomposed as a linear combination of Arctan of inverses.
A few examples:
The irreductible numbers, those which cannot be decomposed this way are called Störmer numbers.
Thus, the inverse of every number which is not a Störmer number can be expressed as a linear combinaison of Arctan of inverse of smaller irreductible integers.
What makes John Machin's formula particularily efficient is the fact that it contains few terms and the numbers involved are not too small: the convergence of the Arctan serie is then quite fast.
Let's define the cost of a formula as the sum of the inverse of the distinct integers without 1.
Numbers like 21 or 31 are not suitable for a computation of
What are the 10 most efficient numbers less than 20000?
Answer format: Comma separated list in ascending order of cost.
Example: 41,46,75,17,32,7,68,93,43,57 (for numbers less than 100)
[My timing : < 10 sec]