求助~如何计算高精度开方的算法~

xunyuetian

比如根号2，要精确到小数点后1000位以上的算法~

曾经考虑用笔算的算法，但精度不够~郁闷~

etangelsky

楼主是干什么用的阿？
#include <math.h>;
{
……
a=sqrt(2)
}这样还用什么算法吗？有不是在考高等数学。

frog_skywalker

楼主要的是小数点后面1000位，直接sqrt用double也就15位吧

好像有现成的处理大数的程序，基本思想都是用字符串吧
然后自己写一遍计算机计算开方的程序，无穷级数。。。
应该可以实现了

随口说的，可能有未知问题。。。楼主实现了贴心得上来=,=

zerozero_nine

利用 (10*a + b)* (10*a + b) = 100*a*a + 20*a*b + b*b
直接计算, eg. 121 = 100 + 21 = 100*1*1 + 20*1*1 + 1*1
计算机是2进制可改用 (2*a + b)* (2*a + b) = 4*a*a + 4*a*b + b*b

unsigned int SqrRoot(unsigned int n)
{
unsigned int m, q, t;
m=q=0;
    for (int i=0;i<16;i++) {
        q=(q<<2)|((n>;>;(30-(i<<1)))&3);
        m<<=1;
        t=(m<<1)|1;
        if (t<=q) {
            m|=1;
            q-=t;
        }
    }
    return m;
}

上面只用了一个unsigned int去作, 改成数组便成.

john56

无穷级数里面有好方法。

ccf

正数n的平方根可以通过计算一系列近似值来获得,每个近拟值都比前一个更加接近准确值.每一个近似值是1.接下来的近拟值则可以能过下面的公式来获得.
a=(a+n/a)/2

公式在这里，但怎样解决浮点数精度的问题就要问其它各位大大了

songtebo

可以采用大数乘除的办法。有一个MAX???的开源数学软件可以帮你实现。有个encyptlib里面也有大数的乘除。

如果想得到2千位准确，然后将2乘10的4000次方，然后用大数除的办法可以得到你想要的结果，如果没有提供直接的开方的算法，可以用encyptlib编程，不断逼近。

如果找不到大数的算法，那就自己写，不麻烦的：）

(除法也可以用乘法代替)

我脑袋里面一片模糊，但我想能解决问题。那些名字我都忘记了：（