问题描述:
给定一个包含小数点的在0.0到99.999的包含6位的输入数 r(以字符串的形式),和一个在0到25之间的数 n, 求 r 的 n 次方,要求是,输出结果不包括无效的0,这里的0既包括了在所得结果的最后几位的无效的0,也包括了0.345....中的0,即结果应当是.345....的形式,当所得结果是一个整数时,不应当打印小数点。
这是一个高精度计算的问题,思路是通过模拟乘法过程,求得结果。这里我们将 r 的每一位倒序存储到一个数组中,例如:3.0001(包含了6位),指数n 为 2,那么它的计算思路为 10003 *10003,那么 i 和 j 从0 开始,设结果为 result [i + j] = temp[j] *mul[i]
这一步实现了结果的每一位都是正确的。
这里的乘法并不是像下面的
10009
* 10009
= 90081
+10009
= 100180081
而实际上90001*90001 = 8100180001,显然这两个结果倒序过来是不同的,实际的计算方法应当是:
10009
*10009
=10009
+ 90081
= 100081(这里的9+9的1进位应当相关数组元素下标大的方向进位)0081 (这里还是有些问题,等我回头再想一下), 实际上均满足
result [i + j] = temp[j] *mul[i]的表达式,i,j在这里最大都是4 , result[8] = temp[4]*mul[4]=18(8向后进位为result[9],这样就满足结果的正确性了)
程序ac:
具体的程序为:
#include<iostream>
#include<string>
using namespace std;
int main()
{
string r;
int n, pos_of_point;
const int r_len = 160;
int result[r_len], temp[r_len], mul[6];
while (cin >> r >> n) {
for (int i = 0;i < r_len;i++)
temp[i] = result[i] = 0;
for (int i = 0;i < 6;i++)
mul[i] = 0;
pos_of_point = 0;
size_t pos = r.find(".");
if (pos != string::npos)
pos_of_point = (5 - pos) * n;
for (int i = 5, j = 0;i >= 0;i--)
{
if (r[i] != '.')
{
result[j] = temp[j] = mul[j] = r[i] - '0'; //按照逆序将字符转换为数字
j++;
}
}
while (n >= 2)//当指数为1时,只要按照规定的输出格式输出原数即可
{
--n;
for (int i = 0;i < r_len;i++)
result[i] = 0;
for (int i = 0; i < 5;i++)//用乘数的每一位去乘过程temp值
{
int t;
for (int j = 0;j < r_len;j++)
{
if (mul[i] == 0)
break;
t = temp[j] * mul[i];
result[i + j] += t; //这两句话是核心
//处理进位问题,并且是连续处理进位,如果使得下一位也产生进位,仍要进行计算
for (int k = i + j;result[k] > 9; ++k)
{
result[k + 1] += result[k] / 10;//注意一定是+=因为result[k + 1]在自己的循环时还要加上这个数
result[k] = result[k] % 10;
}
}
}
for (int i = 0;i < r_len;i++)//实际这里可以用一些方法精简计算
temp[i] = result[i];//数据存储在temp数组中
}
//从后向前找到第一个不为0的数
int firstIndex = -1;
for (int i = r_len - 1;i >= pos_of_point;i--)
{
if (result[i] > 0)
{
firstIndex = i;
break;
}
}
int lastIndex = -1;
for (int j = 0;j < pos_of_point;j++)
{
if (result[j] != 0)
{
lastIndex = j;
break;
}
}
if (firstIndex != -1)
{
while (firstIndex >= pos_of_point)
{
cout << result[firstIndex];
--firstIndex;
}
}
if (lastIndex != -1)
{
cout << '.';
--pos_of_point;
while (lastIndex <= pos_of_point)
{
cout << result[pos_of_point];
--pos_of_point;
}
}
cout << endl;
}
}