大数阶乘 - C++

大数的阶乘

没办法，今天是贼不爽的一天，就这个大数阶乘搞了我一天，搞得我头昏脑胀，所以今天必须得吐槽一下。不得不说，要做出这题，逻辑思维一定要缜密。废话少说，直接上题。

题目描述
输入一个正整数n，输出n!的值。

其中n!=123*…*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

题目已经提示的很清楚了，要做出这题，就必须得用数组来计数相乘N。否则我也没想到其他的法子。解题思想就是你把每次计算出的结果用数组存储,a[0]存个位，a[1]存十位,a[2]存百位~~~~~~再把存取的结果也就是数组的每一位分别与N相乘，在把每一位相乘结果相加又得到一个新的数组，循环这个过程，直到乘以n;一开始想的就是如何让这些结果相加，然后发现了一个规律，b[i] = b[i] + (a[j]/(int)pow(10,cnt) ) % 10，也就是说，假如你要得到百位结果a[2],就可以用a[2]%10+a[1]/10 + a[0]/100 % 10;通过这种规律就可以写出相加的代码。
c++代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int a[100001];
int b[100001];
int n;
int flag;
int length(int num)
{
   int len = 0;
   if(num == 0)
   return 1;
   while(num)
   {
     num /= 10;
     len ++;
   }
   return len;
}
int pow(int num)
{
   int sta = 1;
   while(num)
   {
    sta = sta * 10;
    num --;
   }
   return sta; 
}
void mult(int N,int len)
{
   if(N == n + 1)
   {
    int temp = 0;
    for(int i = len - 1; i >= 0; i --)
    {
       if(a[i] != 0)
   temp = 1;
    if(temp == 1 || a[i] != 0)
    cout << a[i]; 
 }
    return;
   }
   for(int i = 0; i < len; i ++)
   {
    a[i] = a[i] * N; 
   }
   len = length(a[len - 1]) + len + 1;
   int temp = 0;
   for(int i = 0; i < len ; i ++)
   {
    int cnt = i;
    int flag;
    if(i - 5 >= 0)
    {
      cnt = 5;
      flag = i - 5; 
 }
    else flag = 0;
    for(int j = flag; j <= i; j ++)
    {
      //b[i] = b[i] + (a[j]/(int)pow(10,cnt) + temp) % 10;
      if(flag == 0)
      b[i] = b[i] + (a[j]/pow(cnt)) % 10;
      else   
      b[i] = b[i] + (a[j]/pow(cnt)) % 10; 
   cnt --;  
 }
 b[i] += temp;
 if(b[i] >= 10)
      {
        temp = b[i] / 10;
     b[i] %= 10; 
   }else
   temp = 0; 
   }
   
   for(int i = 0; i < len; i ++)
   {
    a[i] = b[i];
   }
   memset(b,0,sizeof(b));
   mult(N + 1,len);
   return;
}

我这里有些东西做的不够精确，也只是为了过测试点，就比如结果的位数我是没有仔细考虑的，只能让位数增大，去除前面的0，得到结果，还有就是用了两个数组，但毕竟是为了刷题，也没有什么好法子了，直接AC就对了。