高精度乘法详解(阶乘、两数相乘为例)

原理：

我们每个人在小学的时候都学过怎么进行加减乘除。
当我们用编程来模拟这样一个过程的时候，我们只需要把这些步骤分解即可。
比如进行乘法的时候，我们首先都是低位对齐，然后让每一个低位乘上另外一个数的每个位，得到的结果后需要有两个东西，一个本位一个进位，本位实际是结果对10取余，而进位则是与10相除得到结果。
现在我们细分步骤：

1. 低位对齐
2. 以低位为基准对另一个数的各个位相乘
3. 每次与一个位相乘后得到的结果需要有本位(本位<10)和进位
4. 本位是由每个位的乘法对10取余，进位是与10相除取整。

实际上这样理论上的空讲还是无法说清楚，大家可以自己手动计算一个乘法再结合这个过程。比如123*34.

两个大数相乘的代码实现

根据前面分解的步骤，由于第一步是低位对齐然后进行操作(由于是从低位开始操作)。而平时输入的结果实际上前面的位置放的是高位，为了方便低位操作，我们需要先把低位逆序。

实现代码

#include <bits/stdc++.h>
using namespace std;
int main()
{
    //用于逆序存储两个相乘的数
	int a[20000]; 
    int b[20000];
    //用于存储两个数相乘只后的答案
    int res[200000];
    //对这三个数组进行初始化为0
    memset(a, 0, sizeof(a));
    memset(b, 0, sizeof(b));
    memset(res, 0, sizeof(res));
    //下面整段为输入处理
    string s1;
    string s2;
    cin>>s1;
    cin>>s2;
    int n = s1.size();
    int m = s2.size();
    for (int i = n-1,j=0; i >= 0; i--,j++){
        a[j] = s1[i]-'0';
    }
    for (int i = m-1,j=0; i >= 0; i--,j++){
        b[j] = s2[i]-'0';
    }
    //以下整段是对乘法模拟的过程
	int temp;
	for(int i=0;i<n;i++){   int up = 0;
        int j = 0;
		for (; j < m; j++){
        //每次的加法还需要加上res[j+i],这方面大家可自己手算理解为什么要这么做
            temp = a[i]*b[j]+up+res[j+i];
            res[j+i] = temp%10;
            up = temp/10;
        }
        //若该轮乘法进行完后，进位仍然不为0时，故继续将下一位更新
        while(up){
            temp = res[j+i]+up;
            res[j+i] = temp%10;
            up = temp/10;
            j++;
        }
        
	}//倒序输出每一位，从后往前一旦碰到非0，则开始倒序输出。
    bool flag = false;
	for(int i=199999;i>=0;i--){
        if(res[i])
            flag = true;
        if(flag)
            cout<<res[i]; 
    }
        
		
	return 0;
}

求大数阶乘的代码实现

和两个数相乘的代码差不多，只是由于不清楚所要相乘的另一个数的位数，所以用digit变量记录下得到的结果的位数。

#include <bits/stdc++.h>
using namespace std;
int main(){
//创建res数组用于存储阶乘的结果,我这里设置最多存20000位
	int res[20000];
	//初始化为0
	memset(res,0,sizeof(res));
	//输入需要的阶乘n
	int n;
	cin>>n;
//开始更新答案,初始化第一个位的res[0]=1(阶乘的结果至少为1位);
	res[0] = 1;
	int temp,digit=1,up=0;
	for(int i=2;i<=n;i++){
	//计数并每次乘以各个位数
		for(int j=0;j<digit;j++){
			temp = i*res[j]+up;
			res[j] = temp%10;
			up = temp/10;
		}
	//当up进位不为0时说明还能再往前进位，增加进位直到up为0为止
	//这里是和上个不同的地方，因为确定了前面没有够到的位肯定是0，而上一个前面可能为非零,所以还需要与上一位相加。
	while(up){
		res[digit++] = up%10;
		up /= 10;
		}
	}
 //答案更新结束，开始打印答案，从后往前检查非0，然后开始逆序打印。  
 bool flag = false;
 for(int i=19999;i>=0;i--){
	if(res[i])
		flag = true;
	if(flag)
		cout<<res[i];
}   
		
	return 0;
}