8

我如何确定取幂的前 n 位(a b)。

eg: for a = 12, b = 13 & n = 4, the first 4 digits are 1069.
4

5 回答 5

21

通过以下迭代计算 a b :

a 1 = a 1 ,
a 2 = a 2 ,
...
a i = a i ,
...
a b = a b

你有a i+1 = a i ×a。计算每个a i不准确。问题是a b的相对误差小于a 的相对误差n倍。 您希望获得小于10 -n的最终相对误差。因此,每一步的相对误差可能是。删除每一步的最后一位数字。
论坛

例如,a=2b=16n=1。最终相对误差为10 -n = 0.1。每一步的相对误差为0.1/16 > 0.001。因此,3 位数字在每一步都很重要。如果 n = 2,则必须保存 4 位。通用规则:每一步保存 [ n+log 10 b ] 个数字。

2 (1), 4 (2), 8 (3), 16 (4), 32 (5), 64 (6), 128 (7), 256 (8), 512 (9), 1024 (10) → 102、204
(11)、408(12)、816(13)、1632(14)→163、326(15)、652(16)。

答案:6。

该算法的复杂度为O(b)。但是很容易将其更改为O(log b)

于 2010-10-06T15:10:44.187 回答
5

另一种解决方案,使用 log10:

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

int main(int argc, char **argv) {
       int a = 12;
       int b = 13;
       int n = 4;
       double x, y;

       x = b * log10(a);
       y = floor(pow(10, x - floor(x) + n - 1));
       printf("Result: %d\n", (int)y);

       return EXIT_SUCCESS;
}
于 2010-10-08T05:31:07.083 回答
5

n=9 k=3 n^n=387420489 答案应该是 387

在此处输入图像描述

这是同一件事,@RC 在他的代码中所做的。谢谢@RC,我刚刚展示了您的代码的数学表示。

于 2013-08-25T18:15:45.457 回答
1

对于这种情况 - 使用幻数 12,13,4:

#include <sstream>
#include <iomanip>
#include <cmath>

double a = 12;
int b = 13;
double result = std::pow(a,b);

std::stringstream strVal;
strVal.setf( ios::fixed, ios::floatfield );
strVal << result;
std::string output(strVal.str().substr(0,4));

输出=“1069”

std::stringstream intStr(output);
int intVal;
intStr >> intVal;

整数值 = 1069

编辑:这应该适用于结果不溢出的任何组合double

于 2010-10-06T15:00:50.867 回答
0

以编程方式执行此操作的最简单方法是使用字符串流将求幂的结果转换为字符串,然后取 n 个最重要(即左侧)字符。

如果您想要一种没有字符串的方法,那么这将起作用:

#include <iostream>
#include <sstream>
#include <math.h>

using namespace std;

double nDigExp( double a, double b, int n )
{
    stringstream ss;
    ss.setf( ios::fixed, ios::floatfield );
    ss << pow(a,b);
    double ret;
    for ( int i = 0; i < n; ++i) ret = (10 * ret) + (ss.get() - '0'); // Yeuch!!
    return ret;
}

int main( )
{
    double result = nDigExp( 12, 13, 4 );

    cout << result << endl;

    return 0;
}

但这并不是最优雅的代码。我相信你可以改进它。

于 2010-10-06T15:02:04.650 回答