基础练习 阶乘计算——高精度乘法

问题描述

问题描述

输入一个正整数n，输出n!的值。其中n!=1*2*3*…*n。

算法描述

n!可能很大，而计算机能表示的整数范围有限，需要使用高精度计算的方法。
使用一个数组A来表示一个大整数a，A[0]表示a的个位，A[1]表示a的十位，依次类推。

将a乘以一个整数k变为将数组A的每一个元素都乘以k，请注意处理相应的进位。

首先将a设为1，然后乘2，乘3，当乘到n时，即得到了n!的值。

输入格式

输入包含一个正整数n，n<=1000。

输出格式

输出n!的准确值。

样例输入

10

样例输出

3628800

代码如下：

#include<iostream>
#include<sstream>
#include<cstring>
using namespace std;
int main()
{
 char ass[30000],as[30000];
 int a[30000] = { 0 }, b[30000] = { 1 }, c[30000] = { 0 };
 int i = 0, j = 0, k = 0;
 int n,alen = 0,blen = 1,clen = 1;
 cin >> n;
 for (i = 1; i <= n; i++)//阶乘总循环
 {
  stringstream ss;
  alen = 0;
  ss << i;
  ss >> ass;  //当前阶乘数，转换为字符串
  while (ass[alen] != '\0')//获得当前阶乘数的位数
  {
   alen++;
  }
  for (j = 0; j < alen; j++)//将字符串逆序，因为需要和b数组顺序保持一致
  {
   as[j] = ass[alen - j - 1];
  }
  for (j = 0; j < alen; j++)//用数组a装阶乘数
  {
   a[j] = as[j] - '0';
  }
  for(j = 0;j < alen;j++) //利用乘法法则计算
   for (k = 0; k < blen; k++)
   {
    c[j + k] += a[j] * b[k];//用数组c来装结果
   }
  clen = alen + blen - 1; //clen为乘过后结果的“最大”位数
  for (j = 0; j < clen; j++)//进位，保证数组c中每一位都是个位数
  {
    c[j + 1] += c[j] / 10;
    c[j] %= 10;
  }
  if (c[clen] > 0) clen++;
  while (c[clen - 1] == 0 && clen > 1) clen--; //while去0
  for (j = 0; j < clen; j++)//将c赋值给b，同时将c清零
  {
   b[j] = c[j];
   c[j] = 0;
  }
  blen = clen;
  clen = 0;
 }
 for (i = blen - 1; i >= 0; i--)//逆序输出结果
 {
  cout << b[i];
 }
 cout << endl;
 return 0;
}

总体思路：

中心代码是利用乘法法则的那一段循环，其它代码都是由它派生出来的。
乘法法则具体参考小学学习的逐位相乘……
需要注意的是：乘法法则是从个位开始乘，所以要注意数组的排列
另外就是C++数据流将数字转换成字符串
还有就是WHILE去0，这一步很重要，否则会导致输出格式不正确