r - 如何使用 R 中的 mpfr 包检查一个数字是否是一个完美的正方形？

Question

我有一个行列式，我知道它是整数的平方，但是因为它的值比.Machine$integer.max我使用的 mpfr 包大。

但我仍然有问题。

这是算法：

> a<- mpfr(sqrt(det(M)), precBits=512);a
1 'mpfr' number of precision  512   bits 
[1] 430080000000001.1875

请你帮助我好吗？

Answer 1

性能有问题吗？如果没有，那么以下应该工作。

> x<-mpfr(31415926535897932384626433832795, 500)
> is.whole(sqrt(x))
[1] FALSE


> y<-mpfr(31415926535897932384626433832794, 500)^2
> y
1 'mpfr' number of precision  500   bits
[1] 986960440108935918772069008410384076085841574993068761741787136
> is.whole(sqrt(y))
[1] TRUE