[百炼]1001:Exponentiation 的问题

问题描述:


给定一个包含小数点的在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

扫描二维码关注公众号,回复: 3127581 查看本文章

而实际上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;
	}

}

               

猜你喜欢

转载自blog.csdn.net/T2777/article/details/82555962