将无限浮点数转换为Int[副本]

Question

我试图检查一个数字是否是一个完美的平方。然而，我正在处理非常大的数字，所以python认为它的无穷大出于某种原因。有办法绕过这个吗？这是代码，第一个变量只包含一堆非常非常非常非常大的数字

import math
from decimal import Decimal
def main():
for i in lst:
    root = math.sqrt(Decimal(i))
    print(root)
    if int(root + 0.5) ** 2 == i:
       print(str(i) + " True")

Answer 1

@casevh有一个正确的答案——使用一个可以对任意大整数进行数学运算的库。因为您在寻找平方，所以您可能正在使用整数，有人可能会认为使用浮点类型（包括decimal.decimal）在某种意义上是不雅观的。

您绝对不应该使用Python的float类型；它的精度有限（约16位小数）。如果你用十进制。小数，请小心指定精度（这取决于数字的大小）。

由于Python有一个大整数类型，因此可以编写一个相当简单的算法来检查平方度；请参阅我对这种算法的实现，以及浮点问题的说明，以及如何使用十进制。下面是十进制。

import math
import decimal

def makendigit(n):
    """Return an arbitraryish n-digit number"""
    return sum((j%9+1)*10**i for i,j in enumerate(range(n)))  
x=makendigit(30)

# it looks like float will work...
print 'math.sqrt(x*x) - x: %.17g' % (math.sqrt(x*x) - x)
# ...but actually they won't
print 'math.sqrt(x*x+1) - x: %.17g' % (math.sqrt(x*x+1) - x)

# by default Decimal won't be sufficient...
print 'decimal.Decimal(x*x).sqrt() - x:',decimal.Decimal(x*x).sqrt() - x
# ...you need to specify the precision
print 'decimal.Decimal(x*x).sqrt(decimal.Context(prec=30)) - x:',decimal.Decimal(x*x).sqrt(decimal.Context(prec=100)) - x

def issquare_decimal(y,prec=1000):
    x=decimal.Decimal(y).sqrt(decimal.Context(prec=prec))
    return x==x.to_integral_value()

print 'issquare_decimal(x*x):',issquare_decimal(x*x)
print 'issquare_decimal(x*x+1):',issquare_decimal(x*x+1)

# you can check for "squareness" without going to floating point.
# one option is a bisection search; this Newton's method approach
# should be faster.

# For "industrial use" you should use gmpy2 or some similar "big
# integer" library.
def isqrt(y):
    """Find largest integer <= sqrt(y)"""
    if not isinstance(y,(int,long)):
        raise ValueError('arg must be an integer')
    if y<0:
        raise ValueError('arg must be positive')
    if y in (0,1):
        return y
    x0=y//2
    while True:
        # newton's rule
        x1= (x0**2+y)//2//x0
        # we don't always get converge to x0=x1, e.g., for y=3
        if abs(x1-x0)<=1:
            # nearly converged; find biggest
            # integer satisfying our condition
            x=max(x0,x1)
            if x**2>y:
                while x**2>y:
                    x-=1
            else:
                while (x+1)**2<=y:
                    x+=1
            return x
        x0=x1

def issquare(y):
    """Return true if non-negative integer y is a perfect square"""
    return y==isqrt(y)**2

print 'isqrt(x*x)-x:',isqrt(x*x)-x

print 'issquare(x*x):',issquare(x*x)
print 'issquare(x*x+1):',issquare(x*x+1)

Answer 2

我认为你需要看看BigFloat模块，例如：

import bigfloat as bf
b = bf.BigFloat('1e1000', bf.precision(21))
print bf.sqrt(b)

打印BigFloat。精确（'9.9999993810013282e 499'，精度=53）

Answer 3

将math.sqrt（Decimal（i））替换为Decimal（i）. sqrt（），以防止您的Decimal衰变为浮动