C++ 仅需要k位时的快速求幂-续 我需要帮助的地方。。。
我现在要做的是将这个解决方案转换为c++，它计算一个数字的值：C++ 仅需要k位时的快速求幂-续 我需要帮助的地方。。。,c++,algorithm,math,exponentiation,C++,Algorithm,Math,Exponentiation,我现在要做的是将这个解决方案转换为c++，它计算一个数字的值： n^m = exp10(m log10(n)) = exp(q (m log(n)/q)) where q = log(10) 从结果中查找前n个数字的方法如下： "the first K digits of exp10(x) = the first K digits of exp10(frac(x)) where frac(x) = the fractional part of x = x - floor(x)." long
n^m = exp10(m log10(n)) = exp(q (m log(n)/q)) where q = log(10)
"the first K digits of exp10(x) = the first K digits of exp10(frac(x))
where frac(x) = the fractional part of x = x - floor(x)."
long double foo = m * logl(n);
foo = fmodl(foo, logl(10.0)) + some_epsilon;
sprintf(some_string, "%.9Lf", expl(foo));
/* boring string parsing code here */
u l l function getPrefix(long double pow /*exponent*/, long double length /*length of prefix*/)
{
long double dummy; //unused but necessary for modf
long double q = log(10);
u l l temp = floor(pow(10.0, exp(q * modf( (pow * log(2)/q), &dummy) + length - 1));
return temp;
}
编辑 我的尝试的输出示例：
n:2 m:0 n^m:1 计算尾数：1.16334
n:2 m:1 n^m:2 计算尾数：2.32667
n:2 m:2 n^m:4 计算尾数：4.65335
n:2 m:98 n^m:3.16913e+29 计算尾数：8.0022
n:2 m:99 n^m:6.33825e+29
计算尾数：2.16596我会避免这一点。众所周知，很难正确实施。有很多这样的问题，人们被标准库中糟糕的
pow"the first K digits of exp10(x) = the first K digits of exp10(frac(x))
where frac(x) = the fractional part of x = x - floor(x)."
long double foo = m * logl(n);
foo = fmodl(foo, logl(10.0)) + some_epsilon;
sprintf(some_string, "%.9Lf", expl(foo));
/* boring string parsing code here */
m log（n）m*logl（n）2e10foolong double一些ε长双精度expl（foo）1e-9鉴于这里的精度困难，我可能建议使用类似MPFR的库，而不是
long doubledouble-doubleexplogfmod9999999999999991000000000000000999…9