text: Title: Add Support for Arbitrary Precision Arithmetic with BigFloat 
  Body: **Is your feature request related to a problem? Please describe.**
  I tried running 'pysr' on a 1,000 row array with 4 integer input variables and one integer output variable - a Goedel Number.
  
  From Mathematica:
  ```
  GoedelNumber[l_List] := Times @@ MapIndexed[Prime[First[#2]]^#1 &, l]
  ```
  E.g. 
  ```
  Data file:
  # 7	1	5	8	6917761200000
  
  julia> 2^7*3^1*5^5*7^8
  6917761200000
  ```
  
  The model returned:
  ```
  Complexity  Loss       Score     Equation
  1           Inf       NaN       0.22984365
  
  ```
  I am just learning 'pysr' and maybe it's just 'user error'.  However, Inf and Nan suggest that Goedel numbers may exceed 
  Float64.
  
  <img width="1430" alt="Screenshot 2022-12-01 at 8 33 44 AM" src="https://user-images.githubusercontent.com/3105499/205108196-cb0c09d2-ccd6-4a9f-978f-44b3032410ae.png">
  
  **Describe the solution you'd like**
  Not sure what happened, because the largest Goedel number in the input is: 
  1.6679880978201e+23
  
  **Additional context**
  I didn't see any parameters to set 'verbose' mode or 'debugging' information.
  
  [GoedelTableFourParameters.txt](https://github.com/MilesCranmer/PySR/files/10134343/GoedelTableFourParameters.txt)
  
  