高精度阶乘--C++实现

本文主要介绍如何将很大的数相乘并按位存入数组中。

题目描述

计算 1*2..*n， 0<n<1000。

分析：众所周知，计算机能计算的结果是有界的，当超过计算机的计算范围就会溢出无法准确的算出，那么在计算这样的数相乘的时候，我们就需要考虑把计算的结果按位相乘，并按位存入数组，这个时候，就把整体相乘分解成了按位相乘，那么其难点就在于如何处理按位相乘后的进位，下面我举一个例子。如2978*25。

先不考虑进位相乘

然后从右边第一位开始逐位只取其个位。把十位以上向前进一位直到没有进位为止。

这样我们就把两个数的乘积逐位相乘并存入了数组中。

所以对于本题的思路就是，n的阶乘，那就从1开始乘，每乘一次将结果数存入数组，循环乘下一位，并更新数组中的内容，直到乘到n位。

具体代码。

#include<iostream>
#include<math.h>
#include<string.h>
using namespace std;
int main()
{
	int n;
	cin >> n;
	long *dp = new long[5000];
	memset(dp, -1, sizeof(long)* 5000);//所有数置-1。主要作用是判断数组中的结果所在位。
	int j = 0;
	dp[0] = 1;
	if (n == 1)
	{
		return 1;//如果输入1，直接返回
	}
	for (int i=2; i <= n; i++)
	{
		j=0;
		while (dp[j]!=-1)//不考虑进位,直接按位相乘并存入数组
		{
			dp[j] = dp[j] * i;
			j++;
		}
		for (j=0; dp[j]!=-1; j++)//考虑进位
		{
			if (dp[j+1]!=-1)//如果进位时的前一位有值。
			{
				dp[j + 1] = dp[j + 1] + dp[j] / 10;
			}
			else {
				if (dp[j] / 10 != 0)//前一位没有值，且存在进位
				{
					dp[j + 1] = dp[j] / 10;
				}
			}
			dp[j] = dp[j] % 10;//当前位取余，这个放在进位计算后，是进位计算是基于未改变当前为时的计算。
		}
	}
	for (int z =j-1;z>=0; z--)
	{
		cout << dp[z];
	}//循环打印结果。由于存入时是先存低位后存高位，因此需要逆序打印。
	return 0;
}