[seqfan] Re: Worry about old sequence, A030077, paths in K_n

Pierre Abbat phma at bezitopo.org
Sun Apr 3 18:21:50 CEST 2022

On Saturday, April 2, 2022 9:03:43 AM EDT David Applegate wrote:
> I worry about using floating point (extended or not) to check if
> different sums of square roots are equal or not.  Using finite precision
> for this is extremely tricky This is a notoriously hard problem in
> general.  For example, to see that
> sqrt(7)+sqrt(14)+sqrt(39)+sqrt(70)+sqrt(72)+sqrt(76)+sqrt(85) !=
> sqrt(13)+sqrt(16)+sqrt(55)+sqrt(67)+sqrt(73)+sqrt(79) you already need
> more than double-precision floating point (their difference is 10^-19,
> see
> http://citeseerx.ist.psu.edu/viewdoc/download?doi=
> ype=pdf). There is a randomized polynomial time algorithm for testing
> equality, but it isn't just compute the result to enough precision.

The sum is missing sqrt(46). I've thought of solving a similar problem in my 
program Quadlods, since such nearly equal sums of quadratic irrationals would 
result in planes of points when generating n-dimensional Richtmyer sequences.


lo ponse be lo mruli ku po'o cu ga'ezga roda lo ka dinko

More information about the SeqFan mailing list